Displaying 1 result from an estimated 1 matches for "qcirc".
Did you mean:
vcirc
2010 Jul 19
5
par("uin") ?
...y2 = 1, right = TRUE, rad = 0.2)
{
uin <- par("uin")
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...