59 Value determinant_gcd;
75 assert(n > 0 && denominator != 0);
81 result[0] = denominator;
90 result[0] =
value_mult(denominator,denominator);
105 result[0] = denominator;
106 for(i=2; i <= n && result[1]!=0; i++) {
109 determinant_gcd =
pgcd_slow(result[0],result[1]);
127 determinant_gcd =
pgcd_slow(result[0],result[1]);
161 for(r=1; r <= n; r++) {
166 for(c=1; c <= n; c++) {
175 for(c=1; c <= n; c++) {
179 for(r=1; r <= n; r++) {
#define value_minus(v1, v2)
#define value_notone_p(val)
#define value_division(ref, val)
#define value_product(v, w)
#define value_mult(v, w)
whether the default is protected or not this define makes no sense any more...
Value pgcd_slow(Value, Value)
pgcd.c
#define MATRIX_NB_LINES(matrix)
#define MATRIX_DENOMINATOR(matrix)
int MATRIX_DENONIMATOR(matrix): acces au denominateur global d'une matrice matrix
#define matrix_free(m)
Allocation et desallocation d'une matrice.
#define MATRIX_ELEM(matrix, i, j)
Macros d'acces aux elements d'une matrice.
Pmatrix matrix_new(int m, int n)
package matrix
void matrix_sub_determinant(Pmatrix a, int i, int j, result)
void matrix_sub_determinant(Pmatrix a,int i, int j, int result[]) calculate sub determinant of a matr...
void matrix_determinant(Pmatrix a, result)
package matrix
void matrix_hermite(Pmatrix MAT, Pmatrix P, Pmatrix H, Pmatrix Q, Value *det_p, Value *det_q)
package matrix