/* Copyright (C) 2008, 2009, 2014 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 2.1 of the License, or (at your option) any later version. See . */ #include #include "flint.h" #include "fmpz.h" #include "fmpz_poly.h" void fmpz_poly_sub_series(fmpz_poly_t res, const fmpz_poly_t poly1, const fmpz_poly_t poly2, slong n) { slong len1, len2, max = FLINT_MAX(poly1->length, poly2->length); if (n < 0) n = 0; max = FLINT_MIN(max, n); len1 = FLINT_MIN(poly1->length, max); len2 = FLINT_MIN(poly2->length, max); fmpz_poly_fit_length(res, max); _fmpz_poly_sub(res->coeffs, poly1->coeffs, len1, poly2->coeffs, len2); _fmpz_poly_set_length(res, max); _fmpz_poly_normalise(res); }