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