Description: Implementation of a conjugate-gradient type method for solving sparse linear equations:
Solve Ax = b or (A- sI)x = b.
The matrix A- sI must be symmetric and nonsingular, but it may be definite or indefinite. The scalar s is a shifting parameter-- it may be any number. The method is based on Lanczos tridiagonalization. You may provide a preconditioner, but it must be positive definite.
- [lesolv] - This package contains a lot of equations
- [lanczos1] - Lanczos infusion using eigenvalue
- [sor] - Relaxation Factor in the SOR iterative a
File list (Check if you may need any files):
symmlq.f
symmlqblas.f
symmlqtest.f
symmlq_f77.README.txt