Package bemsample


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)


Variable Documentation

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.I = eye(A.shape[0])

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)


Generated on Sat Dec 20 18:43:32 2008 for OpenAcoustics by  doxygen 1.5.5