/*
cwpoly solvpoly(int ynpoly,realmp ywapoly[zmaxsizepad])   //need this
wpoly constructpoly(int,realmp root[zmaxsizepad]);
wpoly multpolys(int ,realmp poly1[zmaxsizepad],int,realmp poly2[zmaxsizepad]);
void operator<(realmp p[zmaxsizepad],wpoly);
*/

extern wpoly constructpoly(int yynpoly,realmp yywwxrot[zmaxsizepad]);
extern wpoly multpolys(int yn1in,realmp ywpoly1in[zmaxsizepad],int yn2in,realmp ywpoly2in[zmaxsizepad]);
extern void operator<(realmp polyleft[zmaxsizepad],wpoly);

realmp ywwxrot[zmaxsizepad],ywwyrot[zmaxsizepad];

void operator<(realmp polyleft[zmaxsizepad],wpoly)
  {int zi,zj;
   for (zi=0;zi<zmaxsize+1;zi++)  
    {for (zj=0;zj<prec+2;zj++) polyleft[zi].sv(zj,wpolyout[zi].gv(zj));
    }
   }

cwpoly solvpoly(int ynpoly,realmp ywapoly[zmaxsizepad])   //need this
  { int zi,zj;
    int yj9,ynnn,yi10;                 //yb is the result
    for (zi=0;zi<ynpoly+1;zi++)
      {for (zj=0;zj<prec+2;zj++)
        {yapoly[zi][zj]=ywapoly[zi].gv(zj);}
      }
    solvpoly(ynpoly,yapoly);

    for (zi=1;zi<ynpasses+1;zi++) //the ywxrot's are used in
      {if (ymess[zi]==2)          //wconstructpoly
        {for (zj=0;zj<prec+2;zj++)
          {ywwxrot[zi].sv(zj,yxrot[zi][zj]);}   //xroot
        }
       if (ymess[zi]==3)
        {for (zj=0;zj<prec+2;zj++)
          {ywwxrot[zi].sv(zj,yxrot[zi][zj]);   //xroot
           ywwyrot[zi].sv(zj,yyrot[zi][zj]);   //yroot
          }
        }
      }                 
     //now, always complexify the roots

   yi10=1;                            //counts
   for (yi=1;yi<ynpasses+1;yi++)      //when yi=1    2
     {if (ymess[yi]==2)
        {for (yj=0;yj<prec+2;yj++)
             {cwpolyout[yi10].sv(0,yj,yxrot[yi][yj]);
              cwpolyout[yi10].sv(1,yj,0);}
         cwpolyout[yi10].sv(1,prec+1,2);
         yi10=yi10+1;
        }
     }
   for (yi=1;yi<ynpasses+1;yi++)
     {if (ymess[yi]==3)
        {for (yj=0;yj<prec+2;yj++)
           {cwpolyout[yi10].sv(0,yj,yxrot[yi][yj]);
            cwpolyout[yi10].sv(1,yj,yyrot[yi][yj]);
            cwpolyout[yi10+1].sv(0,yj,yxrot[yi][yj]);
            cwpolyout[yi10+1].sv(1,yj,yyrot[yi][yj]);}
         cwpolyout[yi10+1].sv(1,prec+1,2-yyrot[yi][prec+1]);
         yi10=yi10+2;
        }
     }
   return cwzpolyout;
  }


wpoly constructpoly(int yynpoly,realmp yywwxrot[zmaxsizepad])
  {int zi,zj;
   for (zi=1;zi<yynpoly+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
        {yxrot[zi][zj]=yywwxrot[zi].gv(zj);}
     }
   ynpoly=yynpoly;
   multiplyout();  //all info still present, just translate output
   for (zi=0;zi<ynpoly+1;zi++)
     {for (zj=0;zj<prec+2;zj++) wpolyout[zi].sv(zj,yb[zi][zj]);}
   return zwpolyout;
  }                                           //done

wpoly multpolys(int yn1in,realmp ywpoly1in[zmaxsizepad],int yn2in,realmp ywpoly2in[zmaxsizepad])
  {int zi,zj;
   for (zi=0;zi<yn1in+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {ypoly1in[zi][zj]=ywpoly1in[zi].gv(zj);}
     }
   for (zi=0;zi<yn2in+1;zi++)
     {for (zj=0;zj<prec+2;zj++)
       {ypoly2in[zi][zj]=ywpoly2in[zi].gv(zj);}
     }
   multpolys(yn1in,ypoly1in,yn2in,ypoly2in);
   for (zi=0;zi<yn1in+yn2in+1;zi++)      //will need oversere.hc
    {for (zj=0;zj<prec+2;zj++) wpolyout[zi].sv(zj,yb[zi][zj]);
    }
   return zwpolyout;
  }                                           //done