/*
 * To change this template, choose Tools | Templates
 * and open the template in the editor.
 */

/**
 *
 * @author nbriere
 */
public class MandelBulb extends Object3D
{
    public boolean inside(double x, double y, double z, boolean transform) // , double d)
    {
        double d = 1; // Oops not so easy
        
        double posX = x;
        double posY = y;
        double posZ = z;

        double dr = 1.0;
        double r = 0.0;

        for (int i = 0; i < 10; i++)
        {
            double dist2xy = x*x + y*y;
            r = Math.sqrt(dist2xy + z*z);
            if (r > 10000)
            {
                break;
            }

            // convert to polar coordinates
            //double theta = Math.acos(z / r);
            double theta = Math.atan2(-z, Math.sqrt(dist2xy));
            double phi = Math.atan2(y, x);
            double zr = Math.pow(r, 8);
            dr = zr / r * 8 * dr + 1.0;

            // scale and rotate the point
            theta = theta * 8;
            phi = phi * 8;
            
            // convert back to cartesian coordinates
            double sintheta = zr * Math.sin(theta);
            x = sintheta * Math.cos(phi);
            y = Math.sin(phi) * sintheta;
            z = zr * Math.cos(theta);

            x += posX;
            y += posY;
            z += posZ;
        }
        
        return ((0.5 * Math.log(r) * r / dr) < d);
    }
}