Beispiel #1
0
 def test_files(self):
     snaps = Snapshots(output_name='Pressure', dformat="point")
     snaps.append("tests/test_datasets/matlab_00.vtk")
     snaps.append("tests/test_datasets/matlab_01.vtk")
     expected_files = [
         "tests/test_datasets/matlab_00.vtk",
         "tests/test_datasets/matlab_01.vtk"
     ]
     assert snaps.files == expected_files
Beispiel #2
0
    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()
Beispiel #3
0
 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)
Beispiel #4
0
 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, )
Beispiel #5
0
 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)
Beispiel #7
0
 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")
Beispiel #8
0
 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)
Beispiel #9
0
    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
Beispiel #10
0
 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")
Beispiel #11
0
 def test_dimension(self):
     snaps = Snapshots(output_name='Pressure', dformat="point")
     snaps.append("tests/test_datasets/matlab_00.vtk")
     snaps.append("tests/test_datasets/matlab_01.vtk")
     assert snaps.dimension == 2500
Beispiel #12
0
 def test_append2(self):
     snaps = Snapshots(output_name='Pressure', dformat="point")
     snaps.append("tests/test_datasets/matlab_00.vtk")
     snaps.append("tests/test_datasets/matlab_01.vtk")
     assert snaps.weighted.shape == (2500, 2)
Beispiel #13
0
 def test_snapshots(self):
     snaps = Snapshots(output_name='Pressure', dformat="point")
Beispiel #14
0
 def test_append_wrongname(self):
     snaps = Snapshots(output_name='Pressure', dformat="point")
     with self.assertRaises(TypeError):
         snaps.append(3.6)
Beispiel #15
0
 def test_append_weights(self):
     snaps = Snapshots(
         output_name='Weights', weight_name='Weights', dformat="cell")
     snaps.append("tests/test_datasets/matlab_00.vtk")
     assert snaps.weights.shape == (2401, 1)