from paraBEM.pan2d import doublet_2_1
import paraBEM
import numpy as np
v1 = paraBEM.PanelVector2(-1, 0)
v2 = paraBEM.PanelVector2(1, 0)
panel = paraBEM.Panel2([v2, v1])
vals = ([
doublet_2_1(paraBEM.Vector2(x, 0), panel, True)
for x in np.linspace(-2, 2, 20)
])
print(vals)
import paraBEM
from paraBEM.pan2d import doublet_2_0, source_2_0, doublet_2_0_v, source_2_0_v, doublet_2_1, doublet_2_1_v
from paraBEM.vtk_export import VtkWriter
import numpy as np
from paraBEM.utils import check_path
v1 = paraBEM.PanelVector2(-2, 0)
v2 = paraBEM.PanelVector2(2, 0)
panel = paraBEM.Panel2([v1, v2])
n = 500
space = np.linspace(-5, 5, n)
grid = [paraBEM.Vector2([x, y]) for y in space for x in space]
dub_vals = [doublet_2_0(target, panel) for target in grid]
dub_vals_lin_1 = [doublet_2_1(target, panel, True) for target in grid]
dub_vals_lin_2 = [doublet_2_1(target, panel, False) for target in grid]
src_vals = [source_2_0(target, panel) for target in grid]
dublinv1_vals = [doublet_2_1_v(target, panel, True) for target in grid]
dublinv2_vals = [doublet_2_1_v(target, panel, False) for target in grid]
dubv_vals = [doublet_2_0_v(target, panel) for target in grid]
srcv_vals = [source_2_0_v(target, panel) for target in grid]
writer = VtkWriter()
with open(check_path("results/element_2.vtk"), "w") as _file:
writer.structed_grid(_file, "element_2", [n, n, 1])
writer.points(_file, grid)
writer.data(_file, dub_vals, name="doublet", _type="SCALARS", data_type="POINT_DATA")
writer.data(_file, dub_vals_lin_1, name="doublet_lin_1", _type="SCALARS", data_type="POINT_DATA")
writer.data(_file, dub_vals_lin_2, name="doublet_lin_2", _type="SCALARS", data_type="POINT_DATA")
from paraBEM.pan2d import doublet_2_1 import paraBEM import numpy as np v1 = paraBEM.PanelVector2(-1, 0) v2 = paraBEM.PanelVector2(1, 0) panel = paraBEM.Panel2([v2, v1]) vals = ([doublet_2_1(paraBEM.Vector2(x, 0), panel, True) for x in np.linspace(-2, 2, 20)]) print(vals)
