Rational simplification bug ("Quotient by a polynomial of higher degree" error message)



On 2012-11-28 10:37-0500 Karl-Dieter Crisman wrote:

>>
>> Anyway, change that flag.  How to do it from Sage, I don't know.
>>
>
> sage: maxima_calculus('gcd:subres')
> subres
>
> should work, I *think*.

Thanks for that guidance. I got that to work (a zero result for
Resid[0]) for nbody=2 for my present sage platform that I built
myself.  But it was about two orders of magnitude slower than for the
case when I don't run maxima_calculus('gcd:subres').
(Resid[0].simplify_rational() took minutes rather than seconds to
simplify the nobody=2 expression to zero).

Anyhow, from those really slow results I am expecting the
simplification with maxima_calculus('gcd:subres') will take hours for
the nbody=3 case (if it works at all for that case) for my present
platform.  Therefore, I will report later about those results as well
as results for a downloaded binary version of sage-5.4.1 which
includes and is a front-end for a binary version of maxima-5.26.0.

Alan
__________________________
Alan W. Irwin

Astronomical research affiliation with Department of Physics and Astronomy,
University of Victoria (astrowww.phys.uvic.ca).

Programming affiliations with the FreeEOS equation-of-state
implementation for stellar interiors (freeeos.sf.net); the Time
Ephemerides project (timeephem.sf.net); PLplot scientific plotting
software package (plplot.sf.net); the libLASi project
(unifont.org/lasi); the Loads of Linux Links project (loll.sf.net);
and the Linux Brochure Project (lbproject.sf.net).
__________________________

Linux-powered Science
__________________________