Contains status value of the last Poisson Solver routine called. In general, it should be 0 on exit from a routine to proceed with the Fast Helmholtz Solver. The element has no predefined values. This element can also be used to inform the

?_commit_Helmholtz_2D/?_commit_Helmholtz_3D

/

?_commit_sph_p/?_commit_sph_np

routines of how the Commit step of the computation should be carried out (see Figure

"Typical Order of Invoking Poisson Solver Routines"

). A non-zero value of

ipar[0]

with decimal representation =100

a

+10

b

+

c

, where each of

a

,

b

, and

c

is equal to 0 or 9, indicates that some parts of the Commit step should be omitted.

If
c
=9, the routine omits checking of parameters and initialization of the data structures.
If
b
=9,
In the Cartesian case, the routine omits the adjustment of the right-hand side vector
f
to the Neumann boundary condition (multiplication of boundary values by 0.5 as well as incorporation of the boundary function
g
) and/or the Dirichlet boundary condition (setting boundary values to 0 as well as incorporation of the boundary function
G
).
For the Helmholtz solver on a sphere, the routine omits computation of the spherical weights for the
dpar
/
spar
array.
If
a
=9, the routine omits the normalization of the right-hand side vector
f
. Depending on the solver, the normalization means:
2D Cartesian case: multiplication by
h
y
2
, where
h
y
is the mesh size in the
y
direction (for details, see Poisson Solver Implementation).
3D (Cartesian) case: multiplication by
h
z
2
, where
h
z
is the mesh size in the
z
direction.
Helmholtz solver on a sphere: multiplication by
h
θ
2
, where
h
θ
is the mesh size in the θ direction (for details, see Poisson Solver Implementation).

Using

ipar[0]

you can adjust the routine to your needs and improve efficiency in solving multiple Helmholtz problems that differ only in the right-hand side. You must be cautious when using this method, because any misunderstanding of the commit process may cause incorrect results or program failure (see also Caveat on Parameter Modifications).

Contains error messaging options:

In case of errors, the

stat

parameter contains a non-zero value on exit from a routine regardless of the

ipar[1]

setting.

Contains warning messaging options:

In case of warnings, the

stat

parameter contains a non-zero value on exit from a routine regardless of the

ipar[2]

setting.

Internal parameters.

Parameters

6 through 11

are used only in the Cartesian case.

Takes this value:

Takes this value:

Takes this value:

Takes this value:

Takes this value:

Takes this value:

Takes the value of

Takes the value of

Takes the value of

nz

, the number of intervals along the

z

-axis. This parameter is used only in the 3D case (Cartesian).

Internal parameters which define the internal partitioning of the

dpar

/

spar

array.

The values of

ipar[21]

-

ipar[119]

are assigned regardless of the dimension of the problem for the Cartesian solver or of whether the solver on a sphere is periodic.

Contains message style options

. Specifically

:

Contains the number of OpenMP threads to be used for computations in a multithreaded environment. The default value is 1 in the serial mode, and the result returned by the

mkl_get_max_threads

function otherwise.

Internal parameters which define the internal partitioning of the

dpar

/

spar

array.

Takes the value of

ipar[18]

+1, which specifies the internal partitioning of the

dpar/spar

array in the periodic Cartesian case.

Takes the value of

ipar[23]

+3*

ipar[12]

/4, which specifies the internal partitioning of the

dpar/spar

array in the periodic Cartesian case.

Takes the value of

ipar[20]

+1, which specifies the internal partitioning of the

dpar/spar

array in the periodic 3D Cartesian case.

Takes the value of

ipar[25]

+3*

ipar[13]

/4, which specifies the internal partitioning of the

dpar/spar

array in the periodic 3D Cartesian case.

Unused.

Contain the first twenty elements of the

ipar

array of the first Trigonometric Transform that the solver uses. (For details, see Common Parameters in the "Trigonometric Transform Routines" section.)

Contain the first twenty elements of the

ipar

array of the second Trigonometric Transform that the 3D Cartesian and periodic spherical solvers use. (For details, see Common Parameters in the "Trigonometric Transform Routines" section.)

Contain the first twenty elements of the

ipar

array of the third Trigonometric Transform that the solver uses in case of periodic boundary conditions along the

x

-axis. (For details, see Common Parameters in the "Trigonometric Transform Routines" section.)

Contain the first twenty elements of the

ipar

array of the fourth Trigonometric Transform used by periodic spherical solvers and 3D Cartesian solvers with periodic boundary conditions along the

y

-axis. (For details, see Common Parameters in the "Trigonometric Transform Routines" section.)

Internal parameters used by nonuniform 3D solvers.