/*
void operator<(cpademp cwpout,cpademp czz2p)
cpademp cpadeapprox(int cpn,int cpm,complexmp cwps[cmaxsizepad])
cwpoly cextractnum(int cpn,cpademp cwpademp)     //was cpoly
cwpoly cextractden(int cpm,cpademp cz2pmp)
*/

extern cpademp cpadeapprox(int cpn,int cpm,complexmp cwps[cmaxsizepad]);

cpademp cwpademp;  
cwpoly cpadepolyoutmp;       //changed

//pademp wpademp;          
//wpoly padepolyoutmp;


cpademp cpadeapprox(int cpn,int cpm,complexmp cwps[cmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<cpn+cpm+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
        {cps[zi][0][zj]=cwps[zi].gv(0,zj);
         cps[zi][1][zj]=cwps[zi].gv(1,zj);
        }
     }
   cpadeapprox(cpn,cpm,cps);
   for (zi=0;zi<cpn+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {cwpademp.cpnummp[zi].sv(0,zj,cpnum[zi][0][zj]);
        cwpademp.cpnummp[zi].sv(1,zj,cpnum[zi][1][zj]);
       }
     }
   for (zi=0;zi<cpm+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {cwpademp.cpdenmp[zi].sv(0,zj,cpden[zi][0][zj]);
        cwpademp.cpdenmp[zi].sv(1,zj,cpden[zi][1][zj]);
       }
     }
   return cwpademp;
  }    

void operator<(cpademp cwpout,cpademp czz2p)
  {int zi;
   for (zi=0;zi<cmaxsize+1;zi++)
    {cwpout.cpnummp[zi]=czz2p.cpnummp[zi];
     cwpout.cpdenmp[zi]=czz2p.cpdenmp[zi];}
  }

cwpoly cextractnum(int cpn,cpademp cwpademp)     //was cpoly
{unsigned int zi;
 int zj;
 for (zi=0;zi<cpn+1;zi++)
   {for (zj=0;zj<prec+2;zj++)
     {cwpolyout[zi].sv(0,zj,cwpademp.cpnummp[zi].gv(0,zj));
      cwpolyout[zi].sv(1,zj,cwpademp.cpnummp[zi].gv(1,zj));
     }
   }
 return cpadepolyoutmp;
}

cwpoly cextractden(int cpm,cpademp cz2pmp)
{int zi,zj;
 for (zi=0;zi<cpm+1;zi++)
   {for (zj=0;zj<prec+2;zj++)
     {cwpolyout[zi].sv(0,zj,cwpademp.cpdenmp[zi].gv(0,zj));
      cwpolyout[zi].sv(1,zj,cwpademp.cpdenmp[zi].gv(1,zj));
     }
   }
 return cpadepolyoutmp;
}