nag_sparse_nherm_matvec (f11xnc) computes either the matrix-vector product
, or the conjugate transposed matrix-vector product
, according to the value of the argument
trans, where
is a complex
by
sparse non-Hermitian matrix, of arbitrary sparsity pattern. The matrix
is stored in coordinate storage (CS) format (see
Section 2.1.1 in the f11 Chapter Introduction). The array
a stores all the nonzero elements of
, while arrays
irow and
icol store the corresponding row and column indices respectively.
It is envisaged that a common use of
nag_sparse_nherm_matvec (f11xnc) will be to compute the matrix-vector product required in the application of
nag_sparse_nherm_basic_solver (f11bsc) to sparse complex linear systems. This is illustrated in
Section 10 in
nag_sparse_nherm_precon_ssor_solve (f11drc).
None.
The computed vector
satisfies the error bound:
- , if , or
- ,
if ,
where
is a modest linear function of
, and
is the
machine precision.
Please consult the
x06 Chapter Introduction for information on how to control and interrogate the OpenMP environment used within this function. Please also consult the
Users' Note for your implementation for any additional implementation-specific information.
The time taken for a call to
nag_sparse_nherm_matvec (f11xnc) is proportional to
nnz.
It is expected that a common use of
nag_sparse_nherm_matvec (f11xnc) will be to compute the matrix-vector product required in the application of
nag_sparse_nherm_basic_solver (f11bsc) to sparse complex linear systems. In this situation
nag_sparse_nherm_matvec (f11xnc) is likely to be called many times with the same matrix
. In the interests of both reliability and efficiency you are recommended to set
for the first of such calls, and to set
for all subsequent calls.