Skip to content

Commit

Permalink
Correct bugs and link to the right library for testing Kummer.
Browse files Browse the repository at this point in the history
  • Loading branch information
astamm authored and ocommowi committed Mar 24, 2020
1 parent 61ce58f commit 3efa583
Show file tree
Hide file tree
Showing 2 changed files with 7 additions and 6 deletions.
11 changes: 6 additions & 5 deletions Anima/math-tools/special_functions/animaKummerFunctions.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -129,22 +129,23 @@ double OneHalfLaguerreFunction(const double &x)
double bessel0Value = arf_get_d(arb_midref(bessel0Ball), MPFR_RNDN);

precision = 64;
arb_hypgeom_bessel_i_scaled(bessel1Ball, oneBall, inputBall, precision)
arb_hypgeom_bessel_i_scaled(bessel1Ball, oneBall, inputBall, precision);

while (arb_can_round_arf(bessel1Ball, 53, MPFR_RNDN) == 0)
{
precision *= 2;
arb_hypgeom_bessel_i_scaled(bessel1Ball, oneBall, inputBall, precision)
arb_hypgeom_bessel_i_scaled(bessel1Ball, oneBall, inputBall, precision);
}

double bessel1Value = arf_get_d(arb_midref(bessel1Ball), MPFR_RNDN);

double resVal = (1.0 - x) * bessel0Value - x * bessel1Value;

arb_clear(inputBall);
arb_clear(outputBall);
arb_clear(aBall);
arb_clear(bBall);
arb_clear(bessel0Ball);
arb_clear(bessel1Ball);
arb_clear(zeroBall);
arb_clear(oneBall);

return resVal;

Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,7 @@ add_executable(${PROJECT_NAME}
## #############################################################################

target_link_libraries(${PROJECT_NAME}
AnimaOptimizers
AnimaSpecialFunctions
)


Expand Down

0 comments on commit 3efa583

Please sign in to comment.