[dune-pdelab] NonoverlappingConformingDirichletConstraints

[Tatiana Kim]

Hello,


I am having difficulties parallelising my code and would appreciate any help.


My grid is on a cube defined via YaspGrid:

YaspGrid<dim> grid(L,N,periodic,overlap,helper.getCommunicator());


I chose Nonoverlapping DD and defined the NonoverlappingConformingDirichletConstraints. I am using a vector grid function space:


   typedef PDELab::NonoverlappingConformingDirichletConstraints<GridView> CON;

    CON con(grid.leafGridView());


    typedef PDELab::VectorGridFunctionSpace<GridView,

    FEM,

    dim,

    PDELab::ISTLVectorBackend<PDELab::ISTLParameters::static_blocking>,

    PDELab::ISTLVectorBackend<>,

    CON,

    PDELab::EntityBlockedOrderingTag,

    PDELab::DefaultLeafOrderingTag> GFSF;


Up to this line, everything compiles without errors. When I add


    GFSF gfsf(grid.leafGridView(),fem,con);


I get this error:


/Users/tatianakim/Documents/DUNE/master/testproject/src/testproject.cc:330:38:   required from here

/opt/local/include/gcc49/c++/ext/new_allocator.h:120:4: error: no matching function for call to 'Dune::PDELab::NonoverlappingConformingDirichletConstraints<Dune::GridView<Dune::DefaultLeafGridViewTraits<const Dune::YaspGrid<3>, (Dune::PartitionIteratorType)4u> > >::NonoverlappingConformingDirichletConstraints()'

  { ::new((void *)__p) _Up(std::forward<_Args>(__args)...); }



Is it because NonoverlappingConformingDirichletConstraints are not compatible with Vector Grid Function Spaces?


Many thanks,

Tatiana.


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20150906/30867c5f/attachment.htm>