* Makes the interface to the underlying bignum package available. */
/*
- * GiNaC Copyright (C) 1999-2003 Johannes Gutenberg University Mainz, Germany
+ * GiNaC Copyright (C) 1999-2006 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__
#include "ex.h"
#include <stdexcept>
+#include <vector>
-#include <cln/number.h>
-// forward decln of cln::cl_N, since cln/complex_class.h is not included:
-namespace cln { class cl_N; }
+#include <cln/complex.h>
#if defined(G__CINTVERSION) && !defined(__MAKECINT__)
// Cint @$#$! doesn't like forward declaring classes used for casting operators
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
_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:
- void print(const print_context & c, unsigned level = 0) const;
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 lst & ls, const lst & lr, unsigned options = 0) const { return subs_one_level(ls, lr, options); } // overwrites basic::subs() for performance reasons
- ex normal(lst &sym_lst, lst &repl_lst, int level = 0) const;
- ex to_rational(lst &repl_lst) const;
- ex to_polynomial(lst &repl_lst) 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 */
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;
// 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;
};
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(); }
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__