synced to 1.0

'x+I' was missing the plus sign in the C source output

synced to 1.0

synced to 1.0

let_op() is no longer the default implementation of op(). Rationale: let_op()
needs to clear status_flags::hash_calculated, otherwise we could end up with
equal objects that have different hash values. But clearing hash_calculated
on a read-only op() access would be inefficient. let_op() is now only
implemented by lst and matrix.

- added non-const operator[]
- removed let_op() for power and relational
- removed redundant matrix::op()

synced to 1.0

*** empty log message ***

synced to 1.0

dimension in dirac gamma contractions was miscalculated

synced to 1.0

documentation update

dummy index renamer uses smaller set of global indices

dummy index renamer ignores index dimensions

documentation update

added to_polynomial(), to complement to_rational()

delta tensor makes itself square

synced to 1.0

synced to 1.0

added option to print index dimensions

better error reporting for minimal_dim()

*** empty log message ***

2002 -> 2003

fixed bugs in dummy index symmetrization and clifford contraction

* documentation update.

* Remove hack obsoleted with GCC 3.x.

GiNaC 1.1 requires GCC 3.x

implemented 'smartsubs' algebraic substitution, donated by Chris Dams
(tutorial not yet in sync, will update later)

cleaned up some is_a<> vs. is_exactly_a<> stuff

synced to 1.0

synced to 1.0

- synced to 1.0 (gammaL/R)
- removed dirac_gamma6/7()
- renamed simplify_ncmul() to eval_ncmul(), as it should have been called from
  the beginning

synced to 1.0

documentation update

synced to 1.0

synced to 1.0

synced to 1.0

synced to 1.0

synced to 1.0

synced to 1.0

synced to 1.0

synced to 1.0

synced to 1.0

Redone the way in which function serial numbers are stored. This allows for
using namespace qualifiers in the is_ex_the_function() macro like this

  is_ex_the_function(e, GiNaC::sin);

and it also works with the (overloaded) zeta() and psi() functions.
Unfortunately, it's not possible to

  using GiNaC::sin;
  is_ex_the_function(e, sin);

and the is_the_function<> template is not useful outside the is_ex_the_function()
macro. Even renaming the 'sin_SERIAL' class to 'sin' (which doesn't seem to
work with gcc3) you would still have to write 'is_the_function<class sin>(e)'
instead of 'is_the_function<sin>(e)'. I currently see no better solution,
though...

synced to 1.0

scalar_products allows the specification of an index dimension, e.g. to
give different values to 4- and D-dimensional products

documentation update

* Documentation update.

* lsolve(): Fix type of last argument (oops).

* irem(), iquo(): throw an exception, when second argument vanishes.
* inline some trivial calls to cln::instanceof().

* Update for Doxygen 1.2.18 (includes .gif to .png conversion).

* numeric::print(): increase precision for C source output.

synced to 1.0

* Fix comment.

* Fix CLN output of a-2.  (Chris Dams)

* Fix citation.  (The arctan formula for Pi is due to John Machin and
  not Pierre Me'chain.  This confusion came from Davenport's book.)

* Fix broken C output.  (Chris Dams)

* Fix CLN output of -2*a-2*b.  (Chris Dams)

synced to 1.0

* Add virtual method remove_all() to the containers (by Roberto Bagnara).

* Synced to 1.0.

* Synced to 1.0.

synced to 1.0

synced to 1.0

synced to 1.0

synced to 1.0 (tens4delta and mink4metric replaced by the new index dimension
scheme)

synced to 1.0

synced to 1.0 (for real this time...)

moved even more duplicated code to the 'tensor' class

synced to 1.0 with respect to dummy index repositioning

* Include the operators.h header file in the source files, not in ex.h.
  Rationale: It really gets in the way there, when trying to strip down
  the source files for debugging purposes.

* Enable tests N, Q and Q' by default.
  Rationale: Hubert finally is buying new cool hardware.   ;-)
* New test for expanding x^(x-1)*x.  Complain to JvdH.   ;-)

* Include the operators.h header file in the source files, not in ex.h.
  Rationale: It really gets in the way there, when trying to strip down
  the source files for debugging purposes.
* are_ex_trivially_equal(): use this only when reasonable, i.e. when
  an operation on all children was a no-op such that the this pointer is
  returned.

synced to 1.0

implemented 4-dimensional delta, metric and epsilon tensors

synced to 1.0

mentioned the imaginary unit 'I'

* Sync 1.1 branch to 1.0 mainline.

* rotate_left(): rewrite to help compiler make use of rol insns.
* golden_ratio_hash(): Fix comments.

* Oops, a FIXME has crept into the last commit.

* Globalize the hash-value-space.

* basic::compare(const basic&): Rearrange in a more compact fashion.
  (Rationale: produces better code.)

* expairseq::calchash(): Remove one superfluous rotation.

* numeric::calchash(): Apply golden_ratio_hash() to the value returned
  by CLN.  (Rationale: this populates more bits for the frequent small
  integers resulting in much less opportunity for funny cancellations in
  hash-values of monomials.)

* Added documentation for the add_dyn(), mul_dyn()... methods.

* Bump up repetitions since these tests have turned out to be very useful
  for detecting errors in the polynomial code.

* beta(1,x) evaluates to 1/x.

* mul::expand(): Fix bug in aggressive version of expansion of two sums
  and enable it.

* Clean up the output of all timings to stderr.

* mul::expand() (mul.cpp): When multiplying two sums, be careful not to
  allocate too much.

* expairseq::combine_same_terms_sorted_seq() (expairseq.cpp): make this
  ugh, ohhh, ahh, more... legible.  ;-)

* (time_fateman_expand.cpp): new file.
* (Makefile.am): compile it.
* (times.h, times.cpp): declare and call it.

* expairseq::expairseq(const epvector&, const ex& oc) (expairseq.cpp):
  check assertion: GINAC_ASSERT(is_a<numeric>(oc)).
* expairseq::expairseq(epvector*, const ex& oc) (expairseq.cpp): Likewise.
* expair_rest_is_less (expair.h): New functor...
* expairseq::canonicalize() (expairseq.cpp): ...used here.

* power::expand_add(): allocate the precise amount needed for representing
  the result.  And some minor readability cleanups...

* Don't push void to static Bernoulli vector, it breaks with GCC-3.1.

synced with main branch

* Add namespace declarators to is_ex_the_function().

* Remove " syntax for the execution stack.
  (You'll need a bison 1.28 to regenerate your files, since newer versions
  are all broken.  Dammit.)

* Adjust to up-to-date compiler landscape.

* Removed deprecated macros is_ex_a, is_ex_exactly_a and friends.