extern wpoly solvlin(unsigned int wnin,realmp wwmin[zmaxsizepad][zmaxsizepad],realmp wrin[zmaxsizepad]);
extern int wgenlin(unsigned int wnin,unsigned int * * *wmin,unsigned int * *wrin);

int wnin;
int wninmax;

int wn;
int zrr;
int precp2=prec+2;

realmatrix wm_ram;unsigned int * * *wm=wm_ram.m;
realvector wr_ram;unsigned int * * wr=wr_ram.v;
realvector wx_ram;unsigned int * * wx=wx_ram.v;
realscalar wrecip_ram;unsigned int * wrecip=wrecip_ram.s;
wpoly zwpolyout;
realmp wpolyout[zmaxsizepad];

int wnp1;

wpoly solvlin(unsigned int wnin,realmp wwmin[zmaxsizepad][zmaxsizepad],realmp wwrin[zmaxsizepad])
{ int wi,wj,wk,wl;
  wn=wnin;
  wnp1=wn+1;
  for (wk=0;wk<precp2;wk++)
    {




     for (wi=1;wi<wnp1;wi++)
      {wr[wi][wk]=wwrin[wi].gv(wk);
         for (wj=1;wj<wnp1;wj++)
           {wm[wi][wj][wk]=wwmin[wi][wj].gv(wk);}
      }
   }
  for (wk=2;wk<wnp1;wk++)
     {
      wl=wk-1;
      zconinttomp(1);
      zdiv(z,wm[wl][wl]);
      zstore(wrecip);
      for (wi=wk;wi<wnp1;wi++)
        {
         zmult(wm[wi][wl],wr[wl]);
         zdiv(z,wm[wl][wl]);   //zmult(z,wrecip);
         zsub(wr[wi],z);
	 zstore(wr[wi]);
         for (wj=wk;wj<wnp1;wj++)
	   {zmult(wm[wi][wl],wm[wl][wj]);
            zdiv(z,wm[wl][wl]);   //zmult(z,wrecip);
            zsub(wm[wi][wj],z);     
	    zstore(wm[wi][wj]);        
	   }
        }          //end of wi loop
     }       //end of wkloop
  zdiv(wr[wn],wm[wn][wn]);
  zstore(wx[wn]);
  for (wi=2;wi<wnp1;wi++)
     {wk=wn-wi+1;
      zrecall(wr[wk]);
      zstore(wx[wk]);
      for (wj=1;wj<wi;wj++)
        {wl=wn-wj+1;
	 zmult(wm[wk][wl],wx[wl]);
         zsub(wx[wk],z);
	 zstore(wx[wk]);
	}
      zdiv(wx[wk],wm[wk][wk]);
      zstore(wx[wk]);
     }
  for (wk=0;wk<precp2;wk++) 
    {for (wi=1;wi<wnp1;wi++) wpolyout[wi].sv(wk,wx[wi][wk]);}
  return zwpolyout;
}
   
int wgenlin(unsigned int wnin,unsigned int * * *wmin,unsigned int * *wrin)
{ int wi,wj,wk,wl,wi2,zi;
  int wik,wjk;
  wn=wnin;
  wnp1=wn+1;

for (wk=0;wk<precp2;wk++) 
 {for (wi=1;wi<wnp1;wi++)
    {wr[wi][wk]=wrin[wi][wk];
       for (wj=1;wj<wnp1;wj++)
	 {wm[wi][wj][wk]=wmin[wi][wj][wk];}
    }
 }

  for (wk=2;wk<wnp1;wk++)
     {wl=wk-1;
      zconinttomp(1);
      zdiv(z,wm[wl][wl]);
      zstore(wrecip);
      for (wi=wk;wi<wnp1;wi++)
        {zmult(wm[wi][wl],wr[wl]);
	 zdiv(z,wm[wl][wl]);   //zmult(z,wrecip);
         zsub(wr[wi],z);
	 zstore(wr[wi]);
         for (wj=wk;wj<wnp1;wj++)
	   {zmult(wm[wi][wl],wm[wl][wj]);
            zdiv(z,wm[wl][wl]);   //zmult(z,wrecip);
            zsub(wm[wi][wj],z);     
	    zstore(wm[wi][wj]);        
	   }
        }          //end of wi loop
     }       //end of wkloop
  zdiv(wr[wn],wm[wn][wn]);
  zstore(wx[wn]);
  for (wi=2;wi<wnp1;wi++)
     {wk=wn-wi+1;
      zrecall(wr[wk]);
      zstore(wx[wk]);
      for (wj=1;wj<wi;wj++)
        {wl=wn-wj+1;
	 zmult(wm[wk][wl],wx[wl]);
         zsub(wx[wk],z);
	 zstore(wx[wk]);
	}
      zdiv(wx[wk],wm[wk][wk]);
      zstore(wx[wk]);
     }
  return 0;
}