def sanity_ibtrapz(self):
"""
mtools: Checking ibtrapz
"""
from PyAstronomy.pyaC import mtools
import numpy as np
x = np.arange(-2., 2.01, 0.1)
y = 2. * x
x0 = -1.375
x1 = +1.943
# Analytical value of integral
analyt = x1**2 - x0**2
self.assertAlmostEqual(analyt,
mtools.ibtrapz(x, y, x0, x1),
delta=1e-10,
msg="ibtrapz incorrect for linear function.")
self.assertAlmostEqual(
(-1.9)**2 - (-2.0)**2,
mtools.ibtrapz(x, y, -2.0, -2.0 + 0.1),
delta=1e-10,
msg="ibtrapz incorrect for linear function (-2,-1.9).")
self.assertAlmostEqual(
0.0,
mtools.ibtrapz(x, y, -2.0, +2.0),
delta=1e-10,
msg="ibtrapz incorrect for linear function (-2,+2).")
def sanity_ibtrapz(self):
"""
mtools: Checking ibtrapz
"""
from PyAstronomy.pyaC import mtools
import numpy as np
x = np.arange(-2.0, 2.01, 0.1)
y = 2.0 * x
x0 = -1.375
x1 = +1.943
# Analytical value of integral
analyt = x1 ** 2 - x0 ** 2
self.assertAlmostEqual(
analyt, mtools.ibtrapz(x, y, x0, x1), delta=1e-10, msg="ibtrapz incorrect for linear function."
)
self.assertAlmostEqual(
(-1.9) ** 2 - (-2.0) ** 2,
mtools.ibtrapz(x, y, -2.0, -2.0 + 0.1),
delta=1e-10,
msg="ibtrapz incorrect for linear function (-2,-1.9).",
)
self.assertAlmostEqual(
0.0, mtools.ibtrapz(x, y, -2.0, +2.0), delta=1e-10, msg="ibtrapz incorrect for linear function (-2,+2)."
)
def sanity_ibtrapzExample(self):
"""
mtools: Checking example of ibtrapz
"""
from PyAstronomy.pyaC import mtools
import numpy as np
x = np.arange(-2.,2.01,0.1)
y = x**3 + 1.7
x0 = -1.375
x1 = +1.943
# Analytical value of integral
analyt = 0.25*(x1**4 - x0**4) + 1.7*(x1-x0)
print("Analytical value: ", analyt)
print("ibtrapz: ", mtools.ibtrapz(x, y, x0, x1))
def sanity_ibtrapzExample(self):
"""
mtools: Checking example of ibtrapz
"""
from PyAstronomy.pyaC import mtools
import numpy as np
x = np.arange(-2., 2.01, 0.1)
y = x**3 + 1.7
x0 = -1.375
x1 = +1.943
# Analytical value of integral
analyt = 0.25 * (x1**4 - x0**4) + 1.7 * (x1 - x0)
print("Analytical value: ", analyt)
print("ibtrapz: ", mtools.ibtrapz(x, y, x0, x1))