lsolve(1) - Linux man page
Name
lsolve - linear solver for real sparse matrices
Synopsis
lsolve matrix_filename rhs_setting solution_filename residual_filename [options]
Description
This program inputs the data of the coefficient matrix from matrix_filename and solves the linear equation Ax = b with the solver specified by options. It outputs the solution to solution_filename in the extended Matrix Market format and the residual_history to residual_filename in the PLAIN format (see Appendix of the Lis User Manual). The Matrix Market and extended Matrix Market formats are supported for matrix_filename.
One of the following values can be specified by rhs_setting:
- 0
- Use the right hand side vector b included in matrix_filename
- 1
- Use b = (1, ..., 1)^T
- 2
- Use b = A x (1, ..., 1)^T
- rhs_filename
- The filename for the right-hand side vector
The PLAIN and Matrix Market formats are supported for rhs_filename.
Options
The following options are supported:
- -i linear solver
- The following options are supported for linear solver:
- -i {cg|1}
- CG
- -i {bicg|2}
- BiCG
- -i {cgs|3}
- CGS
- -i {bicgstab|4}
- BiCGSTAB
- -i {bicgstabl|5}
- BiCGSTAB(l)
- -ell [2]
- -i {gpbicg|6}
- GPBiCG
- -i {tfqmr|7}
- TFQMR
- -i {orthomin|8}
- Orthomin(m)
- -restart [40]
- -i {gmres|9}
- GMRES(m)
- -restart [40]
- -i {jacobi|10}
- Jacobi
- -i {gs|11}
- Gauss-Seidel
- -i {sor|12}
- SOR
- -omega [1.9]
- The relaxation coefficient omega (0<omega<2)
- -i {bicgsafe|13}
- BiCGSafe
- -i {cr|14}
- CR
- -i {bicr|15}
- BiCR
- -i {crs|16}
- CRS
- -i {bicrstab|17}
- BiCRSTAB
- -i {gpbicr|18}
- GPBiCR
- -i {bicrsafe|19}
- BiCRSafe
- -i {fgmres|20}
- FGMRES(m)
- -restart [40]
- -i {idrs|21}
- IDR(s)
- -irestart [2]
- -i {minres|22}
- -p preconditioner
- The following options are supported for preconditioner:
- -p {none|0}
- None
- -p {jacobi|1}
- Jacobi
- -p {ilu|2}
- ILU(k)
- -ilu_fill [0]
- -p {ssor|3}
- SSOR
- -ssor_w [1.0]
- The relaxation coefficient omega (0<omega<2)
- -p {hybrid|4}
- Hybrid
- -hybrid_i [sor]
The maximum number of the iterations
- -hybrid_maxiter [25]
The convergence criterion
- -hybrid_tol [1.0e-3]
- -hybrid_w [1.5]
- The relaxation coefficient omega of the SOR (0<omega<2)
- -hybrid_ell [2]
- The degree l of the BiCGSTAB(l)
The restart values of the GMRES and Orthomin
- -hybrid_restart [40]
- -p {is|5}
- I+S
- -is_alpha [1.0]
- The parameter alpha of I+alpha*S(m)
- -is_m [3]
- The parameter m of I+alpha*S(m)
- -p {sainv|6}
- SAINV
- -sainv_drop [0.05]
- -p {saamg|7}
- SA-AMG
- -saamg_unsym [false]
- Select the unsymmetric version (The matrix structure must be symmetric)
The drop criterion
- -saamg_theta [0.05|0.12]
- -p {iluc|8}
- Crout ILU
- -iluc_drop [0.05]
The ration of maximum fill-in
- -iluc_rate [5.0]
- -p {ilut|9}
- ILUT
- -ilut_drop [0.05]
The ration of maximum fill-in
- -ilut_rate [5.0]
- -adds true
- Additive Schwarz
- -adds_iter [1]
Other Options:
- -maxiter [1000]
- The maximum number of the iterations
- -tol [1.0e-12]
- The convergence criterion
- -print [0]
- The display of the residual
- -print {none|0}
Save the residual history
- -print {mem|1}
Display the residual history
- -print {out|2}
Save the residual history and display it on the screen
- -print {all|3}
- -scale [0]
- The scaling
- -scale {none|0}
The Jacobi scaling
- -scale {jacobi|1}
The diagonal scaling
- -scale {symm_diag|2}
- -initx_zeros [true]
- The behavior of the initial vector x_0
- -initx_zero {false|0}
All values are set to 0
- -initx_zero {true|1}
- -omp_num_threads [t]
- The number of the threads (t represents the maximum number of the threads)
- -storage [0]
- The matrix storage format
- -storage_block [2]
- The block size of the BSR and BSC formats
- -f [0]
- The precision of the linear solver
- -f {double|0}
Quadruple precision
- -f {quad|1}
See Lis User Manual for full description.
Exit Status
The following exit values are returned:
- 0
- The process is normally terminated
- unspecified
- An
error occurred
See Also
lis(3), esolve(1), spmvtest1(1), spmvtest2(1), spmvtest3(1), spmvtest4(1), spmvtest5(1)
http://www.ssisc.org/lis/
http://math.nist.gov/MatrixMarket/