15template 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);
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;
static int si_min(const int a, const int b)
std::string toString(const gfan::ZCone *const c)
static const ampf getAlgoPascalEpsilon()
void setbounds(int iLow1, int iHigh1, int iLow2, int iHigh2)
INLINE_THIS void Init(int l=0)
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....
const CanonicalForm int s
matrix mpNew(int r, int c)
create a r x c zero-matrix
#define MATELEM(mat, i, j)
1-based access to matrix
static number & pGetCoeff(poly p)
return an alias to the leading coefficient of p assumes that p != NULL NOTE: not copy
VAR ring currRing
Widely used global variable which specifies the current polynomial ring for Singular interpreter and ...
Compatibility layer for legacy polynomial operations (over currRing)
#define pInit()
allocates a new monomial and initializes everything to 0