[Dune] Problem with GmshReader

[Oliver Sander]

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:
> 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
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune
> 


-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 551 bytes
Desc: OpenPGP digital signature
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20131108/ce20de77/attachment.sig>