Updated 5 months ago
https://github.com/cvanaret/symmetricblockmatrix
Form symmetric block matrices and condense them (eliminating blocks) at compile time using Schur complements. The goal is to manipulate KKT/augmented/saddle-point systems symbolically to obtain equivalent systems (typically: unsymmetric, symmetrized, and normal equations) that can be tackled by various numerical solvers