BesselJ 0.1.0.0 → 0.1.0.1

raw patch · 4 files changed

+49/−8 lines, 4 filesPVP ok

version bump matches the API change (PVP)

API changes (from Hackage documentation)

Files

• BesselJ.cabal +1/−1
• CHANGELOG.md +5/−0
• README.md +2/−0
• tests/Main.hs +41/−7

@@ -1,6 +1,6 @@ cabal-version:       2.2
 name:                BesselJ
-version:             0.1.0.0
+version:             0.1.0.1
 synopsis:            Bessel J-function
 description:         Computation of the Bessel J-function of a complex variable. 
 homepage:            https://github.com/stla/BesselJ#readme

@@ -4,3 +4,8 @@ ## 0.1.0.0 - 2023-09-21  First release.+++## 0.1.0.1 - 2023-09-22++More unit tests.

@@ -8,3 +8,5 @@ ___
 
 ![](https://raw.githubusercontent.com/stla/BesselJ/main/images/BesselJ-nu3.png)
+
+![](https://raw.githubusercontent.com/stla/BesselJ/main/images/BesselJ_nu3-3i.png)

@@ -1,10 +1,9 @@ module Main where
-import           Approx               ( assertAreClose )
-import           Data.Complex         ( Complex(..) )
-import           Math.BesselJ
--- import           Math.Gamma           ( gamma )
-import           Test.Tasty           ( defaultMain, testGroup )
-import           Test.Tasty.HUnit     ( testCase )
+import           Approx           ( assertAreClose )
+import           Data.Complex     ( Complex(..), conjugate )
+import           Test.Tasty       ( defaultMain, testGroup )
+import           Test.Tasty.HUnit ( testCase )
+import           Math.BesselJ     ( BesselResult(_result), besselJ )
 
 i_ :: Complex Double
 i_ = 0.0 :+ 1.0
@@ -42,6 +41,41 @@       let z = 3 :+ 4
           s = sqrt(2 / pi / z) * sin z
       x <- _result <$> besselJ 0.5 z 1e-5 5000
-      assertAreClose "" 1e-9 x s
+      assertAreClose "" 1e-9 x s,
+
+    testCase "remove square root" $ do
+      let z = 2 :+ 1
+          nu = (-0.3) :+ 1
+      x <- _result <$> besselJ nu (sqrt (z*z)) 1e-5 5000
+      y <- _result <$> besselJ nu z 1e-5 5000
+      assertAreClose "" 1e-6 x (z**(-nu) * (z*z)**(nu/2) * y),
+
+    testCase "remove square root --- integer nu" $ do
+      let z = 2 :+ 1
+          nu = -4
+      x <- _result <$> besselJ nu (sqrt (z*z)) 1e-5 5000
+      y <- _result <$> besselJ nu z 1e-5 5000
+      assertAreClose "" 1e-6 x (z**(-nu) * (z*z)**(nu/2) * y),
+
+    testCase "remove minus sign" $ do
+      let z = 2 :+ 1
+          nu = (-0.3) :+ 1
+      x <- _result <$> besselJ nu (-z) 1e-5 5000
+      y <- _result <$> besselJ nu z 1e-5 5000
+      assertAreClose "" 1e-6 x ((-z)**nu * z**(-nu) * y), 
+
+    testCase "conjugate" $ do
+      let z = 2 :+ 5
+          nu = 0.3 :+ (-1)
+      x <- _result <$> besselJ (conjugate nu) (conjugate z) 1e-5 5000
+      y <- _result <$> besselJ nu z 1e-5 5000
+      assertAreClose "" 1e-6 x (conjugate y), 
+
+    testCase "conjugate --- integer nu" $ do
+      let z = 2 :+ 5
+          nu = 7
+      x <- _result <$> besselJ nu (conjugate z) 1e-5 5000
+      y <- _result <$> besselJ nu z 1e-5 5000
+      assertAreClose "" 1e-6 x (conjugate y)
 
   ]