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])
import paraBEM
from paraBEM import pan3d
from paraBEM.vtk_export import CaseToVTK
from paraBEM.mesh import mesh_object
import numpy as np
mesh = mesh_object.from_OBJ("../mesh/sphere_low_tri.obj")
case = pan3d.DirichletDoublet0Case3(mesh.panels)
case.v_inf = paraBEM.Vector3(1, 0, 0.)
case.farfield = 100
case.run()
lin = np.linspace(-0.5, 0.5, 5)
grid = [[-2, k, j] for j in lin for k in lin]
vtk_writer = CaseToVTK(case, "results/vtk_test_case")
vtk_writer.write_panels(data_type="point")
# vtk_writer.write_wake_panels()
vtk_writer.write_body_stream(mesh.panels, 50)
vtk_writer.write_field([-2, 2, 100], [-2, 2, 100], [-1, 1, 3])
