コード例 #1
0
ファイル: ZOOM.C プロジェクト: mmammel/projects
static void _fastcall zmo_calcbf(bf_t bfdx, bf_t bfdy,
    bf_t bfnewx, bf_t bfnewy,bf_t bfplotmx1, bf_t bfplotmx2, bf_t bfplotmy1,
    bf_t bfplotmy2, bf_t bfftemp)
{
    bf_t btmp1, btmp2, btmp3, btmp4, btempx, btempy ;
    bf_t btmp2a, btmp4a;
    int saved; saved = save_stack();

    btmp1  = alloc_stack(rbflength+2);
    btmp2  = alloc_stack(rbflength+2);
    btmp3  = alloc_stack(rbflength+2);
    btmp4  = alloc_stack(rbflength+2);
    btmp2a = alloc_stack(rbflength+2);
    btmp4a = alloc_stack(rbflength+2);
    btempx = alloc_stack(rbflength+2);
    btempy = alloc_stack(rbflength+2);

    /* calc cur screen corner relative to zoombox, when zoombox co-ords
       are taken as (0,0) topleft thru (1,1) bottom right */

    /* tempx = dy * plotmx1 - dx * plotmx2; */
    mult_bf(btmp1,bfdy,bfplotmx1);
    mult_bf(btmp2,bfdx,bfplotmx2);
    sub_bf(btempx,btmp1,btmp2);

    /* tempy = dx * plotmy1 - dy * plotmy2; */
    mult_bf(btmp1,bfdx,bfplotmy1);
    mult_bf(btmp2,bfdy,bfplotmy2);
    sub_bf(btempy,btmp1,btmp2);

    /* calc new corner by extending from current screen corners */
    /* *newx = sxmin + tempx*(sxmax-sx3rd)/ftemp + tempy*(sx3rd-sxmin)/ftemp; */
    sub_bf(btmp1,bfsxmax,bfsx3rd);
    mult_bf(btmp2,btempx,btmp1);
    /* show_three_bf("fact1",btempx,"fact2",btmp1,"prod ",btmp2,70); */
    div_bf(btmp2a,btmp2,bfftemp);
    /* show_three_bf("num  ",btmp2,"denom",bfftemp,"quot ",btmp2a,70); */
    sub_bf(btmp3,bfsx3rd,bfsxmin);
    mult_bf(btmp4,btempy,btmp3);
    div_bf(btmp4a,btmp4,bfftemp);
    add_bf(bfnewx,bfsxmin,btmp2a);
    add_a_bf(bfnewx,btmp4a);

    /* *newy = symax + tempy*(sy3rd-symax)/ftemp + tempx*(symin-sy3rd)/ftemp; */
    sub_bf(btmp1,bfsy3rd,bfsymax);
    mult_bf(btmp2,btempy,btmp1);
    div_bf(btmp2a,btmp2,bfftemp);
    sub_bf(btmp3,bfsymin,bfsy3rd);
    mult_bf(btmp4,btempx,btmp3);
    div_bf(btmp4a,btmp4,bfftemp);
    add_bf(bfnewy,bfsymax,btmp2a);
    add_a_bf(bfnewy,btmp4a);
    restore_stack(saved);
}
コード例 #2
0
int
JuliabfFractal()
{
    // new.x = tmpsqrx - tmpsqry + parm.x;
    sub_a_bf(bftmpsqrx, bftmpsqry);
    add_bf(bfnew.x, bftmpsqrx, bfparm.x);

    // new.y = 2 * bfold.x * bfold.y + parm.y;
    mult_bf(bftmp, bfold.x, bfold.y); // ok to use unsafe here
    double_a_bf(bftmp);
    add_bf(bfnew.y, bftmp, bfparm.y);
    return bigfltbailout();
}
コード例 #3
0
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();
}
コード例 #4
0
int  bfMODbailout()
{
    long longmagnitude;

    square_bf(bftmpsqrx, bfnew.x);
    square_bf(bftmpsqry, bfnew.y);
    add_bf(bftmp, bftmpsqrx, bftmpsqry);

    longmagnitude = bftoint(bftmp);
    if (longmagnitude >= (long)rqlim)
        return 1;
    copy_bf(bfold.x, bfnew.x);
    copy_bf(bfold.y, bfnew.y);
    return 0;
}
コード例 #5
0
BFComplex *cplxmul_bf(BFComplex *t, BFComplex *x, BFComplex *y)
{
    bf_t tmp1;
    int saved;
    saved = save_stack();
    tmp1 = alloc_stack(rbflength+2);
    mult_bf(t->x, x->x, y->x);
    mult_bf(t->y, x->y, y->y);
    sub_bf(t->x, t->x, t->y);

    mult_bf(tmp1, x->x, y->y);
    mult_bf(t->y, x->y, y->x);
    add_bf(t->y, tmp1, t->y);
    restore_stack(saved);
    return (t);
}
コード例 #6
0
int  bfMANRbailout()
{
    long longtempmag;

    square_bf(bftmpsqrx, bfnew.x);
    square_bf(bftmpsqry, bfnew.y);
    add_bf(bftmp, bfnew.x, bfnew.y); // don't need abs since we square it next
    // note: in next two lines, bfold is just used as a temporary variable
    square_bf(bfold.x, bftmp);
    longtempmag = bftoint(bfold.x);
    if (longtempmag >= (long)rqlim)
        return 1;
    copy_bf(bfold.x, bfnew.x);
    copy_bf(bfold.y, bfnew.y);
    return (0);
}
コード例 #7
0
ファイル: ZOOM.C プロジェクト: mmammel/projects
/* 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);
}
コード例 #8
0
int  bfMANHbailout()
{
    long longtempmag;

    square_bf(bftmpsqrx, bfnew.x);
    square_bf(bftmpsqry, bfnew.y);
    // note: in next five lines, bfold is just used as a temporary variable
    abs_bf(bfold.x, bfnew.x);
    abs_bf(bfold.y, bfnew.y);
    add_bf(bftmp, bfold.x, bfold.y);
    square_bf(bfold.x, bftmp);
    longtempmag = bftoint(bfold.x);
    if (longtempmag >= (long)rqlim)
        return 1;
    copy_bf(bfold.x, bfnew.x);
    copy_bf(bfold.y, bfnew.y);
    return (0);
}