/* biginit.c - C routines for bignumbers */ /* Note: This is NOT the biginit.c file that come the standard BigNum library, but is a customized version specific to Fractint. The biggest difference is in the allocations of memory for the big numbers. */ #include #ifdef __APPLE__ #include #elif !defined(BIG_ANSI_C) #include #endif /* see Fractint.c for a description of the "include" hierarchy */ #include "port.h" #include "prototyp.h" #include "fractype.h" /* appears to me that avoiding the start of extraseg is unnecessary. If correct, later we can eliminate ENDVID here. */ #ifdef ENDVID #undef ENDVID #endif #define ENDVID 0 /* globals */ #ifdef BIG_BASED _segment bignum_seg; #endif int bnstep, bnlength, intlength, rlength, padding, shiftfactor, decimals; int bflength, rbflength, bfdecimals; /* used internally by bignum.c routines */ static bn_t bnroot=BIG_NULL; static bn_t stack_ptr; /* memory allocator base after global variables */ bn_t bntmp1, bntmp2, bntmp3, bntmp4, bntmp5, bntmp6; /* rlength */ bn_t bntmpcpy1, bntmpcpy2; /* bnlength */ /* used by other routines */ bn_t bnxmin, bnxmax, bnymin, bnymax, bnx3rd, bny3rd; /* bnlength */ bn_t bnxdel, bnydel, bnxdel2, bnydel2, bnclosenuff; /* bnlength */ bn_t bntmpsqrx, bntmpsqry, bntmp; /* rlength */ _BNCMPLX /* bnold, bnnew, */ bnparm, bnsaved; /* bnlength */ _BNCMPLX bnold, bnnew; /* rlength */ bn_t bn_pi; /* TAKES NO SPACE */ bf_t bftmp1, bftmp2, bftmp3, bftmp4, bftmp5, bftmp6; /* rbflength+2 */ bf_t bftmpcpy1, bftmpcpy2; /* rbflength+2 */ bf_t bfxdel, bfydel, bfxdel2, bfydel2, bfclosenuff; /* rbflength+2 */ bf_t bftmpsqrx, bftmpsqry; /* rbflength+2 */ _BFCMPLX /* bfold, bfnew, */ bfparm, bfsaved; /* bflength+2 */ _BFCMPLX bfold, bfnew; /* rbflength+2 */ bf_t bf_pi; /* TAKES NO SPACE */ bf_t big_pi; /* bflength+2 */ /* for testing only */ /* used by other routines */ bf_t bfxmin, bfxmax, bfymin, bfymax, bfx3rd, bfy3rd; /* bflength+2 */ bf_t bfsxmin, bfsxmax, bfsymin, bfsymax, bfsx3rd, bfsy3rd;/* bflength+2 */ bf_t bfparms[10]; /* (bflength+2)*10 */ bf_t bftmp; bf_t bf10tmp; /* dec+4 */ #define LOG10_256 2.4082399653118 #define LOG_256 5.5451774444795 static int save_bf_vars(void); static int restore_bf_vars(void); /*********************************************************************/ /* given bnlength, calc_lengths will calculate all the other lengths */ void calc_lengths(void) { #if 0 #ifdef USE_BIGNUM_C_CODE bnstep = 2; #else /* use 80x86 asm code */ if (cpu >= 386) bnstep = 4; else /* cpu <= 286 */ bnstep = 2; #endif #else bnstep = 4; /* use 4 in all cases */ #endif if (bnlength % bnstep != 0) bnlength = (bnlength / bnstep + 1) * bnstep; if (bnlength == bnstep) padding = bnlength; else padding = 2*bnstep; rlength = bnlength + padding; /* This shiftfactor assumes non-full multiplications will be performed.*/ /* Change to bnlength-intlength for full multiplications. */ shiftfactor = padding - intlength; bflength = bnlength+bnstep; /* one extra step for added precision */ rbflength = bflength + padding; bfdecimals = (int)((bflength-2)*LOG10_256); } /************************************************************************/ /* intended only to be called from init_bf_dec() or init_bf_length(). */ /* initialize bignumber global variables */ long maxptr = 0; long startstack = 0; long maxstack = 0; int bf_save_len = 0; /* ??? for some strange reason, msc 7.0 hangs here without this pragma. ??? */ #ifndef XFRACT #pragma optimize( "", off ) #endif static void init_bf_2(void) { int i; long ptr; save_bf_vars(); /* copy corners values for conversion */ calc_lengths(); /* allocate all the memory at once within the same segment (DOS) */ #if defined(BIG_FAR) || defined(BIG_ANSI_C) bnroot = (bf_t)MK_FP(extraseg,ENDVID); /* ENDVID is to avoid videotable */ #else /* BASED or NEAR */ bnroot = (bf_t)ENDVID; /* ENDVID is to avoid videotable */ #endif #ifdef BIG_BASED bignum_seg = (_segment)extraseg; #endif /* at present time one call would suffice, but this logic allows multiple kinds of alternate math eg long double */ if((i = find_alternate_math(fractype, BIGNUM)) > -1) bf_math = alternatemath[i].math; else if((i = find_alternate_math(fractype, BIGFLT)) > -1) bf_math = alternatemath[i].math; else bf_math = 1; /* maybe called from cmdfiles.c and fractype not set */ floatflag=1; /* Now split up the memory among the pointers */ /* internal pointers */ ptr = 0; bntmp1 = bnroot+ptr; ptr += rlength; bntmp2 = bnroot+ptr; ptr += rlength; bntmp3 = bnroot+ptr; ptr += rlength; bntmp4 = bnroot+ptr; ptr += rlength; bntmp5 = bnroot+ptr; ptr += rlength; bntmp6 = bnroot+ptr; ptr += rlength; bftmp1 = bnroot+ptr; ptr += rbflength+2; bftmp2 = bnroot+ptr; ptr += rbflength+2; bftmp3 = bnroot+ptr; ptr += rbflength+2; bftmp4 = bnroot+ptr; ptr += rbflength+2; bftmp5 = bnroot+ptr; ptr += rbflength+2; bftmp6 = bnroot+ptr; ptr += rbflength+2; bftmpcpy1 = bnroot+ptr; ptr += (rbflength+2)*2; bftmpcpy2 = bnroot+ptr; ptr += (rbflength+2)*2; bntmpcpy1 = bnroot+ptr; ptr += (rlength*2); bntmpcpy2 = bnroot+ptr; ptr += (rlength*2); if (bf_math == BIGNUM) { bnxmin = bnroot+ptr; ptr += bnlength; bnxmax = bnroot+ptr; ptr += bnlength; bnymin = bnroot+ptr; ptr += bnlength; bnymax = bnroot+ptr; ptr += bnlength; bnx3rd = bnroot+ptr; ptr += bnlength; bny3rd = bnroot+ptr; ptr += bnlength; bnxdel = bnroot+ptr; ptr += bnlength; bnydel = bnroot+ptr; ptr += bnlength; bnxdel2 = bnroot+ptr; ptr += bnlength; bnydel2 = bnroot+ptr; ptr += bnlength; bnold.x = bnroot+ptr; ptr += rlength; bnold.y = bnroot+ptr; ptr += rlength; bnnew.x = bnroot+ptr; ptr += rlength; bnnew.y = bnroot+ptr; ptr += rlength; bnsaved.x = bnroot+ptr; ptr += bnlength; bnsaved.y = bnroot+ptr; ptr += bnlength; bnclosenuff= bnroot+ptr; ptr += bnlength; bnparm.x = bnroot+ptr; ptr += bnlength; bnparm.y = bnroot+ptr; ptr += bnlength; bntmpsqrx = bnroot+ptr; ptr += rlength; bntmpsqry = bnroot+ptr; ptr += rlength; bntmp = bnroot+ptr; ptr += rlength; } if (bf_math == BIGFLT) { bfxdel = bnroot+ptr; ptr += bflength+2; bfydel = bnroot+ptr; ptr += bflength+2; bfxdel2 = bnroot+ptr; ptr += bflength+2; bfydel2 = bnroot+ptr; ptr += bflength+2; bfold.x = bnroot+ptr; ptr += rbflength+2; bfold.y = bnroot+ptr; ptr += rbflength+2; bfnew.x = bnroot+ptr; ptr += rbflength+2; bfnew.y = bnroot+ptr; ptr += rbflength+2; bfsaved.x = bnroot+ptr; ptr += bflength+2; bfsaved.y = bnroot+ptr; ptr += bflength+2; bfclosenuff= bnroot+ptr; ptr += bflength+2; bfparm.x = bnroot+ptr; ptr += bflength+2; bfparm.y = bnroot+ptr; ptr += bflength+2; bftmpsqrx = bnroot+ptr; ptr += rbflength+2; bftmpsqry = bnroot+ptr; ptr += rbflength+2; bftmp = bnroot+ptr; ptr += rbflength+2; } big_pi = bnroot+ptr; ptr += bflength+2; /* needed by both BIGNUM & BIGFLT */ bf10tmp = bnroot+ptr; ptr += bfdecimals+4; /* ptr needs to be 16-bit aligned on some systems */ ptr = (ptr+1)&~1; stack_ptr = bnroot + ptr; startstack = ptr; /* max stack offset from bnroot */ maxstack = (long)0x10000l-(bflength+2)*22-ENDVID; /* sanity check */ /* leave room for NUMVARS variables allocated from stack */ /* also leave room for the safe area at top of segment */ if(ptr + NUMVARS*(bflength+2) > maxstack) { char msg[80]; char nmsg[80]; static FCODE fmsg[] = {"Requested precision of %d too high, aborting"}; far_strcpy(nmsg,fmsg); sprintf(msg,nmsg,decimals); stopmsg(0,msg); goodbye(); } /* room for 6 corners + 6 save corners + 10 params at top of extraseg */ /* this area is safe - use for variables that are used outside fractal*/ /* generation - e.g. zoom box variables */ ptr = maxstack; bfxmin = bnroot+ptr; ptr += bflength+2; bfxmax = bnroot+ptr; ptr += bflength+2; bfymin = bnroot+ptr; ptr += bflength+2; bfymax = bnroot+ptr; ptr += bflength+2; bfx3rd = bnroot+ptr; ptr += bflength+2; bfy3rd = bnroot+ptr; ptr += bflength+2; for(i=0;i<10;i++) { bfparms[i] = bnroot+ptr; ptr += bflength+2; } bfsxmin = bnroot+ptr; ptr += bflength+2; bfsxmax = bnroot+ptr; ptr += bflength+2; bfsymin = bnroot+ptr; ptr += bflength+2; bfsymax = bnroot+ptr; ptr += bflength+2; bfsx3rd = bnroot+ptr; ptr += bflength+2; bfsy3rd = bnroot+ptr; ptr += bflength+2; /* end safe vars */ /* good citizens initialize variables */ if(bf_save_len) /* leave save area */ far_memset(bnroot+(bf_save_len+2)*22,0,(unsigned)(startstack-(bf_save_len+2)*22)); else /* first time through - nothing saved */ { /* high variables */ far_memset(bnroot+maxstack,0,(bflength+2)*22); /* low variables */ far_memset(bnroot,0,(unsigned)startstack); } restore_bf_vars(); /* Initialize the value of pi. Needed for trig functions. */ /* init_big_pi(); */ /* call to init_big_pi() has been moved to fractal setup routine */ /* so as to use only when necessary. */ } /**********************************************************/ /* save current corners and parameters to start of bnroot */ /* to preserve values across calls to init_bf() */ static int save_bf_vars(void) { int ret; unsigned int mem; if(bnroot != BIG_NULL) { mem = (bflength+2)*22; /* 6 corners + 6 save corners + 10 params */ bf_save_len = bflength; far_memcpy(bnroot,bfxmin,mem); /* scrub old high area */ far_memset(bfxmin,0,mem); ret = 0; } else { bf_save_len = 0; ret = -1; } return(ret); } /************************************************************************/ /* copy current corners and parameters from save location */ static int restore_bf_vars(void) { bf_t ptr; int i; if(bf_save_len == 0) return(-1); ptr = bnroot; convert_bf(bfxmin,ptr,bflength,bf_save_len); ptr += bf_save_len+2; convert_bf(bfxmax,ptr,bflength,bf_save_len); ptr += bf_save_len+2; convert_bf(bfymin,ptr,bflength,bf_save_len); ptr += bf_save_len+2; convert_bf(bfymax,ptr,bflength,bf_save_len); ptr += bf_save_len+2; convert_bf(bfx3rd,ptr,bflength,bf_save_len); ptr += bf_save_len+2; convert_bf(bfy3rd,ptr,bflength,bf_save_len); ptr += bf_save_len+2; for(i=0;i<10;i++) { convert_bf(bfparms[i],ptr,bflength,bf_save_len); ptr += bf_save_len+2; } convert_bf(bfsxmin,ptr,bflength,bf_save_len); ptr += bf_save_len+2; convert_bf(bfsxmax,ptr,bflength,bf_save_len); ptr += bf_save_len+2; convert_bf(bfsymin,ptr,bflength,bf_save_len); ptr += bf_save_len+2; convert_bf(bfsymax,ptr,bflength,bf_save_len); ptr += bf_save_len+2; convert_bf(bfsx3rd,ptr,bflength,bf_save_len); ptr += bf_save_len+2; convert_bf(bfsy3rd,ptr,bflength,bf_save_len); ptr += bf_save_len+2; /* scrub save area */ far_memset(bnroot,0,(bf_save_len+2)*22); return(0); } /*******************************************/ /* free corners and parameters save memory */ void free_bf_vars() { bf_save_len = bf_math = 0; bnstep=bnlength=intlength=rlength=padding=shiftfactor=decimals=0; bflength=rbflength=bfdecimals=0; maxptr = startstack = maxstack = 0; } #ifndef XFRACT #pragma optimize( "", on ) #endif /************************************************************************/ /* Memory allocator routines start here. */ /************************************************************************/ /* Allocates a bn_t variable on stack */ bn_t alloc_stack(size_t size) { long stack_addr; if(bf_math == 0) { static FCODE msg[] = {"alloc_stack called with bf_math==0"}; stopmsg(0,msg); return(0); } stack_addr = (long)(stack_ptr-bnroot)+size; /* +ENDVID, part of bnroot */ if(stack_addr > maxstack) { static FCODE msg[] = {"Aborting, Out of Bignum Stack Space"}; stopmsg(0,msg); goodbye(); } /* keep track of max ptr */ if(stack_addr > maxptr) maxptr = stack_addr; stack_ptr += size; /* increment stack pointer */ return(stack_ptr - size); } /************************************************************************/ /* Returns stack pointer offset so it can be saved. */ int save_stack(void) { return(stack_ptr - bnroot); } /************************************************************************/ /* Restores stack pointer, effectively freeing local variables */ /* allocated since save_stack() */ void restore_stack(int old_offset) { stack_ptr = bnroot+old_offset; } /************************************************************************/ /* Memory allocator routines end here. */ /************************************************************************/ /************************************************************************/ /* initialize bignumber global variables */ /* dec = decimal places after decimal point */ /* intl = bytes for integer part (1, 2, or 4) */ void init_bf_dec(int dec) { if(bfdigits) decimals=bfdigits; /* blindly force */ else decimals = dec; if(bailout > 10) /* arbitrary value */ /* using 2 doesn't gain much and requires another test */ intlength = 4; else if (fractype == FPMANDELZPOWER || fractype == FPJULIAZPOWER || fractype == DIVIDEBROT5) intlength = 4; /* 2 leaves artifacts in the center of the lakes */ /* the bailout tests need greater dynamic range */ else if(bailoutest == Real || bailoutest == Imag || bailoutest == And || bailoutest == Manr) intlength = 2; else intlength = 1; /* conservative estimate */ bnlength = intlength + (int)(decimals/LOG10_256) + 1; /* round up */ init_bf_2(); } /************************************************************************/ /* initialize bignumber global variables */ /* bnl = bignumber length */ /* intl = bytes for integer part (1, 2, or 4) */ void init_bf_length(int bnl) { bnlength = bnl; if(bailout > 10) /* arbitrary value */ /* using 2 doesn't gain much and requires another test */ intlength = 4; else if (fractype == FPMANDELZPOWER || fractype == FPJULIAZPOWER || fractype == DIVIDEBROT5) intlength = 4; /* 2 leaves artifacts in the center of the lakes */ /* the bailout tests need greater dynamic range */ else if(bailoutest == Real || bailoutest == Imag || bailoutest == And || bailoutest == Manr) intlength = 2; else intlength = 1; /* conservative estimate */ decimals = (int)((bnlength-intlength)*LOG10_256); init_bf_2(); } void init_big_pi(void) { /* What, don't you recognize the first 700 digits of pi, */ /* in base 256, in reverse order? */ int length, pi_offset; static BFCODE pi_table[] = { 0x44, 0xD5, 0xDB, 0x69, 0x17, 0xDF, 0x2E, 0x56, 0x87, 0x1A, 0xA0, 0x8C, 0x6F, 0xCA, 0xBB, 0x57, 0x5C, 0x9E, 0x82, 0xDF, 0x00, 0x3E, 0x48, 0x7B, 0x31, 0x53, 0x60, 0x87, 0x23, 0xFD, 0xFA, 0xB5, 0x3D, 0x32, 0xAB, 0x52, 0x05, 0xAD, 0xC8, 0x1E, 0x50, 0x2F, 0x15, 0x6B, 0x61, 0xFD, 0xDF, 0x16, 0x75, 0x3C, 0xF8, 0x22, 0x32, 0xDB, 0xF8, 0xE9, 0xA5, 0x8E, 0xCC, 0xA3, 0x1F, 0xFB, 0xFE, 0x25, 0x9F, 0x67, 0x79, 0x72, 0x2C, 0x40, 0xC6, 0x00, 0xA1, 0xD6, 0x0A, 0x32, 0x60, 0x1A, 0xBD, 0xC0, 0x79, 0x55, 0xDB, 0xFB, 0xD3, 0xB9, 0x39, 0x5F, 0x0B, 0xD2, 0x0F, 0x74, 0xC8, 0x45, 0x57, 0xA8, 0xCB, 0xC0, 0xB3, 0x4B, 0x2E, 0x19, 0x07, 0x28, 0x0F, 0x66, 0xFD, 0x4A, 0x33, 0xDE, 0x04, 0xD0, 0xE3, 0xBE, 0x09, 0xBD, 0x5E, 0xAF, 0x44, 0x45, 0x81, 0xCC, 0x2C, 0x95, 0x30, 0x9B, 0x1F, 0x51, 0xFC, 0x6D, 0x6F, 0xEC, 0x52, 0x3B, 0xEB, 0xB2, 0x39, 0x13, 0xB5, 0x53, 0x6C, 0x3E, 0xAF, 0x6F, 0xFB, 0x68, 0x63, 0x24, 0x6A, 0x19, 0xC2, 0x9E, 0x5C, 0x5E, 0xC4, 0x60, 0x9F, 0x40, 0xB6, 0x4F, 0xA9, 0xC1, 0xBA, 0x06, 0xC0, 0x04, 0xBD, 0xE0, 0x6C, 0x97, 0x3B, 0x4C, 0x79, 0xB6, 0x1A, 0x50, 0xFE, 0xE3, 0xF7, 0xDE, 0xE8, 0xF6, 0xD8, 0x79, 0xD4, 0x25, 0x7B, 0x1B, 0x99, 0x80, 0xC9, 0x72, 0x53, 0x07, 0x9B, 0xC0, 0xF1, 0x49, 0xD3, 0xEA, 0x0F, 0xDB, 0x48, 0x12, 0x0A, 0xD0, 0x24, 0xD7, 0xD0, 0x37, 0x3D, 0x02, 0x9B, 0x42, 0x72, 0xDF, 0xFE, 0x1B, 0x06, 0x77, 0x3F, 0x36, 0x62, 0xAA, 0xD3, 0x4E, 0xA6, 0x6A, 0xC1, 0x56, 0x9F, 0x44, 0x1A, 0x40, 0x73, 0x20, 0xC1, 0x85, 0xD8, 0x75, 0x6F, 0xE0, 0xBE, 0x5E, 0x8B, 0x3B, 0xC3, 0xA5, 0x84, 0x7D, 0xB4, 0x9F, 0x6F, 0x45, 0x19, 0x86, 0xEE, 0x8C, 0x88, 0x0E, 0x43, 0x82, 0x3E, 0x59, 0xCA, 0x66, 0x76, 0x01, 0xAF, 0x39, 0x1D, 0x65, 0xF1, 0xA1, 0x98, 0x2A, 0xFB, 0x7E, 0x50, 0xF0, 0x3B, 0xBA, 0xE4, 0x3B, 0x7A, 0x13, 0x6C, 0x0B, 0xEF, 0x6E, 0xA3, 0x33, 0x51, 0xAB, 0x28, 0xA7, 0x0F, 0x96, 0x68, 0x2F, 0x54, 0xD8, 0xD2, 0xA0, 0x51, 0x6A, 0xF0, 0x88, 0xD3, 0xAB, 0x61, 0x9C, 0x0C, 0x67, 0x9A, 0x6C, 0xE9, 0xF6, 0x42, 0x68, 0xC6, 0x21, 0x5E, 0x9B, 0x1F, 0x9E, 0x4A, 0xF0, 0xC8, 0x69, 0x04, 0x20, 0x84, 0xA4, 0x82, 0x44, 0x0B, 0x2E, 0x39, 0x42, 0xF4, 0x83, 0xF3, 0x6F, 0x6D, 0x0F, 0xC5, 0xAC, 0x96, 0xD3, 0x81, 0x3E, 0x89, 0x23, 0x88, 0x1B, 0x65, 0xEB, 0x02, 0x23, 0x26, 0xDC, 0xB1, 0x75, 0x85, 0xE9, 0x5D, 0x5D, 0x84, 0xEF, 0x32, 0x80, 0xEC, 0x5D, 0x60, 0xAC, 0x7C, 0x48, 0x91, 0xA9, 0x21, 0xFB, 0xCC, 0x09, 0xD8, 0x61, 0x93, 0x21, 0x28, 0x66, 0x1B, 0xE8, 0xBF, 0xC4, 0xAF, 0xB9, 0x4B, 0x6B, 0x98, 0x48, 0x8F, 0x3B, 0x77, 0x86, 0x95, 0x28, 0x81, 0x53, 0x32, 0x7A, 0x5C, 0xCF, 0x24, 0x6C, 0x33, 0xBA, 0xD6, 0xAF, 0x1E, 0x93, 0x87, 0x9B, 0x16, 0x3E, 0x5C, 0xCE, 0xF6, 0x31, 0x18, 0x74, 0x5D, 0xC5, 0xA9, 0x2B, 0x2A, 0xBC, 0x6F, 0x63, 0x11, 0x14, 0xEE, 0xB3, 0x93, 0xE9, 0x72, 0x7C, 0xAF, 0x86, 0x54, 0xA1, 0xCE, 0xE8, 0x41, 0x11, 0x34, 0x5C, 0xCC, 0xB4, 0xB6, 0x10, 0xAB, 0x2A, 0x6A, 0x39, 0xCA, 0x55, 0x40, 0x14, 0xE8, 0x63, 0x62, 0x98, 0x48, 0x57, 0x94, 0xAB, 0x55, 0xAA, 0xF3, 0x25, 0x55, 0xE6, 0x60, 0x5C, 0x60, 0x55, 0xDA, 0x2F, 0xAF, 0x78, 0x27, 0x4B, 0x31, 0xBD, 0xC1, 0x77, 0x15, 0xD7, 0x3E, 0x8A, 0x1E, 0xB0, 0x8B, 0x0E, 0x9E, 0x6C, 0x0E, 0x18, 0x3A, 0x60, 0xB0, 0xDC, 0x79, 0x8E, 0xEF, 0x38, 0xDB, 0xB8, 0x18, 0x79, 0x41, 0xCA, 0xF0, 0x85, 0x60, 0x28, 0x23, 0xB0, 0xD1, 0xC5, 0x13, 0x60, 0xF2, 0x2A, 0x39, 0xD5, 0x30, 0x9C, 0xB5, 0x59, 0x5A, 0xC2, 0x1D, 0xA4, 0x54, 0x7B, 0xEE, 0x4A, 0x15, 0x82, 0x58, 0xCD, 0x8B, 0x71, 0x58, 0xB6, 0x8E, 0x72, 0x8F, 0x74, 0x95, 0x0D, 0x7E, 0x3D, 0x93, 0xF4, 0xA3, 0xFE, 0x58, 0xA4, 0x69, 0x4E, 0x57, 0x71, 0xD8, 0x20, 0x69, 0x63, 0x16, 0xFC, 0x8E, 0x85, 0xE2, 0xF2, 0x01, 0x08, 0xF7, 0x6C, 0x91, 0xB3, 0x47, 0x99, 0xA1, 0x24, 0x99, 0x7F, 0x2C, 0xF1, 0x45, 0x90, 0x7C, 0xBA, 0x96, 0x7E, 0x26, 0x6A, 0xED, 0xAF, 0xE1, 0xB8, 0xB7, 0xDF, 0x1A, 0xD0, 0xDB, 0x72, 0xFD, 0x2F, 0xAC, 0xB5, 0xDF, 0x98, 0xA6, 0x0B, 0x31, 0xD1, 0x1B, 0xFB, 0x79, 0x89, 0xD9, 0xD5, 0x16, 0x92, 0x17, 0x09, 0x47, 0xB5, 0xB5, 0xD5, 0x84, 0x3F, 0xDD, 0x50, 0x7C, 0xC9, 0xB7, 0x29, 0xAC, 0xC0, 0x6C, 0x0C, 0xE9, 0x34, 0xCF, 0x66, 0x54, 0xBE, 0x77, 0x13, 0xD0, 0x38, 0xE6, 0x21, 0x28, 0x45, 0x89, 0x6C, 0x4E, 0xEC, 0x98, 0xFA, 0x2E, 0x08, 0xD0, 0x31, 0x9F, 0x29, 0x22, 0x38, 0x09, 0xA4, 0x44, 0x73, 0x70, 0x03, 0x2E, 0x8A, 0x19, 0x13, 0xD3, 0x08, 0xA3, 0x85, 0x88, 0x6A, 0x3F, 0x24, /* . */ 0x03, 0x00, 0x00, 0x00 /* <- up to intlength 4 -> */ /* or bf_t int length of 2 + 2 byte exp */ }; length = bflength+2; /* 2 byte exp */ pi_offset = sizeof pi_table - length; _fmemcpy(big_pi, pi_table + pi_offset, length); /* notice that bf_pi and bn_pi can share the same memory space */ bf_pi = big_pi; bn_pi = big_pi + (bflength-2) - (bnlength-intlength); return; }