[dune-pdelab] parallel pdelab for simple problem

Markus Blatt Markus.Blatt at iwr.uni-heidelberg.de
Tue Sep 27 15:39:20 CEST 2011


Hey Eike,

actually this should work out of the box.

I think it would make sense to post a bug report with a small sample
program to demonstrate the problem at
http://users.dune-project.org/projects/dune-pdelab 

Markus

On Tue, Sep 27, 2011 at 09:39:51AM +0100, Eike Mueller wrote:
> Dear dune-pdelab list,
> 
> I'm trying to adapt my serial code to run in parallel. Basically, I
> took  example 4 from the pdelab howto (I modified the local operator
> slightly) and use a 3d YaspGrid. I then modified my code according
> to section 3.3, i.e. I changed the constraints, use a parallel grid
> (again, YaspGrid with overlap 0) and changed to solver backend to
> one of the non-overlapping backends. However, the results are wrong,
> if I run with 8 cores then the solution is discontinuous across the
> faces that separate the eight domains. Also, alpha_skeleton of the
> local operator does not seem to get called for the faces between the
> domains. Is there an example of how to modify the local operator to
> make it parallel? Or do I have to choose a non-zero overlap in this
> case?
> 
> Thank you very much,
> 
> Eike
> 
> _______________________________________________
> dune-pdelab mailing list
> dune-pdelab at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune-pdelab
> 

-- 
Markus Blatt, Interdisciplinary Center for Scientific Computing,
Heidelberg University, Im Neuenheimer Feld 368, D-69120 Heidelberg
Tel: +49 (0) 6221 548881 Fax: +49 (0) 6221 548884

'Mathematics knows no races or geographic boundaries; for
mathematics, the cultural world is one country' - David Hilbert




More information about the dune-pdelab mailing list