/* List of functions in this file:
complexmp operator+(complexmp,complexmp)          //1
complexmp operator+(realmp,complexmp)             //2
complexmp operator+(complexmp,realmp)             //3
complexmp operator+(complexmp,int)
complexmp operator+(complexmp,unsigned int)
complexmp operator+(int,complexmp)
complexmp operator+(unsigned int,complexmp)
complexmp operator-(complexmp,complexmp)          //4
complexmp operator-(realmp,complexmp)             //5
complexmp operator-(complexmp,realmp)             //6
complexmp operator-(complexmp)                    //7 unary -
complexmp operator-(complexmp,int)
complexmp operator-(complexmp,unsigned int)
complexmp operator-(int,complexmp)
complexmp operator-(unsigned int,complexmp)
complexmp operator*(complexmp,complexmp)          //8
complexmp operator*(complexmp,int)                //9
complexmp operator*(complexmp,unsigned int)       //9.5
complexmp operator*(int,complexmp)                //10
complexmp operator*(unsigned int,complexmp)       //10.5
complexmp operator*(complexmp,realmp)             //25
complexmp operator*(realmp,complexmp)             //26
complexmp operator/(complexmp,unsigned int)       //11
complexmp operator/(complexmp,int)                //11.5
complexmp operator/(int,complexmp)                //28
complexmp operator/(unsigned int,complexmp)       //28.5
complexmp operator/(complexmp,complexmp)          //12
complexmp operator/(complexmp,realmp)             //27
complexmp operator/(realmp,complexmp)             //29
void operator<(complexmp,complexmp)               //13
printmp(complexmp)                                //14
complexmp cinttomp(int,int)                       //15
complexmp cfloattomp(float,float)                 //16
complexmp sqrt(complexmp)                         //17
complexmp sin(complexmp)                          //18
complexmp cos(complexmp)                          //19
realmp sinh(realmp)                               //20
realmp cosh(realmp)                               //21             
complexmp log(complexmp)                          //22
complexmp exp(complexmp)                          //23
complexmp arctan(complexmp)                       //24              
complexmp crectopol(complexmp)                    //30
complexmp cpoltorec(complexmp)                    //31
complexmp compconj(complexmp)                     //32
realmp abs(complexmp)                             //33
realmp realpartof(complexmp)                      //34
realmp imaginarypartof(complexmp)                 //35
complexmp cexpitheta(realmp)                      //36
void cwritetofile(complexmp,char *)               //37
void creadfromfile(complexmp,char *)              //38
void cwrite1Darraytofile(complexmp z[cmaxsizepad], char *)    //39
void cread1Darrayfromfile(complexmp z[cmaxsizepad], char *)   //40
void cwrite2Darraytofile(complexmp z[cmaxsizepad][cmaxsizepad],char *) //41
void cread2Darrayfromfile(complexmp z[cmaxsizepad][cmaxsizepad],char *) //42
*/

#include "cclasses.new"

complexmp cwreturn[depthcountmax];
realmp cxcxcx;
realmp xcxcxc;

complexmp operator+(complexmp cw1,complexmp cw2)            //1
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cw1.gv(0,zi);
      cz1[1][zi]=cw1.gv(1,zi);
      cz2[0][zi]=cw2.gv(0,zi);
      cz2[1][zi]=cw2.gv(1,zi);
     }
   cadd(cz1,cz2);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator+(realmp cwrfac,complexmp cw2)    //2
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwrfac.gv(zi);
      cz1[1][zi]=0;
      cz2[0][zi]=cw2.gv(0,zi);
      cz2[1][zi]=cw2.gv(1,zi);
     }
   cadd(cz1,cz2);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }
                                                  
complexmp operator+(complexmp cw2,realmp cwrfac)      //3
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwrfac.gv(zi);
      cz1[1][zi]=0;
      cz2[0][zi]=cw2.gv(0,zi);
      cz2[1][zi]=cw2.gv(1,zi);
     }
   cadd(cz1,cz2);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator+(complexmp cw2, int cr1)
  {int zi;
   for (zi=0;zi<prec+2;zi++)
      {cz2[0][zi]=cw2.gv(0,zi);
       cz2[1][zi]=cw2.gv(1,zi);
      }
   cconinttomp(cr1,0);      
   cadd(cz2,cz);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }
   
complexmp operator+(complexmp cw2,unsigned int cr1)
  {int zi;
   unsigned int zthatzero;
   for (zi=0;zi<prec+2;zi++)
      {cz2[0][zi]=cw2.gv(0,zi);
       cz2[1][zi]=cw2.gv(1,zi);
      }
   cconinttomp(cr1,zthatzero);      
   cadd(cz2,cz);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator+(int cr1,complexmp cw2)
  {int zi;
   for (zi=0;zi<prec+2;zi++)
      {cz2[0][zi]=cw2.gv(0,zi);
       cz2[1][zi]=cw2.gv(1,zi);
      }
   cconinttomp(cr1,0);
   cadd(cz2,cz);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator+(unsigned int cr1,complexmp cw2)
  {int zi;
   unsigned int zzxxzz=0;
   for (zi=0;zi<prec+2;zi++)
      {cz2[0][zi]=cw2.gv(0,zi);
       cz2[1][zi]=cw2.gv(1,zi);
      }
   cconinttomp(cr1,zzxxzz);
   cadd(cz2,cz);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator-(complexmp cw1,complexmp cw2)      //4
  {
   int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cw1.gv(0,zi);
      cz1[1][zi]=cw1.gv(1,zi);
      cz2[0][zi]=cw2.gv(0,zi);
      cz2[1][zi]=cw2.gv(1,zi);
     }
   csub(cz1,cz2);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator-(realmp cwrfac,complexmp cw2)       //5
  {
   int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwrfac.gv(zi);
      cz1[1][zi]=0;
      cz2[0][zi]=cw2.gv(0,zi);
      cz2[1][zi]=cw2.gv(1,zi);
     }
   csub(cz1,cz2);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator-(complexmp cw2,realmp cwrfac)         //6
  {
   int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz2[0][zi]=cwrfac.gv(zi);
      cz2[1][zi]=0;
      cz1[0][zi]=cw2.gv(0,zi);
      cz1[1][zi]=cw2.gv(1,zi);
     }
   csub(cz1,cz2);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }


complexmp operator-(complexmp cw2)  //a=-b does not change sign of b
 {int zi;                                   //a=-a changes sign of a
  depthcount++;                             //7
  for (zi=0;zi<prec+1;zi++) cwreturn[depthcount].sv(0,zi,cw2.gv(0,zi));
  for (zi=0;zi<prec+1;zi++) cwreturn[depthcount].sv(1,zi,cw2.gv(1,zi));
  cwreturn[depthcount].sv(0,prec+1,2-cw2.gv(0,prec+1));
  cwreturn[depthcount].sv(1,prec+1,2-cw2.gv(1,prec+1));
  return cwreturn[depthcount];
}

complexmp operator-(complexmp cw2, int cr1)
  {int zi;
   int zthatzero=0;
   for (zi=0;zi<prec+2;zi++)
      {cz2[0][zi]=cw2.gv(0,zi);
       cz2[1][zi]=cw2.gv(1,zi);
      }
   cconinttomp(cr1,zthatzero);      
   csub(cz2,cz);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }
   
complexmp operator-(complexmp cw2,unsigned int cr1)
  {int zi;
   unsigned int zj=0;
   for (zi=0;zi<prec+2;zi++)
      {cz2[0][zi]=cw2.gv(0,zi);
       cz2[1][zi]=cw2.gv(1,zi);
      }
   cconinttomp(cr1,zj);      
   csub(cz2,cz);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator-(int cr1,complexmp cw2)
  {int zi;
   for (zi=0;zi<prec+2;zi++)
      {cz2[0][zi]=cw2.gv(0,zi);
       cz2[1][zi]=cw2.gv(1,zi);
      }
   cconinttomp(cr1,0);
   csub(cz2,cz);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator-(unsigned int cr1,complexmp cw2)
  {int zi;
   unsigned int zthatzero;
   for (zi=0;zi<prec+2;zi++)
      {cz2[0][zi]=cw2.gv(0,zi);
       cz2[1][zi]=cw2.gv(1,zi);
      }
   cconinttomp(cr1,zthatzero);
   csub(cz2,cz);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator*(complexmp cw1,complexmp cw2)    //8
  {
   int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cw1.gv(0,zi);
      cz1[1][zi]=cw1.gv(1,zi);
      cz2[0][zi]=cw2.gv(0,zi);
      cz2[1][zi]=cw2.gv(1,zi);
     }
   cmult(cz1,cz2);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator*(complexmp cwz1,int cr1)       //9
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwz1.gv(0,zi);
      cz1[1][zi]=cwz1.gv(1,zi);
     }
   cmultint(cz1,cr1);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator*(complexmp cwz1,unsigned int cr1)       //9.5
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwz1.gv(0,zi);
      cz1[1][zi]=cwz1.gv(1,zi);
     }
   cmultint(cz1,cr1);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator*(int cr1, complexmp cwz1)       //10
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwz1.gv(0,zi);
      cz1[1][zi]=cwz1.gv(1,zi);
     }
   cmultint(cz1,cr1);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator*(unsigned int cr1, complexmp cwz1)       //10.5
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwz1.gv(0,zi);
      cz1[1][zi]=cwz1.gv(1,zi);
     }
   cmultint(cz1,cr1);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator/(complexmp cwz1,unsigned int cr1)   //11
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwz1.gv(0,zi);
      cz1[1][zi]=cwz1.gv(1,zi);
     }
   cdivint(cz1,cr1);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator/(complexmp cwz1,int cr1)   //11.5
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwz1.gv(0,zi);
      cz1[1][zi]=cwz1.gv(1,zi);
     }
   cdivint(cz1,cr1);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator/(complexmp cw1,complexmp cw2)      //12
  {
   int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cw1.gv(0,zi);
      cz1[1][zi]=cw1.gv(1,zi);
      cz2[0][zi]=cw2.gv(0,zi);
      cz2[1][zi]=cw2.gv(1,zi);
     }
   cdiv(cz1,cz2);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

void operator<(complexmp cw1,complexmp cw2)         //13
{ int zi;
  for (zi=0;zi<prec+2;zi++)
    {cw1.sv(0,zi,cw2.gv(0,zi));
     cw1.sv(1,zi,cw2.gv(1,zi));}
  depthcount=2;
}

int printmp(complexmp cwp)                           //14
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwp.gv(0,zi);
      cz1[1][zi]=cwp.gv(1,zi);}
   cprnt(cz1);
   depthcount=2;
   return 0;
  }

complexmp cinttomp(int cwr1,int cwr2)        //15
  {int zi;
   cconinttomp(cwr1,cwr2);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   if (cwr1<0)
     {cwreturn[depthcount].sv(0,prec+1,0);}
   if (cwr2<0)
     {cwreturn[depthcount].sv(1,prec+1,0);}
   return cwreturn[depthcount];
  }

complexmp cfloattomp(float cwr1,float cwr2)   //16
  {int zi;
   xcxcxc=conftomp(cwr1);
   cxcxcx=conftomp(cwr2);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,xcxcxc.gv(zi));
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cxcxcx.gv(zi));
   return cwreturn[depthcount];
  }

complexmp sqrt(complexmp cw0)  //Nov25,99     //17
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz0[0][zi]=cw0.gv(0,zi);
      cz0[1][zi]=cw0.gv(1,zi);
     }
   csqrt(cz0);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp sin(complexmp cw8)        //Nov26,99  //18
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz8[0][zi]=cw8.gv(0,zi);
      cz8[1][zi]=cw8.gv(1,zi);
     }
   csin(cz8);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }
  
complexmp cos(complexmp cw8)      //Nov26,99    //19
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz8[0][zi]=cw8.gv(0,zi);
      cz8[1][zi]=cw8.gv(1,zi);
     }
   ccos(cz8);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

realmp sinh(realmp zzwyy1)      //Nov26,99     //20
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {zzyy1[zi]=zzwyy1.gv(zi);
     }
   zsinh(zzyy1);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) wwreturn[depthcount].sv(zi,z[zi]);
   return wwreturn[depthcount];
  }

realmp cosh(realmp zwy2)         //Nov26,99   //21             
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {zy2[zi]=zwy2.gv(zi);
     }
   zcosh(zy2);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) wwreturn[depthcount].sv(zi,z[zi]);
   return wwreturn[depthcount];
  }

complexmp log(complexmp cwz9)    //Nov26,99   //22
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz9[0][zi]=cwz9.gv(0,zi);
      cz9[1][zi]=cwz9.gv(1,zi);
     }
   cln(cz9);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp exp(complexmp cwz9)    //Nov26,99      //23
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz9[0][zi]=cwz9.gv(0,zi);
      cz9[1][zi]=cwz9.gv(1,zi);
     }
   cexp(cz9);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp arctan(complexmp cwz0)   //Nov26,99    //24              
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz0[0][zi]=cwz0.gv(0,zi);
      cz0[1][zi]=cwz0.gv(1,zi);
     }
   carctan(cz0);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator*(complexmp cwz1,realmp cwrfac)    //25
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwz1.gv(0,zi);
      cz1[1][zi]=cwz1.gv(1,zi);
      crfac[zi]=cwrfac.gv(zi);
      crfac[zi]=cwrfac.gv(zi);
     }
   cmultcbyr(cz1,crfac);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }


complexmp operator*(realmp cwrfac,complexmp cwz1)    //26
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwz1.gv(0,zi);
      cz1[1][zi]=cwz1.gv(1,zi);
      crfac[zi]=cwrfac.gv(zi);
     }
   cmultcbyr(cz1,crfac);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator/(complexmp cwz1,realmp cwrfac)    //27
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwz1.gv(0,zi);
      cz1[1][zi]=cwz1.gv(1,zi);
      crfac[zi]=cwrfac.gv(zi);
     }
   cdivcbyr(cz1,crfac);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }


complexmp operator/(int cr1,complexmp cwz2)         //28
  {cconinttomp(cr1,0);
   int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwz2.gv(0,zi);
      cz1[1][zi]=cwz2.gv(1,zi);
     }
   cconinttomp(cr1,0);
   cdiv(cz,cz1);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp operator/(realmp cwrfac,complexmp cwz3)   //29
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz3[0][zi]=cwz3.gv(0,zi);
      cz3[1][zi]=cwz3.gv(1,zi);
      crfac[zi]=cwrfac.gv(zi);
     }
   zmult(cz3[0],cz3[0]);
   zstore(ctempreal);
   zmult(cz3[1],cz3[1]);
   zadd(z,ctempreal);
   zdiv(crfac,z);
   zstore(ctempreal);
   zmult(ctempreal,cz3[0]);
   zstore(ctempreal2);
   zmult(ctempreal,cz3[1]);
   zchs(z);
   zstore(ctempreal);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,ctempreal2[zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,ctempreal[zi]);
   return cwreturn[depthcount];
  }

complexmp crectopol(complexmp cwz26)    //Nov26,99   //30
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz26[0][zi]=cwz26.gv(0,zi);
      cz26[1][zi]=cwz26.gv(1,zi);
     }
   crectopol(cz26);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp cpoltorec(complexmp cwz26)     //Nov26,99   //31
  {int zi;
   for (zi=0;zi<prec+2;zi++)
     {cz26[0][zi]=cwz26.gv(0,zi);
      cz26[1][zi]=cwz26.gv(1,zi);
     }
   cpoltorec(cz26);
   depthcount++;
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(0,zi,cz[0][zi]);
   for (zi=0;zi<prec+2;zi++) cwreturn[depthcount].sv(1,zi,cz[1][zi]);
   return cwreturn[depthcount];
  }

complexmp compconj(complexmp cwz1)    //Nov26,99      //32
  {int zi;
   depthcount++;
   for (zi=0;zi<prec+2;zi++)
     {cwreturn[depthcount].sv(0,zi,cwz1.gv(0,zi));}
   for (zi=0;zi<prec+1;zi++)
     {cwreturn[depthcount].sv(1,zi,cwz1.gv(1,zi));}
   cwreturn[depthcount].sv(1,prec+1,2-cwz1.gv(1,prec+1)); //chs of imag
   return cwreturn[depthcount];                             //part
  }

realmp abs(complexmp cwz1)           //Nov26,99     //33
  {int zi;
   depthcount++;
   for (zi=0;zi<prec+2;zi++)
     {cz1[0][zi]=cwz1.gv(0,zi);
      cz1[1][zi]=cwz1.gv(1,zi);
     }
   cabs(cz1);
   depthcount++;
   for (zi=0;zi<prec+2;zi++)
     {wwreturn[depthcount].sv(zi,z[zi]);}
   return wwreturn[depthcount];  
  }

realmp realpartof(complexmp cwz1)             //34
  {int zi;
   depthcount++;
   for (zi=0;zi<prec+2;zi++)
     {wwreturn[depthcount].sv(zi,cwz1.gv(0,zi));}
   return wwreturn[depthcount];
  }

realmp imaginarypartof(complexmp cwz1)       //35
  {int zi;
   depthcount++;
   for (zi=0;zi<prec+2;zi++)
     {wwreturn[depthcount].sv(zi,cwz1.gv(1,zi));}
   return wwreturn[depthcount];
  }  

complexmp cexpitheta(realmp cwrfac)    //Nov26,99  //36
  {int zi;
   depthcount++;
   for (zi=0;zi<prec+2;zi++)
     {crfac[zi]=cwrfac.gv(zi);}
   cexpitheta(crfac);
   for (zi=0;zi<prec+2;zi++)
     {cwreturn[depthcount].sv(0,zi,cz[0][zi]);
      cwreturn[depthcount].sv(1,zi,cz[1][zi]);}
   return cwreturn[depthcount];
  }

void cwritetofile(complexmp cvariablename,char *ztoit)  //37
  {int zi3;
   zfileout.open(ztoit, ios::out);
   for (zi3=0;zi3<prec+2;zi3++)
     {zfileout<<cvariablename.gv(0,zi3)<<'\n';
      zfileout<<cvariablename.gv(1,zi3)<<'\n';
     }
   zfileout.close();
  }

void creadfromfile(complexmp cvariablename,char *zfromit) //38
  {int zi3;
   zfilein.open(zfromit, ios::in);
   for (zi3=0;zi3<prec+2;zi3++)
    {zfilein>>cz[0][zi3];
     zfilein>>cz[1][zi3];
    cvariablename.sv(0,zi3,cz[0][zi3]);
    cvariablename.sv(1,zi3,cz[1][zi3]);
    }
   zfilein.close();
  }

void cwrite1Darraytofile(complexmp cvariablename[cmaxsizepad], char *ztoit) //39
  {int zi3,zj3;
   zfileout.open(ztoit, ios::out);
   for (zi3=0;zi3<cmaxsize+1;zi3++)
     {for (zj3=0;zj3<prec+2;zj3++)
       {zfileout<<cvariablename[zi3].gv(0,zj3)<<'\n';
        zfileout<<cvariablename[zi3].gv(1,zj3)<<'\n';
       }
     }
   zfileout.close();
  }

void cread1Darrayfromfile(complexmp cvariablename[cmaxsizepad], char *zfromit) //40
  {int zi3,zj3;
   zfilein.open(zfromit, ios::in);
   for (zi3=0;zi3<cmaxsize+1;zi3++)
     {for (zj3=0;zj3<prec+2;zj3++)
       {zfilein>>cz[0][zj3];
        zfilein>>cz[1][zj3];
        cvariablename[zi3].sv(0,zj3,cz[0][zj3]);
        cvariablename[zi3].sv(1,zj3,cz[1][zj3]);
       }
     }
   zfilein.close();
  }

void cwrite2Darraytofile(complexmp cvariablename[cmaxsizepad][cmaxsizepad],char *ztoit) //41
  {int zi3,zj3,zk3;
   zfileout.open(ztoit, ios::out);
   for (zi3=0;zi3<cmaxsize+1;zi3++)
    {for (zj3=0;zj3<cmaxsize+1;zj3++)
{for (zk3=0;zk3<prec+2;zk3++)
  {zfileout<<cvariablename[zi3][zj3].gv(0,zk3)<<'\n';
   zfileout<<cvariablename[zi3][zj3].gv(1,zk3)<<'\n';}
}
    }
   zfileout.close();

  }
  
void cread2Darrayfromfile(complexmp cvariablename[cmaxsizepad][cmaxsizepad],char *zfromit) //42
  {int zi3,zj3,zk3;
   zfilein.open(zfromit, ios::in);
   for (zi3=0;zi3<cmaxsize+1;zi3++)
    {for (zj3=0;zj3<cmaxsize+1;zj3++)
{for (zk3=0;zk3<prec+2;zk3++)
  {zfilein>>cz[0][zk3];
   zfilein>>cz[1][zk3];
   cvariablename[zi3][zj3].sv(0,zk3,cz[0][zk3]);
   cvariablename[zi3][zj3].sv(1,zk3,cz[1][zk3]);
  }
}
    }
   zfilein.close();
  }