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