torme, toryou, torus

torus.R

# torus

ratio <- sort(c((seq(1, 10, 0.101)) / 10))#$, (seq(1, 10, 2))))

plt <- function(i){
  d <- expand.grid(data.frame(
    R = 1,
    r = 1 * i,
    theta = seq(0, 2 * pi, 0.25),
    phi = seq(0, 2 * pi, 0.25)
  )) |> within({
    x <- (R + r * cos(theta)) * cos(phi)
    y <- (R + r * cos(theta)) * sin(phi)
    z <- r * sin(theta)
  })
  
  rgl::plot3d(d$x, d$y, d$z, box=FALSE, axes=FALSE)
  rgl::snapshot3d(sprintf("foo/foo%s.png", i))
}
sapply(ratio, plt)  

rgl::plot3d(d$x, d$y, d$z, box=FALSE, axes=FALSE, xlab="", ylab="", zlab="")
rgl::snapshot3d(sprintf("foo/foo%s.png", i))

M <- par3d("userMatrix")

movie3d(par3dinterp(time = (0:4)*0.75, userMatrix = list(M, 
                                                         rotate3d(M, pi, 1, 0, 0),
                                                         rotate3d(M, pi, 0, 1, 0),
                                                         rotate3d(M, pi, 1, 0, 1),
                                                         rotate3d(M, pi, 1, 1, 0) ) ),
        duration = 5, movie = "foo", dir = "foo")