diff --git a/dune/contact/assemblers/largedeformationcoupling.hh b/dune/contact/assemblers/largedeformationcoupling.hh index c62e4d3d15a0ffefb4c92b42dd1f66639344daed..dc56e1f5d2e39672fad454ed4dac288791095a22 100644 --- a/dune/contact/assemblers/largedeformationcoupling.hh +++ b/dune/contact/assemblers/largedeformationcoupling.hh @@ -119,18 +119,18 @@ private: constexpr int jDim = Gradient::dimension; FieldMatrix<field_type,dim,jDim> mass(0); - for (size_t j=0; j<corners.size(); j++) - for (size_t i=0; i<corners.size(); i++) + for (size_t j=0; j < corners.size(); ++j) + for (size_t i=0; i < corners.size(); ++i) mass[j].axpy(corners[i][j], jacobians[i]); - FieldMatrix<field_type,jDim,dim> trans; - Dune::MatrixVector::transpose(mass,trans); + FieldMatrix<field_type, jDim, dim> trans; + Dune::MatrixVector::transpose(mass, trans); auto inverse = mass.leftmultiplyany(trans); try { - inverse.invert(); - } catch(Dune::FMatrixError e) { - DUNE_THROW(Dune::Exception, "Computing dual covariant basis failed, metric not invertible: " << inverse.determinant() << " \n"); + inverse.invert(); + } catch(Dune::FMatrixError& e) { + DUNE_THROW(Dune::FMatrixError, "Computing dual covariant basis failed, metric not invertible: " << inverse.determinant() << " \n"); } auto dual = inverse.rightmultiplyany(trans);