void nf_elem_scalar_div_fmpq(nf_elem_t a, const nf_elem_t b, const fmpq_t c, const nf_t nf) { if (nf->flag & NF_LINEAR) { fmpz * den = LNF_ELEM_DENREF(a); fmpz * num = LNF_ELEM_NUMREF(a); const fmpz * const den2 = LNF_ELEM_DENREF(b); const fmpz * const num2 = LNF_ELEM_NUMREF(b); _fmpq_mul(num, den, num2, den2, fmpq_denref(c), fmpq_numref(c)); } else if (nf->flag & NF_QUADRATIC) { fmpz * den = QNF_ELEM_DENREF(a); fmpz * num = QNF_ELEM_NUMREF(a); const fmpz * const den2 = QNF_ELEM_DENREF(b); const fmpz * const num2 = QNF_ELEM_NUMREF(b); _fmpq_poly_scalar_div_fmpq(num, den, num2, den2, 2, fmpq_numref(c), fmpq_denref(c)); } else { fmpq_poly_scalar_div_fmpq(NF_ELEM(a), NF_ELEM(b), c); } }
void _fmpq_submul(fmpz_t rnum, fmpz_t rden, const fmpz_t op1num, const fmpz_t op1den, const fmpz_t op2num, const fmpz_t op2den) { fmpz_t u, v; fmpz_init(u); fmpz_init(v); _fmpq_mul(u, v, op1num, op1den, op2num, op2den); _fmpq_sub(rnum, rden, rnum, rden, u, v); fmpz_clear(u); fmpz_clear(v); }
void fmpq_div_fmpz(fmpq_t res, const fmpq_t op, const fmpz_t x) { fmpz_t y; fmpz_init(y); fmpz_one(y); _fmpq_mul(fmpq_numref(res), fmpq_denref(res), fmpq_numref(op), fmpq_denref(op), y, x); fmpz_clear(y); if (fmpz_sgn(fmpq_denref(res)) < 0) { fmpz_neg(fmpq_numref(res), fmpq_numref(res)); fmpz_neg(fmpq_denref(res), fmpq_denref(res)); } }
void _fmpq_div(fmpz_t rnum, fmpz_t rden, const fmpz_t op1num, const fmpz_t op1den, const fmpz_t op2num, const fmpz_t op2den) { fmpz_t t, u; fmpz_init(t); fmpz_init(u); fmpz_set(t, op2den); fmpz_set(u, op2num); _fmpq_mul(rnum, rden, op1num, op1den, t, u); fmpz_clear(t); fmpz_clear(u); if (fmpz_sgn(rden) < 0) { fmpz_neg(rnum, rnum); fmpz_neg(rden, rden); } }
void fmpq_mul(fmpq_t res, const fmpq_t op1, const fmpq_t op2) { _fmpq_mul(fmpq_numref(res), fmpq_denref(res), fmpq_numref(op1), fmpq_denref(op1), fmpq_numref(op2), fmpq_denref(op2)); }