def test_vol_prefactor():
assert nestle.vol_prefactor(1) == 2.
assert nestle.vol_prefactor(2) == math.pi
assert nestle.vol_prefactor(3) == 4. / 3. * math.pi
assert nestle.vol_prefactor(4) == 1. / 2. * math.pi**2
assert nestle.vol_prefactor(5) == 8. / 15. * math.pi**2
assert nestle.vol_prefactor(9) == 32. / 945. * math.pi**4
def test_vol_prefactor():
assert nestle.vol_prefactor(1) == 2.
assert nestle.vol_prefactor(2) == math.pi
assert nestle.vol_prefactor(3) == 4./3. * math.pi
assert nestle.vol_prefactor(4) == 1./2. * math.pi**2
assert nestle.vol_prefactor(5) == 8./15. * math.pi**2
assert nestle.vol_prefactor(9) == 32./945. * math.pi**4
def test_ellipsoid_sphere():
"""Test that Ellipsoid works like a sphere when ``a`` is proportional to
the identity matrix."""
scale = 5.
for n in range(1, NMAX + 1):
ctr = 2.0 * scale * np.ones(n) # arbitrary non-zero center
a = 1.0 / scale**2 * np.identity(n)
ell = nestle.Ellipsoid(ctr, a)
assert_allclose(ell.vol, nestle.vol_prefactor(n) * scale**n)
assert_allclose(ell.axlens, scale * np.ones(n))
assert_allclose(ell.axes, scale * np.identity(n))
def test_ellipsoid_sphere():
"""Test that Ellipsoid works like a sphere when ``a`` is proportional to
the identity matrix."""
scale = 5.
for n in range(1, NMAX+1):
ctr = 2.0 * scale * np.ones(n) # arbitrary non-zero center
a = 1.0 / scale**2 * np.identity(n)
ell = nestle.Ellipsoid(ctr, a)
assert_allclose(ell.vol, nestle.vol_prefactor(n) * scale**n)
assert_allclose(ell.axlens, scale * np.ones(n))
assert_allclose(ell.axes, scale * np.identity(n))