* Makes the interface to the underlying bignum package available. */
/*
- * GiNaC Copyright (C) 1999-2004 Johannes Gutenberg University Mainz, Germany
+ * GiNaC Copyright (C) 1999-2016 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
*
* 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 <stdexcept>
+#include "archive.h"
#include <cln/complex.h>
-
-#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 <cln/complex_class.h>
-#endif
+#include <stdexcept>
+#include <vector>
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
* cl_default_float_format when it gets changed. The only other
* meaningful thing to do with it is converting it to an unsigned,
- * for temprary storing its value e.g. The user must not create an
+ * for temporarily storing its value e.g. The user must not create an
* own working object of this class! Since C++ forces us to make the
* class definition visible in order to use an object we put in a
* flag which prevents other objects of that class to be created. */
_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<digits_changed_callback> callbacklist;
};
// functions overriding virtual functions from base classes
public:
- unsigned precedence() const {return 30;}
- bool info(unsigned inf) 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;
- 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(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;
+ unsigned precedence() const override {return 30;}
+ bool info(unsigned inf) const override;
+ bool is_polynomial(const ex & var) const override;
+ int degree(const ex & s) const override;
+ int ldegree(const ex & s) const override;
+ ex coeff(const ex & s, int n = 1) const override;
+ bool has(const ex &other, unsigned options = 0) const override;
+ ex eval() const override;
+ ex evalf() const override;
+ ex subs(const exmap & m, unsigned options = 0) const override { return subs_one_level(m, options); } // overwrites basic::subs() for performance reasons
+ ex normal(exmap & repl, exmap & rev_lookup) const override;
+ ex to_rational(exmap & repl) const override;
+ ex to_polynomial(exmap & repl) const override;
+ numeric integer_content() const override;
+ ex smod(const numeric &xi) const override;
+ numeric max_coefficient() const override;
+ ex conjugate() const override;
+ ex real_part() const override;
+ ex imag_part() const override;
+ /** Save (a.k.a. serialize) object into archive. */
+ void archive(archive_node& n) const override;
+ /** Read (a.k.a. deserialize) object from archive. */
+ void read_archive(const archive_node& n, lst& syms) override;
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;
+ ex derivative(const symbol &s) const override { return 0; }
+ bool is_equal_same_type(const basic &other) const override;
+ unsigned calchash() const override;
// new virtual functions which can be overridden by derived classes
// (none)
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;
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;
protected:
cln::cl_N value;
};
+GINAC_DECLARE_UNARCHIVER(numeric);
// global constants
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_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(); }
ex CatalanEvalf();
-// utility functions
-
-/** Specialization of is_exactly_a<numeric>(obj) for numeric objects. */
-template<> inline bool is_exactly_a<numeric>(const basic & obj)
-{
- return obj.tinfo()==TINFO_numeric;
-}
-
} // namespace GiNaC
-#ifdef __MAKECINT__
-#pragma link off defined_in cln/number.h;
-#pragma link off defined_in cln/complex_class.h;
-#endif
-
-#endif // ndef __GINAC_NUMERIC_H__
+#endif // ndef GINAC_NUMERIC_H