Merge pull request #6 from Cactus-proj/dev

test: init all tests
Cactus-proj · Dec 9, 2024 · 8e5e64f · 8e5e64f
2 parents b7b7754 + 412cd82
commit 8e5e64f
Show file tree

Hide file tree

Showing 27 changed files with 184 additions and 12 deletions.
diff --git a/README.md b/README.md
@@ -15,12 +15,18 @@ Therefore, contributions regarding precision and error-related contributions are
 > which is a header-only library containing a modernized (C++17) refactoring of
 > the classic `AMOS/SPECFUN/Cephes/...` special function implementations.
 
+
 ## Docs
 
 [List of Cephes functions](doc/markdown/index.md)
 
-- [bessel](doc/markdown/bessel.md)
-- [exp_int](doc/markdown/exp_int.md)
+- [Exponential integral](doc/markdown/exp_int.md)
+- [Gamma functions](doc/markdown/gamma.md)
+- [Error functions](doc/markdown/gamma.md)
+- [Bessel functions](doc/markdown/bessel.md)
+- [Hypergeometric functions](doc/markdown/hyper.md)
+- [Elliptic functions](doc/markdown/elliptic.md)
+- [Miscellaneous functions](doc/markdown/misc.md)
 
 
 ## Test
@@ -31,6 +37,7 @@ cmake -DCMAKE_BUILD_TYPE=Coverage -S . -B build && cmake --build build
 cd build/ && ctest && make coverage_html
 ```
 
+
 ## License
 
 Using a `BSD-3-Clause` like [LICENSE](LICENSE.txt).

diff --git a/doc/markdown/bessel.md b/doc/markdown/bessel.md
@@ -1,4 +1,4 @@
-# cephes/bessel
+# Bessel functions
 
 ## Airy Functions
 

diff --git a/doc/markdown/exp_int.md b/doc/markdown/exp_int.md
@@ -1,4 +1,4 @@
-# cephes/Exponential integral
+# Exponential integral
 
 * **expn**, [Exponential integral En](doubldoc.md#expn)
 

diff --git a/doc/markdown/hyper.md b/doc/markdown/hyper.md
@@ -1,6 +1,4 @@
-# cephes/hyper
-
-## Hypergeometric Functions
+# Hypergeometric Functions
 
 * **hyp2f1**, [Gauss hypergeometric function](doubldoc.md#hyp2f1)
 * **hyperg**, [Confluent hypergeometric function](doubldoc.md#hyperg)
diff --git a/doc/markdown/misc.md b/doc/markdown/misc.md
@@ -1,4 +1,4 @@
-# cephes/misc
+# Miscellaneous functions
 
 * **spence**, [Dilogarithm](doubldoc.md#spence)
 * **zeta**, [Zeta function of two arguments](doubldoc.md#zeta)

diff --git a/tests/CMakeLists.txt b/tests/CMakeLists.txt
@@ -24,8 +24,8 @@ endfunction()
 
 add_subdirectory(exp_int)
 add_subdirectory(gamma)
-# add_subdirectory(error)
+add_subdirectory(error)
 add_subdirectory(bessel)
 add_subdirectory(hyper)
-# add_subdirectory(elliptic)
+add_subdirectory(elliptic)
 add_subdirectory(misc)
diff --git a/tests/elliptic/CMakeLists.txt b/tests/elliptic/CMakeLists.txt
@@ -0,0 +1,7 @@
+
+add_gtest(ellpk)
+add_gtest(ellik)
+add_gtest(ellpe)
+add_gtest(ellie)
+
+add_gtest(ellpj)
diff --git a/tests/elliptic/ellie.cpp b/tests/elliptic/ellie.cpp
@@ -0,0 +1,8 @@
+#include <xtest.hpp>
+#include <cephes/elliptic.h>
+
+
+TEST(EllipticEInc, Branches) {
+    // m == 0.0
+    EXPECT_REL_NEAR_F64(cephes::ellie(3.0, 0.0), 3.0);
+}
diff --git a/tests/elliptic/ellik.cpp b/tests/elliptic/ellik.cpp
@@ -0,0 +1,8 @@
+#include <xtest.hpp>
+#include <cephes/elliptic.h>
+
+
+TEST(EllipticKInc, Branches) {
+    // m == 0.0
+    EXPECT_REL_NEAR_F64(cephes::ellik(3.0, 0.0), 3.0);
+}
diff --git a/tests/elliptic/ellpe.cpp b/tests/elliptic/ellpe.cpp
@@ -0,0 +1,8 @@
+#include <xtest.hpp>
+#include <cephes/elliptic.h>
+
+
+TEST(EllipticE, Branches) {
+    // x == 0.0
+    EXPECT_REL_NEAR_F64(cephes::ellpe(0.0), 1.0);
+}
diff --git a/tests/elliptic/ellpj.cpp b/tests/elliptic/ellpj.cpp
@@ -0,0 +1,11 @@
+#include <xtest.hpp>
+#include <cephes/elliptic.h>
+
+
+TEST(EllipticJaco, Branches) {
+    double u, m, sn, cn, dn, ph;
+    // m < 0.0 || m > 1.0
+    u = 1.0;
+    EXPECT_REL_NEAR_F64(cephes::ellpj(u, -1.0, &sn, &cn, &dn, &ph), -1);
+    EXPECT_REL_NEAR_F64(cephes::ellpj(u, 2.0, &sn, &cn, &dn, &ph), -1);
+}
diff --git a/tests/elliptic/ellpk.cpp b/tests/elliptic/ellpk.cpp
@@ -0,0 +1,9 @@
+#include <xtest.hpp>
+#include <cephes/elliptic.h>
+
+
+TEST(EllipticK, Branches) {
+    // (x < 0.0) || (x > 1.0)
+    EXPECT_REL_NEAR_F64(cephes::ellpk(-1.0), 0.0);
+    EXPECT_REL_NEAR_F64(cephes::ellpk(2.0), 0.0);
+}
diff --git a/tests/error/CMakeLists.txt b/tests/error/CMakeLists.txt
@@ -0,0 +1,7 @@
+
+add_gtest(erf)
+add_gtest(erfc)
+add_gtest(ndtr)
+
+add_gtest(dawsn)
+add_gtest(fresnl)
diff --git a/tests/error/dawsn.cpp b/tests/error/dawsn.cpp
@@ -0,0 +1,7 @@
+#include <xtest.hpp>
+#include <cephes/error.h>
+
+
+TEST(Dawson, Branches) {
+    EXPECT_GT(cephes::dawsn(1.0), 0.0);
+}
diff --git a/tests/error/erf.cpp b/tests/error/erf.cpp
@@ -0,0 +1,8 @@
+#include <xtest.hpp>
+#include <cephes/error.h>
+
+
+TEST(Erf, Branches) {
+    // fabs(x) > 1.0
+    EXPECT_REL_NEAR_F64(cephes::erf(20.0), 1.0);
+}
diff --git a/tests/error/erfc.cpp b/tests/error/erfc.cpp
@@ -0,0 +1,8 @@
+#include <xtest.hpp>
+#include <cephes/error.h>
+
+
+TEST(ErfC, Branches) {
+    // z < -MAXLOG && a < 0
+    EXPECT_REL_NEAR_F64(cephes::erfc(-20.0), 2.0);
+}
diff --git a/tests/error/fresnl.cpp b/tests/error/fresnl.cpp
@@ -0,0 +1,9 @@
+#include <xtest.hpp>
+#include <cephes/error.h>
+
+
+TEST(Fresnel, Branches) {
+    double x, s, c;
+
+    EXPECT_REL_NEAR_F64(cephes::fresnl(1.0, &s, &c), 0.0);
+}
diff --git a/tests/error/ndtr.cpp b/tests/error/ndtr.cpp
@@ -0,0 +1,7 @@
+#include <xtest.hpp>
+#include <cephes/error.h>
+
+
+TEST(ndtr, Branches) {
+    EXPECT_GT(cephes::ndtr(0.1), 0.0);
+}
diff --git a/tests/gamma/CMakeLists.txt b/tests/gamma/CMakeLists.txt
@@ -6,8 +6,8 @@ add_gtest(lgam)
 add_gtest(igam)
 add_gtest(igamc)
 add_gtest(igami)
-# add_gtest(incbet)
-# add_gtest(incbi)
+add_gtest(incbet)
+add_gtest(incbi)
 
 add_gtest(psi)
 add_gtest(rgamma)

diff --git a/tests/gamma/incbet.cpp b/tests/gamma/incbet.cpp
@@ -0,0 +1,12 @@
+#include <xtest.hpp>
+#include <cephes/gamma.h>
+
+
+TEST(BetaInc, BasicAssertions) {
+
+}
+TEST(BetaInc, Branches) {
+    // aa <= 0.0 || bb <= 0.0
+    EXPECT_REL_NEAR_F64(cephes::incbet(0.0, 1.0, 1.0), 0.0);
+    EXPECT_REL_NEAR_F64(cephes::incbet(1.0, 0.0, 1.0), 0.0);
+}
diff --git a/tests/gamma/incbi.cpp b/tests/gamma/incbi.cpp
@@ -0,0 +1,11 @@
+#include <xtest.hpp>
+#include <cephes/gamma.h>
+
+
+TEST(BetaIncInv, BasicAssertions) {
+
+}
+TEST(BetaIncInv, Branches) {
+    // yy0 <= 0
+    EXPECT_REL_NEAR_F64(cephes::incbi(1.0, 1.0, -1.0), 0.0);
+}
diff --git a/tests/hyper/CMakeLists.txt b/tests/hyper/CMakeLists.txt
@@ -1,2 +1,3 @@
 
 add_gtest(hyp2f1)
+add_gtest(hyperg)
diff --git a/tests/hyper/hyperg.cpp b/tests/hyper/hyperg.cpp
@@ -0,0 +1,7 @@
+#include <xtest.hpp>
+#include <cephes/hyper.h>
+
+
+TEST(Hyp2f1, Branches) {
+    EXPECT_GT(cephes::hyperg(1.0, 2.0, 3.0), 0.0);
+}
diff --git a/tests/misc/CMakeLists.txt b/tests/misc/CMakeLists.txt
@@ -1,2 +1,7 @@
 
+add_gtest(spence)
+
+add_gtest(zeta)
+add_gtest(zetac)
+
 add_gtest(struve)
diff --git a/tests/misc/spence.cpp b/tests/misc/spence.cpp
@@ -0,0 +1,14 @@
+#include <xtest.hpp>
+#include <cephes/misc.h>
+
+/*
+
+*/
+TEST(Spence, Branches) {
+    // x < 0
+    EXPECT_REL_NEAR_F64(cephes::spence(-1.0), 0.0);
+    // x == 1.0
+    EXPECT_REL_NEAR_F64(cephes::spence(1.0), 0.0);
+    // x == 0.0
+    // EXPECT_REL_NEAR_F64(cephes::spence(0.0), pi*pi/6.0);
+}
diff --git a/tests/misc/zeta.cpp b/tests/misc/zeta.cpp
@@ -0,0 +1,10 @@
+#include <xtest.hpp>
+#include <cephes/misc.h>
+
+/*
+
+*/
+TEST(Zeta, Branches) {
+    // x == 1.0
+    EXPECT_GT(cephes::zeta(1.0, 1.0), 1e308);
+}
diff --git a/tests/misc/zetac.cpp b/tests/misc/zetac.cpp
@@ -0,0 +1,10 @@
+#include <xtest.hpp>
+#include <cephes/misc.h>
+
+/*
+
+*/
+TEST(ZetaC, Branches) {
+    // x < -170.6243
+    EXPECT_REL_NEAR_F64(cephes::zetac(-175.0), 0.0);
+}