[GiNaC-list] Bug in numer_denom?
Alexei Sheplyakov
alexei.sheplyakov at gmail.com
Tue Nov 23 17:32:12 CET 2010
Hello,
On Thu, Nov 18, 2010 at 10:50:23PM +0100, Richard B. Kreckel wrote:
> Sorry, I see now that you are right. The old and new computed
> numerators differ by a factor 282901891126422365/highA. That factor
> certainly shouldn't contain a highA.
I think the new one is correct (also I've cross-checked it with Maxima).
> Also, computing the GCD of the original A and B segfaults over here.
Got a simpler example:
gcd(1/282901891126422365*(x + y), 165888/282901891126422365*(x - y));
I've posted the patch to ginac-devel.
Best regards,
Alexei.
More information about the GiNaC-list
mailing list