Parent Directory
| 
Revision Log
| Links to HEAD: | (view) (annotate) |
| Links to release_1-0-10: | (view) (annotate) |
| Sticky Tag: |
powers of indexed objects are now parenthesized correctly in LaTeX output
* power::expand_add(): allocate the precise amount needed for representing the result. And some minor readability cleanups...
- (l)degree(s), coeff(s, n) and collect(s) were extended to accept expressions of any class (except add/mul/ncmul/numeric) for "s". They should even work if "s" is a "power" object, as long as the exponent is non-integer, but with some limitations. For example, you can "collect(a*2^x+b*2^x, 2^x)" to get "(a+b)*2^x", but "degree(2^(3*x), 2^x)" yields 0 instead of 3).
* Finilize version 1.0.4 (version numbers, copyrights and such rubbish).
* Fix thinkos in print_csrc context. (Cebix: read the implementations of ex::compare() and basic::compare()!!!)
* Added output-support for Python bindings and LaTeX printing for indexed objects (Pearu Peterson). * relational.*: Removed superfluous is_equal_same_type()... * inifcns_trans.cpp: ...added extern declaration as workaround for compiler-bug triggered by changes to relational.*, * input_parser.yy: added missing namespace specification.
- input parser recognizes "sqrt()", which is also used in the output - divide(a,b,q) only modifies q if the division succeeds; also, divide(a,b,a) works now - fixed small bug in dummy index renaming which could cause it to not recognize renamable indices in some cases - power::degree() and power::ldegree() throw an exception when encountering a non-integer exponent
* As agreed upon at lunch, remove debugmsg since it is hardly of any use.
* Some internal reorganization WRT flyweight handling and initialization, resulting in a general speed-up. Wherever thou hast written _ex7() thou shalt write _ex7 from now on.
- 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.
- carried on with felonious plot about making ex::bp private.
- Fix the evalf(2^Pi) bug reported by Pearu.
- carried on with felonious plot about making ex::bp private.
fixed potential crash in power:(l)degree() when used with non-integer exponents
- Move several comments into doxygen comment blocks.
- power::eval(): for the case (b_n/b_d)^(e_n/e_d) we now check separately if we can compute numerator and denominator. This allows us to eval (3/8)^(1/3) to 1/2*3^(1/3) instead of holding it. For square roots this is still less clever than what MapleV does, but then again that system has the funny goof not to touch 8^(1/3) where we immediately eval() to plain 2; fun, fun, fun... Oh, and a little memory hole has been squished along the way, too.
- replaced the Derivative() function by a more resonable fderivative class; to see it in action, try "series(abs(x),x==0,3)" and "series(abs(x),x==y,3)" in ginsh with previous GiNaC releases and with this one - add::eval() throws an exception when the first term is non-commutative and the overall_coeff is non-zero - small fixes to print_tree output of indexed objects - ncmul::expand() now actually does expand something like "(dirac_gamma(mu)+1) *dirac_gamma(nu)" correctly (but the change to add::eval() should make this a "can't happen" case) - moved ToString() to tostring.h (fewer header dependencies) - ALL YOUR "BASES CLASSES" ARE BELONG TO US
* Ctors of class numeric are not explicit any more. All built-in callers for pseudofunctions are now templated and default to ex arguments which relaxes the need for explicit ctors.
some cleanups
- added find() (like has(), but returns list of all occurrences)
- added lst::sort() and lst::unique()
- status_flags::expanded is only taken into account when no expand options
are specified. This makes it possible to re-expand with other options.
- added expand_options::expand_function_args
- collect(foo, {}) doesn't crash any more
- introduced match_same_type() method which is slightly similar to
is_equal_same_type() but doesn't check subexpressions. Fixed a number of
match() bugs with this (e.g. a==b matched a!=b, and matrices with different
dimensions but the same number of elements could match). Only container
classes with additional member variables that have to be equal for a match
need to implement match_same_type().
some cleanups
removed a lot of superfluous const_cast<>()s
* Fix a bitch of a bug where 1<I (and all other relationals) returned true. This did not happen for the binary relational operators between numerics. Now the behaviour is more consistent (throwing exceptions when the difference becomes complex). * Some regression tests were actually wrong in the light of this bug. :-(
* Supplement some (now deprecated) macros by inlined template functions: - ex_to_foobar(baz) -> ex_to<foobar>(baz).
* 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)
- Instead of just totally symmetric or antisymmetric, complex symmetries can now be defined for indexed objects. Symmetries are described by a tree of "symmetry" objects that is constructed with the sy_none(), sy_symm(), sy_anti() and sy_cycl() functions. The symmetry of a function with respect to its arguments can also be defined (this is currently only used for the Beta function). - color_trace() and dirac_trace() can be applied to a more general class of expressions, e.g. using it on a relation will take the trace on both sides etc. - Generalized map() to take a function object instead of a function pointer. This allows passing an arbitrary number of additional state to the function being called. - The unarchiving functinos find_bool(), find_unsigned() and find_string() can take an additional "index" argument.
- ++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.
- 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()
- renamed ncpower() to ncpow() - canonicalize_clifford() uses to_rational() now - minor cleanups
- fix LaTeX-output bug reported by Stefan, remove obsolete has(matrix,ex).
- first implementation of pattern matching
replaced "precedence" static member variable by virtual precedence() function
- 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
* Methods of class ex which do absolutely nothing than type dispatch should be inlined for obvious reasons.
* basic::collec() never worked correctly on non-polynomials till now.
We simply correct for the missing terms, which makes it at least
algebraically correct (suggested by Pearu Peterson).
* Added support to print out LaTeX-style, based on some work by Stefan
Weinzierl:
- print.h: added a class print_latex().
- function.pl: added function_options::latex_name() because there was
no good way of setting that name before.
- inifncs*.cpp: use it.
- *.cpp: switch to see if print_context is actually print_latex and
then change the style.
* Fixed some stupid bugs in sections that were #define'd away for gcc-2.95.
- 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
subs() performs "syntactic substitution" as in Maple; you can substitute any object by any other object (even numerics) but only when the object subs() is called on exactly matches an object on the left-hand side of the substitution list; expairseq still doesn't substitute its numeric coefficients so substituting numerics will not replace them in all occurences (this shouldn't matter, though)
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
- added skeleton implementation of color and clifford classes (don't bother trying it, it doesn't work yet) - some macros in utils.h (DEFAULT_*) provide shorthands for implementing many required class member functions
- revamped indexed objects - subs() works on matrices
* 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
expand() always expands the exponent and transforms x^(a+b) -> x^a*x^b
- 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. :-)
- power::derivative(): slight optimization.
- non-integer powers are treated as 0 by (l)degree() and coeff()
- Inserted some more std:: to make it compile under GCC2.96.
- more indentation fixes
- indentation is now done with tabs - printtree() prints the class names as known to the registrar, not the (mangled) class names from typeid().name()
- Complete revamp of methods in class matrix. Some redundant (and poor) implementations of elimination schemes were thrown out. The code is now highly orthogonal, more flexible and much more efficient. The function lsolve() now falls back to matrix::solve(), which in turn calls different elimination strategies. matrix::fraction_free_elim() was thrown out for good. matrix::pivot() has now two integer arguments (ro and co) to allow pivoting when eliminating to echelon form. - Some new checks to guard against regressions in the new matrix code. - Some minor bug-fixes in documentation. - Bumped up version from 0.6.3 to 0.6.4. - ...and the usual changes triggered by other changes.
- 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.
- introduced info_flag::algebraic. - matrix::fraction_free_elimination(bool) finished and used by matrix::determinant(). - introduced some comments.
bug fixed in an incomplete assertion in power.cpp
- 0^a -> 0^a instead of 0.
- According to CLTL 0^I is undefined, 0^(I+epsilon) is 0 and 0^(I-epsilon) is an overflow. Now, this should be honored.
- The status_flags::expanded is now used on some occasions.
- Doxygen'ed a bunch of comments.
- More drastic performance boost on matrix stuff. - Addition of some comments.
- Changed behaviour of .evalf() to not evalf any exponents.
- cosmetic fix
- bugfix: 2^(3/2) does not throw any more. This was now secured in... - ...a new function paranoia_check10.
- 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.
- changed all instances of "foo const &/*" to "const foo &/*" - changed function arguments like "int const" to a simple "int"
- 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
- made nops() return unsigned instead of int - changed dates in copyright notices
- more logic on the trigonometric function stuff. - changed several occurences of numeric const & to const numeric &, which is the same, but doxygen kept being confused because declarations differed from implementations.
- Banned exZERO(), exONE(), exMINUSHALF() and all this from the interface. There is now _ex1() meaning 1, _ex_1() meaning -1, _ex1_2() meaning 1/2 and so on defined in utils.h and implemented in utils.cpp. Feel free to extend them as it pleases you but use them inside the library only. - Added more evaluations of trigonometric functions - Added series expansions for psi(x) and psi(n,x) at all their poles
- introduced info_flags::cinteger, info_flags::crational, info_flags::cinteger_polynomial, info_flags::crational_polynomial with intuitive behaviour. - extended documentation - made things like evlaf(zeta(3)) work - killed several bugs
#ifndef around namespace GiNaC { }
- ASSERT macro renamed to GINAC_ASSERT
- put everything in "GiNaC" namespace - converted some macros into inline functions - split debugmsg.h into assertion.h (public) and debugmsg.h (private)
- modified the comment blocks so the copyright message no longer appears in the reference manual (made by Doxygen)
- modified GiNaC headers to Alexander's liking
- enforced GiNaC coding standards :-)
- switched to automake build environment
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 |