<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Dune team,<div><br></div><div>In my work I need to have access to pointers to a face on a boundary. Do you know how exactly I con do this? </div><div><br></div><div>My problem is mainly with the "intersect" iterator. Can I somehow produce an entityPointer to an intersect in which ( intersect->boundary() = true) ?</div><div><br></div><div>I used the following lines but it seems intersect iterator can not provide me with a pointer.</div><div><br></div><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "><span style="color: #911c67"><br></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "><span style="color: #911c67"><br></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "><span style="color: #911c67">for</span> (<span style="color: #016142">E<span class="Apple-style-span" style="color: rgb(0, 0, 0); ">Grid::<span style="color: #911c67">template</span> Codim<Grid::dimension-3>::<span style="color: #911c67">template</span> Partition<InteriorBorder_Partition>::LeafIterator</span></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "><span style="color: #016142"><span class="Apple-style-span" style="color: rgb(0, 0, 0); "></span><span class="Apple-style-span" style="color: rgb(0, 0, 0); "> elemit = elembit; elemit != elemeit; ++elemit)</span></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "> {</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; min-height: 15px; "><br></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; color: rgb(80, 143, 115); "><span style="color: #000000"> </span>/** Loop over intersections of the current element. */</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "> Grid::LeafIntersectionIterator begIntersectionIter = elemit->ileafbegin();</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "> Grid::LeafIntersectionIterator endIntersectionIter = elemit->ileafend();</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; min-height: 15px; "><br></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "> <span style="color: #911c67">for</span> (Grid::LeafIntersectionIterator intersect = begIntersectionIter; intersect != endIntersectionIter; ++intersect)</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "> {</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "> <span style="color: #911c67">if</span> (intersect->boundary() == <span style="color: #911c67">true</span>)</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "> {</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "> <span style="text-decoration: underline"><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; display: inline !important; "><span style="color: #775843">Grid</span>::<span style="color: #911c67">template</span> <span style="color: #016142">Codim</span><<span style="color: #775843">Grid</span>::dimension-2>::<span style="color: #016142">EntityPointer</span></div></span><span style="text-decoration: underline"> facepointer (*intersect)</span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "> }</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "> }</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "> }</div></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "><br></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 11px/normal Monaco; "><br></div><div><br></div><div>Another solution is to loop over face iterators and find out if they belong to the boundary. I could not find a related function for this. Is there any possibility in Dune to do that?</div><div><br></div><div>Thank you in advance for your assistance</div><div><br></div><div>Best regards</div><div><br></div><div>Arya Fallahi</div><div>
<br></div></body></html>