//finally 4560 ticks

extern int zlongdiv();

unsigned int zugh1lb,zugh1rb,zugh1rbrb,zugh1rblb,zugh1lbrb;
unsigned zproduct,zeax,zebx,zedx,zedxeax;

unsigned int y2lb,y2rb,zcarrylb,zcarryrb,rightright,rightleft,leftright;
unsigned int leftleft,lbofrightleft;
unsigned int rbofrightleft,lbofleftright,rbofleftright;
unsigned int slot1,slot3;
unsigned int center,lbofcenter,rbofcenter;
unsigned int zdivisorlb,zdivisorrb;
unsigned int zlint2lb,zlint2rb;

unsigned int znumerator1,znumerator2;
unsigned int znumerator1lb,znumerator1rb;  
unsigned int znumerator2lb,znumerator2rb; //numerator=znumerator1,znumerator2
unsigned int zdenominator=0xffffffff;     //denominator zdenominator
unsigned int zdenominatorlb,zdenominatorrb;
unsigned int zborrowfromleft;

unsigned int crossl,crossr,rbofcrossr,lbofcrossr;
unsigned int reml,remr,lbofremr;
unsigned int remainderl,remainderr,lbofremainderr;
unsigned int zquotient1,zquotient2;
unsigned int zquotient,zremainder;

unsigned int zdlb,zdrb;

void modex()
 {int zecx;
  for (zecx=1;zecx<prec+2;zecx++) z[zecx-1]=z[zecx];
 }

void modey()             //appears unused.
 {int zecx;
  for (zecx=2;zecx<prec+2;zecx++) z[zecx-1]=z[zecx];
 }

void moxy1()
  {zugh1lb=zugh1>>16;     //left bits of zugh1=leftpart
    zugh1rb=zugh1<<16;
    zugh1rb=zugh1rb>>16;   //right bits of zugh1=rightpart
    zugh1rb=10*zugh1rb+zcarryprnt;   //10*rightpart+zcarryprnt;
    zugh1rbrb=zugh1rb<<16;           
    zugh1rbrb=zugh1rbrb>>16;  //right bits of that=rightbits
    zugh1rblb=zugh1rb-zugh1rbrb; 
    zugh1rblb=zugh1rblb>>16;     //left bits of that=leftbits
    zugh1lb=10*zugh1lb+zugh1rblb;  //10*leftpart+leftbits=left
    zcarryprnt=zugh1lb>>16;   //pick off the carry
    zugh1=zugh1lb<<16; //should have been ffff<<16
    zugh1=zugh1+zugh1rbrb;
}


void mode4()
 {int zecx,zedx;
  z[1]=x[1]+y[1];
  if (z[1]<x[1]) zedx=1; else zedx=0;
  for (zecx=2;zecx<prec+2;zecx++)
    {z[zecx]=x[zecx]+zedx;
     if (z[zecx]<zedx) zedx=1; else zedx=0;
     z[zecx]=z[zecx]+y[zecx];
     if (z[zecx]<y[zecx]) zedx=zedx+1;  //hmm
    }
 }         // x is in edi,y is in esi, z is in ebp



//through this point, up to 1455 ticks and running ok.
//goes up to 2505 when mode6() redone.
//Did move some << and >> outside of loops they were unnecessarily in.
//But this is what I expected.


void mode6()                //zlint28 set in genprecd.h2
  {int zebx,zecx;

   zlint2lb=zlint2>>16;
   zlint2rb=zlint2<<16;
   zlint2rb=zlint2rb>>16;
   zugh1=zlint28;
   zugh1lb=zugh1>>16;
   zugh1rb=zugh1<<16;
   zugh1rb=zugh1rb>>16;
   leftleft=zugh1lb*zlint2lb;
   leftright=zugh1lb*zlint2rb;
   rightleft=zugh1rb*zlint2lb;
   rightright=zugh1rb*zlint2rb;
   slot1=rightright+zcarry;
   if (slot1<rightright) slot3=0x00000001; else slot3=0;
   center=leftright+rightleft;
   if (center<leftright) slot3=slot3+0x00010000;
   lbofcenter=center>>16;   //break them into 16 bit lb and rb
   slot3=slot3+lbofcenter+leftleft;
   rbofcenter=center<<16;
   slot1=slot1+rbofcenter;
   if (slot1<rbofcenter) slot3=slot3+0x00000001;
   zcarry=slot3;  //only after zcarry in initialization
   zebx=0;
   for (zecx=precp1mzkisky+1;zecx<prec+1;zecx++)
     {zebx=zebx++;
      zugh1=z11[zecx];
      zugh1lb=zugh1>>16;
      zugh1rb=zugh1<<16;
      zugh1rb=zugh1rb>>16;
      leftleft=zugh1lb*zlint2lb;
      leftright=zugh1lb*zlint2rb;
      rightleft=zugh1rb*zlint2lb;
      rightright=zugh1rb*zlint2rb;
      slot1=rightright+zcarry;
      if (slot1<rightright) slot3=0x00000001; else slot3=0;
      center=leftright+rightleft;
      if (center<leftright) slot3=slot3+0x00010000;
      lbofcenter=center>>16;   //break them into 16 bit lb and rb
      slot3=slot3+lbofcenter+leftleft;
      rbofcenter=center<<16;
      slot1=slot1+rbofcenter;
      if (slot1<rbofcenter) slot3=slot3+0x00000001;
      x[zebx]=slot1;
      zcarry=slot3;
     }
  }

//at 2500 after mode66 redone

void mode66()
  {int zecx,zebx;
   zlint2lb=zlint2>>16;       //zlint2 has been set
   zlint2rb=zlint2<<16;
   zlint2rb=zlint2rb>>16;
   zebx=0;                    //zcarry has been set
   for (zecx=precp1mzkisky+1;zecx<prec+1;zecx++)
     {zebx=zebx++;
      zugh1=z11[zecx];
      zugh1lb=zugh1>>16;
      zugh1rb=zugh1<<16;
      zugh1rb=zugh1rb>>16;
      leftleft=zugh1lb*zlint2lb;
      leftright=zugh1lb*zlint2rb;
      rightleft=zugh1rb*zlint2lb;
      rightright=zugh1rb*zlint2rb;
      slot1=rightright+zcarry;
      if (slot1<rightright) slot3=0x00000001; else slot3=0;
      center=leftright+rightleft;
      if (center<leftright) slot3=slot3+0x00010000;
      lbofcenter=center>>16;   //break them into 16 bit lb and rb
      slot3=slot3+lbofcenter+leftleft;
      rbofcenter=center<<16;
      slot1=slot1+rbofcenter;
      if (slot1<rbofcenter) slot3=slot3+0x00000001;
      x[zebx]=slot1;
      zcarry=slot3;
     }
   zebx++;
   x[zebx]=zcarry;
  }

void mode77()
 {int zecx;
  unsigned int zedx;
  zedx=0;
  for (zecx=1;zecx<maxinmode7+1;zecx++)
   {zugh1=x[zecx]+zedx;
    if (zugh1<zedx) zedx=1; else zedx=0;
    z[zecx]=zugh1+z[zecx];
    if (z[zecx]<zugh1) zedx=zedx+1;
   }
//  zecx++;  //don't need this outside of assembly
  z[zecx]=x[zecx]+zedx;
 }

//up to 2860 ticks now--but class takes 13800 ticks

void mode777()
 {int zecx;
  unsigned int zedx=0;
  for (zecx=1;zecx<prec+1;zecx++)
    {zugh1=x[zecx]+zedx;
     if (zugh1<zedx) zedx=1;else zedx=0;
     z[zecx]=zugh1+z[zecx];
     if (z[zecx]<zugh1) zedx=zedx+1;
    }
  z[zecx]=x[zecx]+zedx;
 }

//about same number of ticks

int mode8()
  {int zecx;
   znumerator1=0;                   //starts with 0000,z11[prec]
   for (zecx=prec;zecx>0;zecx--)    //           =znumerator1,znumerator2
     {znumerator2=z11[zecx];        //next look will be 
      zdenominator=zkk2;
      zlongdiv();
      z[zecx]=zquotient;            //quotient into z
      znumerator1=zremainder;      //remainder is carried
     }
   return 0;
  }

//mode1(*(zden+precmzknpr));   //cross multiplies zd*zden to get zL
//moved to genprecd.h2
//4135 ticks after mode1 (main divider routine) redone

//up to 3065 ticks after mode2 replaced

int mode22()
  {znumerator1=zftemp1;
   znumerator2=0;
   zdenominator=za;
   zlongdiv();
   zugh1=zremainder;
   return 0;
  }

int mode3()
 {int zecx,zedx,zebx;
  zedx=0;
  for (zecx=0;zecx<zknp2;zecx++)
    {zugh1=zden[precp1mzkn+zecx]+zrem[zecx];
     if (zugh1<zrem[zecx]) zebx=1; else zebx=0;
     zugh1=zugh1+zedx;
     if (zugh1<zedx) zebx=zebx+1;
     zrem[zecx]=zugh1;
     zedx=zebx;
    }
  return 0;
 }

//up to 4245 ticks

int mode33()
  {int zebx,zecx,zedx;
   zebx=0;
   for (zecx=0;zecx<zknp2;zecx++)
     {zugh1=zhex99-zL[zecx];
      zrem[zecx]=zugh1+zrem[zecx];
      if (zrem[zecx]<zugh1) zedx=1; else zedx=0;
      zrem[zecx]=zrem[zecx]+zebx;
      if (zrem[zecx]<zebx) zedx=zedx+1;
      zebx=zedx;
     }
   return 0;
  }

/*
#define mode33()                                                     \
__asm__ (                                                           \
   "movl $0x0,%%ebx\n\t"     0 to ebx                                          \
   "movl $0x0,%%ecx\n\t"     0 to ecs                                          \
   "pushl %%ebp\n\t"  \
   "movl _zLorig,%%ebp\n\t"  \    zL is ebp
   "pushl %%esi\n\t"  \
   "movl _zhex99,%%esi\n\t"  \    zhex99 is esi
   "pushl %%edi\n\t"  \
   "movl _zremorig,%%edi\n\t"  \   zrem is edi
   "9:\n\t"                                         \
   "movl %%esi,%%eax\n\t"            zhex99                            \
   "subl (%%ebp,%%ecx,4),%%eax\n\t"    -zL[ecx]                                  \
   "addl (%%edi,%%ecx,4),%%eax\n\t"    +zrem[ecx]                                 \
   "movl $0x0,%%edx\n\t"                                               \
   "adcl $0x0,%%edx\n\t"            //carry into edx                                   \
   "addl %%ebx,%%eax\n\t"           ebx+mainpart                                    \
   "adcl $0x0,%%edx\n\t"            //add carry to edx                                   \
   "movl %%eax,(%%edi,%%ecx,4)\n\t"  main part to zrem[ecx]                                  \
   "movl %%edx,%%ebx\n\t"           edx to ebx                                   \
   "incl %%ecx\n\t"                                               \
   "cmpl _zknp2,%%ecx\n\t"         up to zknp2-1                               \
   "jne 9b\n\t"                                                         \
   "popl %%edi\n\t"  \
   "popl %%esi\n\t"  \
   "popl %%ebp\n\t"  \
   :                                                                \
   :                                                                \
   : "eax", "ebx", "ecx", "edx", "memory", "cc");
*/

int zlongdiv()
 {int zshiftmarker;
  zshiftmarker=0;
zverytop:
  znumerator1lb=znumerator1>>16;
  znumerator1rb=znumerator1<<16;
  znumerator1rb=znumerator1rb>>16;
  znumerator2lb=znumerator2>>16;
  znumerator2rb=znumerator2<<16;
  znumerator2rb=znumerator2rb>>16;
  zdenominatorlb=zdenominator>>16;  //when zdenominator=6 this is?
  zdenominatorrb=zdenominator<<16;
  zdenominatorrb=zdenominatorrb>>16;
if (zdenominatorlb==0)   //shift whole thing over
  {zdenominator=zdenominatorrb<<16;         //no more reference to zdenominator
   znumerator1=znumerator1rb<<16;
   znumerator1=znumerator1+znumerator2lb;
   znumerator2=znumerator2<<16;
   zshiftmarker=1;                     //remember to shift remainder
   goto zverytop;
  }
  zquotient1=znumerator1/zdenominatorlb;    //first divide::should be 16 bits
  if (zquotient1>0xffff) zquotient1=0xffff;
tryanotherzquotient1:
  crossr=zquotient1*zdenominatorrb;   //first part of cross will be 32 bits
  lbofcrossr=crossr>>16;
  rbofcrossr=crossr<<16;      
  crossr=rbofcrossr>>16;      //crossr willnow be 16 bits
  crossl=zquotient1*zdenominatorlb+lbofcrossr;
crossr=crossr<<16;  //adds a 0 to crossr
if (crossr>znumerator2)   //borrow 1 0000 0000
  {  remr=crossr-znumerator2-1;    //this will be +
     remr=0xffffffff-remr;  //this will be +
     zborrowfromleft=1;
  }
else
  {remr=znumerator2-crossr;
   zborrowfromleft=0;
  }
if (crossl+zborrowfromleft>znumerator1)
  {zquotient1=zquotient1-1;
   goto tryanotherzquotient1;}

reml=znumerator1-zborrowfromleft-crossl;

if (reml>zdenominatorlb)
  {zquotient1=zquotient1+1;
   goto tryanotherzquotient1;
  }

lbofremr=remr>>16;        //shift rem left
reml=reml<<16;
reml=reml+lbofremr;
remr=remr<<16;
zquotient2=reml/zdenominatorlb;

tryanotherzquotient2:

if (zquotient2>0xffff) zquotient2=0xffff;
  crossr=zquotient2*zdenominatorrb;   //first part of cross will be 32 bits
  lbofcrossr=crossr>>16;
  rbofcrossr=crossr<<16;      
  crossr=rbofcrossr>>16;      //crossr willnow be 16 bits
crossl=zquotient2*zdenominatorlb+lbofcrossr;
crossr=crossr<<16;  //adds a 0 to crossr
if (crossr>remr)   //borrow 1 0000 0000      //ok
  {  remainderr=crossr-remr-1;    //this will be +
     remainderr=0xffffffff-remainderr;  //this will be +
     zborrowfromleft=1;
  }
else
  {remainderr=remr-crossr;
   zborrowfromleft=0;
  }
if (crossl+zborrowfromleft>reml)
  {zquotient2=zquotient2-1;
   goto tryanotherzquotient2;
  }

remainderl=reml-zborrowfromleft-crossl;

if (remainderl>zdenominatorlb)
  {zquotient2=zquotient2+1;
   goto tryanotherzquotient2;
  }

lbofremainderr=remainderr>>16;        //shift rem left
remainderl=remainderl<<16;
zremainder=remainderl+lbofremainderr;
//remainderl=remainderl+lbofremainderr;
//remainderr=remainderr<<16;     //does nothing, should be removed
zquotient1=zquotient1<<16;
zquotient=zquotient1+zquotient2;    //output is zquotient,zremainder
if (zshiftmarker==1) zremainder=zremainder>>16;
return 0;
}