/usr/include/pynac/numeric.h is in libpynac-dev 0.3.2+dfsg-1.
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 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 | /** @file numeric.h
*
* Makes the interface to the underlying bignum package available. */
/*
* This is a modified version of code included with Ginac.
* The modifications and modified version is:
*
* GiNaC-SAGE Copyright (C) 2008 William Stein
*
* 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
*/
/* The original copyright:
* 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_NUMERIC_H__
#define __GINAC_NUMERIC_H__
#include "basic.h"
#include "constant.h"
#include "ex.h"
#include "py_funcs.h"
#include <stdexcept>
#include <vector>
#include <iostream>
void ginac_pyinit_Integer(PyObject*);
void ginac_pyinit_Float(PyObject*);
void ginac_pyinit_I(PyObject*);
namespace GiNaC {
enum Type {
LONG,
DOUBLE,
PYOBJECT,
MPFR,
MPZ,
MPQ,
MPFC,
MPQC
};
union Value {
double _double;
long int _long;
PyObject* _pyobject;
};
class Number_T { // abstract base class
protected:
Type t;
Value v;
public:
Number_T();
Number_T(const int& x);
Number_T(const long int& x);
Number_T(const unsigned int& x);
Number_T(const unsigned long& x);
Number_T(const double& x);
Number_T(const Number_T& x);
Number_T(const char* s);
Number_T(PyObject* o); // *STEALS a REFERENCE*
// unarchive
Number_T(const archive_node &n, lst &sym_lst);
// Destructor
~Number_T();
// archive
void archive(archive_node &n) const;
Number_T operator+(Number_T x) const;
Number_T operator*(Number_T x) const;
Number_T operator/(Number_T x) const;
Number_T operator-(Number_T x) const;
Number_T& operator=(const Number_T& x);
Number_T operator()(const int& x);
Number_T operator-();
operator double() const;
operator int() const;
operator long int() const;
unsigned hash() const;
int compare_same_type(const Number_T& right) const;
bool operator==(const Number_T& right) const;
bool operator!=(const Number_T& right) const;
bool operator<=(const Number_T& right) const;
bool operator>=(const Number_T& right) const;
bool operator<(const Number_T& right) const;
bool operator>(const Number_T& right) const;
int csgn() const;
int get_parent_char() const;
bool is_zero() const;
bool is_positive() const;
bool is_negative() const;
bool is_integer() const;
bool is_cinteger() const;
bool is_pos_integer() const;
bool is_nonneg_integer() const;
bool is_even() const;
bool is_odd() const;
bool is_prime() const;
bool is_rational() const;
bool is_crational() const;
bool is_real() const;
Number_T numer() const;
Number_T denom() const;
Number_T lcm(Number_T b) const;
Number_T gcd(Number_T b) const;
Number_T conjugate() const;
Number_T sin() const;
Number_T cos() const;
Number_T exp() const;
Number_T log() const;
Number_T tan() const;
Number_T asin() const;
Number_T acos() const;
Number_T atan() const;
Number_T atan(const Number_T& y) const;
Number_T sinh() const;
Number_T cosh() const;
Number_T tanh() const;
Number_T asinh() const;
Number_T acosh() const;
Number_T atanh() const;
Number_T Li(const Number_T &n, PyObject* parent) const;
Number_T Li2() const;
Number_T zeta() const;
Number_T lgamma() const;
Number_T tgamma() const;
Number_T psi() const;
Number_T psi(const Number_T& y) const;
Number_T factorial() const;
Number_T doublefactorial() const;
Number_T fibonacci() const;
Number_T evalf(PyObject* parent) const;
Number_T step() const;
Number_T isqrt() const;
Number_T sqrt() const;
Number_T abs() const;
Number_T mod(const Number_T &b) const;
Number_T smod(const Number_T &b) const;
Number_T irem(const Number_T &b) const;
Number_T irem(const Number_T &b, Number_T& q) const;
Number_T iquo(const Number_T &b) const;
Number_T iquo(const Number_T &b, Number_T& q) const;
int int_length() const;
friend std::ostream& operator << (std::ostream& os, const Number_T& s);
friend Number_T pow(const Number_T& base, const Number_T& exp);
friend void coerce(Number_T& new_left, Number_T& new_right, const Number_T& left, const Number_T& right);
friend PyObject* Number_T_to_pyobject(const Number_T& x);
};
std::ostream& operator << (std::ostream& os, const Number_T& s);
Number_T pow(const Number_T& base, const Number_T& exp);
/** Exception class thrown when a singularity is encountered. */
class pole_error : public std::domain_error {
public:
explicit pole_error(const std::string& what_arg, int degree);
int degree() const;
private:
int deg;
};
/** This class is a wrapper around CLN-numbers within the GiNaC class
* hierarchy. Objects of this type may directly be created by the user.*/
class numeric : public basic
{
GINAC_DECLARE_REGISTERED_CLASS(numeric, basic)
// member functions
// other constructors
public:
numeric(int i);
numeric(unsigned int i);
numeric(long i);
numeric(unsigned long i);
numeric(long numer, long denom);
numeric(double d);
numeric(const char *);
numeric(const Number_T& x);
numeric(PyObject*);
numeric(archive_node& n, const lst &sym_lst);
// functions overriding virtual functions from base classes
public:
unsigned precedence() const {return 30;}
bool info(unsigned inf) const;
bool is_polynomial(const ex & var) const;
int degree(const ex & s) const;
int ldegree(const ex & s) const;
ex coeff(const ex & s, int n = 1) const;
bool has(const ex &other, unsigned options = 0) const;
ex eval(int level = 0) const;
ex evalf(int level = 0, PyObject* parent=NULL) const;
ex subs(const exmap & m, unsigned options = 0) const { return subs_one_level(m, options); } // overwrites basic::subs() for performance reasons
ex normal(exmap & repl, exmap & rev_lookup, int level = 0) const;
ex to_rational(exmap & repl) const;
ex to_polynomial(exmap & repl) const;
numeric integer_content() const;
ex smod(const numeric &xi) const;
numeric max_coefficient() const;
ex conjugate() const;
ex real_part() const;
ex imag_part() const;
protected:
/** Implementation of ex::diff for a numeric always returns 0.
* @see ex::diff */
ex derivative(const symbol &s) const { return 0; }
bool is_equal_same_type(const basic &other) const;
unsigned calchash() const;
// new virtual functions which can be overridden by derived classes
// (none)
// non-virtual functions in this class
public:
const numeric add(const numeric &other) const;
const numeric sub(const numeric &other) const;
const numeric mul(const numeric &other) const;
const numeric div(const numeric &other) const;
const numeric power(const numeric &other) const;
const numeric & add_dyn(const numeric &other) const;
const numeric & sub_dyn(const numeric &other) const;
const numeric & mul_dyn(const numeric &other) const;
const numeric & div_dyn(const numeric &other) const;
const numeric & power_dyn(const numeric &other) const;
const numeric & operator=(int i);
const numeric & operator=(unsigned int i);
const numeric & operator=(long i);
const numeric & operator=(unsigned long i);
const numeric & operator=(double d);
const numeric & operator=(const char *s);
const numeric inverse() const;
numeric step() const;
int csgn() const;
int compare(const numeric &other) const;
bool is_equal(const numeric &other) const;
bool is_zero() const;
bool is_positive() const;
bool is_negative() const;
bool is_integer() const;
bool is_pos_integer() const;
bool is_nonneg_integer() const;
bool is_even() const;
bool is_odd() const;
bool is_prime() const;
bool is_rational() const;
bool is_real() const;
bool is_cinteger() const;
bool is_crational() const;
bool operator==(const numeric &other) const;
bool operator!=(const numeric &other) const;
bool operator<(const numeric &other) const;
bool operator<=(const numeric &other) const;
bool operator>(const numeric &other) const;
bool operator>=(const numeric &other) const;
bool is_parent_pos_char() const;
int get_parent_char() const;
int to_int() const;
long to_long() const;
double to_double() const;
PyObject* to_pyobject() const;
const numeric real() const;
const numeric imag() const;
const numeric numer() const;
const numeric denom() const;
int int_length() const;
//numeric(const Number_T &z);
protected:
void print_numeric(const print_context & c, const char *par_open,
const char *par_close, const char *imag_sym,
const char *mul_sym, unsigned level, bool latex) const;
void do_print(const print_context & c, unsigned level) const;
void do_print_latex(const print_latex & c, unsigned level) const;
void do_print_csrc(const print_csrc & c, unsigned level) const;
void do_print_csrc_cl_N(const print_csrc_cl_N & c, unsigned level) const;
void do_print_tree(const print_tree & c, unsigned level) const;
void do_print_python_repr(const print_python_repr & c, unsigned level) const;
// member variables
//protected:
public:
Number_T value;
};
// global constants
extern numeric I;
// global functions
const numeric exp(const numeric &x);
const numeric log(const numeric &x);
const numeric sin(const numeric &x);
const numeric cos(const numeric &x);
const numeric tan(const numeric &x);
const numeric asin(const numeric &x);
const numeric acos(const numeric &x);
const numeric atan(const numeric &x);
const numeric atan(const numeric &y, const numeric &x);
const numeric sinh(const numeric &x);
const numeric cosh(const numeric &x);
const numeric tanh(const numeric &x);
const numeric asinh(const numeric &x);
const numeric acosh(const numeric &x);
const numeric atanh(const numeric &x);
const numeric Li2(const numeric &x);
const numeric zeta(const numeric &x);
const numeric lgamma(const numeric &x);
const numeric tgamma(const numeric &x);
const numeric psi(const numeric &x);
const numeric psi(const numeric &n, const numeric &x);
const numeric factorial(const numeric &n);
const numeric doublefactorial(const numeric &n);
const numeric binomial(const numeric &n, const numeric &k);
const numeric bernoulli(const numeric &n);
const numeric fibonacci(const numeric &n);
const numeric isqrt(const numeric &x);
const numeric sqrt(const numeric &x);
const numeric abs(const numeric &x);
const numeric mod(const numeric &a, const numeric &b);
const numeric smod(const numeric &a, const numeric &b);
const numeric irem(const numeric &a, const numeric &b);
const numeric iquo(const numeric &a, const numeric &b);
const numeric iquo(const numeric &a, const numeric &b, numeric &r);
const numeric gcd(const numeric &a, const numeric &b);
const numeric lcm(const numeric &a, const numeric &b);
// wrapper functions around member functions
inline const numeric pow(const numeric &x, const numeric &y)
{ return x.power(y); }
inline const numeric inverse(const numeric &x)
{ return x.inverse(); }
inline numeric step(const numeric &x)
{ return x.step(); }
inline int csgn(const numeric &x)
{ return x.csgn(); }
inline bool is_zero(const numeric &x)
{ return x.is_zero(); }
inline bool is_positive(const numeric &x)
{ return x.is_positive(); }
inline bool is_negative(const numeric &x)
{ return x.is_negative(); }
inline bool is_integer(const numeric &x)
{ return x.is_integer(); }
inline bool is_pos_integer(const numeric &x)
{ return x.is_pos_integer(); }
inline bool is_nonneg_integer(const numeric &x)
{ return x.is_nonneg_integer(); }
inline bool is_even(const numeric &x)
{ return x.is_even(); }
inline bool is_odd(const numeric &x)
{ return x.is_odd(); }
inline bool is_prime(const numeric &x)
{ return x.is_prime(); }
inline bool is_rational(const numeric &x)
{ return x.is_rational(); }
inline bool is_real(const numeric &x)
{ return x.is_real(); }
inline bool is_cinteger(const numeric &x)
{ return x.is_cinteger(); }
inline bool is_crational(const numeric &x)
{ return x.is_crational(); }
inline int to_int(const numeric &x)
{ return x.to_int(); }
inline long to_long(const numeric &x)
{ return x.to_long(); }
inline double to_double(const numeric &x)
{ return x.to_double(); }
inline const numeric real(const numeric &x)
{ return x.real(); }
inline const numeric imag(const numeric &x)
{ return x.imag(); }
inline const numeric numer(const numeric &x)
{ return x.numer(); }
inline const numeric denom(const numeric &x)
{ return x.denom(); }
// numeric evaluation functions for class constant objects:
ex ConstantEvalf(unsigned serial, PyObject* parent=NULL);
ex UnsignedInfinityEvalf(unsigned serial, PyObject* parent=NULL);
ex InfinityEvalf(unsigned serial, PyObject* parent=NULL);
ex NegInfinityEvalf(unsigned serial, PyObject* parent=NULL);
} // namespace GiNaC
#endif // ndef __GINAC_NUMERIC_H__
|