/usr/include/pynac/pseries.h is in libpynac-dev 0.7.12-2build1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 | /** @file pseries.h
*
* Interface to class for extended truncated power series. */
/*
* GiNaC Copyright (C) 1999-2008 Johannes Gutenberg University Mainz, Germany
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the Free Software
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
*/
#ifndef __GINAC_SERIES_H__
#define __GINAC_SERIES_H__
#include "basic.h"
#include "expairseq.h"
namespace GiNaC {
/** This class holds a extended truncated power series (positive and negative
* integer powers). It consists of expression coefficients (only non-zero
* coefficients are stored), an expansion variable and an expansion point.
* Other classes must provide members to convert into this type. */
class pseries : public basic
{
GINAC_DECLARE_REGISTERED_CLASS(pseries, basic)
friend class ex;
// other constructors
public:
pseries(const ex &rel_, epvector ops_);
// functions overriding virtual functions from base classes
public:
unsigned precedence() const override {return 38;} // for clarity just below add::precedence
size_t nops() const override;
const ex op(size_t i) const override;
numeric degree(const ex &s) const override;
numeric ldegree(const ex &s) const override;
ex coeff(const ex &s, const ex & n) const override;
ex collect(const ex &s, bool distributed = false) const override;
ex eval(int level=0) const override;
ex evalf(int level=0, PyObject* parent=nullptr) const override;
ex series(const relational & r, int order, unsigned options = 0) const override;
ex subs(const exmap & m, unsigned options = 0) const override;
ex normal(exmap & repl, exmap & rev_lookup, int level = 0, unsigned options = 0) const override;
ex expand(unsigned options = 0) const override;
ex conjugate() const override;
ex real_part() const override;
ex imag_part() const override;
protected:
ex derivative(const symbol & s) const override;
// non-virtual functions in this class
public:
/** Get the expansion variable. */
ex get_var() const {return var;}
/** Get the expansion point. */
ex get_point() const {return point;}
/** Convert the pseries object to an ordinary polynomial.
*
* @param no_order flag: discard higher order terms */
ex convert_to_poly(bool no_order = false) const;
/** Check whether series is compatible to another series (expansion
* variable and point are the same. */
bool is_compatible_to(const pseries &other) const {return var.is_equal(other.var) && point.is_equal(other.point);}
/** Check whether series has the value zero. */
bool is_zero() const {return seq.size() == 0;}
/** Returns true if there is no order term, i.e. the series terminates and
* false otherwise. */
bool is_terminating() const;
/** Get coefficients and exponents. */
ex coeffop(size_t i) const;
ex exponop(size_t i) const;
ex add_series(const pseries &other) const;
ex mul_const(const numeric &other) const;
ex mul_series(const pseries &other) const;
ex power_const(const numeric &p, int deg) const;
pseries shift_exponents(int deg) const;
protected:
void print_series(const print_context & c, const char *openbrace, const char *closebrace, const char *mul_sym, const char *pow_sym, unsigned level) const;
void do_print(const print_context & c, unsigned level) const override;
void do_print_latex(const print_latex & c, unsigned level) const;
void do_print_tree(const print_tree & c, unsigned level) const override;
void do_print_python(const print_python & c, unsigned level) const;
void do_print_python_repr(const print_python_repr & c, unsigned level) const override;
protected:
/** Vector of {coefficient, power} pairs */
epvector seq;
/** Series variable (holds a symbol) */
ex var;
/** Expansion point */
ex point;
};
// utility functions
/** Convert the pseries object embedded in an expression to an ordinary
* polynomial in the expansion variable. The result is undefined if the
* expression does not contain a pseries object at its top level.
*
* @param e expression
* @return polynomial expression
* @see is_a<>
* @see pseries::convert_to_poly */
inline ex series_to_poly(const ex &e)
{
return (ex_to<pseries>(e).convert_to_poly(true));
}
inline bool is_terminating(const pseries & s)
{
return s.is_terminating();
}
} // namespace GiNaC
#endif // ndef __GINAC_SERIES_H__
|