Module SOMN in SLATECPreconditioned Orthomin Sparse Iterative Ax=b Solver. Routine to solve a general linear system Ax = b using the Preconditioned Orthomin method. (SLAP Package).
CALL SOMN (N, B, X, NELT, IA, JA, A, ISYM, MATVEC, MSOLVE, NSAVE, ITOL, TOL, ITMAX, ITER, ERR, IERR, IUNIT, R, Z, P, AP, EMAP, DZ, CSAV, RWORK, IWORK)| D2a4 | Solution of systems of linear equations, including matrix inverses, and LU and related factorizations (real nonsymmetric sparse matrices) |