END OF AN ERA, FRACTALFORUMS.COM IS CONTINUED ON FRACTALFORUMS.ORG

//MsltoeRiemannSphere     Variation1----------------------------------
//  http://www.fractalforums.com/new-theories-and-research/revisiting-the-riemann-sphere-%28again%29/
void RiemannSphereMsltoeV1Iteration(CVector3 &z, const cFractal *fractal)
{
  double r = z.Length();
  if (r < 1e-21) r = 1e-21;
  z *= fractal->transformCommon.scale/r;
  double w = 1.0 - z.z;
  if (w > -1e-21 && w < 1e-21) w = (w > 0) ? 1e-21 : -1e-21;
  w = 1.0/w;
  CVector3 t3;
  t3.x = z.x * w;
  t3.y = z.y * w;

  t3.z = (r - 1.5) * (1.0 + t3.x * t3.x + t3.y * t3.y);

  if (fractal->transformCommon.rotationEnabled)
  {
   t3 = fractal->transformCommon.rotationMatrix.RotateVector(t3);
  }

  t3.x = t3.x - floor(t3.x + 0.5);
  t3.y = t3.y - floor(t3.y + 0.5);

  z = t3 * fractal->transformCommon.constantMultiplier441;

  z += fractal->transformCommon.additionConstant000;
}