Input

• A - S3L internal array handle for the global general sparse matrix. It is produced by a prior call to one of the following sparse routines:

  • S3L_declare_sparse

  • S3L_read_sparse

  • S3L_rand_sparse

• b - Global array of rank 1, with the same data type and precision as A and x and a length equal to the order of the sparse matrix. b contains the right-hand side vector of the linear problem.

• x - Global array of rank 1, with the same data type and precision as A and b and a length equal to the order of the sparse matrix. On input, x contains the initial guess for the solution to the linear system. Upon successful completion, x contains the converged solution (see the Output section).

• iparm - Integer local array of rank 1 and length s3l_iter_iparm_size, where:

  • iparm[S3l_iter_solver] - Specifies the iterative algorithm to be used. Set it to 0 to use the default solver GMRES. See the Desctription sectino for details. 

  • iparm[S3l_iter_pc] - Specifies the preconditioner to be used. Set it to 0 to use the default option, S3L_none.

  • iparm[S3l_iter_conv] - Selects the criterion to be used for stopping the computation.

• rparm - Specifies options for computing all or part of the matrix U.