Displaying 1 result from an estimated 1 matches for "qcircr".
Did you mean:
qcirc
2010 Jul 19
5
par("uin") ?
...quot;)
ux <- uin[1]
uy <- uin[2]
dx <- x2 - x1
dy <- y2 - y1
alpha <- atan(ux * dx, uy * dy)
scale <- sqrt((ux * dx)^2 + (uy * dy)^2)
if (scale > 5 * rad)
rad <- rad/scale
qcirc <- cbind(cos((0:10) * pi/20), sin((0:10) * pi/20))
qcircr <- cbind(cos((10:0) * pi/20), sin((10:0) * pi/20))
rot <- function(theta) t(cbind(c(cos(theta), sin(theta)),
c(-sin(theta), cos(theta))))
seg1 <- t(t(rad * qcirc %*% rot(-pi/2)) + c(0, rad))
seg4 <- t(t(rad * qcirc) + c(0, 1 - rad))
seg3 <- t(t((rad * qcircr)...