Matrix-free methods
In
- the power method,
- the Lanczos algorithm,[2]
- Locally Optimal Block Preconditioned Conjugate Gradient Method (LOBPCG),[3]
- Wiedemann's coordinate recurrence algorithm,[4] and
- the conjugate gradient method.[5]
- Krylov subspace methods
Distributed solutions have also been explored using coarse-grain parallel software systems to achieve homogeneous solutions of linear systems.[6]
It is generally used in solving non-linear equations like Euler's equations in computational fluid dynamics. Matrix-free conjugate gradient method has been applied in the non-linear elasto-plastic finite element solver.[7] Solving these equations requires the calculation of the Jacobian which is costly in terms of CPU time and storage. To avoid this expense, matrix-free methods are employed. In order to remove the need to calculate the Jacobian, the Jacobian vector product is formed instead, which is in fact a vector itself. Manipulating and calculating this vector is easier than working with a large matrix or linear system.
References
- ISBN 978-0-691-12202-1
- .
- ISBN 978-3-540-54508-8
- S2CID 13305650
- .