f11xac computes either the matrix-vector product
, or the transposed matrix-vector product
, according to the value of the argument
trans, where
is an
by
sparse nonsymmetric 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 nonzero elements of
, while arrays
irow and
icol store the corresponding row and column indices respectively.
It is envisaged that a common use of
f11xac will be to compute the matrix-vector product required in the application of
f11bec to sparse linear systems. An illustration of this usage appears in
Section 10 in
f11ddc.
None.
- A nonzero element has been supplied which does not lie within the matrix , is out of order, or has duplicate row and column indices. Consider calling f11zac to reorder and sum or remove duplicates.
- NE_ALLOC_FAIL
-
Dynamic memory allocation failed.
See
Section 3.1.2 in the Introduction to the NAG Library CL Interface for further information.
- NE_BAD_PARAM
-
On entry, argument had an illegal value.
- NE_INT
-
On entry, .
Constraint: .
On entry, .
Constraint: .
- NE_INT_2
-
On entry, and .
Constraint: .
- NE_INTERNAL_ERROR
-
An internal error has occurred in this function. Check the function call and any array sizes. If the call is correct then please contact
NAG for assistance.
See
Section 7.5 in the Introduction to the NAG Library CL Interface for further information.
- NE_INVALID_CS
-
On entry, , and .
Constraint: and .
On entry, , and .
Constraint: and .
- NE_NO_LICENCE
-
Your licence key may have expired or may not have been installed correctly.
See
Section 8 in the Introduction to the NAG Library CL Interface for further information.
- NE_NOT_STRICTLY_INCREASING
-
On entry, is out of order: .
On entry, the location () is a duplicate: .
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
f11xac is proportional to
nnz.
It is expected that a common use of
f11xac will be to compute the matrix-vector product required in the application of
f11bec to sparse linear systems. In this situation
f11xac 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.