<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    Hi,<br>
    here is the mesh.<br>
    Mladen<br>
    <br>
    <div class="moz-cite-prefix">On 11/08/2013 10:52 AM, Oliver Sander
      wrote:<br>
    </div>
    <blockquote cite="mid:%3C527CB465.7020806@igpm.rwth-aachen.de%3E"
      type="cite">
      <pre wrap="">Hi Mladen,
the code looks okay to me.  Can you please send the simple.msh file you test with?
Thanks,
Oliver

Am 08.11.2013 10:37, schrieb Mladen Jurak:
</pre>
      <blockquote type="cite">
        <pre wrap="">Dear Dune,

I found a strange problem with GmshReader<>::read() function.
Sample code that I send you constructs well a UGGrid grid from "simple.msh" file.
But, if I uncomment the code that iterates over all vertices the
GmshReader<>::read() functions finishes with a segmentation fault.
In the other hand, iteration over elements does not affect the mesh reading.
Can you reproduce this error?

Best regards,
Mladen


#ifdef HAVE_CONFIG_H
# include "config.h"
#endif
#include <iostream>
#include <vector>

#include <dune/common/mpihelper.hh>
#include <dune/grid/uggrid.hh>
#include <dune/grid/io/file/gmshreader.hh>

int main(int argc, char** argv)
{
    Dune::MPIHelper& helper = Dune::MPIHelper::instance(argc, argv);
    // read grid
    const int dim =2;
    typedef Dune::UGGrid<dim> GridType;
    GridType * p_grid = Dune::GmshReader<GridType>::read("simple.msh");

    // iterate over all vertices
//    typedef typename GridType::LeafGridView LeafGridView;
//    LeafGridView leafView = p_grid->leafView();
//
//    typedef typename LeafGridView::template Codim<dim>::Iterator VertexLeafIterator;
//    for (VertexLeafIterator it  = leafView.template begin<dim>();
//                            it != leafView.template end<dim>(); ++it) {
//      //
//
//    }
    return 0;
}







_______________________________________________
Dune mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Dune@dune-project.org">Dune@dune-project.org</a>
<a class="moz-txt-link-freetext" href="http://lists.dune-project.org/mailman/listinfo/dune">http://lists.dune-project.org/mailman/listinfo/dune</a>

</pre>
      </blockquote>
      <pre wrap="">

</pre>
    </blockquote>
    <br>
    <div class="moz-signature">-- <br>
      <small>
        Mladen Jurak <br>
        Department of Mathematics <br>
        Faculty of Science <br>
        University of Zagreb <br>
        Bijenicka 30 <br>
        10000 Zagreb <br>
        Croatia <br>
        <br>
        Tel: +385 (0)1 460 5738 <br>
        Fax: +385 (0)1 468 0335 <br>
        E-mail: <a href="mailto:jurak@math.hr">jurak@math.hr</a> <br>
        Web: <a href="http://web.math.hr/%7Ejurak">http://web.math.hr/~jurak</a>
      </small>
    </div>
  </body>
</html>