[GiNaC-list] Sqrfree_parfrac Issue GiNaC

[Konstantin Leyde]

Hello,

My name is Konstantin Leyde, and I am a master student in physics at the ETH Zurich.

For my thesis project, I have been working with the GiNaC library, and I have encountered an issue with the partial fractioning (sqrfree_parfrac) routine:
Unless I misunderstood the way the sqrfree_parfrac routine works, it appears that for special cases, the output of the sqrfree_parfrac routine misses a numerical prefactor for one term, and thus returns the wrong answer.

Please find attached a pdf, that explains the problem in more detail.
In this pdf (as well as separately) we have attached the code that shows the issue.


Many thanks in advance,

Best wishes,
Konstantin

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.ginac.de/pipermail/ginac-list/attachments/20200722/2a5cab0a/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sqrfree_parfrac_GiNaC.pdf
Type: application/pdf
Size: 107504 bytes
Desc: sqrfree_parfrac_GiNaC.pdf
URL: <http://www.ginac.de/pipermail/ginac-list/attachments/20200722/2a5cab0a/attachment-0001.pdf>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: issue_GiNaC.cpp
URL: <http://www.ginac.de/pipermail/ginac-list/attachments/20200722/2a5cab0a/attachment-0001.ksh>