]> www.ginac.de Git - ginac.git/blob - ginac/polynomial/interpolate_padic_uvar.h
Fix power::subs() in some special cases.
[ginac.git] / ginac / polynomial / interpolate_padic_uvar.h
1 /** @file interpolate_padic_uvar.h
2  *
3  *  Utility function for interpolation. */
4
5 /*
6  *  GiNaC Copyright (C) 1999-2022 Johannes Gutenberg University Mainz, Germany
7  *
8  *  This program is free software; you can redistribute it and/or modify
9  *  it under the terms of the GNU General Public License as published by
10  *  the Free Software Foundation; either version 2 of the License, or
11  *  (at your option) any later version.
12  *
13  *  This program is distributed in the hope that it will be useful,
14  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
15  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
16  *  GNU General Public License for more details.
17  *
18  *  You should have received a copy of the GNU General Public License
19  *  along with this program; if not, write to the Free Software
20  *  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
21  */
22
23 #ifndef GINAC_UPOLY_INTERPOLATE_PADIC_H
24 #define GINAC_UPOLY_INTERPOLATE_PADIC_H
25
26 #include "ring_traits.h"
27
28 namespace cln {
29
30 static inline cln::cl_I smod(const cln::cl_I& x, const cln::cl_I y)
31 {
32         cln::cl_I r = mod(x, y);
33         const cln::cl_I y_2 = y >> 1;
34         if (r > y_2)
35                 r = r - y;
36         return r;
37 }
38
39 } // namespace cln
40
41 namespace GiNaC {
42
43 template<typename T> static void
44 interpolate(T& g, const typename T::value_type& gamma,
45             const typename T::value_type& modulus,
46             const std::size_t degree_hint = 1)
47 {
48         typedef typename T::value_type ring_t;
49         g.clear();
50         g.reserve(degree_hint);
51         ring_t e = gamma;
52         while (!zerop(e)) {
53                 const ring_t gi = smod(e, modulus);
54                 g.push_back(gi);
55                 e = exquo(e - gi, modulus);
56         }
57 }
58
59 } // namespace GiNaC
60
61 #endif // ndef GINAC_UPOLY_INTERPOLATE_PADIC_H