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
41
42
43
|
/*
Copyright (C) 2021 Fredrik Johansson
This file is part of Calcium.
Calcium 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 2.1 of the License, or
(at your option) any later version. See <http://www.gnu.org/licenses/>.
*/
#include "fexpr.h"
int
fexpr_expanded_normal_form(fexpr_t res, const fexpr_t expr, ulong flags)
{
fexpr_vec_t args;
fmpz_mpoly_ctx_t ctx;
fmpz_mpoly_q_t frac;
int success;
fexpr_vec_init(args, 0);
fexpr_arithmetic_nodes(args, expr);
_fexpr_vec_sort_fast(args->entries, args->length);
/* todo: when length == 0, use fmpq arithmetic instead */
fmpz_mpoly_ctx_init(ctx, FLINT_MAX(args->length, 1), ORD_LEX);
fmpz_mpoly_q_init(frac, ctx);
success = fexpr_get_fmpz_mpoly_q(frac, expr, args, ctx);
if (success)
fexpr_set_fmpz_mpoly_q(res, frac, args, ctx);
else
fexpr_set(res, expr);
fmpz_mpoly_q_clear(frac, ctx);
fmpz_mpoly_ctx_clear(ctx);
fexpr_vec_clear(args);
return success;
}
|
