comment { Carr3400.ufm Version 1.6 Ultra Fractal conversions of Robert Carr's formulas 3400 - 3499 by Erik Reckase, February 24, 2000 History: Version 1.6 - minor optimization tweaks + cleanup Version 1.5 - better hints + fix over-optimization optimization tweeks + long line trunc Version 1.4 - major optimization overhaul method = multipass Version 1.3 - Removed unused whitesq code Optimized bailout condition Reduced iteration tests Version 1.2 - Made the following optimizations Carr3417,3459-3460,3464 3465-3466,3469-3471 : flip(flip) removed Carr3401,3412,3417-3418 3426,3428-3430,3452 3459-3460,3464-3466 3486-3487,3491-3492 : conj(conj) removed Carr3407,3413,3421-3424 3427,3431,3446,3472 3473-3478,3481-3482 3491,3497-3499 : abs(conj) -> abs() Version 1.1 - Modified all whitesq code to eliminate warns Version 1.0 - Initial Release Original formulas pulled from 43hCarr.frm 43iCarr.frm 43jCarr.frm Carr3400-3499 : pixel fix Carr3400-3499 : removed unnecessary code Carr3431 : variables not used, removed Carr3436 : | | fix Carr3400,3412,3447,3490 : cosxx fix } Carr3400 (YAXIS) {; Updated for UF2 by Erik Reckase, February 2000 ; Modified Sylvie Gallet frm. init: pixel2=-abs(real(pixel))+flip(imag(pixel)) b6=conj(.025/pixel2), b7=flip(.01/pixel2) b8=cabs(.002/pixel2+cabs(.002/pixel2)) b3=exp(.0045/pixel2)/conj(cos(.1-pixel2))-conj(.01/pixel2) b2=(b8-b6-b7)/.65-.002/pixel2 b10a=pixel2-flip(.01/pixel2)-tanh(.01/pixel2)-conj(.01/pixel2) b10=conj(pixel2/3)*sinh(pixel2)^2*tanh(pixel2) \ *conj(abs(pixel2/2)) z=conj(b10a*b10)^2-b2*.5+b6+b7-.02*b3+.075 IF (@p3==0) z = z^.8456 ELSEIF (@p3==2) IF (whitesq) z = z^.8456 ENDIF ENDIF l1=round(real(p1)),l2=round(imag(p1)) l3=round(real(p2)),l4=round(imag(p2)) float bailout=16,int iter=0, float zro=0 loop: IF ((iter==l1) || (iter==l2) || (iter==l3) || (iter==l4)) z = z*zro ENDIF z=z^2+(-.7602073214,.080996112) iter=iter+1 bailout: |z|<=bailout default: title = "Carr 3400" periodicity = 0 maxiter = 700 magn = 1 center = (0,0) method = multipass param p1 caption = "Iter Limits 1,2" default = (200,300) hint = "0< Re(P1) < Im(P1) < Re(P2) < Im(P2)