forked from tthsqe12/flint2
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathget_str.c
58 lines (48 loc) · 1.37 KB
/
get_str.c
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
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
/*
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 2.1 of the License, or
(at your option) any later version. See <http://www.gnu.org/licenses/>.
*/
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <gmp.h>
#include "flint.h"
#include "fmpz.h"
#include "fmpz_poly.h"
char *
_fmpz_poly_get_str(const fmpz * poly, slong len)
{
slong i, bound;
char *str, *strbase;
if (len == 0)
{
str = (char *) flint_malloc(2 * sizeof(char));
str[0] = '0';
str[1] = '\0';
return str;
}
bound = (slong) (ceil(log10((double) (len + 1))));
for (i = 0; i < len; i++)
bound += fmpz_sizeinbase(poly + i, 10) + 1;
bound += len + 2;
strbase = (char *) flint_malloc(bound * sizeof(char));
str = strbase;
str += flint_sprintf(str, "%wd ", len);
do
{
if (!COEFF_IS_MPZ(*poly))
str += flint_sprintf(str, " %wd", *poly);
else
str += gmp_sprintf(str, " %Zd", COEFF_TO_PTR(*poly));
} while (poly++, --len);
return strbase;
}
char *
fmpz_poly_get_str(const fmpz_poly_t poly)
{
return _fmpz_poly_get_str(poly->coeffs, poly->length);
}