<div dir="ltr"><div><div><div>Hi Jo<br><br><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote">I'm sorry, you appear to have included the wrong version of MeshTraversal in<br>
your mail.  At least I cannot see how it could have produced the output above,<br>
since it does not contain the string "my face".<br></blockquote><div> </div></div>My apologies for the confusion. I should have pasted the precise output from my code, here it is:<br><br></div><br><div>// -------------------------------------------------------------------------------------------------<br>dimension of the grid is: 2<br>visiting element 1(cube, 2)<br>visiting face 1(simplex, 1)<br>node 1 for face 1 is 1 0<br>node 2 for face 1 is 1 1<br>visiting face 2(simplex, 1)<br>node 1 for face 2 is 0 1<br>node 2 for face 2 is 1 1<br>visiting face 3(simplex, 1)<br>node 1 for face 3 is 0 0<br>node 2 for face 3 is 0 1<br>visiting face 4(simplex, 1)<br>node 1 for face 4 is 0 0<br>node 2 for face 4 is 1 0<br></div>// -------------------------------------------------------------------------------------------------<br>// end of output from MeshTraversal<br><br></div><div>I hope it is clear now, please write if it is not. Also please see the mesh data in gmsh format which I am using as input to the code above-<br></div><br>$MeshFormat<br>2.2 0 8<br>$EndMeshFormat<br>$Nodes<br>4<br>1 0 0 0<br>2 1 0 0<br>3 1 1 0<br>4 0 1 0<br>$EndNodes<br>$Elements<br>9<br>1 15 2 0 1 1<br>2 15 2 0 2 2<br>3 15 2 0 3 3<br>4 15 2 0 4 4<br>5 1 2 0 1 1 2<br>6 1 2 0 1 2 3<br>7 1 2 0 1 3 4<br>8 1 2 0 1 4 1<br>9 3 2 0 6 1 2 3 4<br>$EndElements<br><br><br><div>Many thanks for your help,<br></div>Ganesh<br><br><div><br><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Apr 10, 2015 at 9:27 PM, Jö Fahlke <span dir="ltr"><<a href="mailto:jorrit@jorrit.de" target="_blank">jorrit@jorrit.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Am Fri, 10. Apr 2015, 14:09:09 +0100 schrieb Ganesh Diwan:<br>
> Date: Fri, 10 Apr 2015 14:09:09 +0100<br>
> From: Ganesh Diwan <<a href="mailto:gcdiwan83@gmail.com">gcdiwan83@gmail.com</a>><br>
> To: Oliver Sander <<a href="mailto:sander@igpm.rwth-aachen.de">sander@igpm.rwth-aachen.de</a>>, DUNE Liste<br>
>  <<a href="mailto:dune@dune-project.org">dune@dune-project.org</a>><br>
> Subject: Re: [Dune] Local node numbering in Dune<br>
> X-No-Auth: unauthenticated sender<br>
> X-No-Relay: not in my network<br>
> X-Envelope-From: <<a href="mailto:gcdiwan83@gmail.com">gcdiwan83@gmail.com</a>><br>
<span class="">><br>
> Hi Oliver<br>
><br>
> When I loop through the element(s), I see that Dune returns element faces<br>
> that are in conflict with what the link above says. To be specific I get<br>
> my face 2 as Dune reference element face 0,<br>
> my face 0 as Dune reference element face 1,<br>
> my face 3 as Dune reference element face 2,  and<br>
> my face 1 as Dune reference element face 3,<br>
><br>
> Clearly I am doing something wrong!<br>
<br>
</span>I'm sorry, you appear to have included the wrong version of MeshTraversal in<br>
your mail.  At least I cannot see how it could have produced the output above,<br>
since it does not contain the string "my face".<br>
<br>
(In case I misunderstood and MeshTraversal was not meant to produce that<br>
output, please clarify how you obtained the numbers after "my face" and "Dune<br>
reference element face".)<br>
<br>
Regards,<br>
Jö.<br>
<div class="HOEnZb"><div class="h5"><br>
> I have:<br>
><br>
> typedef Dune::ALUGrid<2, 2, Dune::cube, Dune::nonconforming> GridTypeALU;<br>
> std::auto_ptr<GridTypeALU> mygrid( GmshReader<GridTypeALU>::read(<br>
> "./grids/hhgex.msh", true, true ) );<br>
> MeshTraversal (*mygrid); //<br>
><br>
> and then I have MeshTraversal as:<br>
><br>
> template<class G><br>
> void MeshTraversal(G& grid) {<br>
>     //first we extract the dimensions of the grid<br>
><br>
>     const int dim = G::dimension;<br>
>     std::cout << "dimension of the grid is: " << dim << std::endl;<br>
>     const int dimworld = G::dimensionworld;<br>
>     typedef typename G::ctype ct;<br>
>     typedef typename G::LeafGridView GridView;<br>
><br>
>     typedef typename GridView::template Codim<0>::Iterator LeafIterator;<br>
>     typedef typename GridView::IntersectionIterator IntersectionIterator;<br>
><br>
>     GridView gridView = grid.leafView();<br>
>     // get grid view on leaf part<br>
><br>
>     LeafIterator endit = gridView.template end<0>();<br>
><br>
>     int count = 0;<br>
>     //<br>
>     for (LeafIterator it = gridView.template begin<0>(); it != endit; ++it)<br>
> {<br>
><br>
>         count++;<br>
>         Dune::GeometryType gt = it->type();<br>
>         //<br>
>         std::cout << "visiting element " << count << gt << std::endl;<br>
><br>
>         IntersectionIterator isend = gridView.iend(*it);<br>
><br>
>         int countface = 0;<br>
>         for (IntersectionIterator is = gridView.ibegin(*it); is != isend;<br>
>                 ++is) {<br>
>             countface++;<br>
>             Dune::GeometryType gfacet = is->type();<br>
>             std::cout << "visiting face " << countface << gfacet <<<br>
> std::endl;<br>
>             std::cout << "node 1 for face" << countface << " is "<br>
>                     << is->geometry().corner(0) << std::endl;<br>
>             std::cout << "node 2 for face" << countface << " is "<br>
>                     << is->geometry().corner(1) << std::endl;<br>
>             //<br>
>         }<br>
>     }<br>
>     std::cin.get();<br>
> }<br>
><br>
><br>
><br>
> On Wed, Apr 8, 2015 at 3:57 PM, Ganesh Diwan <<a href="mailto:gcdiwan83@gmail.com">gcdiwan83@gmail.com</a>> wrote:<br>
><br>
> > Thankyou!<br>
> ><br>
> > Ganesh<br>
> ><br>
> > On Wed, Apr 8, 2015 at 3:56 PM, Oliver Sander <<a href="mailto:sander@igpm.rwth-aachen.de">sander@igpm.rwth-aachen.de</a>><br>
> > wrote:<br>
> ><br>
> >> Am 08.04.2015 um 16:50 schrieb Ganesh Diwan:<br>
> >> > Dear List<br>
> >> ><br>
> >> > It appears that after reading the grid from a Gmsh format file inside<br>
> >> Dune,<br>
> >> > the local dofs 2 and 3 i.e. node number 3 and 4 are swapped. This is<br>
> >> inside<br>
> >> > dune-grid//dune/grid/io/file/gmshreader.hh<br>
> >> ><br>
> >> > Is there a specific reason this is done? Is there a way to tell Dune to<br>
> >> use<br>
> >> > the same node numbering as in Gmsh file and not to alter it?<br>
> >><br>
> >> Hi Ganesh,<br>
> >> Dune has one fixed numbering for the local node numbers of all reference<br>
> >> elements<br>
> >> (see<br>
> >> <a href="http://www.dune-project.org/doc/doxygen/dune-geometry-html/group___geometry_reference_elements.html" target="_blank">http://www.dune-project.org/doc/doxygen/dune-geometry-html/group___geometry_reference_elements.html</a><br>
> >> )<br>
> >><br>
> >> These numberings differ from the ones used by gmsh, and hence the gmsh<br>
> >> reader<br>
> >> needs to do some node swapping.<br>
> >><br>
> >> Best,<br>
> >> Oliver<br>
> >><br>
> >> ><br>
> >> > Thankyou<br>
> >> > Ganesh<br>
> >> ><br>
> >> ><br>
> >> ><br>
> >> > _______________________________________________<br>
> >> > Dune mailing list<br>
> >> > <a href="mailto:Dune@dune-project.org">Dune@dune-project.org</a><br>
> >> > <a href="http://lists.dune-project.org/mailman/listinfo/dune" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
> >> ><br>
> >><br>
> >><br>
> >><br>
> >> _______________________________________________<br>
> >> Dune mailing list<br>
> >> <a href="mailto:Dune@dune-project.org">Dune@dune-project.org</a><br>
> >> <a href="http://lists.dune-project.org/mailman/listinfo/dune" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
> >><br>
> >><br>
> ><br>
<br>
> _______________________________________________<br>
> Dune mailing list<br>
> <a href="mailto:Dune@dune-project.org">Dune@dune-project.org</a><br>
> <a href="http://lists.dune-project.org/mailman/listinfo/dune" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
<br>
<br>
</div></div><span class="HOEnZb"><font color="#888888">--<br>
Jorrit (Jö) Fahlke, Institute for Computational und Applied Mathematics,<br>
University of Münster, Orleans-Ring 10, D-48149 Münster<br>
Tel: <a href="tel:%2B49%20251%2083%2035146" value="+492518335146">+49 251 83 35146</a> Fax: <a href="tel:%2B49%20251%2083%2032729" value="+492518332729">+49 251 83 32729</a><br>
<br>
Interpunktion, Orthographie und Grammatik der Email ist frei erfunden.<br>
Eine Übereinstimmung mit aktuellen oder ehemaligen Regeln wäre rein<br>
zufällig und ist nicht beabsichtigt.<br>
</font></span></blockquote></div><br></div>