X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Fnumeric.h;h=30091760507018a5943819d9a1e35253f9c3c32c;hp=30d4edf34116f5bc9b675845a9d64604b0dd4b45;hb=02a1c658c196cdf0090fdd134a9aeb884a671922;hpb=4afb5bbe2c0b0a60928120a042997ba7d89e8f5c diff --git a/ginac/numeric.h b/ginac/numeric.h index 30d4edf3..30091760 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 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2000 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 @@ -24,15 +24,15 @@ #define __GINAC_NUMERIC_H__ #include -#include -#include +#include "basic.h" +#include "ex.h" class cl_N; // We want to include cln.h only in numeric.cpp in order to // avoid namespace pollution and keep compile-time low. -#ifndef NO_GINAC_NAMESPACE +#ifndef NO_NAMESPACE_GINAC namespace GiNaC { -#endif // ndef NO_GINAC_NAMESPACE +#endif // ndef NO_NAMESPACE_GINAC #define HASHVALUE_NUMERIC 0x80000001U @@ -63,24 +63,27 @@ private: * hierarchy. Objects of this type may directly be created by the user.*/ class numeric : public basic { + GINAC_DECLARE_REGISTERED_CLASS(numeric, basic) + // friends - friend numeric exp(const numeric & x); - friend numeric log(const numeric & x); - friend numeric sin(const numeric & x); - friend numeric cos(const numeric & x); - friend numeric tan(const numeric & x); - friend numeric asin(const numeric & x); - friend numeric acos(const numeric & x); - friend numeric atan(const numeric & x); - friend numeric atan(const numeric & y, const numeric & x); - friend numeric sinh(const numeric & x); - friend numeric cosh(const numeric & x); - friend numeric tanh(const numeric & x); - friend numeric asinh(const numeric & x); - friend numeric acosh(const numeric & x); - friend numeric atanh(const numeric & x); - friend numeric zeta(const numeric & x); - friend numeric bernoulli(const numeric & n); + friend const numeric exp(const numeric & x); + friend const numeric log(const numeric & x); + friend const numeric sin(const numeric & x); + friend const numeric cos(const numeric & x); + friend const numeric tan(const numeric & x); + friend const numeric asin(const numeric & x); + friend const numeric acos(const numeric & x); + friend const numeric atan(const numeric & x); + friend const numeric atan(const numeric & y, const numeric & x); + friend const numeric sinh(const numeric & x); + friend const numeric cosh(const numeric & x); + friend const numeric tanh(const numeric & x); + friend const numeric asinh(const numeric & x); + friend const numeric acosh(const numeric & x); + friend const numeric atanh(const numeric & x); + friend const numeric zeta(const numeric & x); + friend const numeric bernoulli(const numeric & n); + friend const numeric fibonacci(const numeric & n); friend numeric abs(const numeric & x); friend numeric mod(const numeric & a, const numeric & b); friend numeric smod(const numeric & a, const numeric & b); @@ -114,9 +117,9 @@ public: explicit numeric(unsigned long i); explicit numeric(long numer, long denom); explicit numeric(double d); - explicit numeric(char const *); - numeric(cl_N const & z); - + explicit numeric(const char *); + numeric(const cl_N & z); + // functions overriding virtual functions from bases classes public: basic * duplicate() const; @@ -125,19 +128,18 @@ public: void printtree(ostream & os, unsigned indent) const; void printcsrc(ostream & os, unsigned type, unsigned precedence=0) const; bool info(unsigned inf) const; + bool has(const ex & other) const; + ex eval(int level=0) const; ex evalf(int level=0) const; - ex diff(symbol const & s) const; ex normal(lst &sym_lst, lst &repl_lst, int level=0) const; numeric integer_content(void) const; ex smod(const numeric &xi) const; numeric max_coefficient(void) const; protected: - int compare_same_type(basic const & other) const; - bool is_equal_same_type(basic const & other) const; - unsigned calchash(void) const { - hashvalue=HASHVALUE_NUMERIC; - return HASHVALUE_NUMERIC; - } + ex derivative(const symbol & s) const; + int compare_same_type(const basic & other) const; + bool is_equal_same_type(const basic & other) const; + unsigned calchash(void) const; // new virtual functions which can be overridden by derived classes // (none) @@ -159,7 +161,7 @@ public: const numeric & operator=(long i); const numeric & operator=(unsigned long i); const numeric & operator=(double d); - const numeric & operator=(char const * s); + const numeric & operator=(const char * s); numeric inverse(void) const; int csgn(void) const; int compare(const numeric & other) const; @@ -184,6 +186,7 @@ public: 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; numeric real(void) const; numeric imag(void) const; @@ -202,37 +205,41 @@ protected: extern const numeric some_numeric; extern const numeric I; -extern type_info const & typeid_numeric; +extern const type_info & typeid_numeric; extern _numeric_digits Digits; -#define is_a_numeric_hash(x) ((x)==HASHVALUE_NUMERIC) +//#define is_a_numeric_hash(x) ((x)==HASHVALUE_NUMERIC) // may have to be changed to ((x)>=0x80000000U) +// has been changed +//#define is_a_numeric_hash(x) ((x)&0x80000000U) + // global functions -numeric exp(const numeric & x); -numeric log(const numeric & x); -numeric sin(const numeric & x); -numeric cos(const numeric & x); -numeric tan(const numeric & x); -numeric asin(const numeric & x); -numeric acos(const numeric & x); -numeric atan(const numeric & x); -numeric atan(const numeric & y, const numeric & x); -numeric sinh(const numeric & x); -numeric cosh(const numeric & x); -numeric tanh(const numeric & x); -numeric asinh(const numeric & x); -numeric acosh(const numeric & x); -numeric atanh(const numeric & x); -numeric zeta(const numeric & x); -numeric gamma(const numeric & x); -numeric psi(const numeric & x); -numeric psi(const numeric & n, const numeric & x); -numeric factorial(const numeric & n); -numeric doublefactorial(const numeric & n); -numeric binomial(const numeric & n, const numeric & k); -numeric bernoulli(const numeric & n); +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 zeta(const numeric & x); +const numeric gamma(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); numeric abs(const numeric & x); numeric mod(const numeric & a, const numeric & b); @@ -318,14 +325,16 @@ ex PiEvalf(void); ex EulerGammaEvalf(void); ex CatalanEvalf(void); + // utility functions inline const numeric &ex_to_numeric(const ex &e) { return static_cast(*e.bp); } -#ifndef NO_GINAC_NAMESPACE + +#ifndef NO_NAMESPACE_GINAC } // namespace GiNaC -#endif // ndef NO_GINAC_NAMESPACE +#endif // ndef NO_NAMESPACE_GINAC #endif // ndef __GINAC_NUMERIC_H__