/* Copyright (C) 2015 Vladimir Glazachev 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 "aprcl.h" void unity_zp_coeff_add_fmpz(unity_zp f, ulong ind, const fmpz_t x) { fmpz_t coeff; fmpz_init(coeff); fmpz_mod_poly_get_coeff_fmpz(coeff, f->poly, ind, f->ctx); if (fmpz_is_zero(coeff)) { fmpz_clear(coeff); fmpz_mod_poly_set_coeff_fmpz(f->poly, ind, x, f->ctx); return; } fmpz_clear(coeff); fmpz_add(f->poly->coeffs + ind, f->poly->coeffs + ind, x); if (fmpz_cmp(f->poly->coeffs + ind, fmpz_mod_ctx_modulus(f->ctx)) >= 0) fmpz_sub(f->poly->coeffs + ind, f->poly->coeffs + ind, fmpz_mod_ctx_modulus(f->ctx)); } void unity_zp_coeff_add_ui(unity_zp f, ulong ind, ulong x) { fmpz_t coeff; fmpz_init(coeff); fmpz_mod_poly_get_coeff_fmpz(coeff, f->poly, ind, f->ctx); if (fmpz_is_zero(coeff)) { fmpz_clear(coeff); fmpz_mod_poly_set_coeff_ui(f->poly, ind, x, f->ctx); return; } fmpz_clear(coeff); fmpz_add_ui(f->poly->coeffs + ind, f->poly->coeffs + ind, x); if (fmpz_cmp(f->poly->coeffs + ind, fmpz_mod_ctx_modulus(f->ctx)) >= 0) fmpz_sub(f->poly->coeffs + ind, f->poly->coeffs + ind, fmpz_mod_ctx_modulus(f->ctx)); }