diff --git a/src/Li5.cpp b/src/Li5.cpp index 91a8e540..713c9412 100644 --- a/src/Li5.cpp +++ b/src/Li5.cpp @@ -118,9 +118,9 @@ std::complex Li5(const std::complex& z) if (az <= 1.) { u = -clog(1. - z); } else { // az > 1. - const std::complex lnz = clog(-z); - const std::complex lnz2 = pow2(lnz); - const std::complex lnz4 = pow2(lnz2); + const auto lnz = clog(-z); + const auto lnz2 = pow2(lnz); + const auto lnz4 = pow2(lnz2); u = -clog(1. - 1./z); rest = -1./360.*lnz*(7*PI4 + 10.*PI2*lnz2 + 3.*lnz4); } diff --git a/src/Li6.cpp b/src/Li6.cpp index 8b2c7824..4ac948e5 100644 --- a/src/Li6.cpp +++ b/src/Li6.cpp @@ -122,19 +122,19 @@ std::complex Li6(const std::complex& z) std::complex u, r; double sgn = 1; - if (std::abs(z) <= 1.) { + if (az <= 1.) { u = -clog(1. - z); } else { // az > 1. - const std::complex lnz = clog(-z); - const std::complex lnz2 = pow2(lnz); - const std::complex lnz4 = pow2(lnz2); - const std::complex lnz6 = lnz2*lnz4; + const auto lnz = clog(-z); + const auto lnz2 = pow2(lnz); + const auto lnz4 = pow2(lnz2); + const auto lnz6 = lnz2*lnz4; u = -clog(1. - 1./z); r = -31.*PI6/15120. - 7./720.*PI4*lnz2 - 1./144.*PI2*lnz4 - 1./720.*lnz6; sgn = -1; } - const std::complex sum = + const auto sum = u * (bf[0] + u * (bf[1] + u * (bf[2] +