<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;" dir="ltr">
<p>Hi,</p>
<p>I have some problem using (or maybe misusing :) ) a FilteredGridPart so I need your help.</p>
<p><br>
</p>
<p>I have a HostGridPart and a FilterdGridPart which uses a DomainFilter to select only a portion of the HostGridPart. Moreover <span>useFilteredIndexSet=true therefore the FilteredGridPart has its own IndexSet.</span></p>
<p><span><br>
</span></p>
<p><span>Now I have a P1 continuous Lagrange space for each grid part, let's call it HostSpace and FilteredSpace.</span></p>
<p><span>I want to assemble a LinearOperator which has as domain the filteredSpace and as range the hostSpace.</span></p>
<p><span><br>
</span></p>
<p><span>When I try to fill the stencil I do something like this:</span></p>
<p><span><br>
</span></p>
<p><span></p>
<div>Stencil< FilteredSpace, HostSpace > stencil( filteredSpace, hostSpace );</div>
<div>for( const auto& entity : filteredSpace )</div>
<div> stencil.fill( entity, filteredSpace.gridPart().convert( entity ) );</div>
<div><br>
</div>
<div>but I obtain the following runtime error</div>
<div><br>
</div>
<div><span>/usr/local_machine/dune-trunk/dune-fem/dune/fem/gridpart/codimindexset.hh:422: bool Dune::Fem::CodimIndexSet<GridImp>::checkValidIndex(const IndexType&) const [with GridImp = Dune::GeometryGrid<Dune::AlbertaGrid<2, 2>, Dune::Fem::VertexFunction<Dune::AlbertaGrid<2,
2> >, std::allocator<void> >; Dune::Fem::CodimIndexSet<GridImp>::IndexType = int]: Assertion `idx != invalidIndex()' failed.</span><br>
</div>
<div><br>
</div>
<div>I assume that there is some issue with the IndexSet. What am I doing wrong?<span style="font-family: Calibri, Arial, Helvetica, sans-serif, "Apple Color Emoji", "Segoe UI Emoji", NotoColorEmoji, "Segoe UI Symbol", "Android Emoji", EmojiSymbols; font-size: 16px;"></span></div>
<div><br>
</div>
<div>Thank you very much,</div>
<div>Cheers,</div>
<div>Marco</div>
<br>
</span>
<p></p>
</div>
</body>
</html>