Skip to content
Snippets Groups Projects
Commit b49506c0 authored by Elias Pipping's avatar Elias Pipping
Browse files

[Extern] Entities are no longer returned as pointer

parent 2d9e986e
No related branches found
No related tags found
No related merge requests found
...@@ -39,7 +39,7 @@ class PointTractionBoundaryAssembler ...@@ -39,7 +39,7 @@ class PointTractionBoundaryAssembler
template <class TrialLocalFE, class BoundaryIterator> template <class TrialLocalFE, class BoundaryIterator>
void assemble(const BoundaryIterator& it, LocalVector& localVector, void assemble(const BoundaryIterator& it, LocalVector& localVector,
const TrialLocalFE& tFE) { const TrialLocalFE& tFE) {
auto const inside = it->inside();
localVector = 0.0; localVector = 0.0;
// geometry of the boundary face // geometry of the boundary face
...@@ -47,9 +47,8 @@ class PointTractionBoundaryAssembler ...@@ -47,9 +47,8 @@ class PointTractionBoundaryAssembler
it->geometry(); it->geometry();
// get quadrature rule // get quadrature rule
const int order = (it->inside()->type().isSimplex()) const int order =
? 2 * (order_ - 1) (inside.type().isSimplex()) ? 2 * (order_ - 1) : 2 * (order_ * dim - 1);
: 2 * (order_ * dim - 1);
// get quadrature rule // get quadrature rule
const Dune::QuadratureRule<ctype, dim - 1>& quad = const Dune::QuadratureRule<ctype, dim - 1>& quad =
...@@ -75,8 +74,8 @@ class PointTractionBoundaryAssembler ...@@ -75,8 +74,8 @@ class PointTractionBoundaryAssembler
// evaluate the displacement gradient at quad point of the element // evaluate the displacement gradient at quad point of the element
typename GridFunction::DerivativeType localDispGradient; typename GridFunction::DerivativeType localDispGradient;
if (displacement_->isDefinedOn(*it->inside())) if (displacement_->isDefinedOn(inside))
displacement_->evaluateDerivativeLocal(*it->inside(), elementQuadPos, displacement_->evaluateDerivativeLocal(inside, elementQuadPos,
localDispGradient); localDispGradient);
else else
displacement_->evaluateDerivative(segmentGeometry.global(quadPos), displacement_->evaluateDerivative(segmentGeometry.global(quadPos),
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment