|
Namespaces |
namespace | bemsample |
Variables |
string | bemsample.geo = 'square.geo' |
tuple | bemsample.gmsh = GMsh2D() |
tuple | bemsample.nodes = gmsh.getNodes() |
tuple | bemsample.sf = gmsh.getLines("") |
list | bemsample.nodenum = nodes.shape[0] |
list | bemsample.sfnum = sf.shape[0] |
| bemsample.M = sfnum |
tuple | bemsample.S = zeros(M, 'float64') |
tuple | bemsample.r = zeros((M, 2), 'float64') |
tuple | bemsample.nv = zeros((M, 2), 'float64') |
list | bemsample.node1 = nodes[sf[ii, 0], :] |
list | bemsample.node2 = nodes[sf[ii, 1], :] |
| bemsample.v1 = node2-node1 |
tuple | bemsample.normal = array([v1[1], -v1[0]]) |
list | bemsample.mp = r[ii, :] |
int | bemsample.beta = 0 |
int | bemsample.w = 2 |
int | bemsample.c = 343 |
float | bemsample.ro0 = 1.21 |
| bemsample.k = w/c |
tuple | bemsample.A = zeros((M, M), 'complex64') |
list | bemsample.en = nv[n, :] |
list | bemsample.rs = [0., 0.] |
| node1 = nodes[sf[ii, 0], :] node2 = nodes[sf[ii, 1], :] f = lambda y, x:G(k, r[m, :], array([x, y])) g = lambda y, x:Gn(k, r[m, :], array([x, y]), en) Gi, err = dblquad(f, node1[1], node2[1], lambda x:node1[0], lambda x:node2[0]) Gni, err = dblquad(g, node1[1], node2[1], lambda x:node1[0], lambda x:node2[0]) A[m, n] = -1j * w * ro0 * beta * Gi A[m, n] -= Gni
|
tuple | bemsample.Pd = zeros(M, "complex64") |
tuple | bemsample.I = eye(A.shape[0]) |
tuple | bemsample.Phi = linalg.solve(0.5 * I - A, Pd) |
float | bemsample.dx = 0.1 |
float | bemsample.dy = 0.1 |
tuple | bemsample.xx = arange(dx, 1, dx) |
tuple | bemsample.yy = arange(dy, 1, dy) |
tuple | bemsample.An = zeros(M,'complex64') |
tuple | bemsample.Pr = zeros((len(xx), len(yy)), 'complex64') |
tuple | bemsample.rr = array([xx[ii], yy[jj]]) |
tuple | bemsample.pd = G(k, rs, rr) |