Source code for sfepy.tests.test_mesh_expand

[docs] def test_mesh_expand(): import numpy as nm from sfepy.mesh.mesh_tools import expand2d from sfepy.discrete.fem.mesh import Mesh from sfepy import data_dir mesh2d = Mesh.from_file(data_dir + '/meshes/2d/square_quad.mesh') mesh3d_ref = Mesh.from_file(data_dir +\ '/meshes/3d/cube_medium_hexa.mesh') mesh3d_gen = expand2d(mesh2d, 0.1, 10) mesh3d_gen.coors[:,2] += -0.5 d0 = nm.sort(nm.linalg.norm(mesh3d_ref.coors, axis=1)) d1 = nm.sort(nm.linalg.norm(mesh3d_gen.coors, axis=1)) ok = nm.linalg.norm(d0 - d1) < 1e-6 assert ok