30 gaussElem(
const fglmVector newv,
const fglmVector newp, number & newpdenom, number & newfac ) : v( newv ), p( newp ), pdenom( newpdenom ), fac( newfac )
36 #ifndef HAVE_EXPLICIT_CONSTR 39 void mac_gaussElem(
const fglmVector newv,
const fglmVector newp, number & newpdenom, number & newfac )
62 #ifndef HAVE_EXPLICIT_CONSTR 68 for ( k=
max; k > 0; k-- )
75 #ifndef HAVE_EXPLICIT_CONSTR 79 for ( k=
size; k > 0; k-- )
115 for ( k= 1; k <=
size; k++ ) {
121 temp=
nMult( fac2, pdenom );
125 temp=
nMult( pdenom, elems[k].pdenom );
131 number gcd = v.gcd();
135 number temp =
nMult( pdenom, gcd );
147 temp=
nDiv( pdenom, gcd );
184 isPivot[ pivotcol ]=
TRUE;
185 perm[
size]= pivotcol;
188 #ifndef HAVE_EXPLICIT_CONSTR
number getconstelem(int i) const
#define omFreeSize(addr, size)
bool pivot(const matrix aMat, const int r1, const int r2, const int c1, const int c2, int *bestR, int *bestC, const ring R)
This code computes a score for each non-zero matrix entry in aMat[r1..r2, c1..c2].
gaussElem(const fglmVector newv, const fglmVector newp, number &newpdenom, number &newfac)
fglmVector getDependence()
void nihilate(const number fac1, const number fac2, const fglmVector v)
static int max(int a, int b)
BOOLEAN reduce(fglmVector v)
void setelem(int i, number &n)
ring currRing
Widely used global variable which specifies the current polynomial ring for Singular interpreter and ...
static FORCE_INLINE number n_SubringGcd(number a, number b, const coeffs r)