C:/Documents and Settings/hisaha/My Documents/MyProjects/SourceForge/Foac/trunk/openacoustics/examples/bemsample.py File Reference


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)

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