extern int cprnt(unsigned int * *czp);                       
extern int cadd(unsigned int * *cz1,unsigned int * *cz2);
extern int csub(unsigned int * *cz1,unsigned int * *cz2);
extern int cmult(unsigned int * *cz3,unsigned int * *cz4);
extern int cmultpol(unsigned int * *cz3,unsigned int * *cz4);  
extern int cdiv(unsigned int * *cz5,unsigned int * *cz6);
extern int cdivpol(unsigned int * *cz5,unsigned int * *cz6);   
extern int cabs(unsigned int * *cz1);
extern int ccompconj(unsigned int * *cz1);
extern int ccompconjpol(unsigned int * *cz1);                 
extern int crecall(unsigned int * *cstr);                     
extern int cstore(unsigned int * *cstr);                      
extern int crectopol(unsigned int * *cz26);
extern int cpoltorec(unsigned int * *cz26);
extern int cchs(unsigned int * *cz1);                         
extern int cchspol(unsigned int * *cz1);                      
extern int cxtonthpol(unsigned int * *cz7,unsigned int cn);   
extern int zsinh(unsigned int *zzyy1);                       
extern int zcosh(unsigned int *zy2);                         
extern int csin(unsigned int * *cz8);                        
extern int ccos(unsigned int * *cz9);                        
extern int clogpolinrec(unsigned int * *cz8);                  
extern int cln(unsigned int * *cz9);
extern int cexp(unsigned int * *cz9);
extern int carctan(unsigned int * *cz0);                     
extern int cconinttomp(unsigned int cr1,unsigned int cr2);
extern int cconinttomp(int cr1,int cr2);            
extern int czero(unsigned int * *cz1);                         
extern int cmultint(unsigned int * *cz1,unsigned int cr1);
extern int cmultint(unsigned int * *cz1,int cr1);    
extern int cdivint(unsigned int * *cz1,unsigned int cr1);
extern int cdivint(unsigned int * *cz1,int cr1);   
extern int cmultcbyr(unsigned int * *cz1,unsigned int * crfac);  
extern int cdivcbyr(unsigned int * *cz1,unsigned int * crfac);
extern int caddcplusr(unsigned int * * cz1,unsigned int * crfac);
extern int cexpitheta(unsigned int *crfac);
extern int crecip(unsigned int * *cz1);                        
extern int csqrt(unsigned int * *cz0);
extern int csqrtpol(unsigned int * *cz9);                       
extern int cgenprecmem();

#include "cclasses.h"

int czi;                             
unsigned int cn;                     
unsigned int cr1,cr2;

complexscalar cz_ram;unsigned int * *cz=cz_ram.cs;      
complexscalar cz0_ram;unsigned int * *cz0=cz0_ram.cs;
complexscalar cz1_ram;unsigned int * *cz1=cz1_ram.cs;
complexscalar cz2_ram;unsigned int * *cz2=cz2_ram.cs;
complexscalar cz26_ram;unsigned int * *cz26=cz26_ram.cs;
complexscalar cz3_ram;unsigned int * *cz3=cz3_ram.cs;
complexscalar cz4_ram;unsigned int * *cz4=cz4_ram.cs;
complexscalar cz5_ram;unsigned int * *cz5=cz5_ram.cs;
complexscalar cz6_ram;unsigned int * *cz6=cz6_ram.cs;
complexscalar cz7_ram;unsigned int * *cz7=cz7_ram.cs;
complexscalar cz8_ram;unsigned int * *cz8=cz8_ram.cs;
complexscalar cz88_ram;unsigned int * *cz88=cz88_ram.cs;
complexscalar cz9_ram;unsigned int * *cz9=cz9_ram.cs;
complexscalar cz99_ram;unsigned int * *cz99=cz99_ram.cs;   
complexscalar czp_ram;unsigned int * *czp=czp_ram.cs;
complexscalar cstr_ram;unsigned int * *cstr=cstr_ram.cs;    
complexscalar ctempcomp_ram;unsigned int * *ctempcomp=ctempcomp_ram.cs;
complexscalar ctempcomp1_ram;unsigned int * *ctempcomp1=ctempcomp1_ram.cs;
complexscalar ctempcomp3_ram;unsigned int * *ctempcomp3=ctempcomp3_ram.cs;  
realscalar ctempreal_ram;unsigned int *ctempreal=ctempreal_ram.s;
realscalar ctempreal2_ram;unsigned int *ctempreal2=ctempreal2_ram.s;
realscalar ctempreal3_ram;unsigned int *ctempreal3=ctempreal3_ram.s;
realscalar ctempreal4_ram;unsigned int *ctempreal4=ctempreal4_ram.s;
realscalar ctempreal5_ram;unsigned int *ctempreal5=ctempreal5_ram.s;
realscalar ctempreal6_ram;unsigned int *ctempreal6=ctempreal6_ram.s;   
realscalar cpie_ram;unsigned int *cpie=cpie_ram.s;
realscalar zzyy1_ram;unsigned int *zzyy1=zzyy1_ram.s;
realscalar zy2_ram;unsigned int *zy2=zy2_ram.s;
realscalar crfac_ram;unsigned int *crfac=crfac_ram.s;
realscalar z3boy_ram;unsigned int *z3boy=z3boy_ram.s;   

int cigenprecmem=cgenprecmem();
int cgenprecmem()
 {int zi;
  unsigned int zj;
  zmem=zmem+sizeof(zi)+3*sizeof(zj);
  zmemforreal();
  zmem=zmem+11*zreturnmem;
  zmemforcomplex();
  zmem=zmem+19*zreturnmem;
  return 0;
 }

int cprnt(unsigned int * *czp)
  {zprnt(czp[0]);
   zprnt(czp[1]);
   return 0;
  }

int cstore(unsigned int * *cstr)
  {zrecall(cz[0]);
   zstore(cstr[0]);
   zrecall(cz[1]);
   zstore(cstr[1]);
   return 0;
  }

int crecall(unsigned int * *cstr)
  {zrecall(cstr[0]); 
   zstore(cz[0]);
   zrecall(cstr[1]);
   zstore(cz[1]);
   return 0;
  }

int crecip(unsigned int * *cz1)
  {crecall(cz1);
   cstore(ctempcomp3);
   cconinttomp(1,0);
   cdiv(cz,ctempcomp3);
   return 0;
  }


int cchs(unsigned int * *cz1)
  {zrecall(cz1[0]);
   zchs(z);
   zstore(cz[0]);
   zrecall(cz1[1]);
   zchs(z);
   zstore(cz[1]);
   return 0;
  }

int cchspol(unsigned int * *cz1)
  {zpi();
   zadd(z,cz1[1]);
   zstore(cz[1]);
   zrecall(cz1[0]);
   zstore(cz[0]);
   return 0;
  }

int cadd(unsigned int * *cz1, unsigned int * *cz2) 
  {zadd(cz1[0],cz2[0]);
   zstore(cz[0]);
   zadd(cz1[1],cz2[1]);
   zstore(cz[1]);
   return 0;
  }

int csub(unsigned int * *cz1, unsigned int * *cz2)
  {zsub(cz1[0],cz2[0]);
   zstore(cz[0]);
   zsub(cz1[1],cz2[1]);
   zstore(cz[1]);
   return 0;
  }
 
int cmult(unsigned int * *cz3, unsigned int * *cz4)
  {zmult(cz3[0],cz4[0]);
   zstore(ctempreal);
   zmult(cz3[1],cz4[1]);
   zsub(ctempreal,z);
   zstore(ctempreal6);
   zmult(cz3[0],cz4[1]);
   zstore(ctempreal);
   zmult(cz3[1],cz4[0]);
   zadd(ctempreal,z);
   zstore(cz[1]);
   zrecall(ctempreal6);
   zstore(cz[0]);
   return 0;
  }

int cmultpol(unsigned int * *cz3, unsigned int * *cz4)
  {zmult(cz3[0],cz4[0]);
   zstore(cz[0]);
   zadd(cz3[1],cz4[1]);
   zstore(cz[1]);
   return 0;
  }

int cdiv(unsigned int * *cz5,unsigned int * *cz6)
  {int zij;
   for (zij=0;zij<prec+2;zij++)
     {ctempcomp[0][zij]=cz5[0][zij];
      ctempcomp[1][zij]=cz5[1][zij];
      ctempcomp1[0][zij]=cz6[0][zij];
      ctempcomp1[1][zij]=cz6[1][zij];}
   zmult(ctempcomp1[0],ctempcomp1[0]);
   zstore(ctempreal2);
   zmult(ctempcomp1[1],ctempcomp1[1]);
   zadd(z,ctempreal2);
   zstore(ctempreal2);         
   ccompconj(ctempcomp1);
   cmult(ctempcomp,cz);
   zdiv(cz[0],ctempreal2);
   zstore(cz[0]);               
   zdiv(cz[1],ctempreal2);
   zstore(cz[1]);

   return 0;
  }

int cdivpol(unsigned int * *cz5,unsigned int * *cz6)
  {
   zdiv(cz5[0],cz6[0]);
   zstore(cz[0]);
   zsub(cz5[1],cz6[1]);
   zstore(cz[1]);
   return 0;
  }

int cabs(unsigned int * *cz1)
  {zmult(cz1[0],cz1[0]);
   zstore(ctempreal);
   zmult(cz1[1],cz1[1]);
   zadd(ctempreal,z);
   zsqrt(z);
   return 0;
  }

int ccompconj(unsigned int * *cz1)
  {zchs(z);             
   zrecall(cz1[0]);
   zstore(cz[0]);
   zrecall(cz1[1]);
   zchs(z);
   zstore(cz[1]);
   return 0;
  }

int ccompconjpol(unsigned int * *cz1)
  {zrecall(cz1[0]);
   zstore(cz[0]);
   zrecall(cz1[1]);
   zchs(z);
   zstore(cz[1]);
   return 0;
  }

int crectopol(unsigned int * *cz26)
  {zpi();
   zstore(cpie);
   if (cz26[0][prec]==0)
     {zdivint(cpie,2);      
      goto cyurigg;}
   zdiv(cz26[1],cz26[0]);
   zarctan(z);              
cyurigg: czi=cz26[0][prec+1];
   if (czi==0) 
     {zadd(z,cpie);}     
   zstore(ctempreal3);   
   cabs(cz26);              
   zstore(cz[0]);
   zrecall(ctempreal3);
   zstore(cz[1]);
   return 0;
  }

int cpoltorec(unsigned int * *cz26)
  {zcos(cz26[1]);        
   zmult(z,cz26[0]);     
   zstore(ctempreal3);
   zsin(cz26[1]);        
   zmult(z,cz26[0]);     
   zstore(cz[1]);        
   zrecall(ctempreal3);  
   zstore(cz[0]);
   return 0;
  }

int cxtonthpol(unsigned int * *cz7,unsigned int cn)
  {zxtonth(cz7[0],cn);
   zstore(cz[0]);
   zmultint(cz7[1],cn);
   zstore(cz[1]);
   return 0;
  }

int zsinh(unsigned int *zzyy1)
  {zrecall(zzyy1);
   zstore(z3boy);
   zexp(z3boy);
   zstore(ctempreal);
   zrecall(z3boy);
   zchs(z);
   zexp(z);
   zsub(ctempreal,z);
   zdivint(z,2);
   return 0;
  }

int zcosh(unsigned int *zy2)
  {zrecall(zy2);
   zstore(z3boy);
   zexp(z3boy);
   zstore(ctempreal);
   zrecall(z3boy);
   zchs(z);
   zexp(z);
   zadd(ctempreal,z);
   zdivint(z,2);
   return 0;
  }
  
int csin(unsigned int * *cz8)
  {int zi;
   for (zi=0;zi<prec+2;zi++) 
     {cz88[0][zi]=cz8[0][zi];
      cz88[1][zi]=cz8[1][zi];}
   zsin(cz88[0]);
   zstore(ctempreal4);
   zcosh(cz88[1]);
   zmult(z,ctempreal4);
   zstore(ctempreal4);
   zcos(cz88[0]);
   zstore(ctempreal5);
   zsinh(cz88[1]);
   zmult(z,ctempreal5);
   zstore(cz[1]);
   zrecall(ctempreal4);
   zstore(cz[0]);
   return 0;
  }

int ccos(unsigned int * *cz9)
  {int zi;
   for (zi=0;zi<prec+2;zi++) 
     {cz99[0][zi]=cz9[0][zi];
      cz99[1][zi]=cz9[1][zi];}
   zcos(cz99[0]);
   zstore(ctempreal4);
   zcosh(cz99[1]);
   zmult(z,ctempreal4);
   zstore(ctempreal4);
   zsin(cz99[0]);
   zstore(ctempreal5);
   zsinh(cz99[1]);
   zmult(z,ctempreal5);
   zchs(z);
   zstore(cz[1]);
   zrecall(ctempreal4);
   zstore(cz[0]);
   return 0;
  }

int clogpolinrec(unsigned int * *cz8)
  {zlog(cz8[0]);
   zstore(cz[0]);
   zrecall(cz8[1]);
   zstore(cz[1]);
   return 0;
  }

int cln(unsigned int * *cz9)
  {crectopol(cz9);
   clogpolinrec(cz);
   return 0;
  }

int cexp(unsigned int * *cz9)
  {zexp(cz9[0]);
   zstore(ctempreal5);
   zcos(cz9[1]);
   zmult(z,ctempreal5);
   zstore(cz[0]);
   zsin(cz9[1]);
   zmult(z,ctempreal5);
   zstore(cz[1]);
   return 0;
  }

int carctan(unsigned int * *cz9)
   {int czik,zsignflag;
    zsignflag=0;
    for (czik=0;czik<prec+2;czik++)
      {cz99[0][czik]=cz9[0][czik];
       cz99[1][czik]=cz9[1][czik];}
    zmult(cz99[0],cz99[0]);
    zstore(ctempreal6);
    zmult(cz99[1],cz99[1]);
    zadd(z,ctempreal6);
    zstore(ctempreal);         
    zconinttomp(1);
    zsub(z,ctempreal);
    zstore(ctempreal);
    zmultint(cz99[0],2);
    if (z[prec+1]==0) z[prec+1]=2;
    zstore(ctempreal2);                 
    zrecall(ctempreal);
    if (z[prec+1]==0) z[prec+1]=2;    
    zsub(z,ctempreal2);                 
    if (z[prec+1]==0) goto youch;
    if (ctempreal[prec+1]==0) zsignflag=1;
    zdiv(ctempreal2,ctempreal);
    zarctan(z);
    zstore(ctempreal);
    if (zsignflag==1) 
       {zpi();
        zadd(z,ctempreal);
        zstore(ctempreal);}
    zdivint(ctempreal,2);
    zstore(ctempreal);
    goto youch2;
youch:
    if (ctempreal2[prec+1]==0) zsignflag=1;
    zdiv(ctempreal,ctempreal2);
    zarctan(z);
    zstore(ctempreal);
    if (zsignflag==1) 
      {zpi();
       zadd(z,ctempreal);
       zstore(ctempreal);}
    zpi();
    zdivint(z,2);
    zsub(z,ctempreal);
    zstore(ctempreal);
    zdivint(ctempreal,2);
    zstore(ctempreal);
youch2:
    zconinttomp(1);
    zsub(cz99[1],z);
    zmult(z,z);
    zadd(z,ctempreal6);
    zstore(ctempreal2);
    zconinttomp(1);         
    zadd(cz99[1],z);
    zmult(z,z);
    zadd(z,ctempreal6);
    zdiv(z,ctempreal2);
    zlog(z);
    zdivint(z,4);
    zstore(cz[1]);
    zrecall(ctempreal);
    zstore(cz[0]);
    return 0;
   }

int cconinttomp(unsigned int cr1,unsigned int cr2)
  {zconinttomp(cr1);
   zstore(cz[0]);
   zconinttomp(cr2);
   zstore(cz[1]);
   return 0;
  }

int cconinttomp(int cr1,int cr2)        
  {zconinttomp(cr1);
   zstore(cz[0]);
   zconinttomp(cr2);
   zstore(cz[1]);
   return 0;
  }

int czero(unsigned int * *cz1)
  {zzero(cz1[0]);
   zzero(cz1[1]);
   return 0;
  }

int cmultint(unsigned int * *cz1,unsigned int cr1)
  {zmultint(cz1[0],cr1);
   zstore(cz[0]);
   zmultint(cz1[1],cr1);
   zstore(cz[1]);
   return 0;
  }

int cmultint(unsigned int * *cz1,int cr1)       
  {zmultint(cz1[0],cr1);
   zstore(cz[0]);
   zmultint(cz1[1],cr1);
   zstore(cz[1]);
   return 0;
  }

int cdivint(unsigned int * *cz1,unsigned int cr1)
  {zdivint(cz1[0],cr1);
   zstore(cz[0]);
   zdivint(cz1[1],cr1);
   zstore(cz[1]);
   return 0;
  }

int cdivint(unsigned int * *cz1,int cr1)          
  {zdivint(cz1[0],cr1);
   zstore(cz[0]);
   zdivint(cz1[1],cr1);
   zstore(cz[1]);
   return 0;
  }

int cmultcbyr(unsigned int * *cz1,unsigned int * crfac)
  {
   zrecall(crfac);
   zstore(z3boy);
   zmult(z3boy,cz1[0]);
   zstore(cz[0]);
   zmult(z3boy,cz1[1]);
   zstore(cz[1]);
   return 0;
  }

int cdivcbyr(unsigned int * *cz1,unsigned int * crfac)
  {zrecall(crfac);
   zstore(z3boy);
   zdiv(cz1[0],z3boy);
   zstore(cz[0]);
   zdiv(cz1[1],z3boy);
   zstore(cz[1]);
   return 0;
  }

int caddcplusr(unsigned int * * cz1,unsigned int * crfac)
  {zrecall(crfac);
   zstore(z3boy);
   zadd(cz1[0],z3boy);
   zstore(cz[0]);
   zrecall(cz1[1]);     
   zstore(cz[1]);
   return 0;
  }

int cexpitheta(unsigned int *crfac)
  {zrecall(crfac);
   zstore(z3boy);
   zcos(z3boy);
   zstore(cz[0]);
   zsin(z3boy);
   zstore(cz[1]);
   return 0;
  }

int csqrt(unsigned int * *cz0)    
  {crecall(cz0);
   crectopol(cz);
   csqrtpol(cz);      
   cpoltorec(cz);     
   return 0;
  }

int csqrtpol(unsigned int * *cz9)
  {zrecall(cz9[0]);
   zsqrt(z);
   zstore(cz[0]);
   zrecall(cz9[1]);
   zdivint(z,2);
   zstore(cz[1]);
   return 0;
  }