Some tubular objects with rgl (R)

.gitignore

*.png
.Rproj.user
.Rhistory

cyclideKnot.R

library(rgl)

p <- 2
q <- 3
epsilon <- 0.1

# cyclidoidal knot curve
ccld <- function(t, a, c, d){
  u <- p*t
  v <- q*t
  d <- d + epsilon
  b <- sqrt(a*a-c*c)
  den <- a - c * cos(u) * cos(v)
  cbind(
    (d * (c - a * cos(v) * cos(u)) + b * b * cos(v)) / den,
    (b * sin(v) * (a - d * cos(u))) / den,
    (b * sin(u) * (c * cos(v) - d)) / den
  )
}

t_ <- seq(0, 2*pi, length.out = 200L)
pts <- ccld(t_, a = 1.94, c = 0.34, d = 0.56)
knotMesh <- addNormals(cylinder3d(pts, radius = epsilon, sides = 30L))

# cyclide mesh
ccldMesh <- cgalMeshes::cyclideMesh(1.94, 0.34, 0.56, nu = 150L, nv = 100L)

# plot
open3d(windowRect = c(50, 50, 562, 562))
bg3d(rgb(54, 57, 64, maxColorValue = 255))
view3d(-15, -50, zoom = 0.75)
shade3d(ccldMesh, col = "gold")
shade3d(knotMesh, col = "navy")

# animation ####
movie3d(spin3d(axis = c(0, 0, 1), rpm = 10),
        duration = 6, fps = 20,
        movie = "zzpic", 
        dir = ".",
        convert = FALSE,
        webshot = FALSE,
        startTime = 1/20)

library(gifski)
gifski(
  png_files = Sys.glob("zzpic*.png"),
  gif_file = "cyclideKnot.gif",
  width = 512,
  height = 512,
  delay = 1/15
)

pngfiles <- Sys.glob("zzpic*.png")
file.remove(pngfiles)

cyclideKnot.gif

cyclidoidalRotoid.gif

cyclidoidalRotoid.R

cyclidoidalRotoid_byInversion.gif

cyclidoidalRotoid_byInversion.R

cyclidoidalRotoid_byInversion2.gif

cyclidoidalRotoid_byInversion2.R

cyclidoidalTubularHelix.gif

cyclidoidalTubularHelix.R

gist_tubular.Rproj

helicoidalCurve.gif

helicoidalCurve.R

toroidalTubularHelix.gif

toroidalTubularHelix.R

torusKnotWrappedOnTorus.R

torusWrappedInHelix.gif

torusWrappedInHelix.R

trefoilKnotWrappedOnTorus.gif

tubularKnot.R

tubularKnot_p2_q8.gif

tubularKnot_p3_q5.gif