/*
smult(int,realmp s1[zmaxsizepad],realmp s2[zmaxsizepad]); //s1*s2
stonth(int,realmp s[zmaxsizepad],realmp p);               //s^p
operator<(realmp sstore[zmaxsizepad],fakemp outputmp);    //overloaded <
sinvert(int,realmp s[zmaxsizepad]);                       //invert s
ssubstitute(int,realmp s1[zmaxsizepad],realmp s2[zmaxsizepad]); //s1(s2)
sadd(int,realmp s1[zmaxsizepad],realmp s2[zmaxsizepad]);  //s1+s2
ssub(int,realmp s1[zmaxsizepad],realmp s2[zmaxsizepad]);  //s1-s2
sder(int,realmp s[zmaxsizepad]);                          //ds(x)/dx
slogder(int,realmp s[zmaxsizepad]);                       //dlogs(x)/dx
smultfac(int,realmp s[zmaxsizepad],realmp fac);           //s*fac
seval(int,realmp s[zmaxsizepad], realmp x);               //s(x)
//could add more such as schs (change sign of a series)
*/

extern fakemp smult(int sn1,realmp swin1[zmaxsizepad],realmp swin2[zmaxsizepad]);
extern fakemp stonth(int sn3,realmp swin4[zmaxsizepad],realmp swp);
extern void operator<(realmp swstore[zmaxsizepad],fakemp outputmp);
extern fakemp sinvert(int sn1,realmp swin1[zmaxsizepad]);            
extern fakemp ssubstitute(int sn1, realmp swin1[zmaxsizepad],realmp swin2[zmaxsizepad]);
extern fakemp sadd(int sn1,realmp swin1[zmaxsizepad],realmp swin2[zmaxsizepad]);
extern fakemp ssub(int sn1,realmp swin1[zmaxsizepad],realmp swin2[zmaxsizepad]);
extern fakemp sder(int sn5,realmp swin5[zmaxsizepad]);
extern fakemp slogder(int sn6,realmp swin6[zmaxsizepad]);
extern fakemp smultfac(int sn2,realmp swin3[zmaxsizepad],realmp sfac);
extern realmp seval(int sn6, realmp swin6[zmaxsizepad], realmp sw77);
//could add more such as schs (change sign of a series)

realmp swoutmp[zmaxsizepad];
fakemp outputmp;

fakemp smult(int sn1,realmp swin1[zmaxsizepad],realmp swin2[zmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {sin1[zi][zj]=swin1[zi].gv(zj);
        sin2[zi][zj]=swin2[zi].gv(zj);
       }
     }
   smult(sn1,sin1,sin2);
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {swoutmp[zi].sv(zj,sout[zi][zj]);}
     }
   return outputmp;
  }

fakemp stonth(int sn3,realmp swin4[zmaxsizepad],realmp swp)
  {int zi,zj;
   for (zj=0;zj<prec+2;zj++) sp[zj]=swp.gv(zj);
   for (zi=0;zi<sn3+1;zi++)
     {for (zj=0;zj<prec+2;zj++) sin4[zi][zj]=swin4[zi].gv(zj);}
   stonth(sn3,sin4,sp);
   for (zi=0;zi<sn3+1;zi++)
     {for (zj=0;zj<prec+2;zj++) swoutmp[zi].sv(zj,sout[zi][zj]);}
   return outputmp;
  }

void operator<(realmp swstore[zmaxsizepad],fakemp outputmp)
  {int zi,zj;
   for (zi=0;zi<zmaxsize+1;zi++)
     {for (zj=0;zj<prec+2;zj++) swstore[zi].sv(zj,swoutmp[zi].gv(zj));}
  }

fakemp sinvert(int sn1,realmp swin1[zmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {sin1[zi][zj]=swin1[zi].gv(zj);
       }
     }
   sinvert(sn1,sin1);
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {swoutmp[zi].sv(zj,sout[zi][zj]);}
     }
   return outputmp;
  }

fakemp ssubstitute(int sn1, realmp swin1[zmaxsizepad], realmp swin2[zmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {sin1[zi][zj]=swin1[zi].gv(zj);
        sin2[zi][zj]=swin2[zi].gv(zj);
       }
     }
   ssubstitute(sn1,sin1,sin2);
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {swoutmp[zi].sv(zj,sout[zi][zj]);}
     }
   return outputmp;
  }

fakemp sadd(int sn1,realmp swin1[zmaxsizepad],realmp swin2[zmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {sin1[zi][zj]=swin1[zi].gv(zj);
        sin2[zi][zj]=swin2[zi].gv(zj);
       }
     }
   sadd(sn1,sin1,sin2);
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {swoutmp[zi].sv(zj,sout[zi][zj]);}
     }
   return outputmp;
  }

fakemp ssub(int sn1,realmp swin1[zmaxsizepad],realmp swin2[zmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {sin1[zi][zj]=swin1[zi].gv(zj);
        sin2[zi][zj]=swin2[zi].gv(zj);
       }
     }
   ssub(sn1,sin1,sin2);
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {swoutmp[zi].sv(zj,sout[zi][zj]);}
     }
   return outputmp;
  }

fakemp sder(int sn5,realmp swin5[zmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<sn5+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {sin5[zi][zj]=swin5[zi].gv(zj);
       }
     }
   sder(sn5,sin5);
   for (zi=0;zi<sn5;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {swoutmp[zi].sv(zj,sout[zi][zj]);}
     }
   return outputmp;
  }

fakemp slogder(int sn6,realmp swin6[zmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<sn6+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {sin6[zi][zj]=swin6[zi].gv(zj);
       }
     }
   slogder(sn6,sin6);
   for (zi=0;zi<sn6;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {swoutmp[zi].sv(zj,sout[zi][zj]);}
     }
   return outputmp;
  }

fakemp smultfac(int sn2,realmp swin3[zmaxsizepad],realmp swfac)
  {int zi,zj;
   for (zj=0;zj<prec+2;zj++) sfac[zj]=swfac.gv(zj);
   for (zi=0;zi<sn2+1;zi++)
     {for (zj=0;zj<prec+2;zj++) sin3[zi][zj]=swin3[zi].gv(zj);}
   smultfac(sn2,sin3,sfac);
   for (zi=0;zi<sn2+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {swoutmp[zi].sv(zj,sout[zi][zj]);}
     }
   return outputmp;
  }

extern realmp seval(int sn6, realmp swin6[zmaxsizepad], realmp sw77)
  {int zi,zj;
   for (zj=0;zj<prec+2;zj++) s77[zj]=sw77.gv(zj);
   for (zi=0;zi<sn6+1;zi++)
     {for (zj=0;zj<prec+2;zj++) sin6[zi][zj]=swin6[zi].gv(zj);}
   seval(sn6,sin6,s77);
   depthcount++;
   for (zj=0;zj<prec+2;zj++) wwreturn[depthcount].sv(zj,z[zj]);
   return wwreturn[depthcount];
  }