[Dune] Navier-Stokes parallel with UG and YaspGrid in 2D

manuel.hofmann at iwr.uni-heidelberg.de manuel.hofmann at iwr.uni-heidelberg.de
Tue Mar 20 10:10:47 CET 2012


Hi,

i tried to run the navier-stokes (cgstokes_instat) example in  
pdelab-howto in parallel mode. However i get with YaspGrid and UGGrid  
some interface communication errors. I used UG-3.9.1-patch7.

[findDDDInterfaces_:/home/mhofmann/diplomarbeit/dev/dune/dune-grid/dune/grid/uggrid.hh:968]: Edge communication not supported for interfaces of type  all / all  
interface

[comm:/home/mhofmann/diplomarbeit/dev/dune/dune-grid/dune/grid/yaspgrid.hh:2439]: interface communication not  
implemented

I tried ISTLBackend_OVLP_BCGS_ILU0 and ISTLBackend_NOVLP_CG_NOPREC as  
linear solver and newton solver of pdelab as nonlinear sovler. Both  
with the same error message. I am using the latest dune-grid module

This is how i load the grid with UG8(if thats important?):
Dune::MPIHelper& helper = Dune::MPIHelper::instance(argc, argv);
[...]
typedef Dune::UGGrid<2> GridType;
GridType grid;
std::string grid_file = "grids/cylinder.msh";
Dune::GridFactory<GridType> factory(&grid);
if(helper.rank()==0) {
   Dune::GmshReader<GridType>::read(factory,grid_file,true,false);
}
factory.createGrid();
[...]
grid.globalRefine(parameters.domain_level);
if(!Dune::MPIHelper::isFake)  grid.loadBalance();

---

Dune::FieldVector<double,2> L(1.0);
Dune::FieldVector<int,2> N(1);
Dune::FieldVector<bool,2> B(false);
int overlap=1;
Dune::YaspGrid<2> grid(helper.getCommunicator(),L,N,B,overlap);
typedef Dune::YaspGrid<2>::LeafGridView GV;


Am i doing something wrong or is this still not implemented?

Regards
Manuel




More information about the Dune mailing list