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));
}
