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
/*
Copyright (C) 2008, 2009 William Hart
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 "nmod_poly.h"
int
nmod_poly_equal_trunc(const nmod_poly_t poly1, const nmod_poly_t poly2, slong n)
{
slong i, len1, len2;
if (poly1 == poly2)
return 1; /* same polynomial */
if (n < 0)
n = 0;
len1 = FLINT_MIN(poly1->length, n);
len2 = FLINT_MIN(poly2->length, n);
if (len1 < len2)
{
for (i = len1; i < len2; i++)
{
if (poly2->coeffs[i] != 0)
return 0;
}
} else if (len2 < len1)
{
for (i = len2; i < len1; i++)
{
if (poly1->coeffs[i] != 0)
return 0;
}
}
for (i = 0; i < FLINT_MIN(len1, len2); i++) /* check if coefficients the same */
if (poly1->coeffs[i] != poly2->coeffs[i])
return 0;
return 1;
}