diff --git a/dune/gfe/averageinterface.hh b/dune/gfe/averageinterface.hh index e56b9496290a7c3afe9c8543693aff0a2a7f082d..05afb8e1ddd1ccf010f0e356b63044f74a9b5443 100644 --- a/dune/gfe/averageinterface.hh +++ b/dune/gfe/averageinterface.hh @@ -458,7 +458,7 @@ void computeTotalForceAndTorque(const BoundaryPatchBase<GridView>& interface, // Sum up the total torque \int (x - x_0) \times f dx Dune::FieldVector<double,dim> worldPos = it->geometry().global(quadPos); outputTorque.axpy(quad[qp].weight()*integrationElement, - crossProduct(worldPos - center, localPressure)); + crossProduct(Dune::FieldVector<double,dim>(worldPos - center), localPressure)); } @@ -568,7 +568,7 @@ void computeAveragePressure(const Dune::FieldVector<double,GridType::dimension>& phi_i[j] = shapeFunctionValues[i]; mu_tilde[i][j].axpy(quad[qp].weight() * integrationElement, - crossProduct(worldPos-crossSection.r, phi_i)); + crossProduct(Dune::FieldVector<double,dim>(worldPos-crossSection.r), phi_i)); } diff --git a/dune/gfe/makestraightrod.hh b/dune/gfe/makestraightrod.hh index 4facbc628ddbc2942cf44c81953e5fd22f08e1b0..50df6eab86e8b71c6f2dd996d2abc67c43f376c0 100644 --- a/dune/gfe/makestraightrod.hh +++ b/dune/gfe/makestraightrod.hh @@ -21,7 +21,7 @@ void makeStraightRod(std::vector<RigidBodyMotion<dim> >& rod, int n, Dune::FieldVector<double,3> zAxis(0); zAxis[2] = 1; - Dune::FieldVector<double,3> axis = crossProduct(end-beginning, zAxis); + Dune::FieldVector<double,3> axis = crossProduct(Dune::FieldVector<double,3>(end-beginning), zAxis); if (axis.two_norm() != 0) axis /= -axis.two_norm();