Direct surface.py
From Werner KRAUTH
This page presents the program direct_surface.py, a direct-sampling algorithm for uniform points on the surface of a d-dimensional unit sphere
Contents |
[edit]
Description
[edit]
Program
import random, math dimensions = 5 nsamples = 20 for sample in xrange(nsamples): R = [random.gauss(0.0, 1.0) for d in xrange(dimensions)] radius = math.sqrt(sum(x ** 2 for x in R)) print [x / radius for x in R]
[edit]
Version
See history for version information.