/*
pademp padeapprox(int,int,realmp s[zmaxsizepad])    //padeapprox
void operator<(pademp xpademp,pademp wpademp)       //overload <
wpoly extractnum(int pn,pademp wpademp)             //extractnum poly
wpoly extractden(int pm,pademp wpademp)             //extractden poly
*/

extern pademp padeapprox(int pn,int pm,realmp pws[zmaxsizepad]);

pademp wpademp;          
wpoly padepolyoutmp;

pademp padeapprox(int pn,int pm,realmp pws[zmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<pn+pm+1;zi++)
      {for (zj=0;zj<prec+2;zj++) ps[zi][zj]=pws[zi].gv(zj);}
   padeapprox(pn,pm,ps);
   for (zi=0;zi<pn+1;zi++)
      {for (zj=0;zj<prec+2;zj++)
        {wpademp.nummp[zi].sv(zj,pnum[zi][zj]);}
      }
   for (zi=0;zi<pm+1;zi++)
      {for (zj=0;zj<prec+2;zj++)
        {wpademp.denmp[zi].sv(zj,pden[zi][zj]);}
      }
   return wpademp;
  }

void operator<(pademp xpademp,pademp wpademp)       //one <
{ int zi;
  for (zi=0;zi<zmaxsize+1;zi++)
    {xpademp.nummp[zi]=wpademp.nummp[zi];
     xpademp.denmp[zi]=wpademp.denmp[zi];
    }
  depthcount=2;
}

wpoly extractnum(int pn,pademp wpademp)
{ int zi,zj;
  for (zi=0;zi<pn+1;zi++)                //Nov24
    {for (zj=0;zj<prec+2;zj++) wpolyout[zi].sv(zj,wpademp.nummp[zi].gv(zj));
    }                             //realmp
  return padepolyoutmp;  
}

wpoly extractden(int pm,pademp wpademp)
{ int zi,zj;
  for (zi=0;zi<pm+1;zi++)                //Nov24
    {for (zj=0;zj<prec+2;zj++) wpolyout[zi].sv(zj,wpademp.denmp[zi].gv(zj));
    }                             //realmp
  return padepolyoutmp;  
}