/*
Copyright (C) 2011 Sebastian Pancratz
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 "ulong_extras.h"
#include "nmod_vec.h"
#include "nmod_poly.h"
void nmod_poly_sub_series(nmod_poly_t res,
const nmod_poly_t poly1, const nmod_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);
nmod_poly_fit_length(res, max);
_nmod_poly_sub(res->coeffs, poly1->coeffs, len1,
poly2->coeffs, len2, poly1->mod);
_nmod_poly_set_length(res, max);
_nmod_poly_normalise(res);
}