[Dune] PAMG

[Markus Blatt]

Hi,

You are using OpenMPI, are you not?

The guys from Debian told me that my bug concerning OpenMPI was fixed
and I switched back to supposedly working OpenMPI
version. Unfortunately this lead to the problems you described (no
convergence)!

Surprisingly one level solvers work. Which means that I was tempted to
use sophisticated MPI functionality in my AMG :( I learned my lesson
by now and will never use the newest MPI functionality again.

I will invest this problem again and maybe rewrite some stuff.

Please use lam or mpich in the meantime.

Markus

On Thu, Apr 30, 2009 at 05:37:54PM +0200, Bernd Flemisch wrote:
> Hey Markus,
> 
> danke für deine Antwort.
> 
> Mit SeqSSOR bricht er auch irgendwann ein, z.B. für N=1000, aber mit
> SeqSOR und einer Relaxation von 0.7 geht es recht gut.
> 
> Was mir noch aufgefallen ist: wenn man bei smootherArgs.iterations etwas
> größer 1 reinschreibt wird es schlechter. Wenn man im Konstruktor von
> AMG die entsprechende Anzahl an Iterationen eingibt, macht er das
> Gewünschte.
> 
> Noch eine Frage: Ich würde das gerne mal auf unsere Problemstellungen
> anwenden. Wie haben unsymmetrische, aus konvektionsdominierten Problemen
> entstehende Matrizen. Hast Du da vorab vielleicht Tipps, wie ich das AMG
> dafür am besten einstelle?
> 
> Herzlichen Dank! Viele Grüße
> Bernd
> 
> > Hi,
> >
> > tatsächlich.
> >
> > Zur allgemeinen Verwirrung habe ich anscheinend immer den test ILU
> > gemacht. Bei genauerem hinsehen funktioniert der aber nicht (schon gar
> > nicht parallel). Wäre interessant zu wissen, warum nicht.
> >
> > Einfach SeqSSOR und SeqSOR als Glätter nehmen.
> >
> > BTW: Das ist wirklich ein sehr sehr einfacher Test. Die
> > Daten werden nur in einer Dimension verteilt, was das ganze fur P->N
> > natürlich sehr ineffizient macht.
> >
> >
> > Gruß,
> >
> > Markus
> >
> > On Thu, Apr 30, 2009 at 02:07:56PM +0200, Bernd Flemisch wrote:
> >   
> >> Nachdem es bei uns jetzt mit der Parallelisierung endlich ernst geworden
> >> ist, habe ich mir wieder dein paralleles AMG angeschaut. Der Test
> >> pamgtest komipiliert und läuft prinzipiell! Für die
> >> Standardeinstellungen schaut mit N=100 auch die Konvergenz gut aus, da
> >> brauche ich mit 1, 2, 4 Prozessen 9, 20, 25 Iterationen.
> >>
> >> Leider geht bereits ab N = 120 mit mehr als einem Prozess die Konvergenz
> >> verloren. Hast Du eine Idee, woran das liegen könnte?
> >>
> >> Vielleicht setze ich auch einfach die Einstellungen nicht richtig? Kann
> >> man irgendwo nachlesen, wie man diese am besten wählt?
> >>
> >> Vielen Dank schon mal für deine Antwort! Grüße
> >> Bernd
> >>
> >>     
> >
> >   
> 
>