def test_exp_basics(self):
"""Test exponential tf."""
rtf = ExpRTransform(0.1, 1e1)
gd = np.ones(100) * 0
assert abs(rtf.transform(gd)[0] - 0.1) < 1e-15
gd = np.ones(100) * 99
assert abs(rtf.transform(gd)[0] - 1e1) < 1e-10
self.check_consistency(rtf)
self.check_deriv(rtf)
def setUp(self):
"""Set up radial grid for integral tests."""
pts = HortonLinear(100)
tf = ExpRTransform(1e-3, 1e1)
rad_pts = tf.transform(pts.points)
rad_wts = tf.deriv(pts.points) * pts.weights
self.rgrid = OneDGrid(rad_pts, rad_wts)
def test_basics2():
"""Test basic radial grid transform properties for bigger grid."""
oned = UniformInteger(100)
rtf = ExpRTransform(1e-3, 1e1)
grid = rtf.transform_1d_grid(oned)
assert isinstance(grid, OneDGrid)
assert grid.size == 100
# assert grid.shape == (100,)
# assert grid.rtransform == rtf
assert (grid.weights > 0).all()
assert (grid.points == rtf.transform(oned.points)).all()
def test_basics1():
"""Test basic radial grid transform properties."""
oned = UniformInteger(4)
rtf = ExpRTransform(0.1, 1e1)
grid = rtf.transform_1d_grid(oned)
assert isinstance(grid, OneDGrid)
assert grid.size == 4
# assert grid.shape == (4,)
# assert grid.rtransform == rtf
assert (grid.weights > 0).all()
assert (grid.points == rtf.transform(oned.points)).all()
