class complexmp {
 public:
    complexmp();
    int gv(int zi,int zj) const { return *zelement[zi][zj]; }
    void sv(int zi,int zj,unsigned int zterm) {*zelement[zi][zj]=zterm; }
    complexmp operator=(const complexmp &);
    complexmp operator=(const int &);
 private:
    unsigned int ***zelement;
};

complexmp::complexmp()
  {int zi2;
   zelement=new unsigned int * * [2];
   zelement[0]=new unsigned int * [zprecpad];
   zelement[1]=new unsigned int * [zprecpad];
   for (zi2=0;zi2<zprecpad+2;zi2++)
     {zelement[0][zi2]=new unsigned int [1];
      zelement[1][zi2]=new unsigned int [1];
     }
   for (zi2=0;zi2<prec+2;zi2++)
     {zelement[0][zi2]=new unsigned int;
      zelement[1][zi2]=new unsigned int;
      *zelement[0][zi2]=0;
      *zelement[1][zi2]=0;
     }
  }

complexmp complexmp::operator=(const complexmp & rhs)     //there was an &
  {int zi3;
   if (this==&rhs) return *this;
   depthcount=2;
   for (zi3=0;zi3<prec+2;zi3++)
     {*zelement[0][zi3]=rhs.gv(0,zi3);
      *zelement[1][zi3]=rhs.gv(1,zi3);
     }
   return *this;
  }


complexmp complexmp::operator=(const int & rhs)
  {int zi3;
   unsigned int rhss;
   depthcount=2;
   for (zi3=prec+1;zi3>=0;zi3--)
     {if (zi3==prec+1)
         {*zelement[1][zi3]=2;
          *zelement[0][zi3]=2;
          rhss=rhs;
          if (rhs<0)
            {*zelement[0][zi3]=0;
             rhss=-rhs;
            }
         }
      if (zi3==prec)
         {*zelement[1][zi3]=0;
          *zelement[0][zi3]=rhss;}
      if (zi3<prec)
        {*zelement[0][zi3]=0;
         *zelement[1][zi3]=0;
         if (zi3==0)
           {*zelement[0][zi3]=0x40000001;
            *zelement[1][zi3]=0x40000001;
           }
        }
     }
   return *this;
  }

class cpademp
{public:
 complexmp *cpnummp;
 complexmp *cpdenmp;
 cpademp();
};

cpademp::cpademp()
  {cpnummp=new complexmp[cmaxsizepad];
   cpdenmp=new complexmp[cmaxsizepad];
  }