Parent Directory
| 
Revision Log
| Links to HEAD: | (view) (annotate) |
| Links to release_0-9-0: | (view) (annotate) |
| Sticky Tag: |
normal.cpp, normal.h: Fix Yun's algorithm (there was a mistake in Geddes et al.) for square-free factorization. Also change the syntax to allow for square-free factoring in a recursive way in Q[X]. This is an example how it works in Z[x,y,z] in ginsh: > ?sqrfree sqrfree(expression [, symbol-list]) - square-free factorization of a polynomial > sqrfree(expand((x+y)*(x+y+z)^2*(y-z)^3)); (x+y)*(-x-y-z)^2*(y^3+3*y*z^2-3*y^2*z-z^3) > sqrfree(expand((x+y)*(x+y+z)^2*(y-z)^3),[y]); (x+y)*(y-z)^3*(x+y+z)^2 > sqrfree(expand((x+y)*(x+y+z)^2*(y-z)^3),[z]); (x+y)*(y-z)^3*(x+y+z)^2 > sqrfree(expand((x+y)*(x+y+z)^2*(y-z)^3),[x]); (x+y)*(-x-y-z)^2*(y^3+3*y*z^2-3*y^2*z-z^3) > sqrfree(expand((x+y)*(x+y+z)^2*(y-z)^3),[x,y]); (x+y)*(-x-y-z)^2*(y^3+3*y*z^2-3*y^2*z-z^3) > sqrfree(expand((x+y)*(x+y+z)^2*(y-z)^3),[x,y,z]); (x+y)*(-x-y-z)^2*(y^3+3*y*z^2-3*y^2*z-z^3) > sqrfree(expand((x+y)*(x+y+z)^2*(y-z)^3),[v]); 3*x*y^5+3*x^2*y^4+x^3*y^3+y^6-y*z^5-2*y^4*z^2+y^2*z^4-y^5*z+2*y^3*z^3-3*x^3*y^2*z+3*x^3*y*z^2-5*x*y^4*z-x*y*z^4+3*x^2*y*z^3+6*x*y^2*z^3+3*x^2*y^2*z^2-2*x*y^3*z^2-7*x^2*y^3*z-x^3*z^3-2*x^2*z^4-x*z^5
* 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.
- 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. :-)
- added ex::to_rational() to convert general expression to rational expression
by replacing all non-rational parts with temporary symbols, e.g.:
ex a = pow(sin(x), 2) - pow(cos(x), 2);
ex b = sin(x) + cos(x);
ex d;
lst l;
divide(a.to_rational(l), b.to_rational(l), d);
cout << d.subs(l) << endl;
will print "sin(x)-cos(x)"
- 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.
- made nops() return unsigned instead of int - changed dates in copyright notices
#ifndef around namespace GiNaC { }
- 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 |