flint-sys 0.9.0

Bindings to the FLINT C library
Documentation 
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
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123

/*
    Copyright (C) 2023 Albin Ahlbäck

    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 <string.h>
#include "test_helpers.h"
#include "fmpz.h"
#include "fmpz_factor.h"

#define TMP_FILENAME "tmp"
#define STRLEN(x) (sizeof(x) - sizeof(char))
#define RESULT_STRING   \
    "0\n"               \
    "1\n"               \
    "-1\n"              \
    "3 * 5 * 7\n"       \
    "-1 * 3^2 * 5^3 * 7"

TEST_FUNCTION_START(fmpz_factor_fprint, state)
{
    FILE * fs;
    fmpz_t a;
    fmpz_factor_t f;
    int res, tmp;
    char str[sizeof(RESULT_STRING)];

    fs = fopen(TMP_FILENAME, "w+");
    if (fs == NULL)
        flint_throw(FLINT_TEST_FAIL, "Could not open temporary file \"" TMP_FILENAME "\"\n");

    fmpz_init(a);
    fmpz_factor_init(f);
    res = 0;

    /* Print 0 */
    fmpz_zero(a);
    fmpz_factor(f, a);
    tmp = fmpz_factor_fprint(fs, f);
    if (tmp != STRLEN("0"))
        goto wrong_return_value;
    res += tmp;
    res += (fputc('\n', fs) != EOF);

    /* Print 1 */
    fmpz_one(a);
    fmpz_factor(f, a);
    tmp = fmpz_factor_fprint(fs, f);
    if (tmp != STRLEN("1"))
        goto wrong_return_value;
    res += tmp;
    res += (fputc('\n', fs) != EOF);

    /* Print -1 */
    fmpz_set_si(a, -1);
    fmpz_factor(f, a);
    tmp = fmpz_factor_fprint(fs, f);
    if (tmp != STRLEN("-1"))
        goto wrong_return_value;
    res += tmp;
    res += (fputc('\n', fs) != EOF);

    /* Print 3 * 5 * 7 */
    fmpz_set_si(a, 3 * 5 * 7);
    fmpz_factor(f, a);
    tmp = fmpz_factor_fprint(fs, f);
    if (tmp != STRLEN("3 * 5 * 7"))
        goto wrong_return_value;
    res += tmp;
    res += (fputc('\n', fs) != EOF);

    /* Print -1 * 3^2 * 5^3 * 7 */
    fmpz_set_si(a, -1 * 9 * 125 * 7);
    fmpz_factor(f, a);
    tmp = fmpz_factor_fprint(fs, f);
    if (tmp != STRLEN("-1 * 3^2 * 5^3 * 7"))
    {
wrong_return_value:
        fseek(fs, -tmp, SEEK_CUR);
        fread(str, sizeof(char), tmp, fs);
        str[tmp] = '\0';
        flint_throw(FLINT_TEST_FAIL,
                "Wrong return value for a = %{fmpz}\n"
                "Got string: %s\n",
                a, str);
    }
    res += tmp;
    res += (fputc('\0', fs) != EOF);

    fseek(fs, 0, SEEK_SET);
    if (fread(str, sizeof(char), res, fs) != res)
        flint_throw(FLINT_TEST_FAIL, "Could not read %d bytes from filestream.\n", res);

    if (strcmp(str, RESULT_STRING))
        flint_throw(FLINT_TEST_FAIL,
                "Result differed from expected output.\n"
                "Got:\n"
                "%s\n\n"
                "Expected:\n"
                "%s\n",
                str, RESULT_STRING);

    fmpz_clear(a);
    fmpz_factor_clear(f);

    if (fclose(fs))
        flint_throw(FLINT_TEST_FAIL, "Could not close temporary file \"" TMP_FILENAME "\"\n");

    if (remove(TMP_FILENAME))
        flint_throw(FLINT_TEST_FAIL, "Could not remove temporary file \"" TMP_FILENAME "\"\n");

    TEST_FUNCTION_END(state);
}

#undef TMP_FILENAME
#undef STRLEN
#undef RESULT_STRING