Skip to content
Snippets Groups Projects
Commit 12827a9c authored by Jonathan Youett's avatar Jonathan Youett
Browse files

Minor cleanup

parent c23f576f
No related branches found
No related tags found
No related merge requests found
...@@ -119,18 +119,18 @@ private: ...@@ -119,18 +119,18 @@ private:
constexpr int jDim = Gradient::dimension; constexpr int jDim = Gradient::dimension;
FieldMatrix<field_type,dim,jDim> mass(0); FieldMatrix<field_type,dim,jDim> mass(0);
for (size_t j=0; j<corners.size(); j++) for (size_t j=0; j < corners.size(); ++j)
for (size_t i=0; i<corners.size(); i++) for (size_t i=0; i < corners.size(); ++i)
mass[j].axpy(corners[i][j], jacobians[i]); mass[j].axpy(corners[i][j], jacobians[i]);
FieldMatrix<field_type,jDim,dim> trans; FieldMatrix<field_type, jDim, dim> trans;
Dune::MatrixVector::transpose(mass,trans); Dune::MatrixVector::transpose(mass, trans);
auto inverse = mass.leftmultiplyany(trans); auto inverse = mass.leftmultiplyany(trans);
try { try {
inverse.invert(); inverse.invert();
} catch(Dune::FMatrixError e) { } catch(Dune::FMatrixError& e) {
DUNE_THROW(Dune::Exception, "Computing dual covariant basis failed, metric not invertible: " << inverse.determinant() << " \n"); DUNE_THROW(Dune::FMatrixError, "Computing dual covariant basis failed, metric not invertible: " << inverse.determinant() << " \n");
} }
auto dual = inverse.rightmultiplyany(trans); auto dual = inverse.rightmultiplyany(trans);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment