Albrecht, I may have missed your saying this, but are the images generated by the parameter files you have posted available online somewhere? Also, there are some interesting complexities in the formulas. Could you say anything about your process of generating the formulas? Thanks, Tim On Sat, Jan 16, 2021 at 3:08 PM Albrecht Niekamp <xxfractalist@outlook.com> wrote:
The Splash1 formula possibly is not available everywhere, so I include it now: FRM:splash_1 { ; Basic data : FractMask08 ; Albrecht Niekamp March,03 ; Lake effect Sylvie Gallet ; p1: Julia seed ; real (p2) : water level in % of the screen height ; imag (p2) : amplitude of the wave (try 0.2) ; real (p3) : frequency (try 300) ; imag (p3) : size y (>0 smaller -1 no text) ; real (p4) : horizontal factor ; imag (p4) : vertical factor ; real (p5) : size left (>0 smaller) ; imag (p6) : size right (>0 smaller) in=0 t=0 z=0 w=0 start=1 c=pixel if (whitesq) row=imag(scrnpix) rows=imag(scrnmax) y=row/(rows-1) col=real(scrnpix) cols=real(scrnmax) x=col/(cols-1) size=imag(p3)/1000 hf=real(p4)/100 vf=imag(p4)/1000 cl=real(p5)/100 cr=imag(p5)/100 tx=(imag(p3)!=-1) if (tx) xx=x yy=y f=1-y d=x-0.5 cr=cr+hf cl=cl+hf If (x>0.5) e=d*d-0.25 hd=-e*hf g=(0.25-e)/5 x=x+g*cr+((hf*0.005)*(hf<0)) else e=x*x-0.25 hd=e*hf*2 g=(0.25-e)/10 x=x-g*cl+((hf*0.01)*(hf<0)) endif d=1-xx if (y>0.5) if (xx>0.5) e=d*d b=f*(0.5-d)*vf*40*(0.25-e) else e=xx*xx b=f*(0.5-xx)*vf*40*(0.25-e) endif y=y+f*vf*1.6+size else if (xx>0.5) e=d*d b=-(f*(0.5-d)*vf*40*(0.25-e)) else e=xx*xx b=-(f*(0.5-xx)*vf*40*(0.25-e)) endif y=y-y*vf*1.6-size endif e=exp(hd) b=exp(b) in=in+(x>.06*e&&x<.92*e&&y<.693*b&&y>.68*b) in=in+(x>.07*e&&x<.13*e&&y<.626*b&&y>.613*b) in=in+(x>.19*e&&x<.25*e&&y<.626*b&&y>.613*b) in=in+(x>.32*e&&x<.41*e&&y<.626*b&&y>.6*b) in=in+(x>.45*e&&x<.5*e&&y<.626*b&&y>.613*b) in=in+(x>.56*e&&x<.6*e&&y<.626*b&&y>.6*b) in=in+(x>.63*e&&x<.69*e&&y<.626*b&&y>.613*b) in=in+(x>.7*e&&x<.75*e&&y<.626*b&&y>.613*b) in=in+(x>.76*e&&x<.81*e&&y<.626*b&&y>.613*b) in=in+(x>.85*e&&x<.89*e&&y<.626*b&&y>.6*b) in=in+(x>.08*e&&x<.12*e&&y<.613*b&&y>.266*b) in=in+(x>.2*e&&x<.24*e&&y<.613*b&&y>.36*b) in=in+(x>.31*e&&x<.32*e&&y<.613*b&&y>.48*b) in=in+(x>.44*e&&x<.47*e&&y<.613*b&&y>.373*b) in=in+(x>.49*e&&x<.51*e&&y<.613*b&&y>.586*b) in=in+(x>.55*e&&x<.56*e&&y<.613*b&&y>.346*b) in=in+(x>.6*e&&x<.61*e&&y<.613*b&&y>.52*b) in=in+(x>.64*e&&x<.68*e&&y<.613*b&&y>.36*b) in=in+(x>.71*e&&x<.75*e&&y<.613*b&&y>.36*b) in=in+(x>.76*e&&x<.8*e&&y<.613*b&&y>.373*b) in=in+(x>.84*e&&x<.85*e&&y<.613*b&&y>.346*b) in=in+(x>.89*e&&x<.9*e&&y<.613*b&&y>.52*b) in=in+(x>.32*e&&x<.35*e&&y<.6*b&&y>.573*b) in=in+(x>.36*e&&x<.42*e&&y<.6*b&&y>.586*b) in=in+(x>.56*e&&x<.59*e&&y<.6*b&&y>.306*b) in=in+(x>.61*e&&x<.62*e&&y<.6*b&&y>.52*b) in=in+(x>.85*e&&x<.88*e&&y<.6*b&&y>.306*b) in=in+(x>.9*e&&x<.91*e&&y<.6*b&&y>.52*b) in=in+(x>.36*e&&x<.4*e&&y<.586*b&&y>.373*b) in=in+(x>.5*e&&x<.51*e&&y<.586*b&&y>.572*b) in=in+(x>.32*e&&x<.34*e&&y<.573*b&&y>.466*b) in=in+(x>.43*e&&x<.44*e&&y<.573*b&&y>.413*b) in=in+(x>.14*e&&x<.15*e&&y<.52*b&&y>.373*b) in=in+(x>.34*e&&x<.35*e&&y<.52*b&&y>.466*b) in=in+(x>.13*e&&x<.14*e&&y<.493*b&&y>.4*b) in=in+(x>.35*e&&x<.36*e&&y<.48*b&&y>.453*b) in=in+(x>.12*e&&x<.13*e&&y<.466*b&&y>.426*b) in=in+(x>.5*e&&x<.51*e&&y<.412*b&&y>.36*b) in=in+(x>.25*e&&x<.28*e&&y<.44*b&&y>.36*b) in=in+(x>.31*e&&x<.35*e&&y<.44*b&&y>.386*b) in=in+(x>.28*e&&x<.29*e&&y<.426*b&&y>.373*b) in=in+(x>.24*e&&x<.25*e&&y<.413*b&&y>.36*b) in=in+(x>.49*e&&x<.5*e&&y<.413*b&&y>.373*b) in=in+(x>.16*e&&x<.18*e&&y<.4*b&&y>.266*b) in=in+(x>.32*e&&x<.34*e&&y<.386*b&&y>.373*b) in=in+(x>.48*e&&x<.49*e&&y<.386*b&&y>.36*b) in=in+(x>.75*e&&x<.76*e&&y<.386*b&&y>.36*b) in=in+(x>.19*e&&x<.2*e&&y<.373*b&&y>.36*b) in=in+(x>.33*e&&x<.39*e&&y<.373*b&&y>.36*b) in=in+(x>.45*e&&x<.48*e&&y<.373*b&&y>.36*b) in=in+(x>.54*e&&x<.55*e&&y<.373*b&&y>.36*b) in=in+(x>.59*e&&x<.62*e&&y<.373*b&&y>.36*b) in=in+(x>.63*e&&x<.64*e&&y<.373*b&&y>.36*b) in=in+(x>.7*e&&x<.71*e&&y<.373*b&&y>.36*b) in=in+(x>.76*e&&x<.79*e&&y<.373*b&&y>.36*b) in=in+(x>.83*e&&x<.84*e&&y<.373*b&&y>.36*b) in=in+(x>.88*e&&x<.91*e&&y<.373*b&&y>.36*b) in=in+(x>.15*e&&x<.16*e&&y<.346*b&&y>.266*b) in=in+(x>.14*e&&x<.15*e&&y<.306*b&&y>.266*b) in=in+(x>.57*e&&x<.59*e&&y<.306*b&&y>.266*b) in=in+(x>.64*e&&x<.68*e&&y<.306*b&&y>.266*b) in=in+(x>.86*e&&x<.88*e&&y<.306*b&&y>.266*b) in=in+(x>.07*e&&x<.08*e&&y<.28*b&&y>.266*b) in=in+(x>.12*e&&x<.14*e&&y<.28*b&&y>.266*b) endif if (in) w=1 z=(xx+yy/8)*689 else if (imag(p2)>0 && imag(p2)<=100) level=real(p2)/100 ; water level ampl=imag(p2) ; amplitude of the wave freq=real(p3) ; frequency angle=real(rotskew*pi/180) exp_irot=exp(-flip(angle)) h=1/real(magxmag) q=h/0.75*imag(magxmag) tanskew=tan(imag(rotskew*pi/180)) u=2*q*exp_irot v=2*h*(tanskew+flip(1))*exp_irot z3rd=center+(-q-h*tanskew-flip(h))*exp_irot z=pixel-z3rd s=imag(conj(u)*z)/imag(conj(u)*v) if (s<=level) dy=level-s z=z+2*dy*(1+ampl*sin(freq*dy^0.2))*v endif z=z+z3rd if (ismand) c=z z=0 else c=p1 endif endif endif endif : if (w) z=z-1 phc_bailout=|z|>=1 else z=z^2+c if (start) t=t+1 if (|z|>=128) u=2*(exp(t/768)) z=z^u start=0 endif if (t>513) u=2*(sin(t/768)) z=z*u c=p1*u start=0 endif endif phc_Bailout=|z|<=2048 endif phc_bailout } ------------end------- _______________________________________________ Fractint mailing list Fractint@mailman.xmission.com https://mailman.xmission.com/cgi-bin/mailman/listinfo/fractint