Application of modified homotopy perturbation method for solving the augmented systems (Q2389561)

A new approach is presented for solving the augmented linear system \[ \left(\begin{matrix} A & B \\ B^T & 0 \end{matrix}\right) \left(\begin{matrix} x \\ y \end{matrix}\right) = \left(\begin{matrix} b \\ q \end{matrix}\right), \tag{1} \] where \(A \in \mathbb R^{m \times m}\) is a symmetric and positive definite matrix, \(B \in \mathbb R^{m \times n} (m \geq n)\) is a matrix of full column rank, \(B^T\) is the transposed of \(B\), and \(b \in \mathbb R^m, q \in \mathbb R^n\) are given vectors. Since the second diagonal block matrix is null, a full-rank matrix \(Q\) with small parameter \(p \in [0, 1]\) is introduced and a new iterative method is constructed by using the modified homotopy perturbation method. Sufficient and necessary conditions to guarantee the iteration to converge are derived. Four special choices of the full-rank matrix \(Q\) are considered for solving the problem (1). Numerical experiments show that the method proposed is more simple and effective than the existing methods.