Matrix-Free Conjugate Gradient Implementation of Implicit Schemes

Abstract

Matrix-free conjugate gradient algorithms are described as applied to large systems of algebraic equations arising in the implementation of implicit schemes intended for problems in continuum mechanics. The algorithms are shown to reduce the computer memory requirements and provide a higher computer performance. Additional advantages include the robustness of the algorithms and the simplicity of their implementation and debugging.