|
1 | 1 | #include "../testing_utils.h" |
2 | | - |
3 | 2 | #include <complex> |
4 | 3 | #include <tuple> |
5 | 4 | #include <xsf/bessel.h> |
6 | 5 |
|
7 | | - |
8 | 6 | TEST_CASE("cyl_bessel_k gh-46", "[cyl_bessel_k][xsf_tests]") { |
9 | 7 | using test_case = std::tuple<double, std::complex<double>, std::complex<double>, double>; |
| 8 | + using std::complex; |
10 | 9 | auto [v, z, ref, rtol] = GENERATE( |
11 | | - test_case{0.0, std::complex{680.0, -1000.0}, |
12 | | - std::complex{1.901684871999608e-298, 1.713412341479591e-297}, 1e-13}, |
13 | | - test_case{0.0, std::complex{680.0, -680.0}, |
14 | | - std::complex{-4.553730032944803e-298, 1.878727010109855e-297}, 1e-13}, |
15 | | - test_case(0.0, std::complex{25.0, 100.0}, |
16 | | - std::complex{1.699267100365868e-12, -2.234890030902166e-13}, 5e-16) |
| 10 | + test_case{0.0, complex{680.0, -1000.0}, complex{1.901684871999608e-298, 1.713412341479591e-297}, 1e-13}, |
| 11 | + test_case{0.0, complex{680.0, -680.0}, complex{-4.553730032944803e-298, 1.878727010109855e-297}, 1e-13}, |
| 12 | + test_case{0.0, complex{25.0, 100.0}, complex{1.699267100365868e-12, -2.234890030902166e-13}, 5e-16} |
17 | 13 | ); |
18 | 14 | const std::complex w = xsf::cyl_bessel_k(v, z); |
19 | 15 | const auto rel_error = xsf::extended_relative_error(w, ref); |
|
0 commit comments