X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Fnumeric.h;h=3e0dd18500b7998f08236436224a24f339e46ffb;hp=b7c1e37d0eab208273c4a86c538b94417ea0ae39;hb=9177e7536ea82b739c72d3e41a319af7dbc15661;hpb=b4be7b0f30fbb6178cf4ee83e1b3952e084bd8ca diff --git a/ginac/numeric.h b/ginac/numeric.h index b7c1e37d..3e0dd185 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-2003 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2010 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,30 +17,34 @@ * * 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__ -#define __GINAC_NUMERIC_H__ +#ifndef GINAC_NUMERIC_H +#define GINAC_NUMERIC_H #include "basic.h" #include "ex.h" +#include "archive.h" -#include - -#include -// forward decln of cln::cl_N, since cln/complex_class.h is not included: -namespace cln { class cl_N; } - +#include #if defined(G__CINTVERSION) && !defined(__MAKECINT__) // Cint @$#$! doesn't like forward declaring classes used for casting operators // so we have to include the definition of cln::cl_N here, but it is enough to // do so for the compiler, hence the !defined(__MAKECINT__). #include #endif +#include +#include namespace GiNaC { +/** 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 * than a dumber basic type since as a side-effect we let it change @@ -57,11 +61,14 @@ 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; }; @@ -97,19 +104,27 @@ public: 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) const; + bool has(const ex &other, unsigned options = 0) const; ex eval(int level = 0) const; ex evalf(int level = 0) 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(lst &repl_lst) const; - ex to_polynomial(lst &repl_lst) 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; + /** Save (a.k.a. serialize) object into archive. */ + void archive(archive_node& n) const; + /** Read (a.k.a. deserialize) object from archive. */ + void read_archive(const archive_node& n, lst& syms); protected: /** Implementation of ex::diff for a numeric always returns 0. * @see ex::diff */ @@ -139,6 +154,7 @@ public: 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; @@ -171,7 +187,7 @@ public: const numeric denom() const; int int_length() const; // converting routines for interfacing with CLN: - numeric(const cln::cl_N &z); + explicit 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; @@ -185,8 +201,9 @@ protected: // member variables protected: - cln::cl_number value; + cln::cl_N value; }; +GINAC_DECLARE_UNARCHIVER(numeric); // global constants @@ -241,6 +258,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(); } @@ -250,6 +270,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(); } @@ -308,14 +331,6 @@ ex EulerEvalf(); ex CatalanEvalf(); -// utility functions - -/** Specialization of is_exactly_a(obj) for numeric objects. */ -template<> inline bool is_exactly_a(const basic & obj) -{ - return obj.tinfo()==TINFO_numeric; -} - } // namespace GiNaC #ifdef __MAKECINT__ @@ -323,4 +338,4 @@ template<> inline bool is_exactly_a(const basic & obj) #pragma link off defined_in cln/complex_class.h; #endif -#endif // ndef __GINAC_NUMERIC_H__ +#endif // ndef GINAC_NUMERIC_H