[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