Commit 88f09d23 authored by Dominic Kempf's avatar Dominic Kempf
Browse files

Store fibre intersection information by reference

parent 364f6038
......@@ -159,7 +159,7 @@ class EulerBernoulli2DLocalOperator
for (std::size_t fibindex=0; fibindex<fibre_parametrizations.size(); ++fibindex)
{
// Check whether this fiber actually intersects the cell
auto fibintersection = fibre_intersections[fibindex];
auto& fibintersection = fibre_intersections[fibindex];
auto it = fibintersection.element_fibre_intersections.find(is.index(entity));
if (it == fibintersection.element_fibre_intersections.end())
continue;
......@@ -282,7 +282,7 @@ class EulerBernoulli2DLocalOperator
bool flipped = false;
// Check whether this fiber actually intersects the cell
auto fibintersection = fibre_intersections[fibindex];
auto& fibintersection = fibre_intersections[fibindex];
auto it = fibintersection.facet_fibre_intersections.find(std::make_pair(is.index(ig.inside()), is.index(ig.outside())));
if (it == fibintersection.facet_fibre_intersections.end())
{
......
Markdown is supported
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