[dune-pdelab] nonoverlappingsinglephaseflow

[Felix Heimann]

Hi Bernd,

I think the changes to the boundary problems of problem "C" in commit to 
revision 96 probably broke something. However, problem "D" still works fine. So 
if you just want a parallel example replace

-    std::string problem="C";
+    std::string problem="D";

in nonoverlappingsinglephaseflow.cc . This problem realizes a simple pressure 
gradient  as boundary condition. So it should be a good parallel example to 
start from.

Best Regards,
Felix


On Donnerstag, 13. August 2009 12:10:50 Bernd Flemisch wrote:
> Dear Dune-pdelab,
>
> I tried to run the example "nonoverlappingsinglephaseflow" without
> modifying anything. The CG solver does not converge, I get the following
> output:
> parallel run on 1 process(es)
> GridFunctionSpace(static size version):
>  1 degrees of freedom in codim 2
> codim=2 offset=0 size=66049
> WARNING: cannot handle multiple geometry types in static size grid
> function space
> /0/ constrained dofs=514
> === matrix setup 0.416026 s
> === jacobian assembly 0.372023 s
> === CGSolver
> 25001        0.0174095
> === rate=0.999621, T=124.972, TIT=0.00499867, IT=25001
>
> Can you maybe check that and/or tell me from where to start when I want
> to test dune-pdelab in parallel?
>
> Thank you! Kind regards
> Bernd