<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<p>Dear PDELab users,</p>
<p>I am trying to run my DG code for two phase flow on parallel
processes.</p>
<p>I use UGGrid because I plan to use GMSH for mesh generation. <br>
</p>
<p>I was wondering if anyone can help me with nonoverlapping case.</p>
<p>Some details of what I am doing:</p>
<blockquote><font color="#3333ff" size="-1">typedef
Dune::UGGrid<dim> GridType;<br>
Dune::FieldVector<typename GridType::ctype,dim> ll(0);<br>
Dune::FieldVector<typename GridType::ctype,dim> ur(1);<br>
std::array<unsigned int,dim> elements;<br>
std::fill(elements.begin(), elements.end(), 10);<br>
auto grid =
Dune::StructuredGridFactory<GridType>::createCubeGrid(lowerLeft,
upperRight, elements);<br>
grid->loadBalance();<br>
typedef GridType::LeafGridView GV;<br>
const GV& gv=grid->leafGridView();<br>
<br>
// GFS<br>
typedef Dune::PDELab::P0ParallelGhostConstraints CON0;<br>
typedef Dune::PDELab::istl::VectorBackend<> VBE0;<br>
typedef
Dune::PDELab::QkDGLocalFiniteElementMap<Coord,Real,1,dim>
FEM0;<br>
FEM0 fem0;<br>
typedef Dune::PDELab::GridFunctionSpace<GV, FEM0, CON0,
VBE0> GFS0;<br>
GFS0 gfs0(gv, fem0);<br>
typedef Dune::PDELab::istl::VectorBackend<
Dune::PDELab::istl::Blocking::none, block_size> VBE;<br>
typedef Dune::PDELab::LexicographicOrderingTag OrderingTag;<br>
typedef Dune::PDELab::PowerGridFunctionSpace<
GFS0,block_size,VBE, OrderingTag> GFS;<br>
<br>
// solution vector<br>
--U--<br>
<br>
// Local operator<br>
--LOP--<br>
<br>
// Grid operator<br>
--GO--<br>
<br>
// linear solver backend<br>
typedef Dune::PDELab::ISTLBackend_NOVLP_BCGS_SSORk<GO> LS;<br>
LS ls (go,5000,3,2);<br>
<br>
// Solver for non-linear problem<br>
typedef Dune::PDELab::Newton< GO, LS, U > SOLVER;<br>
SOLVER solver( go, u, ls );<br>
<br>
// Solve problem<br>
solver.apply(u);</font><br>
</blockquote>
<p>The P0ParallelGhostConstraints is what I tried in a 'last-ditch'
attempt to get rid of the seg-fault, but I am not sure if it is
correct to use it here... <br>
</p>
<p><br>
</p>
<p>I will really appreciate any help on this. <br>
</p>
<p>Thanks a lot in advance,</p>
<p>Best wishes,</p>
<p>Shubhangi<br>
</p>
<pre class="moz-signature" cols="72">--
Dr. Shubhangi Gupta
FE Marine Geosystems
Helmholtz Center for Ocean Research, GEOMAR
Wischhofstraße 1-3, D-24148 Kiel
Office:
Room: 12-206
Phone: +49 431 600-1402
Email: sgupta(at)geomar.de</pre>
</body>
</html>