15 #if __FLINT_RELEASE >= 20500 26 int convFlintISingI (fmpz_t
f)
33 void convSingIFlintI(fmpz_t
f,
int p)
40 void convFlintNSingN (mpz_t z, fmpz_t
f)
46 void convSingNFlintN(fmpz_t
f, mpz_t z)
59 fmpz_mat_init(M, r, c);
62 fmpz_mat_init(Transf, T->
rows(), T->
rows());
72 convSingNFlintN(dummy,n);
74 fmpz_set(fmpz_mat_entry(M, i-1, j-1), dummy);
80 for(i=T->
rows();i>0;i--)
82 for(j=T->
rows();j>0;j--)
85 convSingNFlintN(dummy,n);
87 fmpz_set(fmpz_mat_entry(Transf, i-1, j-1), dummy);
93 fmpz_lll_context_init_default(fl);
95 fmpz_lll(M, Transf, fl);
97 fmpz_lll(M,
NULL, fl);
102 convFlintNSingN(n, fmpz_mat_entry(M, i-1, j-1));
110 for(i=T->
rows();i>0;i--)
112 for(j=T->
cols();j>0;j--)
114 convFlintNSingN(n, fmpz_mat_entry(Transf, i-1, j-1));
130 fmpz_mat_init(M, r, c);
132 fmpz_mat_init(Transf, r, r);
139 convSingIFlintI(dummy,
IMATELEM(*m,i,j));
140 fmpz_set(fmpz_mat_entry(M, i-1, j-1), dummy);
146 for(i=T->
rows();i>0;i--)
148 for(j=T->
rows();j>0;j--)
150 convSingIFlintI(dummy,
IMATELEM(*T,i,j));
151 fmpz_set(fmpz_mat_entry(Transf, i-1, j-1), dummy);
157 fmpz_lll_context_init_default(fl);
159 fmpz_lll(M, Transf, fl);
161 fmpz_lll(M,
NULL, fl);
166 IMATELEM(*res,i,j)=convFlintISingI(fmpz_mat_entry(M, i-1, j-1));
171 for(i=Transf->r;i>0;i--)
173 for(j=Transf->r;j>0;j--)
175 IMATELEM(*T,i,j)=convFlintISingI(fmpz_mat_entry(Transf, i-1, j-1));
This file is work in progress and currently not part of the official Singular.
Coefficient rings, fields and other domains suitable for Singular polynomials.
static FORCE_INLINE number n_InitMPZ(mpz_t n, const coeffs r)
conversion of a GMP integer to number
#define BIMATELEM(M, I, J)
coeffs basecoeffs() const
static FORCE_INLINE void n_Delete(number *p, const coeffs r)
delete 'p'
#define IMATELEM(M, I, J)
static FORCE_INLINE void n_MPZ(mpz_t result, number &n, const coeffs r)
conversion of n to a GMP integer; 0 if not possible