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

[Matthew Knepley]

On Thu, Sep 25, 2008 at 12:42 PM, Etienne PERCHAT
<etienne.perchat at transvalor.com> wrote:
> Hi Barry,
>
> yes we are in mpich2 1.0.2p1.

As before, we really need to see the system and rhs, or the whole code.

  Matt

> 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>
>
>
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which
their experiments lead.
-- Norbert Wiener