1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
/*
Copyright (C) 2010 Sebastian Pancratz
This file is part of FLINT.
FLINT is free software: you can redistribute it and/or modify it under
the terms of the GNU Lesser General Public License (LGPL) as published
by the Free Software Foundation; either version 3 of the License, or
(at your option) any later version. See <https://www.gnu.org/licenses/>.
*/
#include <gmp.h>
#include "ulong_extras.h"
#include "fmpz.h"
void fmpz_sqrtrem(fmpz_t f, fmpz_t r, const fmpz_t g)
{
if (fmpz_sgn(g) < 0)
{
flint_throw(FLINT_ERROR, "Exception (fmpz_sqrtrem). g is negative.\n");
}
if (!COEFF_IS_MPZ(*g))
{
if (COEFF_IS_MPZ(*r))
_fmpz_clear_mpz(*r);
fmpz_set_ui(f, n_sqrtrem((ulong *) r, *g));
}
else
{
mpz_ptr r_mpz_ptr, f_mpz_ptr;
_fmpz_promote(f); /* must not hang on to pointer whilst promoting */
r_mpz_ptr = _fmpz_promote(r);
f_mpz_ptr = COEFF_TO_PTR(*f);
mpz_sqrtrem(f_mpz_ptr, r_mpz_ptr, COEFF_TO_PTR(*g));
_fmpz_demote_val(f);
_fmpz_demote_val(r);
}
}