[Dune] Problem with GmshReader

Oliver Sander sander at igpm.rwth-aachen.de
Fri Nov 8 12:47:08 CET 2013


Hi Mladen,
I ran your code using dune from git master, and there are no errors.
What version do you use?  Do you have built Dune with mpi installed?
--
Oliver

Am 08.11.2013 12:18, schrieb Mladen Jurak:
> Hi,
> here is the mesh.
> Mladen
> 
> On 11/08/2013 10:52 AM, Oliver Sander wrote:
>> 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
>>>
>>
> 
> 
> 
> _______________________________________________
> 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/756997ff/attachment.sig>


More information about the Dune mailing list