The ITSOLVER call has the following syntax and arguments:
call ITSOLVER (x, error, iter, method, A, b, precond, tol, maxiter, start, history);
The conjugate gradient and minimum residual algorithms (method = ’CG’ or method = ’MINRES’) require to be symmetric; hence you must specify only the lower triangular part of , while the remaining algorithms require all nonzero coefficients to be listed. The following table lists valid values for the precond parameter for each class of algorithm.
Table 21.1: Subroutine Definitions and Valid Preconditioners
Method Value |
Algorithm |
Preconditioners |
---|---|---|
|
conjugate gradient |
|
|
minimum residual |
|
|
biconjugate gradient |
|
|
conjugate gradient squared |
|
solution vector
final solution error (optional)
resultant number of iterations (optional)
three-column matrix of triplets, where the first column contains the value, the next column contains the row indices, and the third column contains the column indices of the nonzero matrix coefficients. The order in which triplets are listed is insignificant. For symmetric matrices specify only the lower triangular part, including the main diagonal (row indices must be greater than or equal to the corresponding column indices). Zero coefficients should not be included. No missing values or duplicate entries are allowed.
the right-hand-side vector
preconditioner, default value "NONE"
desired tolerance, default value
maximum number of iterations, default value
initial guess
the history of errors for each iteration