1. change FEMBase
2. change SparseMatrix transpose() to const