X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Fnumeric.h;h=f52ffc0bd73e090b8b7f28621bbed14074a86672;hp=4c7fd0b33162ab7572f8de42a75e7ee0d3747e83;hb=6d72e4b12dbd1bb053563bd9fab5fad56e052d1f;hpb=0917777d60656298ba44e092712c2b403bfaf559 diff --git a/ginac/numeric.h b/ginac/numeric.h index 4c7fd0b3..f52ffc0b 100644 --- a/ginac/numeric.h +++ b/ginac/numeric.h @@ -3,7 +3,7 @@ * Makes the interface to the underlying bignum package available. */ /* - * GiNaC Copyright (C) 1999-2001 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2007 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 @@ -17,7 +17,7 @@ * * 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., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ #ifndef __GINAC_NUMERIC_H__ @@ -26,9 +26,10 @@ #include "basic.h" #include "ex.h" -#include -// forward decln of cln::cl_N, since cln/complex_class.h is not included: -namespace cln { class cl_N; } +#include +#include + +#include #if defined(G__CINTVERSION) && !defined(__MAKECINT__) // Cint @$#$! doesn't like forward declaring classes used for casting operators @@ -39,7 +40,11 @@ namespace cln { class cl_N; } namespace GiNaC { -#define HASHVALUE_NUMERIC 0x80000001U +/** Function pointer to implement callbacks in the case 'Digits' gets changed. + * Main purpose of such callbacks is to adjust look-up tables of certain + * functions to the new precision. Parameter contains the signed difference + * between new Digits and old Digits. */ +typedef void (* digits_changed_callback)(long); /** This class is used to instantiate a global singleton object Digits * which behaves just like Maple's Digits. We need an object rather @@ -57,13 +62,27 @@ public: _numeric_digits(); _numeric_digits& operator=(long prec); operator long(); - void print(std::ostream &os) const; + void print(std::ostream& os) const; + void add_callback(digits_changed_callback callback); // member variables private: long digits; ///< Number of decimal digits static bool too_late; ///< Already one object present + // Holds a list of functions that get called when digits is changed. + std::vector callbacklist; +}; + + +/** 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 @@ -72,35 +91,43 @@ class numeric : public basic // member functions - // other ctors + // other constructors public: - explicit numeric(int i); - explicit numeric(unsigned int i); - explicit numeric(long i); - explicit numeric(unsigned long i); - explicit numeric(long numer, long denom); - explicit numeric(double d); - explicit numeric(const char *); + 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 *); - // functions overriding virtual functions from bases classes + // functions overriding virtual functions from base classes public: - void print(const print_context & c, unsigned level = 0) const; - unsigned precedence(void) const {return 30;} + unsigned precedence() const {return 30;} bool info(unsigned inf) const; - bool has(const ex &other) 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) const; - ex normal(lst &sym_lst, lst &repl_lst, int level = 0) const; - ex to_rational(lst &repl_lst) const; - numeric integer_content(void) 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(void) 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 _ex0(); } + ex derivative(const symbol &s) const { return 0; } bool is_equal_same_type(const basic &other) const; - unsigned calchash(void) const; + unsigned calchash() const; // new virtual functions which can be overridden by derived classes // (none) @@ -123,55 +150,63 @@ public: const numeric & operator=(unsigned long i); const numeric & operator=(double d); const numeric & operator=(const char *s); - const numeric inverse(void) const; - int csgn(void) const; + 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(void) const; - bool is_positive(void) const; - bool is_negative(void) const; - bool is_integer(void) const; - bool is_pos_integer(void) const; - bool is_nonneg_integer(void) const; - bool is_even(void) const; - bool is_odd(void) const; - bool is_prime(void) const; - bool is_rational(void) const; - bool is_real(void) const; - bool is_cinteger(void) const; - bool is_crational(void) 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; - int to_int(void) const; - long to_long(void) const; - double to_double(void) const; - cln::cl_N to_cl_N(void) const; - const numeric real(void) const; - const numeric imag(void) const; - const numeric numer(void) const; - const numeric denom(void) const; - int int_length(void) const; + int to_int() const; + long to_long() const; + double to_double() const; + cln::cl_N to_cl_N() const; + const numeric real() const; + const numeric imag() const; + const numeric numer() const; + const numeric denom() const; + int int_length() const; // converting routines for interfacing with CLN: numeric(const cln::cl_N &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) 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: - cln::cl_number value; + cln::cl_N value; }; + // global constants extern const numeric I; extern _numeric_digits Digits; -// deprecated macro, for internal use only -#define is_a_numeric_hash(x) ((x)&0x80000000U) - // global functions const numeric exp(const numeric &x); @@ -219,6 +254,9 @@ inline const numeric pow(const numeric &x, const numeric &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(); } @@ -228,6 +266,9 @@ inline bool is_zero(const numeric &x) 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(); } @@ -281,16 +322,10 @@ inline const numeric denom(const numeric &x) // numeric evaluation functions for class constant objects: -ex PiEvalf(void); -ex EulerEvalf(void); -ex CatalanEvalf(void); - +ex PiEvalf(); +ex EulerEvalf(); +ex CatalanEvalf(); -// utility functions -inline const numeric &ex_to_numeric(const ex &e) -{ - return static_cast(*e.bp); -} } // namespace GiNaC