<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<p>Hi Gregor,</p>
<p>I ran the following test for a yasp grid of 40x40 on two ranks
with overlap=1:</p>
<p>//do node- type report:<br>
uint interior_nodes = 0;<br>
uint border_nodes= 0;<br>
uint overlap_nodes= 0;<br>
uint front_nodes= 0;<br>
uint ghost_nodes= 0;</p>
<p> for (typename GridViewType::template
Codim<DIM>::Iterator vert = (*gv_p).template
begin<DIM>(); vert != (*gv_p).template end<DIM>();
++vert) {<br>
if (vert->partitionType() ==
Dune::PartitionType::OverlapEntity) overlap_nodes++;<br>
if (vert->partitionType() ==
Dune::PartitionType::InteriorEntity) interior_nodes++;<br>
if (vert->partitionType() ==
Dune::PartitionType::BorderEntity) border_nodes++;<br>
if (vert->partitionType() ==
Dune::PartitionType::FrontEntity) front_nodes++;<br>
if (vert->partitionType() ==
Dune::PartitionType::FrontEntity) ghost_nodes++;<br>
}</p>
<p> std::cout <<"--- rank: " << rank << "
VERTICES: interior: "<< interior_nodes<<"; border:
"<< border_nodes<<"; overlap: "<br>
<<overlap_nodes <<"; front: "<<
front_nodes<<"; ghost:
"<<ghost_nodes<<std::endl;</p>
<p> //do cell- type report:<br>
uint interior_cells = 0;<br>
uint border_cells= 0;<br>
uint overlap_cells= 0;<br>
uint front_cells= 0;<br>
uint ghost_cells= 0;</p>
<p> for (typename GridViewType::template
Codim<0>::Iterator cell = (*gv_p).template begin<0>();
cell != (*gv_p).template end<0>(); ++cell) {<br>
if (cell->partitionType() ==
Dune::PartitionType::OverlapEntity) overlap_cells++;<br>
if (cell->partitionType() ==
Dune::PartitionType::InteriorEntity) interior_cells++;<br>
if (cell->partitionType() ==
Dune::PartitionType::BorderEntity) border_cells++;<br>
if (cell->partitionType() ==
Dune::PartitionType::FrontEntity) front_cells++;<br>
if (cell->partitionType() ==
Dune::PartitionType::FrontEntity) ghost_cells++;<br>
}<br>
std::cout <<"--- rank: " << rank << "
CELLS: interior: "<< interior_cells<<"; border:
"<< border_cells<<"; overlap: "<br>
<<overlap_cells <<"; front: "<<
front_cells<<"; ghost:
"<<ghost_cells<<std::endl;</p>
<p><br>
</p>
<p>The output was:</p>
<p>--- rank: 0 VERTICES: interior: 820; border: 41; overlap: 0;
front: 41; ghost: 41<br>
--- rank: 1 VERTICES: interior: 820; border: 41; overlap: 0;
front: 41; ghost: 41</p>
<p>--- rank: 0 CELLS: interior: 800; border: 0; overlap: 40; front:
0; ghost: 0<br>
--- rank: 1 CELLS: interior: 800; border: 0; overlap: 40; front:
0; ghost: 0</p>
<p>so yes, I think that should do it, at least if you iterate over
cells.</p>
<p>lg,</p>
<p>Michael<br>
</p>
<p><br>
</p>
<p><br>
</p>
<p><br>
</p>
<p><br>
</p>
<div class="moz-cite-prefix">
On 11.12.18 14:33, Gregor Corbin wrote:<br>
</div>
<blockquote type="cite"
cite="mid:05e907be-4d20-1bf2-588c-f1a4f92ea1b4@mathematik.uni-kl.de">
<pre class="moz-quote-pre" wrap="">Hi,
i think that if the l2norm is modified with something like this in the
element-loop it will prevent double integration of overlapping parts:
if (eit->partitionType() != Dune::InteriorEntity) {continue;}
@Steffen: is this correct or am i overlooking aome subtleties in the
partitioning?
lg,
Gregor
Am 11.12.18 um 14:22 schrieb Michael Wenske:
</pre>
<blockquote type="cite">
<pre class="moz-quote-pre" wrap="">The l2norm.hh iterates via a templated iterator over the gridview. You
said it is not correct, does that mean
that the PDELab iterators iterate over the overlap-entities doubly?
lg,
Michael
</pre>
</blockquote>
<pre class="moz-quote-pre" wrap="">
</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
dune-pdelab mailing list
<a class="moz-txt-link-abbreviated" href="mailto:dune-pdelab@lists.dune-project.org">dune-pdelab@lists.dune-project.org</a>
<a class="moz-txt-link-freetext" href="https://lists.dune-project.org/mailman/listinfo/dune-pdelab">https://lists.dune-project.org/mailman/listinfo/dune-pdelab</a></pre>
</blockquote>
</body>
</html>