Skip to content

Resolve underflow bug for high-w expansion

Wuttke, Joachim requested to merge higw into main

When kwwc_hig returned a result that was finite in extended precision, but zero in regular double precision, then kwwc fell back to integration, which then failed.

Example reported by Markus Appel:

$ ./runkww 0 c a 1.89 6.3e111
4.940656458412465442e-324 3      1
$ ./runkww 0 c a 1.89 6.4e111
kwwc: numeric integration failed for omega= 6.40000000000000015e+111, beta=       1.8899999999999999; error code -9

Solution: kwwc now analysing return values from series in extended precision.

At this occasion, split source kww_lowlevel from kww.

Merge request reports