def test_save_scalar_field_hexagonal_mesh(tmpdir):
mesh = HexagonalMesh(1, 3, 3)
s = scalar_field(mesh, lambda r: r[0] + r[1])
vtk = VTK(mesh, directory=str(tmpdir), filename="scalar_hexagonal")
vtk.save_scalar(s, name="s")
assert same_as_ref(vtk.write_file(), REF_DIR)
def test_initialise_scalar():
mesh = CuboidMesh(1, 1, 1, 1, 1, 1)
f = scalar_field(mesh, lambda r: r[0] + r[1] + r[2])
assert very_close(f, np.array((1.5)))
def test_save_scalar_field(tmpdir):
mesh = CuboidMesh(4, 3, 2, 4, 3, 2)
s = scalar_field(mesh, lambda r: r[0] + r[1] + r[2])
vtk = VTK(mesh, directory=str(tmpdir), filename="save_scalar")
vtk.save_scalar(s, name="s")
assert same_as_ref(vtk.write_file(), REF_DIR)
def test_save_scalar_field_hexagonal_mesh(tmpdir):
mesh = HexagonalMesh(1, 3, 3)
s = scalar_field(mesh, lambda r: r[0] + r[1])
vtk = VTK(mesh, directory=str(tmpdir), filename="scalar_hexagonal")
vtk.save_scalar(s, name="s")
assert same_as_ref(vtk.write_file(), REF_DIR)
def test_save_scalar_field(tmpdir):
mesh = CuboidMesh(4, 3, 2, 4, 3, 2)
s = scalar_field(mesh, lambda r: r[0] + r[1] + r[2])
vtk = VTK(mesh, directory=str(tmpdir), filename="save_scalar")
vtk.save_scalar(s, name="s")
assert same_as_ref(vtk.write_file(), REF_DIR)
def test_save_scalar_field_hexagonal_mesh():
mesh = HexagonalMesh(1, 3, 3)
s = scalar_field(mesh, lambda r: r[0] + r[1])
vtk = VTK(mesh, directory=OUTPUT_DIR, filename="scalar_hexagonal")
vtk.save_scalar(s, name="s")
def test_save_scalar_field():
mesh = CuboidMesh(4, 3, 2, 4, 3, 2)
s = scalar_field(mesh, lambda r: r[0] + r[1] + r[2])
vtk = VTK(mesh, directory=OUTPUT_DIR, filename="save_scalar")
vtk.save_scalar(s, name="s")
