extern fakemp csmult(int sn1,complexmp cwsin1[cmaxsizepad],complexmp cwsin2[cmaxsizepad]);
extern fakemp cstonth(int csn3,complexmp cwsin4[cmaxsizepad],complexmp csp);
extern fakemp csder(int csn5,complexmp cwsin5[cmaxsizepad]);
extern complexmp cseval(int csn6,complexmp cwsin6[cmaxsizepad],complexmp cwz77); 
extern void operator<(complexmp cwsstore[cmaxsizepad],fakemp coutputmp);
extern fakemp scontoc(int sn6,realmp swin6[cmaxsizepad]); //moved from
                                                             //oversere.h
extern fakemp csadd(int sn1,complexmp cwsin1[cmaxsizepad],complexmp cwsin2[cmaxsizepad]); 
extern fakemp cssub(int sn1,complexmp cwsin1[cmaxsizepad],complexmp cwsin2[cmaxsizepad]);
extern fakemp csmultfac(int csn2,complexmp cwsin3[cmaxsizepad],complexmp cwsfac);
extern fakemp cslogder(int csn6,complexmp cwsin6[cmaxsizepad]);  //csin6?


complexmp cwsoutmp[cmaxsizepad];
fakemp coutputmp;

fakemp csmult(int sn1,complexmp cwsin1[cmaxsizepad],complexmp cwsin2[cmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {csin1[zi][0][zj]=cwsin1[zi].gv(0,zj);
        csin1[zi][1][zj]=cwsin1[zi].gv(1,zj);
        csin2[zi][0][zj]=cwsin2[zi].gv(0,zj);
        csin2[zi][1][zj]=cwsin2[zi].gv(1,zj);
       }
     }
   csmult(sn1,csin1,csin2);
   for (zi=0;zi<sn1+1;zi++)
    {for (zj=0;zj<prec+2;zj++)
      {cwsoutmp[zi].sv(0,zj,csout[zi][0][zj]);
       cwsoutmp[zi].sv(1,zj,csout[zi][1][zj]);
      }
    }
   return coutputmp;
  }

fakemp cstonth(int csn3,complexmp cwsin4[cmaxsizepad],complexmp cwsp)
  {int zi,zj;
   for (zj=0;zj<prec+2;zj++)
    {csp[0][zj]=cwsp.gv(0,zj);
     csp[1][zj]=cwsp.gv(1,zj);
    }
   for (zi=0;zi<csn3+1;zi++)
    {for (zj=0;zj<prec+2;zj++)
      {csin4[zi][0][zj]=cwsin4[zi].gv(0,zj);
       csin4[zi][1][zj]=cwsin4[zi].gv(1,zj);
      }
    }
   cstonth(csn3,csin4,csp);
   for (zi=0;zi<csn3+1;zi++)
    {for (zj=0;zj<prec+2;zj++)
      {cwsoutmp[zi].sv(0,zj,csout[zi][0][zj]);
       cwsoutmp[zi].sv(1,zj,csout[zi][1][zj]);
      }
    }
   return coutputmp;
  }

fakemp csder(int csn5,complexmp cwsin5[cmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<csn5+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
        {csin5[zi][0][zj]=cwsin5[zi].gv(0,zj);
         csin5[zi][1][zj]=cwsin5[zi].gv(1,zj);
        }
     }
   csder(csn5,csin5);
   for (zi=0;zi<csn5;zi++)            //down by 1
     {for (zj=0;zj<prec+2;zj++)
       {cwsoutmp[zi].sv(0,zj,csout[zi][0][zj]);
        cwsoutmp[zi].sv(1,zj,csout[zi][1][zj]);}
     }
   return coutputmp;
  }

complexmp cseval(int csn6,complexmp cwsin6[cmaxsizepad],complexmp cwz77) 
  {
   int zi,zj;
   for (zj=0;zj<prec+2;zj++)
     {cz77[0][zj]=cwz77.gv(0,zj);
      cz77[1][zj]=cwz77.gv(1,zj);
     }
   for (zi=0;zi<csn6+1;zi++)
    {for (zj=0;zj<prec+2;zj++)
      {csin6[zi][0][zj]=cwsin6[zi].gv(0,zj);
       csin6[zi][1][zj]=cwsin6[zi].gv(1,zj);
      }
    }
   cseval(csn6,csin6,cz77);     //returns cz
   depthcount++;
   for (zj=0;zj<prec+2;zj++)
     {cwreturn[depthcount].sv(0,zj,cz[0][zj]);
      cwreturn[depthcount].sv(1,zj,cz[1][zj]);
     }
   return cwreturn[depthcount];
  }

void operator<(complexmp cwsstore[cmaxsizepad],fakemp coutputmp)
  {int zi,zj;
   for (zi=0;zi<cmaxsize+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {cwsstore[zi].sv(0,zj,cwsoutmp[zi].gv(0,zj));
        cwsstore[zi].sv(1,zj,cwsoutmp[zi].gv(1,zj));}
     }
  }

fakemp scontoc(int sn6,realmp swin6[cmaxsizepad])   //soutc
  {int zi,zj;
   for (zi=0;zi<sn6+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {sin6[zi][zj]=swin6[zi].gv(zj);
       }
     }
   scontoc(sn6,sin6);
   for (zi=0;zi<sn6+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {cwsoutmp[zi].sv(0,zj,soutc[zi][0][zj]);
        cwsoutmp[zi].sv(1,zj,soutc[zi][1][zj]);
       }
     }
   return coutputmp;
  }

fakemp csadd(int sn1,complexmp cwsin1[cmaxsizepad],complexmp cwsin2[cmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {csin1[zi][0][zj]=cwsin1[zi].gv(0,zj);
        csin1[zi][1][zj]=cwsin1[zi].gv(1,zj);
        csin2[zi][0][zj]=cwsin2[zi].gv(0,zj);
        csin2[zi][1][zj]=cwsin2[zi].gv(1,zj);
       }
     }
   csadd(sn1,csin1,csin2);
   for (zi=0;zi<sn1+1;zi++)
    {for (zj=0;zj<prec+2;zj++)
      {cwsoutmp[zi].sv(0,zj,csout[zi][0][zj]);
       cwsoutmp[zi].sv(1,zj,csout[zi][1][zj]);
      }
    }
   return coutputmp;
  }

fakemp cssub(int sn1,complexmp cwsin1[cmaxsizepad],complexmp cwsin2[cmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<sn1+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {csin1[zi][0][zj]=cwsin1[zi].gv(0,zj);
        csin1[zi][1][zj]=cwsin1[zi].gv(1,zj);
        csin2[zi][0][zj]=cwsin2[zi].gv(0,zj);
        csin2[zi][1][zj]=cwsin2[zi].gv(1,zj);
       }
     }
   cssub(sn1,csin1,csin2);
   for (zi=0;zi<sn1+1;zi++)
    {for (zj=0;zj<prec+2;zj++)
      {cwsoutmp[zi].sv(0,zj,csout[zi][0][zj]);
       cwsoutmp[zi].sv(1,zj,csout[zi][1][zj]);
      }
    }
   return coutputmp;
  }

fakemp csmultfac(int csn2,complexmp cwsin3[cmaxsizepad],complexmp cwsfac)
  {int zi,zj;
   for (zj=0;zj<prec+2;zj++)
     {csfac[0][zj]=cwsfac.gv(0,zj);
      csfac[1][zj]=cwsfac.gv(1,zj);
     }
   for (zi=0;zi<csn2+1;zi++)
    {for (zj=0;zj<prec+2;zj++)
      {csin3[zi][0][zj]=cwsin3[zi].gv(0,zj);
       csin3[zi][1][zj]=cwsin3[zi].gv(1,zj);
      }
    }
   csmultfac(csn2,csin3,csfac);
   for (zi=0;zi<csn2+1;zi++)
    {for (zj=0;zj<prec+2;zj++)
      {cwsoutmp[zi].sv(0,zj,csout[zi][0][zj]);
       cwsoutmp[zi].sv(1,zj,csout[zi][1][zj]);
      }
    }
   return coutputmp;
  }

fakemp cslogder(int csn6,complexmp cwsin6[cmaxsizepad])  //csin6?
  {int zi,zj;
   for (zi=0;zi<csn6+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
        {csin6[zi][0][zj]=cwsin6[zi].gv(0,zj);
         csin6[zi][1][zj]=cwsin6[zi].gv(1,zj);
        }
     }
   cslogder(csn6,csin6);
   for (zi=0;zi<csn6;zi++)            //down by 1
     {for (zj=0;zj<prec+2;zj++)
       {cwsoutmp[zi].sv(0,zj,csout[zi][0][zj]);
        cwsoutmp[zi].sv(1,zj,csout[zi][1][zj]);}
     }
   return coutputmp;
  }