spqr solver with block matrix
Summary
Attempt to repair the SPQR solver. Although it is specialized for BCRSMatrix<FieldMatrix<T,n,m>> block-sizes n,m > 1 cannot work. It seems that also the T=std::complex overload does not work. This is not fixed in the MR.