<div dir="ltr"><div><div><div>Hi Oliver<br><br></div>When I loop through the element(s), I see that Dune returns element faces 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></div><div><br>Clearly I am doing something wrong!<br><br></div>I have:<br><br>typedef Dune::ALUGrid<2, 2, Dune::cube, Dune::nonconforming> GridTypeALU;<br>std::auto_ptr<GridTypeALU> mygrid( GmshReader<GridTypeALU>::read( "./grids/hhgex.msh", true, true ) );<br>MeshTraversal (*mygrid); //<br><br></div>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>        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 << 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><div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Apr 8, 2015 at 3:57 PM, Ganesh Diwan <span dir="ltr"><<a href="mailto:gcdiwan83@gmail.com" target="_blank">gcdiwan83@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Thankyou!<span class="HOEnZb"><font color="#888888"><br><br></font></span></div><span class="HOEnZb"><font color="#888888">Ganesh<br></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Apr 8, 2015 at 3:56 PM, Oliver Sander <span dir="ltr"><<a href="mailto:sander@igpm.rwth-aachen.de" target="_blank">sander@igpm.rwth-aachen.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>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 Dune,<br>
> the local dofs 2 and 3 i.e. node number 3 and 4 are swapped. This is 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 use<br>
> the same node numbering as in Gmsh file and not to alter it?<br>
<br>
</div></div>Hi Ganesh,<br>
Dune has one fixed numbering for the local node numbers of all reference elements<br>
(see <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>
These numberings differ from the ones used by gmsh, and hence the gmsh 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" target="_blank">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" target="_blank">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></blockquote></div><br></div>
</div></div></blockquote></div><br></div>