Skip to content
Snippets Groups Projects
Commit c73aa08d authored by Oliver Sander's avatar Oliver Sander Committed by sander@PCPOOL.MI.FU-BERLIN.DE
Browse files

fix uses of subIndex and LevelBoundaryPatch 

[[Imported from SVN: r4473]]
parent fa27000d
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 2 additions and 2 deletions
src/rodassembler.cc
+ 2
− 2
View file @ c73aa08d
...@@ -216,7 +216,7 @@ getResultantForce(const LevelBoundaryPatch<GridType>& boundary, ...@@ -216,7 +216,7 @@ getResultantForce(const LevelBoundaryPatch<GridType>& boundary,
{ {
using namespace Dune; using namespace Dune;
if (grid_ != &boundary.getGrid()) if (grid_ != &boundary.gridView().grid())
DUNE_THROW(Dune::Exception, "The boundary patch has to match the grid of the assembler!"); DUNE_THROW(Dune::Exception, "The boundary patch has to match the grid of the assembler!");
const typename GridType::Traits::LeafIndexSet& indexSet = grid_->leafIndexSet(); const typename GridType::Traits::LeafIndexSet& indexSet = grid_->leafIndexSet();
...@@ -274,7 +274,7 @@ getResultantForce(const LevelBoundaryPatch<GridType>& boundary, ...@@ -274,7 +274,7 @@ getResultantForce(const LevelBoundaryPatch<GridType>& boundary,
// the canonical basis of R^3 // the canonical basis of R^3
FieldMatrix<double,3,3> orientationMatrix; FieldMatrix<double,3,3> orientationMatrix;
sol[indexSet.subIndex(*eIt,nIt->indexInInside())].q.matrix(orientationMatrix,1); sol[indexSet.subIndex(*eIt,nIt->indexInInside(),1)].q.matrix(orientationMatrix);
orientationMatrix.umv(localStress, canonicalStress); orientationMatrix.umv(localStress, canonicalStress);
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment