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