public Julia( float global_real, float global_imag, float global_theta )
{
this.global_real = global_real;
this.global_imag = global_imag;
this.global_theta = global_theta;
oc.r = global_real;
oc.i = global_imag;
oc.j = oc.k = 0.0f;
eio.r = (float)System.Math.Cos( global_theta );
eio.i = (float)System.Math.Sin( global_theta );
eio.j = 0.0f;
eio.k = 0.0f;
emio.r = (float)System.Math.Cos( -global_theta );
emio.i = (float)System.Math.Sin( -global_theta );
emio.j = 0.0f;
emio.k = 0.0f;
/***
*** multiply eio*c only once at the beginning of iteration
*** (since q |-> sqrt(eio*(q-eio*c)))
*** q -> e-io*q^2 - eio*c
***/
qmult( ref c, ref eio, ref oc );
}
}