Subroutine ledhmz2b

subroutine ledhmz2b

        ! Variables
    integer, intent(in) :: flag1
    integer, intent(in) :: flag2
    integer, intent(in) :: fdim1
    integer, intent(in) :: fdim2
    integer, intent(in) :: initb
    integer, intent(in) :: init
    integer, intent(in) :: n1
    integer, intent(in) :: m1
    integer, intent(in) :: nu
    integer, intent(in) :: nun
    integer, intent(in) :: nf
    double precision, intent(in), dimension (fdim1,fdim1) :: e1
    double precision, intent(in), dimension (fdim1-flag1,fdim1-flag1) :: e2
    double precision, intent(in), dimension (fdim1) :: w1
    double precision, intent(in), dimension (fdim1-flag1) :: w2
    double precision, intent(in), dimension (0:nun,9) :: un
    double precision, intent(in), dimension (0:n1,0:n1) :: bn
    double precision, intent(in), dimension (0:m1,0:m1) :: bm
    double precision, intent(in) :: rm
    double precision, intent(in) :: beta
    double precision, intent(inout), dimension (0:nu,0:m1) :: u
    double precision, intent(inout), dimension (0:nf,0:m1) :: f
    integer :: i

end subroutine ledhmz2b

subroutine ledhmz2b(n1,m1,rm,beta,un,nun,bn,bm,f,nf,u,nu,e1,e2,w1,w2,flag1,flag2,fdim1,fdim2,initb,init)

Solve the 2D Helmholtz equation: rm*u - u_xx - beta u_yy= f with nonhomogeneous Dirichlet B.C. in IxI.

Input:

f: f(0:nf, 0:m1) coefficients in Legendre of right hand side f
un: contains the physical boundary data:\, nun>=max(n1,m1)
nf: nf>=n1
bn,bm: output from leinit.f
e,w: Output from ledeig.f, contain eigenfunctions and eigenvalues of -u_{xx}= \lambda u.
ws: working array. There are initialized with init=1, they can be repeatedly used with init.ne.1.
rm,beta: constants in the Helmholtz equation
init: initialization flag: init=1, initialize ws.
initb: initialization flag: init=1, initialize un.

Output:

u, coefficients in Legendre of the solution. (nu >= n1)

Author: Jan Ivar Moldekleiv, Jie Shen

Version: 0.8

Description of Variables

flag1

integer, intent(in) :: flag1

flag2

integer, intent(in) :: flag2

fdim1

integer, intent(in) :: fdim1

fdim2

integer, intent(in) :: fdim2

initb

integer, intent(in) :: initb

init

integer, intent(in) :: init

n1

integer, intent(in) :: n1

m1

integer, intent(in) :: m1

nu

integer, intent(in) :: nu

nun

integer, intent(in) :: nun

nf

integer, intent(in) :: nf

e1

double precision, intent(in), dimension (fdim1,fdim1) :: e1

e2

double precision, intent(in), dimension (fdim1-flag1,fdim1-flag1) :: e2

w1

double precision, intent(in), dimension (fdim1) :: w1

w2

double precision, intent(in), dimension (fdim1-flag1) :: w2

un

double precision, intent(in), dimension (0:nun,9) :: un

bn

double precision, intent(in), dimension (0:n1,0:n1) :: bn

bm

double precision, intent(in), dimension (0:m1,0:m1) :: bm

rm

double precision, intent(in) :: rm

beta

double precision, intent(in) :: beta

u

double precision, intent(inout), dimension (0:nu,0:m1) :: u

f

double precision, intent(inout), dimension (0:nf,0:m1) :: f

i

integer :: i