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

Windows



http://www.fractalforums.com/index.php?action=gallery;sa=view;id=15406

Formula:

void AmazingOct(inout dvec4 p, float scale, float Fold, float m2, inout double rr)
{
   Menger(p, 2.0, vec3(1.0,0.5,0.5));
   rr = dot(p.xyz, p.xyz);
   if (rr<double(m2))
   {
      p = p / double(m2);
   }
   else
   {
      if (rr<1.0)
         p = p / rr;
   }
   p *= double(scale);
}

double DE(dvec3 p)
{
   double rr = 1.0;
   for (int i = 0; i < 30; i++)
   {
      AmazingOct(p, 0.8, 1.0, 0.0, rr);
      if (rr>1000.0)
         break;
   }
   return 0.35*length(p.xyz) / p.w;
}

Cool