Skip to content

Commit

Permalink
Update comment
Browse files Browse the repository at this point in the history
  • Loading branch information
kimwalisch committed Feb 10, 2024
1 parent 7b4a394 commit 3973527
Showing 1 changed file with 3 additions and 2 deletions.
5 changes: 3 additions & 2 deletions src/RiemannR.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -202,8 +202,9 @@ long double RiemannR_prime(long double x)

long double epsilon = std::numeric_limits<long double>::epsilon();

// RiemannR_prime(1) = NaN
// Hence we return RiemannR_prime(1.0000000000000001)
// RiemannR_prime(1) = NaN.
// Hence we return RiemannR_prime(1.0000000000000001).
// Required because: sum / log(1) = 0 / 0.
if (std::abs(x - 1.0) < epsilon)
return 0.60792710185402643042L;

Expand Down

0 comments on commit 3973527

Please sign in to comment.