The function torusUnif samples n points from the 3D torus, uniformly with respect to its surface.
Usage
torusUnif(n, a, c)
Value
The function torusUnif returns an n by 3 matrix of coordinates.
Arguments
n
an integer specifying the number of points in the sample.
a
the radius of the torus tube.
c
the radius from the center of the hole to the center of the torus tube.
Author
Fabrizio Lecci
Details
This function torusUnif is an implementation of Algorithm 1 in the reference.
References
Diaconis P, Holmes S, and Shahshahani M (2013). "Sampling from a manifold." Advances in Modern Statistical Theory and Applications: A Festschrift in honor of Morris L. Eaton. Institute of Mathematical Statistics, 102-125.