Using an augmented Lagrangian matrix approach, we analytically solve in this paper a broad class of linear systems that includes symmetric and nonsymmetric problems in saddle point form. To this end, some mild assumptions are made and a preconditioning is specially designed to improve the sensitivity of the systems before the calculation of their solutions. In this way, they can be properly employed for computational purposes. As an example of such a procedure, we develop a direct method for solving the symmetric saddle point problem. Then we show that our method is intrinsically related with the null space method, but it doesn’t depend on finding a basis of the null space of the $(2,1)$ block of the saddle point systems. For this reason, the proposed method can be seen as a basis-free null space method. This method is also extended to nonsymmetric problems. To the best of our knowledge, the analytical solution of the nonsymmetric linear systems considered here is not known. So it is of major interest, since it can lead to the development of efficient algorithms, besides our own.
Citation
Institute of Mathematics, Statistics and Computing Science, University of Campinas, Campinas, 13083-859, SP, Brazil; Department of Mathematics, Federal University of Paraná, Curitiba, 81531-980, PR, Brazil; 08/Set/2015