[petsc-users] Réf. : Re: Unreproducible number of iterations for consecutives solves (GCR + ML)

[Thomas DE-SOZA]

Hi Dave,

Strangely, it only occurs when we use ML. Other preconditioners we have 
tested so far include BoomerAMG from Hypre and a homemade preconditioner 
based on a single precision factorization with MUMPS. Both are OK.

Thomas





dave.mayhem23 at gmail.com
27/06/2013 11:21
 
        Pour :  thomas.de-soza at edf.fr
        cc :    petsc-users at mcs.anl.gov
        Objet : Re: [petsc-users] Unreproducible number of iterations for 
consecutives solves (GCR + ML)


Hi Thomas,
Does this behaviour only occur when using ml, or do you see this with 
other preconditioners as well?



On Thursday, 27 June 2013, Thomas DE-SOZA wrote:

Dear PETSc users, 

We've been using the Krylov solvers in PETSc for a while in our software 
package for implicit structural mechanics and recently added algebraic 
multigrid preconditioning through the ML and Hypre libraries provided as 
part of PETSc. 
Though we're quite happy with this new feature, we've encountered a 
strange behaviour when using KSPGCR in conjunction with the PCML 
preconditioner. 
Indeed, identical and  consecutives solves inside the same run do not 
display the same number of Krylov iterations even in sequential and we 
were wondering why : 

1st solve 
723 KSP unpreconditioned resid norm 2.911385065051e-02 true resid norm 
2.911385065051e-02 ||r(i)||/||b|| 9.979426047969e-09 

2nd solve 
787 KSP unpreconditioned resid norm 2.896035212670e-02 true resid norm 
2.896035212670e-02 ||r(i)||/||b|| 9.926810982197e-09 

3rd solve 
721 KSP unpreconditioned resid norm 2.913123687343e-02 true resid norm 
2.913123687343e-02 ||r(i)||/||b|| 9.985385566274e-09 


Would you say this case requires a great number of iterations and 
therefore reproductibility is not ensured (indeed it is ill-conditioned 
and is not very suited for algebraic multigrid) ? 
Or is there a seed somewhere in ML that would explain this ? I searched 
the ML manual for that and couln't find any. We're using the uncoupled 
coarsening scheme in ML (call PetscOptionsSetValue('-pc_ml_CoarsenScheme', 
'Uncoupled', ierr)). 

Final notes : several consecutive runs of the program do display identical 
behaviour (that is the 1st solve always require 723 iterations, the 2nd 
787, etc). Moreover other cases that are well-conditioned require the same 
number number of iterations for each consecutive solve (though the 
converged residual differs a bit). 

Thanks for any hints, 
Thomas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130627/2bb8ad83/attachment.html>
-------------- next part --------------



Ce message et toutes les pi?ces jointes (ci-apr?s le 'Message') sont ?tablis ? l'intention exclusive des destinataires et les informations qui y figurent sont strictement confidentielles. Toute utilisation de ce Message non conforme ? sa destination, toute diffusion ou toute publication totale ou partielle, est interdite sauf autorisation expresse.

Si vous n'?tes pas le destinataire de ce Message, il vous est interdit de le copier, de le faire suivre, de le divulguer ou d'en utiliser tout ou partie. Si vous avez re?u ce Message par erreur, merci de le supprimer de votre syst?me, ainsi que toutes ses copies, et de n'en garder aucune trace sur quelque support que ce soit. Nous vous remercions ?galement d'en avertir imm?diatement l'exp?diteur par retour du message.

Il est impossible de garantir que les communications par messagerie ?lectronique arrivent en temps utile, sont s?curis?es ou d?nu?es de toute erreur ou virus.
____________________________________________________

This message and any attachments (the 'Message') are intended solely for the addressees. The information contained in this Message is confidential. Any use of information contained in this Message not in accord with its purpose, any dissemination or disclosure, either whole or partial, is prohibited except formal approval.

If you are not the addressee, you may not copy, forward, disclose or use any part of it. If you have received this message in error, please delete it and all copies from your system and notify the sender immediately by return message.

E-mail communication cannot be guaranteed to be timely secure, error or virus-free.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130627/2bb8ad83/attachment-0001.html>