47 template<
unsigned int Precision>
53 template<
unsigned int Precision>
60 template<
unsigned int Precision>
70 template<
unsigned int Precision>
77 template<
unsigned int Precision>
87 template<
unsigned int Precision>
94 template<
unsigned int Precision>
100 template<
unsigned int Precision>
107 template<
unsigned int Precision>
117 template<
unsigned int Precision>
124 template<
unsigned int Precision>
134 template<
unsigned int Precision>
193 template<
unsigned int Precision>
223 tauq.setbounds(0, n-1);
224 taup.setbounds(0, n-1);
228 tauq.setbounds(0,
m-1);
229 taup.setbounds(0,
m-1);
237 for(
i=0;
i<=n-1;
i++)
283 for(
i=0;
i<=
m-1;
i++)
347 template<
unsigned int Precision>
361 if(
m==0 || n==0 || qcolumns==0 )
369 q.setbounds(0,
m-1, 0, qcolumns-1);
370 for(
i=0;
i<=
m-1;
i++)
372 for(
j=0;
j<=qcolumns-1;
j++)
421 template<
unsigned int Precision>
441 if(
m<=0 || n<=0 || zrows<=0 || zcolumns<=0 )
499 while(
i!=i2+istep );
547 while(
i!=i2+istep );
574 template<
unsigned int Precision>
588 if(
m==0 || n==0 || ptrows==0 )
596 pt.setbounds(0, ptrows-1, 0, n-1);
597 for(
i=0;
i<=ptrows-1;
i++)
599 for(
j=0;
j<=n-1;
j++)
648 template<
unsigned int Precision>
668 if(
m<=0 || n<=0 || zrows<=0 || zcolumns<=0 )
730 while(
i!=i2+istep );
777 while(
i!=i2+istep );
804 template<
unsigned int Precision>
824 for(
i=0;
i<=n-2;
i++)
835 for(
i=0;
i<=
m-2;
i++)
849 template<
unsigned int Precision>
873 taup.setbounds(1, minmn);
874 tauq.setbounds(1, minmn);
980 template<
unsigned int Precision>
997 if(
m==0 || n==0 || qcolumns==0 )
1005 q.setbounds(1,
m, 1, qcolumns);
1014 for(
j=1;
j<=qcolumns;
j++)
1042 ap::vmove(
v.getvector(1, vm), qp.getcolumn(
i, ip1,
m));
1054 template<
unsigned int Precision>
1076 if(
m<=0 || n<=0 || zrows<=0 || zcolumns<=0 )
1135 while(
i!=i2+istep );
1173 ap::vmove(
v.getvector(1, vm), qp.getcolumn(
i, ip1,
m));
1185 while(
i!=i2+istep );
1195 template<
unsigned int Precision>
1212 if(
m==0 || n==0 || ptrows==0 )
1220 pt.setbounds(1, ptrows, 1, n);
1227 for(
i=1;
i<=ptrows;
i++)
1247 ap::vmove(
v.getvector(1, vm), qp.getrow(
i, ip1, n));
1269 template<
unsigned int Precision>
1291 if(
m<=0 || n<=0 || zrows<=0 || zcolumns<=0 )
1343 ap::vmove(
v.getvector(1, vm), qp.getrow(
i, ip1, n));
1355 while(
i!=i2+istep );
1403 while(
i!=i2+istep );
1412 template<
unsigned int Precision>
1432 for(
i=1;
i<=n-1;
i++)
1443 for(
i=1;
i<=
m-1;
i++)
static void make_assertion(bool bClause)
raw_vector< T > getvector(int iStart, int iEnd)
void setbounds(int iLow, int iHigh)
const Variable & v
< [in] a sqrfree bivariate poly
int maxint(int m1, int m2)
void vmove(raw_vector< T > vdst, const_raw_vector< T > vsrc)
int minint(int m1, int m2)
void unpackptfrombidiagonal(const ap::template_2d_array< amp::ampf< Precision > > &qp, int m, int n, const ap::template_1d_array< amp::ampf< Precision > > &taup, int ptrows, ap::template_2d_array< amp::ampf< Precision > > &pt)
void multiplybyqfrombidiagonal(const ap::template_2d_array< amp::ampf< Precision > > &qp, int m, int n, const ap::template_1d_array< amp::ampf< Precision > > &tauq, ap::template_2d_array< amp::ampf< Precision > > &z, int zrows, int zcolumns, bool fromtheright, bool dotranspose)
void rmatrixbdunpackpt(const ap::template_2d_array< amp::ampf< Precision > > &qp, int m, int n, const ap::template_1d_array< amp::ampf< Precision > > &taup, int ptrows, ap::template_2d_array< amp::ampf< Precision > > &pt)
void unpackqfrombidiagonal(const ap::template_2d_array< amp::ampf< Precision > > &qp, int m, int n, const ap::template_1d_array< amp::ampf< Precision > > &tauq, int qcolumns, ap::template_2d_array< amp::ampf< Precision > > &q)
void rmatrixbdunpackq(const ap::template_2d_array< amp::ampf< Precision > > &qp, int m, int n, const ap::template_1d_array< amp::ampf< Precision > > &tauq, int qcolumns, ap::template_2d_array< amp::ampf< Precision > > &q)
void rmatrixbdmultiplybyq(const ap::template_2d_array< amp::ampf< Precision > > &qp, int m, int n, const ap::template_1d_array< amp::ampf< Precision > > &tauq, ap::template_2d_array< amp::ampf< Precision > > &z, int zrows, int zcolumns, bool fromtheright, bool dotranspose)
void unpackdiagonalsfrombidiagonal(const ap::template_2d_array< amp::ampf< Precision > > &b, int m, int n, bool &isupper, ap::template_1d_array< amp::ampf< Precision > > &d, ap::template_1d_array< amp::ampf< Precision > > &e)
void rmatrixbdunpackdiagonals(const ap::template_2d_array< amp::ampf< Precision > > &b, int m, int n, bool &isupper, ap::template_1d_array< amp::ampf< Precision > > &d, ap::template_1d_array< amp::ampf< Precision > > &e)
void tobidiagonal(ap::template_2d_array< amp::ampf< Precision > > &a, int m, int n, ap::template_1d_array< amp::ampf< Precision > > &tauq, ap::template_1d_array< amp::ampf< Precision > > &taup)
void rmatrixbd(ap::template_2d_array< amp::ampf< Precision > > &a, int m, int n, ap::template_1d_array< amp::ampf< Precision > > &tauq, ap::template_1d_array< amp::ampf< Precision > > &taup)
void rmatrixbdmultiplybyp(const ap::template_2d_array< amp::ampf< Precision > > &qp, int m, int n, const ap::template_1d_array< amp::ampf< Precision > > &taup, ap::template_2d_array< amp::ampf< Precision > > &z, int zrows, int zcolumns, bool fromtheright, bool dotranspose)
void multiplybypfrombidiagonal(const ap::template_2d_array< amp::ampf< Precision > > &qp, int m, int n, const ap::template_1d_array< amp::ampf< Precision > > &taup, ap::template_2d_array< amp::ampf< Precision > > &z, int zrows, int zcolumns, bool fromtheright, bool dotranspose)
void generatereflection(ap::template_1d_array< amp::ampf< Precision > > &x, int n, amp::ampf< Precision > &tau)
void applyreflectionfromtheright(ap::template_2d_array< amp::ampf< Precision > > &c, amp::ampf< Precision > tau, const ap::template_1d_array< amp::ampf< Precision > > &v, int m1, int m2, int n1, int n2, ap::template_1d_array< amp::ampf< Precision > > &work)
void applyreflectionfromtheleft(ap::template_2d_array< amp::ampf< Precision > > &c, amp::ampf< Precision > tau, const ap::template_1d_array< amp::ampf< Precision > > &v, int m1, int m2, int n1, int n2, ap::template_1d_array< amp::ampf< Precision > > &work)