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()