#include <Eigen/src/SparseCore/BlockSparseMatrix.h>
Proxy returned by adjoint(): solveInPlace solves T^H x = b.