using () operator for setting matrix elements

documentation update

- mentioned that powers of matrices are not automatically expanded
- fixed typos

- Removed obsolete (even for Cint!) malloc_alloc template arguments.

- Removed cint subdir for good because the new separate package seems to
  work properly now (after some longish argument with makecint).

- config.*: updated from FSF.
- Removed cint subdir for good because the new separate package seems to
  work properly now (after some longish argument with makecint).

- inserted a couple of missing namepace std:: resolutions.

- document: collect's argument must be expanded.

- Small LaTeX output fixes as suggested by Stefan.

- fixed typos
- c*M^n with M being a matrix wasn't correctly evalm()'d for n!=1

- ++version_major.
- added matrix::pow() to handle integer exponents with the least amount of
  multiplications possible and...
- ...added power::evalm() to actually use it.
- classhierarchy.fig: add class wildcard.
- some cleanups.

- list delimiters are now { } and matrix delimiters are [ ]
- added evalm() function

- in the output, list delimiters are now { } and matrix delimiters are [ ]
- added evalm() method for evaluating sums and products of matrices
- added map() method for applying a function to subexpressions
- power(noncommut, posint) automatically expands the product, removed ncpow()

- list delimiters are now { }, matrix delimiters are now [ ]
- added an example for evalm()
- removed ncpow()

added test for evalm()

documentation update

- added numer_denom()
- minor cleanups

- added numer_denom() to get numerator and denominator in one pass
- symmetrize() and antisymmetrize() are available as ex methods
- minor cleanups

- mentioned numer_denom()
- symmetrize() and antisymmetrize() are available as ex methods

documentation update

determinant_algo::laplace sometimes had the sign wrong because of the change
to permutation_sign()

- permutation_sign() and shaker_sort() are no longer inlined
- symmetrize() no longer limited to 32-bit integers

- dirac_trace() is twice as fast
- permutation_sign() uses shaker sort
- shaker_sort() doesn't require less-than comparable iterators any more

documented symmetrize() and antisymmetrize()

- added symmetrize() and antisymmetrize() functions
- generalized permutation_sign() template
- moved shaker_sort() template to utils.h

- renamed ncpower() to ncpow()
- canonicalize_clifford() uses to_rational() now
- minor cleanups

mentioned ncpow()

fixed typos

added description of wildcards

- removed debugging code in match()
- added ex::match(const ex & pattern) which doesn't take a list as second
  argument
- expairseq::subs() substitutes in recombined pairs when the pattern is a
  product or a power; this is necessary for things like subs(2*x^2,x^2==3)
  because the first expression is mul(<x,2>,2) and only substituting in the
  "rest" part of the expair is not enough

added section "Pattern matching and advanced substitutions" which covers
match(), has() and subs()

added wildcard and spinormetric classes to check

- fix LaTeX-output bug reported by Stefan, remove obsolete has(matrix,ex).

- first implementation of pattern matching

- Clifford traces of many gammas are a lot faster now (especially with gamma5)
- new simplifications for color classes:
    d.abc T.b T.c = 5/6 T.a
    f.abc T.b T.c = 3/2 I T.a
- added indexed::has_dummy_index_for() member function

- mentioned dummy index renaming
- Clifford and color class description slightly updated

updated for changes in Clifford/color classes

replaced "precedence" static member variable by virtual precedence() function

- adjust preprocessor Cint-checks to the actual situation.

- dummy index renamer didn't account for internal dummy indices of objects
  in products; a~mu.mu-a~nu.nu gets simplified to 0 now
- made a little more use of STL facilities for exvectors, especially in the
  indexed stuff; append_exvector_to_exvector() and index_set_difference() are
  gone and utils.h defines the functors ex_is_less and ex_is_equal

added one more check

documentation update

- simplify_indexed() renames dummy indices so, e.g., "a.i*a.i+a.j*a.j" gets
  simplified to "2*a.i*a.i" (or 2*a.j*a.j, you can't know which one)
- fixed possible crash when calling subs() on expressions with non-commutative
  products
- added canonicalize_clifford()

corrected typo

added checks for dummy index renaming and canonicalize_clifford()

documentation update

- Missed to #include "config.h" for HAVE_SSTREAM.

fixed bug in mul::print with print_csrc

- Fix thinko in csgn_eval().

- Indentation Fix.

- added patch contributors, added remark to mailing lists.

- forgot to close #97136.

- applied `dch -i'.

- Renamed our detected libreadline version in order to avoid conflicts
  with Chat Rameys (still broken) macros and future CPP definitions in
  libreadline's headers.  Oh what a pain in the ass!

- dirac_trace() takes an optional third argument: the expression to use for
  the trace of the unit matrix
- more conventional factor I in gamma5

dirac_trace() now takes three arguments

added missing factor I in gamma5 trace

- Added inclusion of <sstream> or <strstream> since <ginac/ginac.h> now
  does not do so any more.

* Zapped header dependency from <strstream>.

* comparison operators should return const relationals.

* node assignment was not returning *this.

- color and clifford classes are quite functional now
- new "spinidx" class for dotted/undotted indices
- predefined spinor metric tensor (created by spinor_metric())

added documentation for spinidx, spinor_metric() and noncommutative stuff

more checks (more, MORE, BWAHAHHAAA!... uhm, ...) for noncommutative classes

documentation update

- Prepare for landing of 0.8.3.

- minor cleanup.

* Added new operator semantics.

* Added some messy exams for new operator semantics.

* All operators are now strictly following the semantics we all know from
  C/C++.  I.e. e1+e2=42; is now impossible and other such stuff.  For all
  remaining errors (i.e. ++e1 = e1+1;) the compiler alone is to blame...

- Remved obsolete remainders of preprocessor symbol NO_NAMESPACE_GINAC.

* Methods of class ex which do absolutely nothing than type dispatch should
  be inlined for obvious reasons.

* Updated Cint status.

Changes needed for compilation with -DDO_GINAC_ASSERT:
* basic.cpp: if needed, #include <typeinfo>,
* mul.h: mul::print() needed a default for 2nd arg.
Hey, Cebix, could you please have a look at the assertions of type
basis.return_type()==return_types::commutative in power.cpp ctors?
Watch out, though: they are not just bogus because of the rule
(x*y)^2 -> x^2*y^2 in power::eval().

- symbols can have a LaTeX name, e.g. symbol s("s", "\\sigma");
- LaTeX output of indexed objects and tensors is much nicer
- simplify_indexed() can do arbitrary contractions in Dirac gamma strings
- constructing expressions like
    (x * dirac_ONE() + dirac_ONE()) * dirac_gamma(mu)
  resulted in an exception because mul::simplify_ncmul() wasn't implemented
  (maybe afrink didn't anticipate that an implementation would really be
  needed...)

- more clifford checks
- archive checks has a symbol with a LaTeX name

documentation and grammar updates

* add::coeff(): don't epvector::reserve() more terms than we'll actually
  need and don't epvector::push_back() all those _ex0()'s into coeffseq.

* Fix incompatibilities with recent versions of libreadline by having
  configure find out the version and then doing ugly things with #ifdef.

* Fix incompatibilities with recent versions of libreadline by having
  configure find out the version and then doing ugly things with #ifdef.
* config.*: updates from upstream.

* Just for the record: this is what has been used for building.

- dirac_trace() handles arbitrary expressions (even unexpanded), including
  gamma5
- implemented color_trace()
- eps0123() creates a 4-dimensional epsilon tensor without checking the
  indices (so they may be D-dimensional); this may become a different class
  soon, but for now all the implemented properties of the standard epsilon
  tensor also apply to eps0123

added more trace checks

dirac_trace() can handle arbitrary gamma strings without gamma5. The only
remaining unimplemented case is now Tr[gamma5 S_2n]

added more trace checks

* Minor packaging fixes.
* Pushing towards version 0.8.2-1.

* Cint fix.  It looks like the hacks are going to go away really soon,
  after some recent work by Masaharu.

* Example: Don't confuse users with _ex0().

* Push towards 0.8.2.

- added get_representation_label()
- clifford contractions honor repr. label
- first implementation of dirac_trace(), works with 0..3 clifford objects

added first trace check

\18* class matrix: solve() doesn't call algorithms redundantly any more and
  inverse() falls back to solve() which works in more general cases.

* class matrix: solve() doesn't call algorithms redundantly any more and
  inverse() falls back to solve() which works in more general cases.

- degree(), ldegree(), coeff(), lcoeff(), tcoeff() and collect() work with
  non-symbols as the second argument
- collect_distributed() produces result in distributed form (second argument
  should be a list of symbols)

- color/clifford objects have representation label to distinguish elements
  of different algebras; objects with different labels commute with each
  other
- dirac_one() -> dirac_ONE()
- added gamma5 clifford objects which is constructed by dirac_gamma5()
- clifford::simplify_ncmul() anticommutes gamma5's to the front and removes
  squares of gamma5
- the argument to collect() can be a list of objects in which case the
  result is either a recursively collected polynomial, or a polynomial in
  a distributed form with terms like coeff*x1^e1*...*xn^en, as specified by
  the second argument to collect(). For example (ginsh):
    > f=a*x-x+a*x*y+x*y+a*x*y^2-2*x*y^2+y^2;
    a*x+y*a*x+y^2*a*x+y^2-x+y*x-2*y^2*x
    > collect(f,x);
    (-1+y+y^2*a+y*a+a-2*y^2)*x+y^2
    > collect(f,y);
    a*x-x+y^2*(1+a*x-2*x)+y*(a*x+x)
    > collect(f,[x,y]);
    (-1+y*(1+a)+y^2*(-2+a)+a)*x+y^2
    > collect(f,[y,x]);
    (-1+a)*x+y*(1+a)*x+(1+(-2+a)*x)*y^2
    > collect_distributed(f,[x,y]);
    (-1+a)*x+y^2*(-2+a)*x+y*(1+a)*x+y^2

- dirac_one() -> dirac_ONE()
- more clifford checks

documentation update

* fixed some mindboggling tremendous superhuge gigantic #*@$&! memory leaks.

* color.cpp: remove duplicated default args.
* numeric.cpp: insert some missing cln:: namespace decls.
  Why did this ever work?
* structure.cpp, structure.h: adjust to the new print scheme.

* Now this four-package split-up is gonna get in.  Promised.