tests/testthat/test-wsigma.R

In jacobi: Jacobi Theta Functions and Related Functions

test_that("wsigma values", {
  expect_equal(
    wsigma(1, g = c(12, -8)), 
    sin(1i*sqrt(3))/(1i*sqrt(3)*sqrt(exp(1)))
  )
  expect_equal(
    wsigma(2, g = c(1, 2i)), 
    1.8646253716-0.3066001355i
  )
  expect_equal(
    wsigma(1, omega = c(1, 1i)) / 2, 
    as.complex(2^(5/4)*sqrt(pi)*exp(pi/8) / gamma(1/4)^2)
  )
})