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
