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

double PrimitiveSimpleWater(CVector3 point, double height, double amplitude1, double amplitude2, double rotation, int nangwave, double animSpeed, int frame)
{
        CVector3 plane(0,0,-1);
        CVector3 centre(0,0,height);
        plane.Normalize();
        CVector3 ptc = point - center;
        double planeDistance = plane.Dot( ptc );
        if(planeDistance < max(amplitude1,amplitude2) * 10.0)
        {
                CRotationMatrix rotMatrix;
                rotMatrix.RotateZ(rotation/180*M_PI);
                point = rotMatrix.RotateVector(point);

                double phase = animSpeed * frame;
                double dorig = hypot(ptc.x,ptc.y);
                double waveX = sin(phase  +  dorig ); // interesting to add another phase and try to multiply dorig by something
                double waveY = sin(phase  + nangwave*atan2(ptc.y,ptc.x)); // interesting to add another phase
                waveY = waveY*(1.-1./(1.+dorig)); // less perturbed in the center
                planeDistance += amplitude1*waveX + amplitude2*waveY;
                // interesting to add another wave (interference effect) ... etc
        }

        return planeDistance;
}