/* big number declarations */
void calc_corner(bf_t target,bf_t p1,double p2,bf_t p3,double p4,bf_t p5)
{
bf_t btmp1, btmp2 ,btmp3;
int saved; saved = save_stack();
btmp1 = alloc_stack(rbflength+2);
btmp2 = alloc_stack(rbflength+2);
btmp3 = alloc_stack(rbflength+2);
/* use target as temporary variable */
floattobf(btmp3, p2);
mult_bf(btmp1,btmp3,p3);
mult_bf(btmp2,floattobf(target, p4),p5);
add_bf(target,btmp1,btmp2);
add_a_bf(target,p1);
restore_stack(saved);
}
int mandelbf_per_pixel()
{
// parm.x = xxmin + col*delx + row*delx2
mult_bf_int(bfparm.x, bfxdel, (U16)col);
mult_bf_int(bftmp, bfxdel2, (U16)row);
add_a_bf(bfparm.x, bftmp);
add_a_bf(bfparm.x, bfxmin);
// parm.y = yymax - row*dely - col*dely2;
// note: in next four lines, bfold is just used as a temporary variable
mult_bf_int(bfold.x, bfydel, (U16)row);
mult_bf_int(bfold.y, bfydel2, (U16)col);
add_a_bf(bfold.x, bfold.y);
sub_bf(bfparm.y, bfymax, bfold.x);
copy_bf(bfold.x, bfparm.x);
copy_bf(bfold.y, bfparm.y);
if ((inside == BOF60 || inside == BOF61) && !nobof)
{
/* kludge to match "Beauty of Fractals" picture since we start
Mandelbrot iteration with init rather than 0 */
floattobf(bfold.x, param[0]); // initial pertubation of parameters set
floattobf(bfold.y, param[1]);
coloriter = -1;
}
else
{
floattobf(bfnew.x, param[0]);
floattobf(bfnew.y, param[1]);
add_a_bf(bfold.x, bfnew.x);
add_a_bf(bfold.y, bfnew.y);
}
// square has side effect - must copy first
copy_bf(bfnew.x, bfold.x);
copy_bf(bfnew.y, bfold.y);
// Square these to rbflength bytes of precision
square_bf(bftmpsqrx, bfnew.x);
square_bf(bftmpsqry, bfnew.y);
return (1); // 1st iteration has been done
}
int
JuliaZpowerbfFractal()
{
BFComplex parm2;
int saved;
saved = save_stack();
parm2.x = alloc_stack(bflength+2);
parm2.y = alloc_stack(bflength+2);
floattobf(parm2.x, param[2]);
floattobf(parm2.y, param[3]);
ComplexPower_bf(&bfnew, &bfold, &parm2);
add_bf(bfnew.x, bfparm.x, bfnew.x);
add_bf(bfnew.y, bfparm.y, bfnew.y);
restore_stack(saved);
return bigfltbailout();
}
void drawbox(int drawit)
{
struct coords tl,bl,tr,br; /* dot addr of topleft, botleft, etc */
double tmpx,tmpy,dx,dy,rotcos,rotsin,ftemp1,ftemp2;
double fxwidth,fxskew,fydepth,fyskew,fxadj;
bf_t bffxwidth, bffxskew, bffydepth, bffyskew, bffxadj;
int saved;
if (zwidth==0) { /* no box to draw */
if (boxcount!=0) { /* remove the old box from display */
clearbox(); /* asm routine */
boxcount = 0; }
reset_zoom_corners();
return; }
if(bf_math)
{
saved = save_stack();
bffxwidth = alloc_stack(rbflength+2);
bffxskew = alloc_stack(rbflength+2);
bffydepth = alloc_stack(rbflength+2);
bffyskew = alloc_stack(rbflength+2);
bffxadj = alloc_stack(rbflength+2);
}
ftemp1 = PI*zrotate/72; /* convert to radians */
rotcos = cos(ftemp1); /* sin & cos of rotation */
rotsin = sin(ftemp1);
/* do some calcs just once here to reduce fp work a bit */
fxwidth = sxmax-sx3rd;
fxskew = sx3rd-sxmin;
fydepth = sy3rd-symax;
fyskew = symin-sy3rd;
fxadj = zwidth*zskew;
if(bf_math)
{
/* do some calcs just once here to reduce fp work a bit */
sub_bf(bffxwidth,bfsxmax,bfsx3rd);
sub_bf(bffxskew,bfsx3rd,bfsxmin);
sub_bf(bffydepth,bfsy3rd,bfsymax);
sub_bf(bffyskew,bfsymin,bfsy3rd);
floattobf(bffxadj, fxadj);
}
/* calc co-ords of topleft & botright corners of box */
tmpx = zwidth/-2+fxadj; /* from zoombox center as origin, on xdots scale */
tmpy = zdepth*finalaspectratio/2;
dx = (rotcos*tmpx - rotsin*tmpy) - tmpx; /* delta x to rotate topleft */
dy = tmpy - (rotsin*tmpx + rotcos*tmpy); /* delta y to rotate topleft */
/* calc co-ords of topleft */
ftemp1 = zbx + dx + fxadj;
ftemp2 = zby + dy/finalaspectratio;
tl.x = (int)(ftemp1*(dxsize+PIXELROUND)); /* screen co-ords */
tl.y = (int)(ftemp2*(dysize+PIXELROUND));
xxmin = sxmin + ftemp1*fxwidth + ftemp2*fxskew; /* real co-ords */
yymax = symax + ftemp2*fydepth + ftemp1*fyskew;
if(bf_math)
{
calc_corner(bfxmin,bfsxmin,ftemp1,bffxwidth,ftemp2,bffxskew);
calc_corner(bfymax,bfsymax,ftemp2,bffydepth,ftemp1,bffyskew);
}
/* calc co-ords of bottom right */
ftemp1 = zbx + zwidth - dx - fxadj;
ftemp2 = zby - dy/finalaspectratio + zdepth;
br.x = (int)(ftemp1*(dxsize+PIXELROUND));
br.y = (int)(ftemp2*(dysize+PIXELROUND));
xxmax = sxmin + ftemp1*fxwidth + ftemp2*fxskew;
yymin = symax + ftemp2*fydepth + ftemp1*fyskew;
if(bf_math)
{
calc_corner(bfxmax,bfsxmin,ftemp1,bffxwidth,ftemp2,bffxskew);
calc_corner(bfymin,bfsymax,ftemp2,bffydepth,ftemp1,bffyskew);
}
/* do the same for botleft & topright */
tmpx = zwidth/-2 - fxadj;
tmpy = 0.0-tmpy;
dx = (rotcos*tmpx - rotsin*tmpy) - tmpx;
dy = tmpy - (rotsin*tmpx + rotcos*tmpy);
ftemp1 = zbx + dx - fxadj;
ftemp2 = zby + dy/finalaspectratio + zdepth;
bl.x = (int)(ftemp1*(dxsize+PIXELROUND));
bl.y = (int)(ftemp2*(dysize+PIXELROUND));
xx3rd = sxmin + ftemp1*fxwidth + ftemp2*fxskew;
yy3rd = symax + ftemp2*fydepth + ftemp1*fyskew;
if(bf_math)
{
calc_corner(bfx3rd,bfsxmin,ftemp1,bffxwidth,ftemp2,bffxskew);
calc_corner(bfy3rd,bfsymax,ftemp2,bffydepth,ftemp1,bffyskew);
restore_stack(saved);
}
ftemp1 = zbx + zwidth - dx + fxadj;
ftemp2 = zby - dy/finalaspectratio;
tr.x = (int)(ftemp1*(dxsize+PIXELROUND));
tr.y = (int)(ftemp2*(dysize+PIXELROUND));
if (boxcount!=0) { /* remove the old box from display */
clearbox(); /* asm routine */
boxcount = 0; }
if (drawit) { /* caller wants box drawn as well as co-ords calc'd */
#ifndef XFRACT
/* build the list of zoom box pixels */
addbox(tl); addbox(tr); /* corner pixels */
addbox(bl); addbox(br);
drawlines(tl,tr,bl.x-tl.x,bl.y-tl.y); /* top & bottom lines */
drawlines(tl,bl,tr.x-tl.x,tr.y-tl.y); /* left & right lines */
#else
boxx[0] = tl.x + sxoffs;
boxy[0] = tl.y + syoffs;
boxx[1] = tr.x + sxoffs;
boxy[1] = tr.y + syoffs;
boxx[2] = br.x + sxoffs;
boxy[2] = br.y + syoffs;
boxx[3] = bl.x + sxoffs;
boxy[3] = bl.y + syoffs;
boxcount = 1;
#endif
dispbox(); /* asm routine to paint it */
}
}