def test_porous_random_seeding(self):
'''Test the reproducibility of the meshing algorithm.'''
L = np.array([2.0, 2.0])
R = 0.1
dx = 0.1
rs = np.array([
[-0.5, 0.5],
[0.0, 0.0],
[0.1, 0.4],
[0.88, 0.88],
[0.88, 0.0],
])
np.random.seed(2)
mesh_1 = mesh.porous_mesh_factory(rs, R, dx, L)
np.random.seed(3)
mesh_2 = mesh.porous_mesh_factory(rs, R, dx, L)
self.assertTrue(np.allclose(mesh_1.cellCenters.value,
mesh_2.cellCenters.value))
def get_mesh(self, L, dx):
return mesh.porous_mesh_factory(self.rs, self.R, dx, L)
