Vector3 sample_points[32];
for(unsigned i=0; i<n; ++i)
{
- float z = static_cast<float>(i)/n;
- float r = sqrt(1.0f-z*z);
+ float r = static_cast<float>(i)/n;
+ float z = sqrt(1.0f-r*r);
float d = radical_inverse(i);
Geometry::Angle<float> a = Geometry::Angle<float>::from_turns(d);
sample_points[i] = Vector3(cos(a)*r, sin(a)*r, z)*(0.1f+0.9f*d*d);