Parent Directory
| 
Revision Log
| Links to HEAD: | (view) (annotate) |
| Links to release_1-0-10: | (view) (annotate) |
| Sticky Tag: |
* Finilize version 1.0.4 (version numbers, copyrights and such rubbish).
- Cleanups: My evil plot of making ex::bp private may finally be carried out, provided one changes all the is_of_type(obj,type) to is_a<type>(obj) which we don't do right now because it would degrade performance on gcc-2.95.x. Also, ex_to_type(obj) has gone for good now, we have been having ex_to<type>(obj) long enough. ex_to_nonconst_type(obj) have gone as well, since they are almost never used and one can always cast aways the constness explicitly if need should arise.
* Supplement some (now deprecated) macros by inlined template functions: - is_of_type(foo, type) -> is_a<type>(foo) - is_ex_of_type(foo, type) -> is_a<type>(foo) - is_exaclty_of_type(foo, type) -> is_exaclty_a<type>(foo) - is_ex_exaclty_of_type(foo, type) -> is_exaclty_a<type>(foo)
- first implementation of pattern matching
replaced "precedence" static member variable by virtual precedence() function
- 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
- renamed
function::getserial() -> function::get_serial()
symbol::setname() -> symbol::set_name()
symbol::getname() -> symbol::get_name()
for consistency (since the interface has changed anyway)
- added pseries::get_var() and pseries::get_point()
- replaced the various print*() member functions by a single print() that takes a print_context object that determines the output formatting; this should make it easier to add more output types - print_tree output of indexed objects looks better
degree(), ldegree(), coeff(), lcoeff(), tcoeff() and collect() can now be used with constants and indexed expressions as well, so you can use it to collect by powers of Pi or find the coefficient of gamma~0. Limitations: - it only works with symbols, constants and indexed expressions; trying to find the coefficient of, e.g., "x^2" or "x+y" won't work - it does not know about dummy index summations; the coefficient of gamma~0 in p.mu*gamma~mu should be p.0 but is returned as 0 - using the functions on elements of noncommutative products might return wrong or surprising results
* ginac/registrar.h: dtor is inlined now. * ginac/basic.h, ginac/ex.h: removed INLINE_FOO_CONSTRUCTORS and hand-inlined stuff instead, since the macros turned out to get in the way while fine-tuning the inlining. * ginac/utils.h: if 'long long' works, use that one for computing the hash value instead of floating point tricks. 2 Reasons: on Intel the assembler is terrible and slow otherwise, 'long long' will appear in C++ anyways when they adapt to the C99 standard. * Several other small performance tweaks. * ginac/constant.cpp: (constant::calchash()) implement, caring for serial. * ginac/function.pl (function::calchash()): implement, caring for serial. * ginac/expairseq.cpp: honor coeff's hash value. * ginac/pseries.cpp: (pseries::power_const()) fix problems with zero pseries. * Added several pounds of in-source documentation in doxygen style.
GINAC_DECLARE_REGISTERED_CLASS declares duplicate() and compare_same_type(), GINAC_IMPLEMENT_REGISTERED_CLASS implements duplicate()
the destructor, copy constructor, and assignment operator (which were the same for all subclasses of basic, with very few exceptions) are now included in the GINAC_IMPLEMENT_REGISTERED_CLASS macro; the GINAC_DECLARE_REGISTERED_CLASS macro also defines these (and other common) member functions
removed the "some_*" and "typeid_*" definitions since we are using our own RTTI for quite some time now
- See if __GNUC__ < 2.97 before using std::vector<..,malloc_alloc>. Sorry, there is no way having this depend on whether we are setting up Cint or not because the malloc_alloc goes into the mangled signature and the result won't link if one tries to be more clever. It really sucks. Let's all just pray Masaharu throws away that old STL implementation really soon now. (Jeez, he still calls this old junk 'reference'!) - Bumped up a few dates. :-)
- pseries::print(): did not insert parenthesis when needed for precedence. - pseries::precedence: was missing (thus being inherited from basic). - add::print(): was not honoring upper_precedence for rest.
- Inserted some more std:: to make it compile under GCC2.96.
- indentation is now done with tabs - printtree() prints the class names as known to the registrar, not the (mangled) class names from typeid().name()
- Fixed the bug that broke xloop's po_redux. - Interface Change: No argument 'bool branchcut' for series methods, instead 'unsigned options'. Sorry. : ----------------------------------------------------------------------
- Derivatives are now assembled in a slightly different manner (i.e. they 'look' different on first sight). Under certain circumstances this can result in a dramatic speedup because it gives hashing a better chance, especially when computing higher derivatives. - Series expansion accepts an additional bool argument now, telling it whether branch cuts are to be honored (as Maple does it) or not (as Mathematica does it). It defaults to true (i.e. the Maple way). - Many series expansions of built-in functions have been reengineered. - The Dilogarithm (Li2) now has floating point evaluation, derivative and a proper series expansion. - Namespace 'std' cleanly disentangled, as demanded by ISO/EIC 14882-1998(E). - Some helpers updated from upstream.
- pseries::print(): Make it print more nicely. - bernoulli(): A slightly slower but only half as memory consuming algorithm ripped from Pari.
- Overwritten: pseries::printtree().
- We now write f(x).series(x==3,5) instead of f(x).series(x,3,5) and f(x).series(x,4) instead of f(x).series(x,0,4). We also don't allow default arguments any more.
- collect() works better
- series expansion of pseries objects should now work as expected
- diff() is now only defined on classes ex and basic, where it handles higher derivatives. derivative() is now the method to be implemented by the user that does single derivatives, it is protected. - Implemented some very clumpsy automatic generation of ginaccint. It is ugly because we have to compile libginac without namespace GiNaC. - Fixed some wrong logic in numeric::power().
- Renamed flag NO_GINAC_NAMESPACE to NO_NAMESPACE_GINAC because of m4. - Made configure check for cint and makecint and call makecint to dump out cint/Makefile.cint, which will later be called by cint/Makefile. See configure option --with-cint. - Added some documentation about GiNaC-cint (e.g. a manpage). - All dummies in subdir cint/ are created by cint/dummies.pl now, removed older stuff.
- building GiNaC in a separate directory now works
- changed all instances of "foo const &/*" to "const foo &/*" - changed function arguments like "int const" to a simple "int"
- implemented nops() and op() for pseries objects
- implemented global class registry (for class basic and derived classes) - implemented archiving of expressions (except for class idx, and all classes derived from idx and indexed) - added series() wrapper function - class series renamed to pseries to avoid name clashes with global series() wrapper function, ex::series(), and basic::series() - corrected the series expansion of single symbols with respect to the order term - calling subs() on a series didn't work
This form allows you to request diffs between any two revisions of this file. For each of the two "sides" of the diff, select a symbolic revision name using the selection box, or choose 'Use Text Field' and enter a numeric revision.
| Christian Bauer">Christian Bauer | ViewVC Help |
| Powered by ViewVC 1.1.15 |