15 template class std::vector< amp::mpfr_record* >;
30 const unsigned int Precision=300;
68 for(i=1; i<=max_i; i++)
70 for(j=1; j<=max_j; j++)
73 Print(
" to svd:%d,%d=%s\n",i,j,str);
81 svd::svddecomposition<Precision>(a, max_i, max_j, 2, 2, 2,
w, u, vt);
83 Mu=
mpNew(max_i,max_i);
84 Mw=
mpNew(max_i,max_j);
85 Mvt=
mpNew(max_j,max_j);
94 for(i=1;i<=
si_min(max_i,max_j);i++)
101 for(j=1;j<=max_j;j++)
112 L->
m[0].
data=(
char*)Mu;
113 L->
m[1].
data=(
char*)Mw;
114 L->
m[2].
data=(
char*)Mvt;
const CanonicalForm int s
static FORCE_INLINE const char * n_Read(const char *s, number *a, const coeffs r)
!!! Recommendation: This method is too cryptic to be part of the user- !!! interface. As defined here, it is merely a helper !!! method for parsing number input strings.
static int si_min(const int a, const int b)
Compatiblity layer for legacy polynomial operations (over currRing)
static const ampf getAlgoPascalEpsilon()
void setbounds(int iLow1, int iHigh1, int iLow2, int iHigh2)
static number & pGetCoeff(poly p)
return an alias to the leading coefficient of p assumes that p != NULL NOTE: not copy ...
INLINE_THIS void Init(int l=0)
matrix mpNew(int r, int c)
create a r x c zero-matrix
#define pInit()
allocates a new monomial and initializes everything to 0
ring currRing
Widely used global variable which specifies the current polynomial ring for Singular interpreter and ...
std::string toString(const gfan::ZCone *const c)
#define MATELEM(mat, i, j)