can you let me have the result of running the program below (*after* you've done an SVN update, as I've had to update the debugging/tracing code)
Here you go: .../expint.hpp:405 z = 0.5 .../expint.hpp:419 result = 0.028178840877963713110148 .../expint.hpp:422 result = 0.028178840877963713110148 .../expint.hpp:427 result = 0.028178840877963713110148 .../expint.hpp:428 max_iter = 0 0.028178840877963713110148 .../beta.hpp:872 a = 150.599945068359375 .../beta.hpp:873 b = 306722 .../beta.hpp:874 x = 0.91338449716567993164 .../beta.hpp:875 inv = 0 .../beta.hpp:876 normalised = 0 .../beta.hpp:1096 invert = 1 .../beta.hpp:1181 fract = 0 3.6860205160884932633e-565 .../detail/igamma_inverse.hpp:378 a = 16.18250274658203125 .../detail/igamma_inverse.hpp:379 p = 0.007907813991753274790426 .../detail/igamma_inverse.hpp:201 s = -2.41314250292881458062 .../detail/igamma_inverse.hpp:209 ra = 4.02274815848345535584 .../detail/igamma_inverse.hpp:216 w = 8.09144898812073984544 .../detail/igamma_inverse.hpp:278 z = 8.1669017875574152633 .../detail/igamma_inverse.hpp:291 result = 8.08379894264235578944 .../detail/igamma_inverse.hpp:393 guess = 8.08379894264235578944 Halley iteration, limit = 4.5474735088646411895e-13 Halley iteration, delta = -0.007452388058006306973197 Halley iteration, delta = -4.2590653528058823273e-08 Halley iteration, delta = nan Halley iteration, final count = 3 Maximum iterations: 3 .../detail/igamma_inverse.hpp:412 guess = nan nan .../bessel.hpp:340 v = 0 .../bessel.hpp:341 x = 1000 .../bessel.hpp:356 r = 0.004715917977622812778975 0.004715917977622812778975