double*

for

d_init_Helmholtz_2D/d_init_Helmholtz_3D

,

float*

for

s_init_Helmholtz_2D/s_init_Helmholtz_3D

.

The coordinate of the leftmost boundary of the domain along the x-axis.

double*

for

d_init_Helmholtz_2D/d_init_Helmholtz_3D

,

float*

for

s_init_Helmholtz_2D/s_init_Helmholtz_3D

.

The coordinate of the rightmost boundary of the domain along the x-axis.

double*

for

d_init_Helmholtz_2D/d_init_Helmholtz_3D

,

float*

for

s_init_Helmholtz_2D/s_init_Helmholtz_3D

.

The coordinate of the leftmost boundary of the domain along the y-axis.

double*

for

d_init_Helmholtz_2D/d_init_Helmholtz_3D

,

float*

for

s_init_Helmholtz_2D/s_init_Helmholtz_3D

.

The coordinate of the rightmost boundary of the domain along the y-axis.

double*

for

d_init_Helmholtz_3D

,

float*

for

s_init_Helmholtz_3D

.

The coordinate of the leftmost boundary of the domain along the z-axis. This parameter is needed only for the

?_init_Helmholtz_3D

routine.

double*

for

d_init_Helmholtz_3D

,

float*

for

s_init_Helmholtz_3D

.

The coordinate of the rightmost boundary of the domain along the z-axis. This parameter is needed only for the

?_init_Helmholtz_3D

routine.

MKL_INT*

. The number of mesh intervals along the x-axis.

MKL_INT*

. The number of mesh intervals along the y-axis.

MKL_INT*

. The number of mesh intervals along the z-axis. This parameter is needed only for the

?_init_Helmholtz_3D

routine.

char*

. Contains the type of boundary conditions on each boundary. Must contain four characters for

?_init_Helmholtz_2D

and six characters for

?_init_Helmholtz_3D

. Each of the characters can be 'N' (Neumann boundary condition), 'D' (Dirichlet boundary condition), or 'P' (periodic boundary conditions). Specify the types of boundary conditions for the boundaries in the following order:

bd_a

x

,

bd_b

x

,

bd_a

y

,

bd_b

y

,

bd_a

z

, and

bd_b

z

. Specify periodic boundary conditions on the respective boundaries in pairs (for example, 'PPDD' or 'NNPP' in the 2D case). The types of boundary conditions for the last two boundaries are needed only in the 3D case.

double*

for

d_init_Helmholtz_2D/d_init_Helmholtz_3D

,

float*

for

s_init_Helmholtz_2D/s_init_Helmholtz_3D

.

The constant Helmholtz coefficient. Note that to solve Poisson or Laplace problem, you should set the value of

q

to 0.

MKL_INT

array of size 128. Contains integer data to be used by Fast Helmholtz Solver (for details, refer to

ipar

).

double

array of size 5*

nx

/2+7 in the 2D case or 5*(

nx

+

ny

)/2+9 in the 3D case. Contains double-precision data to be used by Fast Helmholtz Solver (for details, refer to

dpar

and

spar

).

float

array of size 5*

nx

/2+7 in the 2D case or 5*(

nx

+

ny

)/2+9 in the 3D case. Contains single-precision data to be used by Fast Helmholtz Solver (for details, refer to

dpar

and

spar

).

MKL_INT*

. Routine completion status, which is also written to

ipar[0]

. Continue to call other Poisson Solver routines only if the status is 0.