def test_dimension_2(self):
pts = Points()
pts.append([0, 1, 2])
pts.append([3, 4, 5])
pts.append([4, 4, 5])
pts.append([6, 4, 5])
pts.append([6, 7, 5])
assert isinstance(pts.triangulation, scipy.spatial.qhull.Delaunay)
def __init__(self,
output_name,
space_type=PODInterpolation,
weight_name=None,
dformat='cell'):
self.mu = Points()
self.snapshots = Snapshots(output_name, weight_name, dformat)
self.space = space_type()
def test_loo_error2(self):
mu = Points()
snap = Snapshots(output_name="Pressure", dformat="point")
space = Pod()
mu.append([-.5, -.5])
mu.append([.5, -.5])
mu.append([.5, .5])
mu.append([-.5, .5])
snap.append("tests/test_datasets/matlab_00.vtk")
snap.append("tests/test_datasets/matlab_01.vtk")
snap.append("tests/test_datasets/matlab_02.vtk")
snap.append("tests/test_datasets/matlab_03.vtk")
error = space.loo_error(mu, snap)
np.testing.assert_almost_equal(max(error), 0.149130165577, decimal=4)
def test_loo_error(self):
mu = Points()
snap = Snapshots(output_name="Pressure", dformat="point")
space = Pod()
mu.append([-.5, -.5])
mu.append([.5, -.5])
mu.append([.5, .5])
mu.append([-.5, .5])
snap.append("tests/test_datasets/matlab_00.vtk")
snap.append("tests/test_datasets/matlab_01.vtk")
snap.append("tests/test_datasets/matlab_02.vtk")
snap.append("tests/test_datasets/matlab_03.vtk")
error = space.loo_error(mu, snap)
assert error.shape == (4, )
def test_generate(self):
mu = Points()
snap = Snapshots(output_name="Pressure", dformat="point")
space = Pod()
mu.append([-.5, -.5])
mu.append([.5, -.5])
mu.append([.5, .5])
mu.append([-.5, .5])
snap.append("tests/test_datasets/matlab_00.vtk")
snap.append("tests/test_datasets/matlab_01.vtk")
snap.append("tests/test_datasets/matlab_02.vtk")
snap.append("tests/test_datasets/matlab_03.vtk")
space.generate(mu, snap)
assert space.state['pod_basis'].shape == (2500, 4)
def test_interpolator(self):
mu = Points()
snap = Snapshots(output_name="Pressure", dformat="point")
space = PODInterpolation()
mu.append([-.5, -.5])
mu.append([.5, -.5])
mu.append([.5, .5])
mu.append([-.5, .5])
snap.append("tests/test_datasets/matlab_00.vtk")
snap.append("tests/test_datasets/matlab_01.vtk")
snap.append("tests/test_datasets/matlab_02.vtk")
snap.append("tests/test_datasets/matlab_03.vtk")
space.generate(mu, snap)
assert isinstance(space.interpolator, LinearInterpolator)
def test_save(self):
mu = Points()
snap = Snapshots(output_name="Pressure", dformat="point")
space = Pod()
mu.append([-.5, -.5])
mu.append([.5, -.5])
mu.append([.5, .5])
mu.append([-.5, .5])
snap.append("tests/test_datasets/matlab_00.vtk")
snap.append("tests/test_datasets/matlab_01.vtk")
snap.append("tests/test_datasets/matlab_02.vtk")
snap.append("tests/test_datasets/matlab_03.vtk")
space.generate(mu, snap)
space.save("tests/test_datasets/podspace")
assert os.path.isfile("tests/test_datasets/podspace")
def test_call(self):
mu = Points()
snap = Snapshots(output_name="Pressure", dformat="point")
space = Pod()
mu.append([-.5, -.5])
mu.append([.5, -.5])
mu.append([.5, .5])
#mu.append([-.5, .5])
snap.append("tests/test_datasets/matlab_00.vtk")
snap.append("tests/test_datasets/matlab_01.vtk")
snap.append("tests/test_datasets/matlab_02.vtk")
#snap.append("tests/test_datasets/matlab_03.vtk")
space.generate(mu, snap)
solution = space([0, 0])
assert solution.shape == (2500, 1)
def __init__(self,
output_name,
space_type=PODInterpolation,
weight_name=None,
dformat='cell',
truncate=None,
interpolator=LinearNDInterpolator,
smoothness=0.0):
self.mu = Points()
self.snapshots = Snapshots(output_name, weight_name, dformat)
self.space = space_type()
self.truncate = truncate
self.interpolator = interpolator
self.smoothness = smoothness
def test_load(self):
mu = Points()
snap = Snapshots(output_name="Pressure", dformat="point")
space = PODInterpolation()
mu.append([-.5, -.5])
mu.append([.5, -.5])
mu.append([.5, .5])
mu.append([-.5, .5])
snap.append("tests/test_datasets/matlab_00.vtk")
snap.append("tests/test_datasets/matlab_01.vtk")
snap.append("tests/test_datasets/matlab_02.vtk")
snap.append("tests/test_datasets/matlab_03.vtk")
space.generate(mu, snap)
space.save("tests/test_datasets/podspace")
another_space = ParametricSpace.load("tests/test_datasets/podspace")
assert another_space.pod_basis.shape == (2500, 4)
os.remove("tests/test_datasets/podspace")
def test_dimension_1(self):
pts = Points()
pts.append([0, 1, 2])
pts.append([3, 4, 5])
assert pts.dimension == 3
def test_size(self):
pts = Points()
pts.append([0, 1, 2])
pts.append([3, 4, 5])
assert pts.size == 2
def test_append2(self):
pts = Points()
pts.append([0, 1, 2])
pts.append([3, 4, 5])
expected_values = np.array([[0, 1, 2], [3, 4, 5]]).T
np.testing.assert_array_almost_equal(pts.values, expected_values)
def test_points(self):
pts = Points()
