BesselJ-0.1.0.1: tests/Main.hs

module Main where
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


main :: IO ()
main = defaultMain $
  testGroup "Tests"
  [ 
    testCase "nu = 1+2i -- z = 3+4i" $ do
      my <- _result <$> besselJ (1 :+ 2) (3 :+ 4) 1e-5 5000
      let wolfram = 0.31925 :+ (-0.66956) 
      assertAreClose "" 1e-5 my wolfram,

    -- testCase "Relation Bessel-J" $ do
    --   let nu = 0.5 :+ 2
    --       z = 3 :+ 4
    --       y = 2 * sin (nu * pi) / (pi * z)
    --   x1 <- _result <$> jnu (nu-1) z 
    --   x2 <- _result <$> jnu (-nu) z
    --   x3 <- _result <$> jnu (1-nu) z
    --   x4 <- _result <$> jnu nu z 
    --   assertAreClose "" 1e-3 (x1*x2 + x3*x4) y

    testCase "recurrence relation" $ do
      let nu = 0.5 :+ 2
          z = 3 :+ 4
      x1 <- _result <$> besselJ nu z 1e-5 5000
      x2 <- _result <$> besselJ (nu+1) z 1e-5 5000
      x3 <- _result <$> besselJ (nu+2) z 1e-5 5000
      let y = 2*(nu+1)/z * x2 - x3
      assertAreClose "" 1e-7 x1 y,

    testCase "elementary equality" $ do
      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,

    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)

  ]