from paraBEM.mesh import mesh_object from paraBEM.pan3d import DirichletDoublet0Source0Case3 as Case from paraBEM.vtk_export import CaseToVTK mesh = mesh_object.from_OBJ("../mesh/sphere_half.obj") for panel in mesh.panels: panel.set_symmetric() case = Case(mesh.panels) case.run() writer = CaseToVTK(case, "results/symmetric_test") writer.write_panels(data_type="point") writer.write_field([-2, 2, 20], [-2, 2, 20], [-2, 2, 20])