****MJ-REJECTED**** RE: Non repetability issue and difference between 2.3.0 and 2.3.3

[Etienne PERCHAT]

Hi Barry, 

yes we are in mpich2 1.0.2p1.

Etienne




-----Message d'origine-----
De : owner-petsc-users at mcs.anl.gov [mailto:owner-petsc-users at mcs.anl.gov] De la part de Barry Smith
Envoyé : jeudi 25 septembre 2008 18:17
À : Etienne PERCHAT
Cc : petsc-users at mcs.anl.gov
Objet : Re: ****MJ-REJECTED**** RE: Non repetability issue and difference between 2.3.0 and 2.3.3



    Are the old and new version of PETSc both using the EXACT SAME MPI?

    Barry

On Sep 25, 2008, at 11:15 AM, Etienne PERCHAT wrote:

>
> Matt,
>
> Thanks a lot for your answers.
>
> Please find attached the asked files zipped (I used an ascii output).
>
> I tested with GMRES but I have PETSC_ERROR s. I don't understand  
> why ...
>
> [0]PETSC ERROR: KSPSolve_GMRES() line 231 in
> src/ksp/ksp/impls/gmres/d:\DEVELO~1\TEST_EP\PETSC-~1.3-P\PETSC-~1.3-P 
> \sr
> c\
> sp\ksp\impls\gmres\gmres.c
> [0]PETSC ERROR: KSPSolve() line 379 in
> src/ksp/ksp/interface/d:\DEVELO~1\TEST_EP\PETSC-~1.3-P\PETSC-~1.3-P 
> \src\
> ksp\ksp\
> NTERF~1\itfunc.c
> [2]PETSC ERROR: Object is in wrong state!
> [2]PETSC ERROR: Currently can use GMRES with only preconditioned
> residual (right preconditioning not coded)!
>
>
> I tried with BICG, BICGSTAB and FGMRS. It works and I notice the same
> behaviour.
>
> I tried using another preconditionner (ie SOR) with again PCR it is  
> even
> worse:
>
> Iteration= 171  residual=   1.164976140862921e+001
> Iteration= 171  residual=   1.171859643971711e+001
>
> I would like to stress also that it was inexact for me to say that the
> residual vary around iteration 40. It happens earlier around  
> increment 6
> (see below)
>
>
>
> *     Initial Residual               : 1.50972105381228e+006
> Iteration= 1    residual=   9.592364159822660e+004
> Iteration= 2    residual=   8.460449880782055e+004
> Iteration= 3    residual=   1.020430403443245e+005
> Iteration= 4    residual=   3.924181572770592e+004
> Iteration= 5    residual=   2.715343433370132e+004
> Iteration= 6    residual=   3.795381174726432e+004
> Iteration= 7    residual=   2.435841639251617e+004
> Iteration= 8    residual=   1.455427563694975e+004
> Iteration= 9    residual=   2.386233690946949e+004
> Iteration= 10   residual=   2.733829425421920e+004
> Iteration= 15   residual=   1.322501801249381e+004
> Iteration= 20   residual=   7.271229329992430e+003
>
>
> *     Initial Residual               : 1.50972105381228e+06
> Iteration= 1    residual=   9.592364159822660e+004
> Iteration= 2    residual=   8.460449880782055e+004
> Iteration= 3    residual=   1.020430403443245e+005
> Iteration= 4    residual=   3.924181572770592e+004
> Iteration= 5    residual=   2.715343433370132e+004
> Iteration= 6    residual=   3.795381174726426e+004
> Iteration= 7    residual=   2.435841639251620e+004
> Iteration= 8    residual=   1.455427563694976e+004
> Iteration= 9    residual=   2.386233690946967e+004
> Iteration= 10   residual=   2.733829425421909e+004
> Iteration= 15   residual=   1.322501801248572e+004
> Iteration= 20   residual=   7.271229329990202e+003
>
>
> Etienne
> <MatAndVec.zip>