pcg_orthogonal_proj​ector.m

Projected Conjugate Gradients Method by orthogonal projector.

Iterative solver is suitable for solving system of linear equations arising in domain decomposition method, especially in the case, FETI is used.
Routine attempts to solve the system of linear equations A*x=b subject to condition B*x=c. The N-by-N coefficient matrix A must be SPD and the right hand side column vector b must have length N. The M-by-N matrix B contains linearly independent rows, column vector c must have length M.