December 2014 - petsc-dev - mailman3.cels.anl.gov

Re: [petsc-dev] [EXTERNAL] [pflotran-dev: 2700] Unintended consequence of CHKERRQ()
by Gautam Bisht 04 Dec '14

04 Dec '14

Hi Jed, Earlier I completely missed your point. I believe the utility of SNESSetFunctionDomainError is that an application can determine if the state (xx) being passed to compute residual function is meaningful (eg. a pair of (P,T) could result in negative density). If the state is not meaningful, the application could simply call SNESSetFunctionDomainError() and not update the residual vector. This would cause SNES to exit with SNES_DIVERGED_FUNCTION_DOMAIN. Looking at PFLOTRAN code more carefully, I now realize that Peter is using SNESSetFunctionDomainError() in FLASH2/MPHASE/MISCIBLE/IMMISCIBLE modes. I will now implement similar check for TH mode. Thanks, -Gautam. On Wed, Dec 3, 2014 at 12:13 PM, Gautam Bisht <gbisht(a)lbl.gov> wrote: > The user function was putting NaNs in the residual vector. I ended up > fixing the user function to avoid putting NaNs in the residual vector. > > -Gautam. > > > On Wed, Dec 3, 2014 at 10:48 AM, Jed Brown <jed(a)jedbrown.org> wrote: > >> "Hammond, Glenn E" <gehammo(a)sandia.gov> writes: >> >> > Can you explain why the NaNs or infs are being produced in the first >> place, or at least pinpoint the location and catch it there? >> >> If the SNESFunction computes NaN because the state is invalid, please >> call SNESSetFunctionDomainError and don't update the residual vector (or >> don't put NaN in there). SNES will clean up and return with >> SNESConvergedReason SNES_DIVERGED_FUNCTION_DOMAIN. >> >> We don't do this by default because detecting NaN collectively requires >> an extra reduction. Perhaps SNES could have a mode that automatically >> checks. >> >> > >> > Glenn >> > >> > From: pflotran-dev(a)googlegroups.com [mailto: >> pflotran-dev(a)googlegroups.com] On Behalf Of Gautam Bisht >> > Sent: Wednesday, December 03, 2014 8:36 AM >> > To: pflotran-dev(a)googlegroups.com >> > Subject: [EXTERNAL] [pflotran-dev: 2695] Unintended consequence of >> CHKERRQ() >> > >> > The issue is that when SNES diverges with a "not-a-number or infinite", >> the code crashes without giving timestepper_BE a chance to reduce the >> time-step. But, if CHKERRQ() is commented out after call to SNESSolve(), >> the code is able to reduce dt sufficiently to converge due to "itol_res" >> (PFLOTRAN's custom convergence reason == inorm_residual < >> solver%newton_inf_res_tol). Attached below is inputdeck (with an outrages >> initial dt = 1yr to illustrate the issue). Apart from commenting out >> CHKERRQ(), is there an alternate way to let timestepper_BE reduce dt? >> > >> > Below are relevant information to reproduce the error. >> > >> >>hg tip >> > changeset: 8541:fa2d4295aa14 >> > tag: tip >> > user: Nathan Collier <nathaniel.collier(a)gmail.com<mailto: >> nathaniel.collier(a)gmail.com>> >> > date: Wed Dec 03 08:52:00 2014 -0500 >> > summary: made output process-model identifier lines to be of >> uniform length (80 characters). Yeah it is a little OCD of me. >> > >> >>./pflotran -pflotranin da.in<http://da.in> >> > >> > == TH FLOW >> ===================================================================== >> > 1 2f: 2.55E+00 2x: 1.63E+07 2s: 1.72E+07 ir: 1.49E+00 iu: 3.10E+06 >> rsn: 0 >> > 2 2f: 1.88E-01 2x: 1.57E+14 2s: 1.57E+14 ir: 5.33E-02 iu: 4.37E+13 >> rsn: 0 >> > 3 2f: 1.07E-01 2x: 1.60E+14 2s: 4.28E+12 ir: 2.27E-02 iu: 9.90E+11 >> rsn: 0 >> > 4 2f: 1.24E+01 2x: 1.60E+14 2s: 9.53E+11 ir: 1.16E+01 iu: 7.15E+11 >> rsn: 0 >> > 5 2f: 1.29E+06 2x: 7.75E+20 2s: 7.75E+20 ir: 1.15E+06 iu: 7.21E+20 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 1[ 0] t= 0.00000E+00 dt= >> 5.00000E-01 >> > 1 2f: 2.22E+00 2x: 1.38E+07 2s: 1.47E+07 ir: 1.47E+00 iu: 3.06E+06 >> rsn: 0 >> > 2 2f: 2.94E-01 2x: 1.11E+14 2s: 1.11E+14 ir: 1.44E-01 iu: 4.05E+13 >> rsn: 0 >> > 3 2f: 6.05E+20 2x: 2.13E+14 2s: 1.79E+14 ir: 4.44E+20 iu: 8.93E+13 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 2[ 0] t= 0.00000E+00 dt= >> 2.50000E-01 >> > 1 2f: 2.17E+00 2x: 1.12E+07 2s: 1.21E+07 ir: 1.44E+00 iu: 3.01E+06 >> rsn: 0 >> > 2 2f: 5.55E-01 2x: 8.08E+13 2s: 8.08E+13 ir: 4.07E-01 iu: 3.66E+13 >> rsn: 0 >> > [0]PETSC ERROR: --------------------- Error Message >> -------------------------------------------------------------- >> > [0]PETSC ERROR: Floating point exception >> > [0]PETSC ERROR: Vec entry at local location 50 is not-a-number or >> infinite at end of function: Parameter number 3 >> > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html >> for trouble shooting. >> > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown >> > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a >> darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< >> http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:16:17 2014 >> > [0]PETSC ERROR: Configure options --download-hdf5=1 >> --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate >> --download-parmetis=yes --download-metis=yes >> --with-cc=/opt/local/bin/mpicc-openmpi-mp >> --with-cxx=/opt/local/bin/mpicxx-openmpi-mp >> --with-fc=/opt/local/bin/mpif90-openmpi-mp >> --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp >> PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 >> > [0]PETSC ERROR: #1 VecValidValues() line 34 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c >> > [0]PETSC ERROR: #2 SNESComputeFunction() line 2043 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > [0]PETSC ERROR: #3 SNESLineSearchApply_Basic() line 40 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c >> > [0]PETSC ERROR: #4 SNESLineSearchApply() line 563 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c >> > [0]PETSC ERROR: #5 SNESSolve_NEWTONLS() line 254 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c >> > [0]PETSC ERROR: #6 SNESSolve() line 3743 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > >> -------------------------------------------------------------------------- >> > MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD >> > with errorcode 72. >> > >> > NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. >> > You may or may not see output from other processes, depending on >> > exactly when Open MPI kills them. >> > >> -------------------------------------------------------------------------- >> > >> > >> >> hg diff >> > diff -r fa2d4295aa14 src/pflotran/timestepper_BE.F90 >> > --- a/src/pflotran/timestepper_BE.F90 Wed Dec 03 08:52:00 >> 2014 -0500 >> > +++ b/src/pflotran/timestepper_BE.F90 Wed Dec 03 08:22:35 2014 >> -0800 >> > @@ -310,7 +310,7 @@ >> > call PetscTime(log_start_time, ierr);CHKERRQ(ierr) >> > >> > call SNESSolve(solver%snes,PETSC_NULL_OBJECT, & >> > - process_model%solution_vec,ierr);CHKERRQ(ierr) >> > + process_model%solution_vec,ierr);!CHKERRQ(ierr) >> > >> > call PetscTime(log_end_time, ierr);CHKERRQ(ierr) >> > >> > >> >>./pflotran -pflotranin da.in<http://da.in> >> > >> > == TH FLOW >> ===================================================================== >> > 1 2f: 2.55E+00 2x: 1.63E+07 2s: 1.72E+07 ir: 1.49E+00 iu: 3.10E+06 >> rsn: 0 >> > 2 2f: 1.88E-01 2x: 1.57E+14 2s: 1.57E+14 ir: 5.33E-02 iu: 4.37E+13 >> rsn: 0 >> > 3 2f: 1.07E-01 2x: 1.60E+14 2s: 4.28E+12 ir: 2.27E-02 iu: 9.90E+11 >> rsn: 0 >> > 4 2f: 1.24E+01 2x: 1.60E+14 2s: 9.53E+11 ir: 1.16E+01 iu: 7.15E+11 >> rsn: 0 >> > 5 2f: 1.29E+06 2x: 7.75E+20 2s: 7.75E+20 ir: 1.15E+06 iu: 7.21E+20 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 1[ 0] t= 0.00000E+00 dt= >> 5.00000E-01 >> > 1 2f: 2.22E+00 2x: 1.38E+07 2s: 1.47E+07 ir: 1.47E+00 iu: 3.06E+06 >> rsn: 0 >> > 2 2f: 2.94E-01 2x: 1.11E+14 2s: 1.11E+14 ir: 1.44E-01 iu: 4.05E+13 >> rsn: 0 >> > 3 2f: 6.05E+20 2x: 2.13E+14 2s: 1.79E+14 ir: 4.44E+20 iu: 8.93E+13 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 2[ 0] t= 0.00000E+00 dt= >> 2.50000E-01 >> > 1 2f: 2.17E+00 2x: 1.12E+07 2s: 1.21E+07 ir: 1.44E+00 iu: 3.01E+06 >> rsn: 0 >> > 2 2f: 5.55E-01 2x: 8.08E+13 2s: 8.08E+13 ir: 4.07E-01 iu: 3.66E+13 >> rsn: 0 >> > [0]PETSC ERROR: --------------------- Error Message >> -------------------------------------------------------------- >> > [0]PETSC ERROR: Floating point exception >> > [0]PETSC ERROR: Vec entry at local location 50 is not-a-number or >> infinite at end of function: Parameter number 3 >> > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html >> for trouble shooting. >> > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown >> > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a >> darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< >> http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 >> > [0]PETSC ERROR: Configure options --download-hdf5=1 >> --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate >> --download-parmetis=yes --download-metis=yes >> --with-cc=/opt/local/bin/mpicc-openmpi-mp >> --with-cxx=/opt/local/bin/mpicxx-openmpi-mp >> --with-fc=/opt/local/bin/mpif90-openmpi-mp >> --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp >> PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 >> > [0]PETSC ERROR: #1 VecValidValues() line 34 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c >> > [0]PETSC ERROR: #2 SNESComputeFunction() line 2043 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > [0]PETSC ERROR: #3 SNESLineSearchApply_Basic() line 40 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c >> > [0]PETSC ERROR: #4 SNESLineSearchApply() line 563 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c >> > [0]PETSC ERROR: #5 SNESSolve_NEWTONLS() line 254 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c >> > [0]PETSC ERROR: #6 SNESSolve() line 3743 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > -> Cut time step: snes= 0 icut= 3[ 0] t= 0.00000E+00 dt= >> 1.25000E-01 >> > 1 2f: 2.12E+00 2x: 9.08E+06 2s: 9.85E+06 ir: 1.41E+00 iu: 2.94E+06 >> rsn: 0 >> > 2 2f: 7.89E-01 2x: 5.95E+13 2s: 5.95E+13 ir: 6.28E-01 iu: 3.20E+13 >> rsn: 0 >> > [0]PETSC ERROR: --------------------- Error Message >> -------------------------------------------------------------- >> > [0]PETSC ERROR: Floating point exception >> > [0]PETSC ERROR: Vec entry at local location 34 is not-a-number or >> infinite at end of function: Parameter number 3 >> > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html >> for trouble shooting. >> > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown >> > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a >> darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< >> http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 >> > [0]PETSC ERROR: Configure options --download-hdf5=1 >> --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate >> --download-parmetis=yes --download-metis=yes >> --with-cc=/opt/local/bin/mpicc-openmpi-mp >> --with-cxx=/opt/local/bin/mpicxx-openmpi-mp >> --with-fc=/opt/local/bin/mpif90-openmpi-mp >> --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp >> PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 >> > [0]PETSC ERROR: #7 VecValidValues() line 34 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c >> > [0]PETSC ERROR: #8 SNESComputeFunction() line 2043 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > [0]PETSC ERROR: #9 SNESLineSearchApply_Basic() line 40 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c >> > [0]PETSC ERROR: #10 SNESLineSearchApply() line 563 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c >> > [0]PETSC ERROR: #11 SNESSolve_NEWTONLS() line 254 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c >> > [0]PETSC ERROR: #12 SNESSolve() line 3743 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > -> Cut time step: snes= 0 icut= 4[ 0] t= 0.00000E+00 dt= >> 6.25000E-02 >> > 1 2f: 2.04E+00 2x: 7.38E+06 2s: 8.06E+06 ir: 1.36E+00 iu: 2.86E+06 >> rsn: 0 >> > 2 2f: 5.60E+04 2x: 4.28E+13 2s: 4.28E+13 ir: 4.47E+04 iu: 2.66E+13 >> rsn: 0 >> > 3 2f: 1.16E+18 2x: 4.91E+13 2s: 2.30E+13 ir: 9.47E+17 iu: 1.69E+13 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 5[ 0] t= 0.00000E+00 dt= >> 3.12500E-02 >> > 1 2f: 1.94E+00 2x: 6.00E+06 2s: 6.57E+06 ir: 1.30E+00 iu: 2.74E+06 >> rsn: 0 >> > 2 2f: 6.76E+04 2x: 2.91E+13 2s: 2.91E+13 ir: 5.18E+04 iu: 2.06E+13 >> rsn: 0 >> > [0]PETSC ERROR: --------------------- Error Message >> -------------------------------------------------------------- >> > [0]PETSC ERROR: Floating point exception >> > [0]PETSC ERROR: Vec entry at local location 14 is not-a-number or >> infinite at end of function: Parameter number 3 >> > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html >> for trouble shooting. >> > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown >> > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a >> darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< >> http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 >> > [0]PETSC ERROR: Configure options --download-hdf5=1 >> --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate >> --download-parmetis=yes --download-metis=yes >> --with-cc=/opt/local/bin/mpicc-openmpi-mp >> --with-cxx=/opt/local/bin/mpicxx-openmpi-mp >> --with-fc=/opt/local/bin/mpif90-openmpi-mp >> --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp >> PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 >> > [0]PETSC ERROR: #13 VecValidValues() line 34 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c >> > [0]PETSC ERROR: #14 SNESComputeFunction() line 2043 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > [0]PETSC ERROR: #15 SNESLineSearchApply_Basic() line 40 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c >> > [0]PETSC ERROR: #16 SNESLineSearchApply() line 563 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c >> > [0]PETSC ERROR: #17 SNESSolve_NEWTONLS() line 254 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c >> > [0]PETSC ERROR: #18 SNESSolve() line 3743 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > -> Cut time step: snes= 0 icut= 6[ 0] t= 0.00000E+00 dt= >> 1.56250E-02 >> > 1 2f: 1.81E+00 2x: 4.84E+06 2s: 5.31E+06 ir: 1.21E+00 iu: 2.58E+06 >> rsn: 0 >> > 2 2f: 1.56E+04 2x: 1.81E+13 2s: 1.81E+13 ir: 1.22E+04 iu: 1.43E+13 >> rsn: 0 >> > 3 2f: 9.46E+15 2x: 1.99E+13 2s: 7.41E+12 ir: 7.34E+15 iu: 5.84E+12 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 7[ 0] t= 0.00000E+00 dt= >> 7.81250E-03 >> > 1 2f: 1.64E+00 2x: 3.85E+06 2s: 4.21E+06 ir: 1.11E+00 iu: 2.37E+06 >> rsn: 0 >> > 2 2f: 1.37E+04 2x: 9.82E+12 2s: 9.82E+12 ir: 1.02E+04 iu: 8.52E+12 >> rsn: 0 >> > 3 2f: 5.35E+14 2x: 1.86E+13 2s: 1.65E+13 ir: 3.78E+14 iu: 1.62E+13 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 8[ 0] t= 0.00000E+00 dt= >> 3.90625E-03 >> > 1 2f: 1.42E+00 2x: 3.01E+06 2s: 3.25E+06 ir: 9.67E-01 iu: 2.09E+06 >> rsn: 0 >> > 2 2f: 6.03E+03 2x: 4.39E+12 2s: 4.39E+12 ir: 4.38E+03 iu: 4.08E+12 >> rsn: 0 >> > [0]PETSC ERROR: --------------------- Error Message >> -------------------------------------------------------------- >> > [0]PETSC ERROR: Floating point exception >> > [0]PETSC ERROR: Vec entry at local location 30 is not-a-number or >> infinite at end of function: Parameter number 3 >> > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html >> for trouble shooting. >> > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown >> > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a >> darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< >> http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 >> > [0]PETSC ERROR: Configure options --download-hdf5=1 >> --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate >> --download-parmetis=yes --download-metis=yes >> --with-cc=/opt/local/bin/mpicc-openmpi-mp >> --with-cxx=/opt/local/bin/mpicxx-openmpi-mp >> --with-fc=/opt/local/bin/mpif90-openmpi-mp >> --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp >> PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 >> > [0]PETSC ERROR: #19 VecValidValues() line 34 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c >> > [0]PETSC ERROR: #20 SNESComputeFunction() line 2043 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > [0]PETSC ERROR: #21 SNESLineSearchApply_Basic() line 40 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c >> > [0]PETSC ERROR: #22 SNESLineSearchApply() line 563 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c >> > [0]PETSC ERROR: #23 SNESSolve_NEWTONLS() line 254 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c >> > [0]PETSC ERROR: #24 SNESSolve() line 3743 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > -> Cut time step: snes= 0 icut= 9[ 0] t= 0.00000E+00 dt= >> 1.95312E-03 >> > 1 2f: 1.16E+00 2x: 2.30E+06 2s: 2.42E+06 ir: 8.00E-01 iu: 1.76E+06 >> rsn: 0 >> > 2 2f: 4.32E+03 2x: 1.47E+12 2s: 1.47E+12 ir: 3.07E+03 iu: 1.42E+12 >> rsn: 0 >> > [0]PETSC ERROR: --------------------- Error Message >> -------------------------------------------------------------- >> > [0]PETSC ERROR: Floating point exception >> > [0]PETSC ERROR: Vec entry at local location 4 is not-a-number or >> infinite at end of function: Parameter number 3 >> > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html >> for trouble shooting. >> > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown >> > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a >> darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< >> http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 >> > [0]PETSC ERROR: Configure options --download-hdf5=1 >> --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate >> --download-parmetis=yes --download-metis=yes >> --with-cc=/opt/local/bin/mpicc-openmpi-mp >> --with-cxx=/opt/local/bin/mpicxx-openmpi-mp >> --with-fc=/opt/local/bin/mpif90-openmpi-mp >> --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp >> PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 >> > [0]PETSC ERROR: #25 VecValidValues() line 34 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c >> > [0]PETSC ERROR: #26 SNESComputeFunction() line 2043 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > [0]PETSC ERROR: #27 SNESLineSearchApply_Basic() line 40 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c >> > [0]PETSC ERROR: #28 SNESLineSearchApply() line 563 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c >> > [0]PETSC ERROR: #29 SNESSolve_NEWTONLS() line 254 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c >> > [0]PETSC ERROR: #30 SNESSolve() line 3743 in >> /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c >> > -> Cut time step: snes= 0 icut= 10[ 0] t= 0.00000E+00 dt= >> 9.76562E-04 >> > 1 2f: 8.79E-01 2x: 1.73E+06 2s: 1.70E+06 ir: 6.11E-01 iu: 1.37E+06 >> rsn: 0 >> > 2 2f: 1.34E+03 2x: 3.17E+11 2s: 3.17E+11 ir: 9.30E+02 iu: 3.14E+11 >> rsn: 0 >> > 3 2f: 3.99E+13 2x: 6.13E+11 2s: 3.04E+11 ir: 3.26E+13 iu: 2.98E+11 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 11[ 0] t= 0.00000E+00 dt= >> 4.88281E-04 >> > 1 2f: 5.96E-01 2x: 1.33E+06 2s: 1.10E+06 ir: 4.19E-01 iu: 9.69E+05 >> rsn: 0 >> > 2 2f: 4.60E+04 2x: 3.54E+10 2s: 3.54E+10 ir: 3.25E+04 iu: 3.54E+10 >> rsn: 0 >> > 3 2f: 9.64E+15 2x: 3.91E+11 2s: 3.55E+11 ir: 7.87E+15 iu: 3.55E+11 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 12[ 0] t= 0.00000E+00 dt= >> 2.44141E-04 >> > 1 2f: 3.54E-01 2x: 1.11E+06 2s: 6.36E+05 ir: 2.50E-01 iu: 5.95E+05 >> rsn: 0 >> > 2 2f: 5.42E+04 2x: 1.33E+09 2s: 1.33E+09 ir: 3.93E+04 iu: 1.33E+09 >> rsn: 0 >> > 3 2f: 3.51E+06 2x: 1.75E+10 2s: 1.75E+10 ir: 2.87E+06 iu: 1.75E+10 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 13[ 0] t= 0.00000E+00 dt= >> 1.22070E-04 >> > 1 2f: 5.74E-02 2x: 1.04E+06 2s: 3.17E+05 ir: 4.02E-02 iu: 3.07E+05 >> rsn: 0 >> > 2 2f: 3.60E-02 2x: 6.66E+09 2s: 6.66E+09 ir: 2.88E-02 iu: 6.66E+09 >> rsn: 0 >> > 3 2f: 5.16E+02 2x: 2.91E+08 2s: 6.95E+09 ir: 3.31E+02 iu: 6.95E+09 >> rsn: 0 >> > 4 2f: 7.20E+04 2x: 7.65E+15 2s: 7.65E+15 ir: 6.73E+04 iu: 7.65E+15 >> rsn: 0 >> > 5 2f: 4.01E+12 2x: 3.15E+16 2s: 3.92E+16 ir: 3.27E+12 iu: 3.92E+16 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 14[ 0] t= 0.00000E+00 dt= >> 6.10352E-05 >> > 1 2f: 8.13E-02 2x: 1.04E+06 2s: 1.33E+05 ir: 8.02E-02 iu: 1.31E+05 >> rsn: 0 >> > 2 2f: 1.94E+02 2x: 1.14E+08 2s: 1.14E+08 ir: 1.33E+02 iu: 1.14E+08 >> rsn: 0 >> > 3 2f: 2.16E-01 2x: 4.22E+08 2s: 5.36E+08 ir: 1.79E-01 iu: 5.36E+08 >> rsn: 0 >> > 4 2f: 3.21E+28 2x: 1.07E+19 2s: 1.07E+19 ir: 2.62E+28 iu: 1.07E+19 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 15[ 0] t= 0.00000E+00 dt= >> 3.05176E-05 >> > 1 2f: 1.13E-05 2x: 1.05E+06 2s: 4.69E+04 ir: 1.13E-05 iu: 4.67E+04 >> rsn: 0 >> > 2 2f: 1.21E-10 2x: 1.05E+06 2s: 6.87E+02 ir: 1.05E-10 iu: 6.87E+02 >> rsn: itol_res >> > >> > Step 1 Time= 3.05176E-05 Dt= 3.05176E-05 [y] snes_conv_reason: >> 10 >> > newton = 47 [ 47] linear = 52 [ 52] cuts = 15 [ 15] >> > --> SNES Linear/Non-Linear Iterations = 2 / 2 >> > --> SNES Residual: 1.209152E-10 1.209152E-12 1.045172E-10 >> > --> max chng: dpmx= 4.5981E+04 dtmpmx= 1.6500E+00 >> > >> > >> > == TH FLOW >> ===================================================================== >> > 1 2f: 1.62E-01 2x: 1.04E+06 2s: 7.34E+04 ir: 1.60E-01 iu: 7.27E+04 >> rsn: 0 >> > 2 2f: 7.18E+01 2x: 6.44E+07 2s: 6.44E+07 ir: 4.68E+01 iu: 6.44E+07 >> rsn: 0 >> > 3 2f: 2.52E-01 2x: 5.48E+07 2s: 1.19E+08 ir: 2.02E-01 iu: 1.19E+08 >> rsn: 0 >> > 4 2f: 1.29E+09 2x: 1.44E+18 2s: 1.44E+18 ir: 1.29E+09 iu: 1.44E+18 >> rsn: max_norm >> > -> Cut time step: snes=-10 icut= 1[ 15] t= 3.05176E-05 dt= >> 1.52588E-05 >> > 1 2f: 3.23E-01 2x: 1.04E+06 2s: 3.15E+04 ir: 3.19E-01 iu: 3.13E+04 >> rsn: 0 >> > 2 2f: 1.17E-03 2x: 3.08E+06 2s: 2.83E+06 ir: 1.09E-03 iu: 2.83E+06 >> rsn: 0 >> > 3 2f: 3.23E-01 2x: 1.04E+06 2s: 2.83E+06 ir: 3.19E-01 iu: 2.83E+06 >> rsn: 0 >> > 4 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 5 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 6 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 7 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 8 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 9 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 10 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 11 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 12 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 13 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 14 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 15 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 16 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 17 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 18 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 19 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 20 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 21 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 22 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 23 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 24 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 25 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 26 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 27 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 28 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 29 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 30 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 31 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 32 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 33 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 34 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 35 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 36 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 37 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 38 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 39 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 40 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 41 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 42 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 43 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 44 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 45 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 46 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 47 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 48 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 49 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 50 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 51 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 52 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 53 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 54 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 55 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 56 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 57 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 58 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 59 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 60 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 61 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 62 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 63 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 64 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 65 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 66 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 67 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 68 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 69 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 70 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 71 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 72 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 73 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 74 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 75 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 76 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 77 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 78 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 79 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 80 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 81 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 82 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 83 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 84 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 85 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 86 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 87 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 88 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 89 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 90 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 91 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 92 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 93 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 94 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 95 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 96 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 97 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 98 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > 99 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 >> rsn: 0 >> > 100 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 >> rsn: 0 >> > -> Cut time step: snes= -5 icut= 2[ 15] t= 3.05176E-05 dt= >> 7.62939E-06 >> > 1 2f: 6.23E-01 2x: 1.04E+06 2s: 1.39E+04 ir: 6.15E-01 iu: 1.38E+04 >> rsn: 0 >> > 2 2f: 6.18E-04 2x: 1.06E+06 2s: 1.15E+05 ir: 5.25E-04 iu: 1.15E+05 >> rsn: 0 >> > 3 2f: 6.23E-01 2x: 1.04E+06 2s: 1.15E+05 ir: 6.15E-01 iu: 1.15E+05 >> rsn: 0 >> > 4 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 5 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 6 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 7 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 8 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 9 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 10 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 11 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 12 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 13 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 14 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 15 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 16 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 17 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 18 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 19 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 20 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 21 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 22 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 23 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 24 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 25 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 26 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 27 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 28 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 29 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 30 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 31 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 32 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 33 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 34 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 35 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 36 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 37 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 38 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 39 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 40 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 41 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 42 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 43 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 44 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 45 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 46 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 47 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 48 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 49 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 50 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 51 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 52 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 53 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 54 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 55 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 56 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 57 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 58 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 59 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 60 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 61 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 62 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 63 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 64 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 65 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 66 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 67 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 68 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 69 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 70 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 71 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 72 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 73 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 74 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 75 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 76 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 77 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 78 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 79 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 80 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 81 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 82 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 83 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 84 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 85 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 86 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 87 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 88 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 89 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 90 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 91 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 92 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 93 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 94 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 95 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 96 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 97 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 98 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > 99 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 >> rsn: 0 >> > 100 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 >> rsn: 0 >> > -> Cut time step: snes= -5 icut= 3[ 15] t= 3.05176E-05 dt= >> 3.81470E-06 >> > 1 2f: 8.15E-01 2x: 1.05E+06 2s: 6.36E+03 ir: 8.05E-01 iu: 6.33E+03 >> rsn: 0 >> > 2 2f: 1.34E-03 2x: 1.05E+06 2s: 4.20E+03 ir: 1.31E-03 iu: 4.19E+03 >> rsn: 0 >> > 3 2f: 4.65E-04 2x: 1.05E+06 2s: 2.99E+02 ir: 4.59E-04 iu: 2.50E+02 >> rsn: 0 >> > 4 2f: 1.43E-04 2x: 1.05E+06 2s: 1.18E+02 ir: 1.41E-04 iu: 1.18E+02 >> rsn: 0 >> > 5 2f: 4.35E-05 2x: 1.05E+06 2s: 8.37E+01 ir: 4.30E-05 iu: 8.37E+01 >> rsn: 0 >> > 6 2f: 1.27E-05 2x: 1.05E+06 2s: 5.80E+01 ir: 1.25E-05 iu: 5.80E+01 >> rsn: 0 >> > 7 2f: 3.26E-06 2x: 1.05E+06 2s: 3.69E+01 ir: 3.22E-06 iu: 3.69E+01 >> rsn: 0 >> > 8 2f: 5.57E-07 2x: 1.05E+06 2s: 1.83E+01 ir: 5.50E-07 iu: 1.83E+01 >> rsn: 0 >> > 9 2f: 2.98E-08 2x: 1.05E+06 2s: 4.65E+00 ir: 2.95E-08 iu: 4.65E+00 >> rsn: 0 >> > 10 2f: 9.92E-11 2x: 1.05E+06 2s: 2.79E-01 ir: 9.78E-11 iu: 2.79E-01 >> rsn: itol_res >> > >> > Step 2 Time= 3.43323E-05 Dt= 3.81470E-06 [y] snes_conv_reason: >> 10 >> > newton = 214 [ 261] linear = 214 [ 266] cuts = 3 [ 18] >> > --> SNES Linear/Non-Linear Iterations = 10 / 10 >> > --> SNES Residual: 9.921423E-11 9.921423E-13 9.778422E-11 >> > --> max chng: dpmx= 1.6597E+03 dtmpmx= 1.9780E-01 >> > >> > >> > -- >> > You received this message because you are subscribed to the Google >> Groups "pflotran-dev" group. >> > To view this discussion on the web visit >> https://groups.google.com/d/msgid/pflotran-dev/CAPz1TnfKQJov5Uxvghp_31uKpxz… >> < >> https://groups.google.com/d/msgid/pflotran-dev/CAPz1TnfKQJov5Uxvghp_31uKpxz… >> >. >> > >> > -- >> > You received this message because you are subscribed to the Google >> Groups "pflotran-dev" group. >> > To view this discussion on the web visit >> https://groups.google.com/d/msgid/pflotran-dev/92b8b510fceb47309cc0b8f7ada4… >> . >> >> -- >> You received this message because you are subscribed to the Google Groups >> "pflotran-dev" group. >> To view this discussion on the web visit >> https://groups.google.com/d/msgid/pflotran-dev/87mw748t8f.fsf%40jedbrown.org >> . >> > >

1 0

Re: [petsc-dev] bitbucket giving wrong information
by Barry Smith 04 Dec '14

04 Dec '14

> On Dec 3, 2014, at 10:43 PM, Patrick Sanan <patrick.sanan(a)gmail.com> wrote: > > Whoops, guess that wasn't the cause then! > > Actually, looking at bitbucket, it looks like the PR is coming from a different repository (https://bitbucket.org/iaannn/petsc/) > > This worked on my machine: > [petsc (master=)]$ git remote add iaannn [email protected]:iaannn/petsc.git > [petsc (master=)]$ git fetch iaannn > [petsc (master=)]$ git checkout -b ianwilliamson/fix-matlab-write-bin-complex iaannn/ianwilliamson/fix-matlab-write-bin-complex > [petsc (ianwilliamson/fix-matlab-write-bin-complex=)]$ Sure thanks, but there is no way I can remember this; hence bitbucket should walk me through it. Barry > > On 04/12/14 12:57 PM, Barry Smith wrote: >> $ git checkout -b ianwilliamson/fix-matlab-write-bin-complex origin/ianwilliamson/fix-matlab-write-bin-complex >> fatal: Cannot update paths and switch to branch 'ianwilliamson/fix-matlab-write-bin-complex' at the same time. >> Did you intend to checkout 'origin/ianwilliamson/fix-matlab-write-bin-complex' which can not be resolved as commit? >> >> >> >>> On Dec 3, 2014, at 9:54 PM, Patrick Sanan <patrick.sanan(a)gmail.com> wrote: >>> >>> Is it possible that you have a second remote with a branch with the same name? In that case, perhaps something like: >>> >>> git checkout -b ianwilliamson/fix-matlab-write-bin-complex origin/ianwilliamson/fix-matlab-write-bin-complex >>> >>> On 04/12/14 12:37 PM, Barry Smith wrote: >>>> So some guy posts a pull request. I, in trying to follow PETSc-rules, want to check out his branch and then merge that rather than use the merge directly from the bitbucket website. So I go to the branch in the bitbucket website and bitbucket claims I can look at it with >>>> >>>> git fetch && git checkout ianwilliamson/fix-matlab-write-bin-complex >>>> >>>> (see attachment). I try that and of course it does not work: >>>> >>>> $ git fetch && git checkout ianwilliamson/fix-matlab-write-bin-complex >>>> remote: Counting objects: 56, done. >>>> remote: Compressing objects: 100% (56/56), done. >>>> remote: Total 56 (delta 48), reused 0 (delta 0) >>>> Unpacking objects: 100% (56/56), done. >>>> From bitbucket.org:petsc/petsc >>>> f685883..cf8fc5e master -> origin/master >>>> 7cadcfe..39ea070 knepley/fix-plex-orient -> origin/knepley/fix-plex-orient >>>> 565f831..ce993ff knepley/pylith -> origin/knepley/pylith >>>> + 3f11d05...21345e8 mlange/refactor-plex-distribute -> origin/mlange/refactor-plex-distribute (forced update) >>>> e9d5acc..916966e next -> origin/next >>>> error: pathspec 'ianwilliamson/fix-matlab-write-bin-complex' did not match any file(s) known to git. >>>> >>>> So is bitbucket lying to me or am I doing something wrong. Should this be a feature request to bitbucket ? Remember I am living in the world where I cannot remember how to do anything and thus expect my tools to tell me how to do something. So why cannot bitbucket tell me the truth of how to get this branch? >>>> >>>> Barry >>>> >>>> >>>> >>>> >>>> <Mail Attachment.png> >

1 0

Re: [petsc-dev] Can ViennaCL coexist with Cusp?
by Karl Rupp 04 Dec '14

04 Dec '14

Hi Jose, > It seems ViennaCL and Cusp are exclusive. Yes, this is the current status. Needs to be resolved with the next bigger refactoring. > So I guess configure should complain when both are enabled at configure time, otherwise the build fails. I fully agree and will take care of this. Best regards, Karli

1 0

Can ViennaCL coexist with Cusp?
by Jose E. Roman 04 Dec '14

04 Dec '14

It seems ViennaCL and Cusp are exclusive. From vecimpl.h: #if defined(PETSC_HAVE_CUSP) PetscCUSPFlag valid_GPU_array; /* indicates where the most recently modified vector data is (GPU or CPU) */ void *spptr; /* if we're using CUSP, then this is the special pointer to the array on the GPU */ #endif #if defined(PETSC_HAVE_VIENNACL) PetscViennaCLFlag valid_GPU_array; /* indicates where the most recently modified vector data is (GPU or CPU) */ void *spptr; /* if we're using ViennaCL, then this is the special pointer to the array on the GPU */ #endif So I guess configure should complain when both are enabled at configure time, otherwise the build fails. Jose

1 0

Re: [petsc-dev] bitbucket giving wrong information
by Patrick Sanan 03 Dec '14

03 Dec '14

Whoops, guess that wasn't the cause then! Actually, looking at bitbucket, it looks like the PR is coming from a different repository (https://bitbucket.org/iaannn/petsc/) This worked on my machine: [petsc (master=)]$ git remote add iaannn [email protected]:iaannn/petsc.git [petsc (master=)]$ git fetch iaannn [petsc (master=)]$ git checkout -b ianwilliamson/fix-matlab-write-bin-complex iaannn/ianwilliamson/fix-matlab-write-bin-complex [petsc (ianwilliamson/fix-matlab-write-bin-complex=)]$ On 04/12/14 12:57 PM, Barry Smith wrote: > $ git checkout -b ianwilliamson/fix-matlab-write-bin-complex origin/ianwilliamson/fix-matlab-write-bin-complex > fatal: Cannot update paths and switch to branch 'ianwilliamson/fix-matlab-write-bin-complex' at the same time. > Did you intend to checkout 'origin/ianwilliamson/fix-matlab-write-bin-complex' which can not be resolved as commit? > > > >> On Dec 3, 2014, at 9:54 PM, Patrick Sanan <patrick.sanan(a)gmail.com> wrote: >> >> Is it possible that you have a second remote with a branch with the same name? In that case, perhaps something like: >> >> git checkout -b ianwilliamson/fix-matlab-write-bin-complex origin/ianwilliamson/fix-matlab-write-bin-complex >> >> On 04/12/14 12:37 PM, Barry Smith wrote: >>> So some guy posts a pull request. I, in trying to follow PETSc-rules, want to check out his branch and then merge that rather than use the merge directly from the bitbucket website. So I go to the branch in the bitbucket website and bitbucket claims I can look at it with >>> >>> git fetch && git checkout ianwilliamson/fix-matlab-write-bin-complex >>> >>> (see attachment). I try that and of course it does not work: >>> >>> $ git fetch && git checkout ianwilliamson/fix-matlab-write-bin-complex >>> remote: Counting objects: 56, done. >>> remote: Compressing objects: 100% (56/56), done. >>> remote: Total 56 (delta 48), reused 0 (delta 0) >>> Unpacking objects: 100% (56/56), done. >>> From bitbucket.org:petsc/petsc >>> f685883..cf8fc5e master -> origin/master >>> 7cadcfe..39ea070 knepley/fix-plex-orient -> origin/knepley/fix-plex-orient >>> 565f831..ce993ff knepley/pylith -> origin/knepley/pylith >>> + 3f11d05...21345e8 mlange/refactor-plex-distribute -> origin/mlange/refactor-plex-distribute (forced update) >>> e9d5acc..916966e next -> origin/next >>> error: pathspec 'ianwilliamson/fix-matlab-write-bin-complex' did not match any file(s) known to git. >>> >>> So is bitbucket lying to me or am I doing something wrong. Should this be a feature request to bitbucket ? Remember I am living in the world where I cannot remember how to do anything and thus expect my tools to tell me how to do something. So why cannot bitbucket tell me the truth of how to get this branch? >>> >>> Barry >>> >>> >>> >>> >>> <Mail Attachment.png>

1 0

bitbucket giving wrong information
by Barry Smith 03 Dec '14

03 Dec '14

So some guy posts a pull request. I, in trying to follow PETSc-rules, want to check out his branch and then merge that rather than use the merge directly from the bitbucket website. So I go to the branch in the bitbucket website and bitbucket claims I can look at it with git fetch && git checkout ianwilliamson/fix-matlab-write-bin-complex (see attachment). I try that and of course it does not work: $ git fetch && git checkout ianwilliamson/fix-matlab-write-bin-complex remote: Counting objects: 56, done. remote: Compressing objects: 100% (56/56), done. remote: Total 56 (delta 48), reused 0 (delta 0) Unpacking objects: 100% (56/56), done. From bitbucket.org:petsc/petsc f685883..cf8fc5e master -> origin/master 7cadcfe..39ea070 knepley/fix-plex-orient -> origin/knepley/fix-plex-orient 565f831..ce993ff knepley/pylith -> origin/knepley/pylith + 3f11d05...21345e8 mlange/refactor-plex-distribute -> origin/mlange/refactor-plex-distribute (forced update) e9d5acc..916966e next -> origin/next error: pathspec 'ianwilliamson/fix-matlab-write-bin-complex' did not match any file(s) known to git. So is bitbucket lying to me or am I doing something wrong. Should this be a feature request to bitbucket ? Remember I am living in the world where I cannot remember how to do anything and thus expect my tools to tell me how to do something. So why cannot bitbucket tell me the truth of how to get this branch? Barry

1 0

Re: [petsc-dev] TS not counting linear iterations properly?
by Barry Smith 03 Dec '14

03 Dec '14

Beats the heck out of me, too hard to run myself. Requires all kinds of data files and other crapola I don't have. > On Dec 3, 2014, at 8:11 PM, Emil Constantinescu <emconsta(a)mcs.anl.gov> wrote: > > On 12/3/14 11:21 AM, Barry Smith wrote: >> >> Jed and Emil, >> >> type: arkimex >> total number of nonlinear solver iterations=1522 >> total number of nonlinear solve failures=158 >> total number of linear solver iterations=4871 >> total number of rejected steps=158 >> >> KSPGMRESOrthog 1601671 1.0 9.8384e+03 1.1 6.96e+12 1.1 0.0e+00 0.0e+00 >> 1.6e+06 8 9 0 0 14 8 9 0 0 14 2777 >> KSPSetUp 6721 1.0 2.0089e-02 1.6 0.00e+00 0.0 0.0e+00 0.0e+00 >> 2.0e+01 0 0 0 0 0 0 0 0 0 0 0 >> KSPSolve 1680 1.0 1.1547e+05 1.0 7.63e+13 1.1 8.6e+07 3.5e+02 >> 1.1e+07 99100100 96100 99100100 96100 2593 >> PCSetUp 5040 1.0 5.0914e+02 1.0 1.01e+09 1.0 8.0e+01 2.8e+04 >> 1.4e+04 0 0 0 0 0 0 0 0 0 0 8 >> PCApply 1584871 1.0 9.0297e+04 1.0 5.83e+13 1.1 7.6e+07 3.8e+02 >> 8.0e+06 77 76 89 93 71 77 76 89 93 71 2535 >> >> Seems to be missing counting some of the linear solves. >> >> What's up with this? > > What are the counts if you run with -ts_dt 1e-12 -ts_max_snes_failures 0 -ts_max_steps 1 ? > > Emil

1 0

Re: [petsc-dev] TS not counting linear iterations properly?
by Emil Constantinescu 03 Dec '14

03 Dec '14

On 12/3/14 11:21 AM, Barry Smith wrote: > > Jed and Emil, > > type: arkimex > total number of nonlinear solver iterations=1522 > total number of nonlinear solve failures=158 > total number of linear solver iterations=4871 > total number of rejected steps=158 > > KSPGMRESOrthog 1601671 1.0 9.8384e+03 1.1 6.96e+12 1.1 0.0e+00 0.0e+00 > 1.6e+06 8 9 0 0 14 8 9 0 0 14 2777 > KSPSetUp 6721 1.0 2.0089e-02 1.6 0.00e+00 0.0 0.0e+00 0.0e+00 > 2.0e+01 0 0 0 0 0 0 0 0 0 0 0 > KSPSolve 1680 1.0 1.1547e+05 1.0 7.63e+13 1.1 8.6e+07 3.5e+02 > 1.1e+07 99100100 96100 99100100 96100 2593 > PCSetUp 5040 1.0 5.0914e+02 1.0 1.01e+09 1.0 8.0e+01 2.8e+04 > 1.4e+04 0 0 0 0 0 0 0 0 0 0 8 > PCApply 1584871 1.0 9.0297e+04 1.0 5.83e+13 1.1 7.6e+07 3.8e+02 > 8.0e+06 77 76 89 93 71 77 76 89 93 71 2535 > > Seems to be missing counting some of the linear solves. > > What's up with this? What are the counts if you run with -ts_dt 1e-12 -ts_max_snes_failures 0 -ts_max_steps 1 ? Emil > >>> >> >> >

1 0

Re: [petsc-dev] [EXTERNAL] [pflotran-dev: 2700] Unintended consequence of CHKERRQ()
by Gautam Bisht 03 Dec '14

03 Dec '14

The user function was putting NaNs in the residual vector. I ended up fixing the user function to avoid putting NaNs in the residual vector. -Gautam. On Wed, Dec 3, 2014 at 10:48 AM, Jed Brown <jed(a)jedbrown.org> wrote: > "Hammond, Glenn E" <gehammo(a)sandia.gov> writes: > > > Can you explain why the NaNs or infs are being produced in the first > place, or at least pinpoint the location and catch it there? > > If the SNESFunction computes NaN because the state is invalid, please > call SNESSetFunctionDomainError and don't update the residual vector (or > don't put NaN in there). SNES will clean up and return with > SNESConvergedReason SNES_DIVERGED_FUNCTION_DOMAIN. > > We don't do this by default because detecting NaN collectively requires > an extra reduction. Perhaps SNES could have a mode that automatically > checks. > > > > > Glenn > > > > From: pflotran-dev(a)googlegroups.com [mailto: > pflotran-dev(a)googlegroups.com] On Behalf Of Gautam Bisht > > Sent: Wednesday, December 03, 2014 8:36 AM > > To: pflotran-dev(a)googlegroups.com > > Subject: [EXTERNAL] [pflotran-dev: 2695] Unintended consequence of > CHKERRQ() > > > > The issue is that when SNES diverges with a "not-a-number or infinite", > the code crashes without giving timestepper_BE a chance to reduce the > time-step. But, if CHKERRQ() is commented out after call to SNESSolve(), > the code is able to reduce dt sufficiently to converge due to "itol_res" > (PFLOTRAN's custom convergence reason == inorm_residual < > solver%newton_inf_res_tol). Attached below is inputdeck (with an outrages > initial dt = 1yr to illustrate the issue). Apart from commenting out > CHKERRQ(), is there an alternate way to let timestepper_BE reduce dt? > > > > Below are relevant information to reproduce the error. > > > >>hg tip > > changeset: 8541:fa2d4295aa14 > > tag: tip > > user: Nathan Collier <nathaniel.collier(a)gmail.com<mailto: > nathaniel.collier(a)gmail.com>> > > date: Wed Dec 03 08:52:00 2014 -0500 > > summary: made output process-model identifier lines to be of uniform > length (80 characters). Yeah it is a little OCD of me. > > > >>./pflotran -pflotranin da.in<http://da.in> > > > > == TH FLOW > ===================================================================== > > 1 2f: 2.55E+00 2x: 1.63E+07 2s: 1.72E+07 ir: 1.49E+00 iu: 3.10E+06 > rsn: 0 > > 2 2f: 1.88E-01 2x: 1.57E+14 2s: 1.57E+14 ir: 5.33E-02 iu: 4.37E+13 > rsn: 0 > > 3 2f: 1.07E-01 2x: 1.60E+14 2s: 4.28E+12 ir: 2.27E-02 iu: 9.90E+11 > rsn: 0 > > 4 2f: 1.24E+01 2x: 1.60E+14 2s: 9.53E+11 ir: 1.16E+01 iu: 7.15E+11 > rsn: 0 > > 5 2f: 1.29E+06 2x: 7.75E+20 2s: 7.75E+20 ir: 1.15E+06 iu: 7.21E+20 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 1[ 0] t= 0.00000E+00 dt= > 5.00000E-01 > > 1 2f: 2.22E+00 2x: 1.38E+07 2s: 1.47E+07 ir: 1.47E+00 iu: 3.06E+06 > rsn: 0 > > 2 2f: 2.94E-01 2x: 1.11E+14 2s: 1.11E+14 ir: 1.44E-01 iu: 4.05E+13 > rsn: 0 > > 3 2f: 6.05E+20 2x: 2.13E+14 2s: 1.79E+14 ir: 4.44E+20 iu: 8.93E+13 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 2[ 0] t= 0.00000E+00 dt= > 2.50000E-01 > > 1 2f: 2.17E+00 2x: 1.12E+07 2s: 1.21E+07 ir: 1.44E+00 iu: 3.01E+06 > rsn: 0 > > 2 2f: 5.55E-01 2x: 8.08E+13 2s: 8.08E+13 ir: 4.07E-01 iu: 3.66E+13 > rsn: 0 > > [0]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > > [0]PETSC ERROR: Floating point exception > > [0]PETSC ERROR: Vec entry at local location 50 is not-a-number or > infinite at end of function: Parameter number 3 > > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html > for trouble shooting. > > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a > darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< > http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:16:17 2014 > > [0]PETSC ERROR: Configure options --download-hdf5=1 > --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate > --download-parmetis=yes --download-metis=yes > --with-cc=/opt/local/bin/mpicc-openmpi-mp > --with-cxx=/opt/local/bin/mpicxx-openmpi-mp > --with-fc=/opt/local/bin/mpif90-openmpi-mp > --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp > PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > > [0]PETSC ERROR: #1 VecValidValues() line 34 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > > [0]PETSC ERROR: #2 SNESComputeFunction() line 2043 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > [0]PETSC ERROR: #3 SNESLineSearchApply_Basic() line 40 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > > [0]PETSC ERROR: #4 SNESLineSearchApply() line 563 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > > [0]PETSC ERROR: #5 SNESSolve_NEWTONLS() line 254 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > > [0]PETSC ERROR: #6 SNESSolve() line 3743 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > > -------------------------------------------------------------------------- > > MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD > > with errorcode 72. > > > > NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. > > You may or may not see output from other processes, depending on > > exactly when Open MPI kills them. > > > -------------------------------------------------------------------------- > > > > > >> hg diff > > diff -r fa2d4295aa14 src/pflotran/timestepper_BE.F90 > > --- a/src/pflotran/timestepper_BE.F90 Wed Dec 03 08:52:00 2014 > -0500 > > +++ b/src/pflotran/timestepper_BE.F90 Wed Dec 03 08:22:35 2014 > -0800 > > @@ -310,7 +310,7 @@ > > call PetscTime(log_start_time, ierr);CHKERRQ(ierr) > > > > call SNESSolve(solver%snes,PETSC_NULL_OBJECT, & > > - process_model%solution_vec,ierr);CHKERRQ(ierr) > > + process_model%solution_vec,ierr);!CHKERRQ(ierr) > > > > call PetscTime(log_end_time, ierr);CHKERRQ(ierr) > > > > > >>./pflotran -pflotranin da.in<http://da.in> > > > > == TH FLOW > ===================================================================== > > 1 2f: 2.55E+00 2x: 1.63E+07 2s: 1.72E+07 ir: 1.49E+00 iu: 3.10E+06 > rsn: 0 > > 2 2f: 1.88E-01 2x: 1.57E+14 2s: 1.57E+14 ir: 5.33E-02 iu: 4.37E+13 > rsn: 0 > > 3 2f: 1.07E-01 2x: 1.60E+14 2s: 4.28E+12 ir: 2.27E-02 iu: 9.90E+11 > rsn: 0 > > 4 2f: 1.24E+01 2x: 1.60E+14 2s: 9.53E+11 ir: 1.16E+01 iu: 7.15E+11 > rsn: 0 > > 5 2f: 1.29E+06 2x: 7.75E+20 2s: 7.75E+20 ir: 1.15E+06 iu: 7.21E+20 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 1[ 0] t= 0.00000E+00 dt= > 5.00000E-01 > > 1 2f: 2.22E+00 2x: 1.38E+07 2s: 1.47E+07 ir: 1.47E+00 iu: 3.06E+06 > rsn: 0 > > 2 2f: 2.94E-01 2x: 1.11E+14 2s: 1.11E+14 ir: 1.44E-01 iu: 4.05E+13 > rsn: 0 > > 3 2f: 6.05E+20 2x: 2.13E+14 2s: 1.79E+14 ir: 4.44E+20 iu: 8.93E+13 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 2[ 0] t= 0.00000E+00 dt= > 2.50000E-01 > > 1 2f: 2.17E+00 2x: 1.12E+07 2s: 1.21E+07 ir: 1.44E+00 iu: 3.01E+06 > rsn: 0 > > 2 2f: 5.55E-01 2x: 8.08E+13 2s: 8.08E+13 ir: 4.07E-01 iu: 3.66E+13 > rsn: 0 > > [0]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > > [0]PETSC ERROR: Floating point exception > > [0]PETSC ERROR: Vec entry at local location 50 is not-a-number or > infinite at end of function: Parameter number 3 > > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html > for trouble shooting. > > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a > darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< > http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 > > [0]PETSC ERROR: Configure options --download-hdf5=1 > --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate > --download-parmetis=yes --download-metis=yes > --with-cc=/opt/local/bin/mpicc-openmpi-mp > --with-cxx=/opt/local/bin/mpicxx-openmpi-mp > --with-fc=/opt/local/bin/mpif90-openmpi-mp > --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp > PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > > [0]PETSC ERROR: #1 VecValidValues() line 34 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > > [0]PETSC ERROR: #2 SNESComputeFunction() line 2043 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > [0]PETSC ERROR: #3 SNESLineSearchApply_Basic() line 40 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > > [0]PETSC ERROR: #4 SNESLineSearchApply() line 563 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > > [0]PETSC ERROR: #5 SNESSolve_NEWTONLS() line 254 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > > [0]PETSC ERROR: #6 SNESSolve() line 3743 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > -> Cut time step: snes= 0 icut= 3[ 0] t= 0.00000E+00 dt= > 1.25000E-01 > > 1 2f: 2.12E+00 2x: 9.08E+06 2s: 9.85E+06 ir: 1.41E+00 iu: 2.94E+06 > rsn: 0 > > 2 2f: 7.89E-01 2x: 5.95E+13 2s: 5.95E+13 ir: 6.28E-01 iu: 3.20E+13 > rsn: 0 > > [0]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > > [0]PETSC ERROR: Floating point exception > > [0]PETSC ERROR: Vec entry at local location 34 is not-a-number or > infinite at end of function: Parameter number 3 > > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html > for trouble shooting. > > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a > darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< > http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 > > [0]PETSC ERROR: Configure options --download-hdf5=1 > --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate > --download-parmetis=yes --download-metis=yes > --with-cc=/opt/local/bin/mpicc-openmpi-mp > --with-cxx=/opt/local/bin/mpicxx-openmpi-mp > --with-fc=/opt/local/bin/mpif90-openmpi-mp > --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp > PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > > [0]PETSC ERROR: #7 VecValidValues() line 34 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > > [0]PETSC ERROR: #8 SNESComputeFunction() line 2043 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > [0]PETSC ERROR: #9 SNESLineSearchApply_Basic() line 40 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > > [0]PETSC ERROR: #10 SNESLineSearchApply() line 563 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > > [0]PETSC ERROR: #11 SNESSolve_NEWTONLS() line 254 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > > [0]PETSC ERROR: #12 SNESSolve() line 3743 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > -> Cut time step: snes= 0 icut= 4[ 0] t= 0.00000E+00 dt= > 6.25000E-02 > > 1 2f: 2.04E+00 2x: 7.38E+06 2s: 8.06E+06 ir: 1.36E+00 iu: 2.86E+06 > rsn: 0 > > 2 2f: 5.60E+04 2x: 4.28E+13 2s: 4.28E+13 ir: 4.47E+04 iu: 2.66E+13 > rsn: 0 > > 3 2f: 1.16E+18 2x: 4.91E+13 2s: 2.30E+13 ir: 9.47E+17 iu: 1.69E+13 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 5[ 0] t= 0.00000E+00 dt= > 3.12500E-02 > > 1 2f: 1.94E+00 2x: 6.00E+06 2s: 6.57E+06 ir: 1.30E+00 iu: 2.74E+06 > rsn: 0 > > 2 2f: 6.76E+04 2x: 2.91E+13 2s: 2.91E+13 ir: 5.18E+04 iu: 2.06E+13 > rsn: 0 > > [0]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > > [0]PETSC ERROR: Floating point exception > > [0]PETSC ERROR: Vec entry at local location 14 is not-a-number or > infinite at end of function: Parameter number 3 > > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html > for trouble shooting. > > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a > darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< > http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 > > [0]PETSC ERROR: Configure options --download-hdf5=1 > --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate > --download-parmetis=yes --download-metis=yes > --with-cc=/opt/local/bin/mpicc-openmpi-mp > --with-cxx=/opt/local/bin/mpicxx-openmpi-mp > --with-fc=/opt/local/bin/mpif90-openmpi-mp > --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp > PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > > [0]PETSC ERROR: #13 VecValidValues() line 34 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > > [0]PETSC ERROR: #14 SNESComputeFunction() line 2043 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > [0]PETSC ERROR: #15 SNESLineSearchApply_Basic() line 40 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > > [0]PETSC ERROR: #16 SNESLineSearchApply() line 563 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > > [0]PETSC ERROR: #17 SNESSolve_NEWTONLS() line 254 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > > [0]PETSC ERROR: #18 SNESSolve() line 3743 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > -> Cut time step: snes= 0 icut= 6[ 0] t= 0.00000E+00 dt= > 1.56250E-02 > > 1 2f: 1.81E+00 2x: 4.84E+06 2s: 5.31E+06 ir: 1.21E+00 iu: 2.58E+06 > rsn: 0 > > 2 2f: 1.56E+04 2x: 1.81E+13 2s: 1.81E+13 ir: 1.22E+04 iu: 1.43E+13 > rsn: 0 > > 3 2f: 9.46E+15 2x: 1.99E+13 2s: 7.41E+12 ir: 7.34E+15 iu: 5.84E+12 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 7[ 0] t= 0.00000E+00 dt= > 7.81250E-03 > > 1 2f: 1.64E+00 2x: 3.85E+06 2s: 4.21E+06 ir: 1.11E+00 iu: 2.37E+06 > rsn: 0 > > 2 2f: 1.37E+04 2x: 9.82E+12 2s: 9.82E+12 ir: 1.02E+04 iu: 8.52E+12 > rsn: 0 > > 3 2f: 5.35E+14 2x: 1.86E+13 2s: 1.65E+13 ir: 3.78E+14 iu: 1.62E+13 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 8[ 0] t= 0.00000E+00 dt= > 3.90625E-03 > > 1 2f: 1.42E+00 2x: 3.01E+06 2s: 3.25E+06 ir: 9.67E-01 iu: 2.09E+06 > rsn: 0 > > 2 2f: 6.03E+03 2x: 4.39E+12 2s: 4.39E+12 ir: 4.38E+03 iu: 4.08E+12 > rsn: 0 > > [0]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > > [0]PETSC ERROR: Floating point exception > > [0]PETSC ERROR: Vec entry at local location 30 is not-a-number or > infinite at end of function: Parameter number 3 > > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html > for trouble shooting. > > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a > darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< > http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 > > [0]PETSC ERROR: Configure options --download-hdf5=1 > --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate > --download-parmetis=yes --download-metis=yes > --with-cc=/opt/local/bin/mpicc-openmpi-mp > --with-cxx=/opt/local/bin/mpicxx-openmpi-mp > --with-fc=/opt/local/bin/mpif90-openmpi-mp > --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp > PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > > [0]PETSC ERROR: #19 VecValidValues() line 34 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > > [0]PETSC ERROR: #20 SNESComputeFunction() line 2043 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > [0]PETSC ERROR: #21 SNESLineSearchApply_Basic() line 40 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > > [0]PETSC ERROR: #22 SNESLineSearchApply() line 563 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > > [0]PETSC ERROR: #23 SNESSolve_NEWTONLS() line 254 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > > [0]PETSC ERROR: #24 SNESSolve() line 3743 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > -> Cut time step: snes= 0 icut= 9[ 0] t= 0.00000E+00 dt= > 1.95312E-03 > > 1 2f: 1.16E+00 2x: 2.30E+06 2s: 2.42E+06 ir: 8.00E-01 iu: 1.76E+06 > rsn: 0 > > 2 2f: 4.32E+03 2x: 1.47E+12 2s: 1.47E+12 ir: 3.07E+03 iu: 1.42E+12 > rsn: 0 > > [0]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > > [0]PETSC ERROR: Floating point exception > > [0]PETSC ERROR: Vec entry at local location 4 is not-a-number or > infinite at end of function: Parameter number 3 > > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html > for trouble shooting. > > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a > darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov< > http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 > > [0]PETSC ERROR: Configure options --download-hdf5=1 > --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate > --download-parmetis=yes --download-metis=yes > --with-cc=/opt/local/bin/mpicc-openmpi-mp > --with-cxx=/opt/local/bin/mpicxx-openmpi-mp > --with-fc=/opt/local/bin/mpif90-openmpi-mp > --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp > PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > > [0]PETSC ERROR: #25 VecValidValues() line 34 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > > [0]PETSC ERROR: #26 SNESComputeFunction() line 2043 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > [0]PETSC ERROR: #27 SNESLineSearchApply_Basic() line 40 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > > [0]PETSC ERROR: #28 SNESLineSearchApply() line 563 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > > [0]PETSC ERROR: #29 SNESSolve_NEWTONLS() line 254 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > > [0]PETSC ERROR: #30 SNESSolve() line 3743 in > /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > > -> Cut time step: snes= 0 icut= 10[ 0] t= 0.00000E+00 dt= > 9.76562E-04 > > 1 2f: 8.79E-01 2x: 1.73E+06 2s: 1.70E+06 ir: 6.11E-01 iu: 1.37E+06 > rsn: 0 > > 2 2f: 1.34E+03 2x: 3.17E+11 2s: 3.17E+11 ir: 9.30E+02 iu: 3.14E+11 > rsn: 0 > > 3 2f: 3.99E+13 2x: 6.13E+11 2s: 3.04E+11 ir: 3.26E+13 iu: 2.98E+11 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 11[ 0] t= 0.00000E+00 dt= > 4.88281E-04 > > 1 2f: 5.96E-01 2x: 1.33E+06 2s: 1.10E+06 ir: 4.19E-01 iu: 9.69E+05 > rsn: 0 > > 2 2f: 4.60E+04 2x: 3.54E+10 2s: 3.54E+10 ir: 3.25E+04 iu: 3.54E+10 > rsn: 0 > > 3 2f: 9.64E+15 2x: 3.91E+11 2s: 3.55E+11 ir: 7.87E+15 iu: 3.55E+11 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 12[ 0] t= 0.00000E+00 dt= > 2.44141E-04 > > 1 2f: 3.54E-01 2x: 1.11E+06 2s: 6.36E+05 ir: 2.50E-01 iu: 5.95E+05 > rsn: 0 > > 2 2f: 5.42E+04 2x: 1.33E+09 2s: 1.33E+09 ir: 3.93E+04 iu: 1.33E+09 > rsn: 0 > > 3 2f: 3.51E+06 2x: 1.75E+10 2s: 1.75E+10 ir: 2.87E+06 iu: 1.75E+10 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 13[ 0] t= 0.00000E+00 dt= > 1.22070E-04 > > 1 2f: 5.74E-02 2x: 1.04E+06 2s: 3.17E+05 ir: 4.02E-02 iu: 3.07E+05 > rsn: 0 > > 2 2f: 3.60E-02 2x: 6.66E+09 2s: 6.66E+09 ir: 2.88E-02 iu: 6.66E+09 > rsn: 0 > > 3 2f: 5.16E+02 2x: 2.91E+08 2s: 6.95E+09 ir: 3.31E+02 iu: 6.95E+09 > rsn: 0 > > 4 2f: 7.20E+04 2x: 7.65E+15 2s: 7.65E+15 ir: 6.73E+04 iu: 7.65E+15 > rsn: 0 > > 5 2f: 4.01E+12 2x: 3.15E+16 2s: 3.92E+16 ir: 3.27E+12 iu: 3.92E+16 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 14[ 0] t= 0.00000E+00 dt= > 6.10352E-05 > > 1 2f: 8.13E-02 2x: 1.04E+06 2s: 1.33E+05 ir: 8.02E-02 iu: 1.31E+05 > rsn: 0 > > 2 2f: 1.94E+02 2x: 1.14E+08 2s: 1.14E+08 ir: 1.33E+02 iu: 1.14E+08 > rsn: 0 > > 3 2f: 2.16E-01 2x: 4.22E+08 2s: 5.36E+08 ir: 1.79E-01 iu: 5.36E+08 > rsn: 0 > > 4 2f: 3.21E+28 2x: 1.07E+19 2s: 1.07E+19 ir: 2.62E+28 iu: 1.07E+19 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 15[ 0] t= 0.00000E+00 dt= > 3.05176E-05 > > 1 2f: 1.13E-05 2x: 1.05E+06 2s: 4.69E+04 ir: 1.13E-05 iu: 4.67E+04 > rsn: 0 > > 2 2f: 1.21E-10 2x: 1.05E+06 2s: 6.87E+02 ir: 1.05E-10 iu: 6.87E+02 > rsn: itol_res > > > > Step 1 Time= 3.05176E-05 Dt= 3.05176E-05 [y] snes_conv_reason: > 10 > > newton = 47 [ 47] linear = 52 [ 52] cuts = 15 [ 15] > > --> SNES Linear/Non-Linear Iterations = 2 / 2 > > --> SNES Residual: 1.209152E-10 1.209152E-12 1.045172E-10 > > --> max chng: dpmx= 4.5981E+04 dtmpmx= 1.6500E+00 > > > > > > == TH FLOW > ===================================================================== > > 1 2f: 1.62E-01 2x: 1.04E+06 2s: 7.34E+04 ir: 1.60E-01 iu: 7.27E+04 > rsn: 0 > > 2 2f: 7.18E+01 2x: 6.44E+07 2s: 6.44E+07 ir: 4.68E+01 iu: 6.44E+07 > rsn: 0 > > 3 2f: 2.52E-01 2x: 5.48E+07 2s: 1.19E+08 ir: 2.02E-01 iu: 1.19E+08 > rsn: 0 > > 4 2f: 1.29E+09 2x: 1.44E+18 2s: 1.44E+18 ir: 1.29E+09 iu: 1.44E+18 > rsn: max_norm > > -> Cut time step: snes=-10 icut= 1[ 15] t= 3.05176E-05 dt= > 1.52588E-05 > > 1 2f: 3.23E-01 2x: 1.04E+06 2s: 3.15E+04 ir: 3.19E-01 iu: 3.13E+04 > rsn: 0 > > 2 2f: 1.17E-03 2x: 3.08E+06 2s: 2.83E+06 ir: 1.09E-03 iu: 2.83E+06 > rsn: 0 > > 3 2f: 3.23E-01 2x: 1.04E+06 2s: 2.83E+06 ir: 3.19E-01 iu: 2.83E+06 > rsn: 0 > > 4 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 5 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 6 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 7 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 8 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 9 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 10 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 11 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 12 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 13 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 14 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 15 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 16 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 17 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 18 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 19 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 20 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 21 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 22 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 23 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 24 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 25 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 26 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 27 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 28 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 29 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 30 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 31 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 32 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 33 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 34 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 35 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 36 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 37 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 38 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 39 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 40 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 41 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 42 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 43 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 44 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 45 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 46 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 47 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 48 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 49 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 50 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 51 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 52 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 53 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 54 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 55 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 56 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 57 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 58 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 59 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 60 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 61 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 62 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 63 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 64 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 65 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 66 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 67 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 68 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 69 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 70 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 71 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 72 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 73 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 74 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 75 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 76 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 77 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 78 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 79 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 80 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 81 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 82 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 83 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 84 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 85 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 86 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 87 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 88 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 89 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 90 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 91 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 92 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 93 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 94 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 95 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 96 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 97 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 98 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > 99 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 > rsn: 0 > > 100 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 > rsn: 0 > > -> Cut time step: snes= -5 icut= 2[ 15] t= 3.05176E-05 dt= > 7.62939E-06 > > 1 2f: 6.23E-01 2x: 1.04E+06 2s: 1.39E+04 ir: 6.15E-01 iu: 1.38E+04 > rsn: 0 > > 2 2f: 6.18E-04 2x: 1.06E+06 2s: 1.15E+05 ir: 5.25E-04 iu: 1.15E+05 > rsn: 0 > > 3 2f: 6.23E-01 2x: 1.04E+06 2s: 1.15E+05 ir: 6.15E-01 iu: 1.15E+05 > rsn: 0 > > 4 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 5 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 6 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 7 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 8 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 9 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 10 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 11 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 12 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 13 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 14 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 15 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 16 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 17 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 18 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 19 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 20 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 21 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 22 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 23 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 24 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 25 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 26 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 27 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 28 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 29 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 30 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 31 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 32 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 33 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 34 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 35 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 36 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 37 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 38 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 39 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 40 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 41 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 42 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 43 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 44 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 45 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 46 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 47 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 48 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 49 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 50 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 51 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 52 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 53 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 54 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 55 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 56 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 57 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 58 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 59 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 60 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 61 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 62 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 63 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 64 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 65 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 66 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 67 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 68 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 69 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 70 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 71 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 72 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 73 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 74 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 75 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 76 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 77 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 78 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 79 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 80 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 81 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 82 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 83 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 84 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 85 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 86 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 87 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 88 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 89 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 90 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 91 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 92 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 93 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 94 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 95 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 96 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 97 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 98 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > 99 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 > rsn: 0 > > 100 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 > rsn: 0 > > -> Cut time step: snes= -5 icut= 3[ 15] t= 3.05176E-05 dt= > 3.81470E-06 > > 1 2f: 8.15E-01 2x: 1.05E+06 2s: 6.36E+03 ir: 8.05E-01 iu: 6.33E+03 > rsn: 0 > > 2 2f: 1.34E-03 2x: 1.05E+06 2s: 4.20E+03 ir: 1.31E-03 iu: 4.19E+03 > rsn: 0 > > 3 2f: 4.65E-04 2x: 1.05E+06 2s: 2.99E+02 ir: 4.59E-04 iu: 2.50E+02 > rsn: 0 > > 4 2f: 1.43E-04 2x: 1.05E+06 2s: 1.18E+02 ir: 1.41E-04 iu: 1.18E+02 > rsn: 0 > > 5 2f: 4.35E-05 2x: 1.05E+06 2s: 8.37E+01 ir: 4.30E-05 iu: 8.37E+01 > rsn: 0 > > 6 2f: 1.27E-05 2x: 1.05E+06 2s: 5.80E+01 ir: 1.25E-05 iu: 5.80E+01 > rsn: 0 > > 7 2f: 3.26E-06 2x: 1.05E+06 2s: 3.69E+01 ir: 3.22E-06 iu: 3.69E+01 > rsn: 0 > > 8 2f: 5.57E-07 2x: 1.05E+06 2s: 1.83E+01 ir: 5.50E-07 iu: 1.83E+01 > rsn: 0 > > 9 2f: 2.98E-08 2x: 1.05E+06 2s: 4.65E+00 ir: 2.95E-08 iu: 4.65E+00 > rsn: 0 > > 10 2f: 9.92E-11 2x: 1.05E+06 2s: 2.79E-01 ir: 9.78E-11 iu: 2.79E-01 > rsn: itol_res > > > > Step 2 Time= 3.43323E-05 Dt= 3.81470E-06 [y] snes_conv_reason: > 10 > > newton = 214 [ 261] linear = 214 [ 266] cuts = 3 [ 18] > > --> SNES Linear/Non-Linear Iterations = 10 / 10 > > --> SNES Residual: 9.921423E-11 9.921423E-13 9.778422E-11 > > --> max chng: dpmx= 1.6597E+03 dtmpmx= 1.9780E-01 > > > > > > -- > > You received this message because you are subscribed to the Google > Groups "pflotran-dev" group. > > To view this discussion on the web visit > https://groups.google.com/d/msgid/pflotran-dev/CAPz1TnfKQJov5Uxvghp_31uKpxz… > < > https://groups.google.com/d/msgid/pflotran-dev/CAPz1TnfKQJov5Uxvghp_31uKpxz… > >. > > > > -- > > You received this message because you are subscribed to the Google > Groups "pflotran-dev" group. > > To view this discussion on the web visit > https://groups.google.com/d/msgid/pflotran-dev/92b8b510fceb47309cc0b8f7ada4… > . > > -- > You received this message because you are subscribed to the Google Groups > "pflotran-dev" group. > To view this discussion on the web visit > https://groups.google.com/d/msgid/pflotran-dev/87mw748t8f.fsf%40jedbrown.org > . >

1 0

Re: [petsc-dev] [EXTERNAL] [pflotran-dev: 2696] Unintended consequence of CHKERRQ()
by Jed Brown 03 Dec '14

03 Dec '14

"Hammond, Glenn E" <gehammo(a)sandia.gov> writes: > Can you explain why the NaNs or infs are being produced in the first place, or at least pinpoint the location and catch it there? If the SNESFunction computes NaN because the state is invalid, please call SNESSetFunctionDomainError and don't update the residual vector (or don't put NaN in there). SNES will clean up and return with SNESConvergedReason SNES_DIVERGED_FUNCTION_DOMAIN. We don't do this by default because detecting NaN collectively requires an extra reduction. Perhaps SNES could have a mode that automatically checks. > > Glenn > > From: pflotran-dev(a)googlegroups.com [mailto:[email protected]] On Behalf Of Gautam Bisht > Sent: Wednesday, December 03, 2014 8:36 AM > To: pflotran-dev(a)googlegroups.com > Subject: [EXTERNAL] [pflotran-dev: 2695] Unintended consequence of CHKERRQ() > > The issue is that when SNES diverges with a "not-a-number or infinite", the code crashes without giving timestepper_BE a chance to reduce the time-step. But, if CHKERRQ() is commented out after call to SNESSolve(), the code is able to reduce dt sufficiently to converge due to "itol_res" (PFLOTRAN's custom convergence reason == inorm_residual < solver%newton_inf_res_tol). Attached below is inputdeck (with an outrages initial dt = 1yr to illustrate the issue). Apart from commenting out CHKERRQ(), is there an alternate way to let timestepper_BE reduce dt? > > Below are relevant information to reproduce the error. > >>hg tip > changeset: 8541:fa2d4295aa14 > tag: tip > user: Nathan Collier <nathaniel.collier(a)gmail.com<mailto:[email protected]>> > date: Wed Dec 03 08:52:00 2014 -0500 > summary: made output process-model identifier lines to be of uniform length (80 characters). Yeah it is a little OCD of me. > >>./pflotran -pflotranin da.in<http://da.in> > > == TH FLOW ===================================================================== > 1 2f: 2.55E+00 2x: 1.63E+07 2s: 1.72E+07 ir: 1.49E+00 iu: 3.10E+06 rsn: 0 > 2 2f: 1.88E-01 2x: 1.57E+14 2s: 1.57E+14 ir: 5.33E-02 iu: 4.37E+13 rsn: 0 > 3 2f: 1.07E-01 2x: 1.60E+14 2s: 4.28E+12 ir: 2.27E-02 iu: 9.90E+11 rsn: 0 > 4 2f: 1.24E+01 2x: 1.60E+14 2s: 9.53E+11 ir: 1.16E+01 iu: 7.15E+11 rsn: 0 > 5 2f: 1.29E+06 2x: 7.75E+20 2s: 7.75E+20 ir: 1.15E+06 iu: 7.21E+20 rsn: max_norm > -> Cut time step: snes=-10 icut= 1[ 0] t= 0.00000E+00 dt= 5.00000E-01 > 1 2f: 2.22E+00 2x: 1.38E+07 2s: 1.47E+07 ir: 1.47E+00 iu: 3.06E+06 rsn: 0 > 2 2f: 2.94E-01 2x: 1.11E+14 2s: 1.11E+14 ir: 1.44E-01 iu: 4.05E+13 rsn: 0 > 3 2f: 6.05E+20 2x: 2.13E+14 2s: 1.79E+14 ir: 4.44E+20 iu: 8.93E+13 rsn: max_norm > -> Cut time step: snes=-10 icut= 2[ 0] t= 0.00000E+00 dt= 2.50000E-01 > 1 2f: 2.17E+00 2x: 1.12E+07 2s: 1.21E+07 ir: 1.44E+00 iu: 3.01E+06 rsn: 0 > 2 2f: 5.55E-01 2x: 8.08E+13 2s: 8.08E+13 ir: 4.07E-01 iu: 3.66E+13 rsn: 0 > [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- > [0]PETSC ERROR: Floating point exception > [0]PETSC ERROR: Vec entry at local location 50 is not-a-number or infinite at end of function: Parameter number 3 > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov<http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:16:17 2014 > [0]PETSC ERROR: Configure options --download-hdf5=1 --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate --download-parmetis=yes --download-metis=yes --with-cc=/opt/local/bin/mpicc-openmpi-mp --with-cxx=/opt/local/bin/mpicxx-openmpi-mp --with-fc=/opt/local/bin/mpif90-openmpi-mp --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > [0]PETSC ERROR: #1 VecValidValues() line 34 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > [0]PETSC ERROR: #2 SNESComputeFunction() line 2043 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > [0]PETSC ERROR: #3 SNESLineSearchApply_Basic() line 40 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > [0]PETSC ERROR: #4 SNESLineSearchApply() line 563 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > [0]PETSC ERROR: #5 SNESSolve_NEWTONLS() line 254 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > [0]PETSC ERROR: #6 SNESSolve() line 3743 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > -------------------------------------------------------------------------- > MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD > with errorcode 72. > > NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. > You may or may not see output from other processes, depending on > exactly when Open MPI kills them. > -------------------------------------------------------------------------- > > >> hg diff > diff -r fa2d4295aa14 src/pflotran/timestepper_BE.F90 > --- a/src/pflotran/timestepper_BE.F90 Wed Dec 03 08:52:00 2014 -0500 > +++ b/src/pflotran/timestepper_BE.F90 Wed Dec 03 08:22:35 2014 -0800 > @@ -310,7 +310,7 @@ > call PetscTime(log_start_time, ierr);CHKERRQ(ierr) > > call SNESSolve(solver%snes,PETSC_NULL_OBJECT, & > - process_model%solution_vec,ierr);CHKERRQ(ierr) > + process_model%solution_vec,ierr);!CHKERRQ(ierr) > > call PetscTime(log_end_time, ierr);CHKERRQ(ierr) > > >>./pflotran -pflotranin da.in<http://da.in> > > == TH FLOW ===================================================================== > 1 2f: 2.55E+00 2x: 1.63E+07 2s: 1.72E+07 ir: 1.49E+00 iu: 3.10E+06 rsn: 0 > 2 2f: 1.88E-01 2x: 1.57E+14 2s: 1.57E+14 ir: 5.33E-02 iu: 4.37E+13 rsn: 0 > 3 2f: 1.07E-01 2x: 1.60E+14 2s: 4.28E+12 ir: 2.27E-02 iu: 9.90E+11 rsn: 0 > 4 2f: 1.24E+01 2x: 1.60E+14 2s: 9.53E+11 ir: 1.16E+01 iu: 7.15E+11 rsn: 0 > 5 2f: 1.29E+06 2x: 7.75E+20 2s: 7.75E+20 ir: 1.15E+06 iu: 7.21E+20 rsn: max_norm > -> Cut time step: snes=-10 icut= 1[ 0] t= 0.00000E+00 dt= 5.00000E-01 > 1 2f: 2.22E+00 2x: 1.38E+07 2s: 1.47E+07 ir: 1.47E+00 iu: 3.06E+06 rsn: 0 > 2 2f: 2.94E-01 2x: 1.11E+14 2s: 1.11E+14 ir: 1.44E-01 iu: 4.05E+13 rsn: 0 > 3 2f: 6.05E+20 2x: 2.13E+14 2s: 1.79E+14 ir: 4.44E+20 iu: 8.93E+13 rsn: max_norm > -> Cut time step: snes=-10 icut= 2[ 0] t= 0.00000E+00 dt= 2.50000E-01 > 1 2f: 2.17E+00 2x: 1.12E+07 2s: 1.21E+07 ir: 1.44E+00 iu: 3.01E+06 rsn: 0 > 2 2f: 5.55E-01 2x: 8.08E+13 2s: 8.08E+13 ir: 4.07E-01 iu: 3.66E+13 rsn: 0 > [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- > [0]PETSC ERROR: Floating point exception > [0]PETSC ERROR: Vec entry at local location 50 is not-a-number or infinite at end of function: Parameter number 3 > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov<http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 > [0]PETSC ERROR: Configure options --download-hdf5=1 --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate --download-parmetis=yes --download-metis=yes --with-cc=/opt/local/bin/mpicc-openmpi-mp --with-cxx=/opt/local/bin/mpicxx-openmpi-mp --with-fc=/opt/local/bin/mpif90-openmpi-mp --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > [0]PETSC ERROR: #1 VecValidValues() line 34 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > [0]PETSC ERROR: #2 SNESComputeFunction() line 2043 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > [0]PETSC ERROR: #3 SNESLineSearchApply_Basic() line 40 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > [0]PETSC ERROR: #4 SNESLineSearchApply() line 563 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > [0]PETSC ERROR: #5 SNESSolve_NEWTONLS() line 254 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > [0]PETSC ERROR: #6 SNESSolve() line 3743 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > -> Cut time step: snes= 0 icut= 3[ 0] t= 0.00000E+00 dt= 1.25000E-01 > 1 2f: 2.12E+00 2x: 9.08E+06 2s: 9.85E+06 ir: 1.41E+00 iu: 2.94E+06 rsn: 0 > 2 2f: 7.89E-01 2x: 5.95E+13 2s: 5.95E+13 ir: 6.28E-01 iu: 3.20E+13 rsn: 0 > [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- > [0]PETSC ERROR: Floating point exception > [0]PETSC ERROR: Vec entry at local location 34 is not-a-number or infinite at end of function: Parameter number 3 > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov<http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 > [0]PETSC ERROR: Configure options --download-hdf5=1 --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate --download-parmetis=yes --download-metis=yes --with-cc=/opt/local/bin/mpicc-openmpi-mp --with-cxx=/opt/local/bin/mpicxx-openmpi-mp --with-fc=/opt/local/bin/mpif90-openmpi-mp --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > [0]PETSC ERROR: #7 VecValidValues() line 34 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > [0]PETSC ERROR: #8 SNESComputeFunction() line 2043 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > [0]PETSC ERROR: #9 SNESLineSearchApply_Basic() line 40 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > [0]PETSC ERROR: #10 SNESLineSearchApply() line 563 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > [0]PETSC ERROR: #11 SNESSolve_NEWTONLS() line 254 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > [0]PETSC ERROR: #12 SNESSolve() line 3743 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > -> Cut time step: snes= 0 icut= 4[ 0] t= 0.00000E+00 dt= 6.25000E-02 > 1 2f: 2.04E+00 2x: 7.38E+06 2s: 8.06E+06 ir: 1.36E+00 iu: 2.86E+06 rsn: 0 > 2 2f: 5.60E+04 2x: 4.28E+13 2s: 4.28E+13 ir: 4.47E+04 iu: 2.66E+13 rsn: 0 > 3 2f: 1.16E+18 2x: 4.91E+13 2s: 2.30E+13 ir: 9.47E+17 iu: 1.69E+13 rsn: max_norm > -> Cut time step: snes=-10 icut= 5[ 0] t= 0.00000E+00 dt= 3.12500E-02 > 1 2f: 1.94E+00 2x: 6.00E+06 2s: 6.57E+06 ir: 1.30E+00 iu: 2.74E+06 rsn: 0 > 2 2f: 6.76E+04 2x: 2.91E+13 2s: 2.91E+13 ir: 5.18E+04 iu: 2.06E+13 rsn: 0 > [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- > [0]PETSC ERROR: Floating point exception > [0]PETSC ERROR: Vec entry at local location 14 is not-a-number or infinite at end of function: Parameter number 3 > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov<http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 > [0]PETSC ERROR: Configure options --download-hdf5=1 --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate --download-parmetis=yes --download-metis=yes --with-cc=/opt/local/bin/mpicc-openmpi-mp --with-cxx=/opt/local/bin/mpicxx-openmpi-mp --with-fc=/opt/local/bin/mpif90-openmpi-mp --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > [0]PETSC ERROR: #13 VecValidValues() line 34 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > [0]PETSC ERROR: #14 SNESComputeFunction() line 2043 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > [0]PETSC ERROR: #15 SNESLineSearchApply_Basic() line 40 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > [0]PETSC ERROR: #16 SNESLineSearchApply() line 563 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > [0]PETSC ERROR: #17 SNESSolve_NEWTONLS() line 254 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > [0]PETSC ERROR: #18 SNESSolve() line 3743 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > -> Cut time step: snes= 0 icut= 6[ 0] t= 0.00000E+00 dt= 1.56250E-02 > 1 2f: 1.81E+00 2x: 4.84E+06 2s: 5.31E+06 ir: 1.21E+00 iu: 2.58E+06 rsn: 0 > 2 2f: 1.56E+04 2x: 1.81E+13 2s: 1.81E+13 ir: 1.22E+04 iu: 1.43E+13 rsn: 0 > 3 2f: 9.46E+15 2x: 1.99E+13 2s: 7.41E+12 ir: 7.34E+15 iu: 5.84E+12 rsn: max_norm > -> Cut time step: snes=-10 icut= 7[ 0] t= 0.00000E+00 dt= 7.81250E-03 > 1 2f: 1.64E+00 2x: 3.85E+06 2s: 4.21E+06 ir: 1.11E+00 iu: 2.37E+06 rsn: 0 > 2 2f: 1.37E+04 2x: 9.82E+12 2s: 9.82E+12 ir: 1.02E+04 iu: 8.52E+12 rsn: 0 > 3 2f: 5.35E+14 2x: 1.86E+13 2s: 1.65E+13 ir: 3.78E+14 iu: 1.62E+13 rsn: max_norm > -> Cut time step: snes=-10 icut= 8[ 0] t= 0.00000E+00 dt= 3.90625E-03 > 1 2f: 1.42E+00 2x: 3.01E+06 2s: 3.25E+06 ir: 9.67E-01 iu: 2.09E+06 rsn: 0 > 2 2f: 6.03E+03 2x: 4.39E+12 2s: 4.39E+12 ir: 4.38E+03 iu: 4.08E+12 rsn: 0 > [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- > [0]PETSC ERROR: Floating point exception > [0]PETSC ERROR: Vec entry at local location 30 is not-a-number or infinite at end of function: Parameter number 3 > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov<http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 > [0]PETSC ERROR: Configure options --download-hdf5=1 --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate --download-parmetis=yes --download-metis=yes --with-cc=/opt/local/bin/mpicc-openmpi-mp --with-cxx=/opt/local/bin/mpicxx-openmpi-mp --with-fc=/opt/local/bin/mpif90-openmpi-mp --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > [0]PETSC ERROR: #19 VecValidValues() line 34 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > [0]PETSC ERROR: #20 SNESComputeFunction() line 2043 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > [0]PETSC ERROR: #21 SNESLineSearchApply_Basic() line 40 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > [0]PETSC ERROR: #22 SNESLineSearchApply() line 563 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > [0]PETSC ERROR: #23 SNESSolve_NEWTONLS() line 254 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > [0]PETSC ERROR: #24 SNESSolve() line 3743 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > -> Cut time step: snes= 0 icut= 9[ 0] t= 0.00000E+00 dt= 1.95312E-03 > 1 2f: 1.16E+00 2x: 2.30E+06 2s: 2.42E+06 ir: 8.00E-01 iu: 1.76E+06 rsn: 0 > 2 2f: 4.32E+03 2x: 1.47E+12 2s: 1.47E+12 ir: 3.07E+03 iu: 1.42E+12 rsn: 0 > [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- > [0]PETSC ERROR: Floating point exception > [0]PETSC ERROR: Vec entry at local location 4 is not-a-number or infinite at end of function: Parameter number 3 > [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. > [0]PETSC ERROR: Petsc Release Version 3.5.1, unknown > [0]PETSC ERROR: ../../../src/pflotran/pflotran on a darwin-gnu-fort-debug named gautams-mac-pro.dhcp.lbl.gov<http://gautams-mac-pro.dhcp.lbl.gov> by gbisht Wed Dec 3 08:17:35 2014 > [0]PETSC ERROR: Configure options --download-hdf5=1 --with-blas-lapack-lib=/System/Library/Frameworks/Accelerate.framework/Versions/Current/Accelerate --download-parmetis=yes --download-metis=yes --with-cc=/opt/local/bin/mpicc-openmpi-mp --with-cxx=/opt/local/bin/mpicxx-openmpi-mp --with-fc=/opt/local/bin/mpif90-openmpi-mp --with-mpiexec=/opt/local/bin/mpiexec-openmpi-mp PETSC_ARCH=darwin-gnu-fort-debug --download-sowing=1 > [0]PETSC ERROR: #25 VecValidValues() line 34 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/vec/vec/interface/rvector.c > [0]PETSC ERROR: #26 SNESComputeFunction() line 2043 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > [0]PETSC ERROR: #27 SNESLineSearchApply_Basic() line 40 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/impls/basic/linesearchbasic.c > [0]PETSC ERROR: #28 SNESLineSearchApply() line 563 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/linesearch/interface/linesearch.c > [0]PETSC ERROR: #29 SNESSolve_NEWTONLS() line 254 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/impls/ls/ls.c > [0]PETSC ERROR: #30 SNESSolve() line 3743 in /Users/gbisht/Research/Models/PETSc/petsc_821a792/src/snes/interface/snes.c > -> Cut time step: snes= 0 icut= 10[ 0] t= 0.00000E+00 dt= 9.76562E-04 > 1 2f: 8.79E-01 2x: 1.73E+06 2s: 1.70E+06 ir: 6.11E-01 iu: 1.37E+06 rsn: 0 > 2 2f: 1.34E+03 2x: 3.17E+11 2s: 3.17E+11 ir: 9.30E+02 iu: 3.14E+11 rsn: 0 > 3 2f: 3.99E+13 2x: 6.13E+11 2s: 3.04E+11 ir: 3.26E+13 iu: 2.98E+11 rsn: max_norm > -> Cut time step: snes=-10 icut= 11[ 0] t= 0.00000E+00 dt= 4.88281E-04 > 1 2f: 5.96E-01 2x: 1.33E+06 2s: 1.10E+06 ir: 4.19E-01 iu: 9.69E+05 rsn: 0 > 2 2f: 4.60E+04 2x: 3.54E+10 2s: 3.54E+10 ir: 3.25E+04 iu: 3.54E+10 rsn: 0 > 3 2f: 9.64E+15 2x: 3.91E+11 2s: 3.55E+11 ir: 7.87E+15 iu: 3.55E+11 rsn: max_norm > -> Cut time step: snes=-10 icut= 12[ 0] t= 0.00000E+00 dt= 2.44141E-04 > 1 2f: 3.54E-01 2x: 1.11E+06 2s: 6.36E+05 ir: 2.50E-01 iu: 5.95E+05 rsn: 0 > 2 2f: 5.42E+04 2x: 1.33E+09 2s: 1.33E+09 ir: 3.93E+04 iu: 1.33E+09 rsn: 0 > 3 2f: 3.51E+06 2x: 1.75E+10 2s: 1.75E+10 ir: 2.87E+06 iu: 1.75E+10 rsn: max_norm > -> Cut time step: snes=-10 icut= 13[ 0] t= 0.00000E+00 dt= 1.22070E-04 > 1 2f: 5.74E-02 2x: 1.04E+06 2s: 3.17E+05 ir: 4.02E-02 iu: 3.07E+05 rsn: 0 > 2 2f: 3.60E-02 2x: 6.66E+09 2s: 6.66E+09 ir: 2.88E-02 iu: 6.66E+09 rsn: 0 > 3 2f: 5.16E+02 2x: 2.91E+08 2s: 6.95E+09 ir: 3.31E+02 iu: 6.95E+09 rsn: 0 > 4 2f: 7.20E+04 2x: 7.65E+15 2s: 7.65E+15 ir: 6.73E+04 iu: 7.65E+15 rsn: 0 > 5 2f: 4.01E+12 2x: 3.15E+16 2s: 3.92E+16 ir: 3.27E+12 iu: 3.92E+16 rsn: max_norm > -> Cut time step: snes=-10 icut= 14[ 0] t= 0.00000E+00 dt= 6.10352E-05 > 1 2f: 8.13E-02 2x: 1.04E+06 2s: 1.33E+05 ir: 8.02E-02 iu: 1.31E+05 rsn: 0 > 2 2f: 1.94E+02 2x: 1.14E+08 2s: 1.14E+08 ir: 1.33E+02 iu: 1.14E+08 rsn: 0 > 3 2f: 2.16E-01 2x: 4.22E+08 2s: 5.36E+08 ir: 1.79E-01 iu: 5.36E+08 rsn: 0 > 4 2f: 3.21E+28 2x: 1.07E+19 2s: 1.07E+19 ir: 2.62E+28 iu: 1.07E+19 rsn: max_norm > -> Cut time step: snes=-10 icut= 15[ 0] t= 0.00000E+00 dt= 3.05176E-05 > 1 2f: 1.13E-05 2x: 1.05E+06 2s: 4.69E+04 ir: 1.13E-05 iu: 4.67E+04 rsn: 0 > 2 2f: 1.21E-10 2x: 1.05E+06 2s: 6.87E+02 ir: 1.05E-10 iu: 6.87E+02 rsn: itol_res > > Step 1 Time= 3.05176E-05 Dt= 3.05176E-05 [y] snes_conv_reason: 10 > newton = 47 [ 47] linear = 52 [ 52] cuts = 15 [ 15] > --> SNES Linear/Non-Linear Iterations = 2 / 2 > --> SNES Residual: 1.209152E-10 1.209152E-12 1.045172E-10 > --> max chng: dpmx= 4.5981E+04 dtmpmx= 1.6500E+00 > > > == TH FLOW ===================================================================== > 1 2f: 1.62E-01 2x: 1.04E+06 2s: 7.34E+04 ir: 1.60E-01 iu: 7.27E+04 rsn: 0 > 2 2f: 7.18E+01 2x: 6.44E+07 2s: 6.44E+07 ir: 4.68E+01 iu: 6.44E+07 rsn: 0 > 3 2f: 2.52E-01 2x: 5.48E+07 2s: 1.19E+08 ir: 2.02E-01 iu: 1.19E+08 rsn: 0 > 4 2f: 1.29E+09 2x: 1.44E+18 2s: 1.44E+18 ir: 1.29E+09 iu: 1.44E+18 rsn: max_norm > -> Cut time step: snes=-10 icut= 1[ 15] t= 3.05176E-05 dt= 1.52588E-05 > 1 2f: 3.23E-01 2x: 1.04E+06 2s: 3.15E+04 ir: 3.19E-01 iu: 3.13E+04 rsn: 0 > 2 2f: 1.17E-03 2x: 3.08E+06 2s: 2.83E+06 ir: 1.09E-03 iu: 2.83E+06 rsn: 0 > 3 2f: 3.23E-01 2x: 1.04E+06 2s: 2.83E+06 ir: 3.19E-01 iu: 2.83E+06 rsn: 0 > 4 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 5 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 6 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 7 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 8 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 9 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 10 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 11 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 12 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 13 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 14 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 15 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 16 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 17 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 18 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 19 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 20 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 21 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 22 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 23 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 24 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 25 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 26 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 27 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 28 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 29 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 30 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 31 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 32 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 33 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 34 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 35 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 36 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 37 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 38 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 39 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 40 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 41 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 42 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 43 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 44 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 45 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 46 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 47 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 48 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 49 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 50 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 51 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 52 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 53 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 54 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 55 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 56 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 57 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 58 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 59 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 60 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 61 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 62 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 63 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 64 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 65 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 66 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 67 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 68 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 69 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 70 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 71 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 72 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 73 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 74 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 75 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 76 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 77 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 78 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 79 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 80 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 81 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 82 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 83 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 84 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 85 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 86 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 87 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 88 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 89 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 90 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 91 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 92 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 93 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 94 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 95 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 96 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 97 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 98 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > 99 2f: 3.23E-01 2x: 1.04E+06 2s: 2.89E+06 ir: 3.19E-01 iu: 2.89E+06 rsn: 0 > 100 2f: 1.18E-03 2x: 3.14E+06 2s: 2.89E+06 ir: 1.10E-03 iu: 2.89E+06 rsn: 0 > -> Cut time step: snes= -5 icut= 2[ 15] t= 3.05176E-05 dt= 7.62939E-06 > 1 2f: 6.23E-01 2x: 1.04E+06 2s: 1.39E+04 ir: 6.15E-01 iu: 1.38E+04 rsn: 0 > 2 2f: 6.18E-04 2x: 1.06E+06 2s: 1.15E+05 ir: 5.25E-04 iu: 1.15E+05 rsn: 0 > 3 2f: 6.23E-01 2x: 1.04E+06 2s: 1.15E+05 ir: 6.15E-01 iu: 1.15E+05 rsn: 0 > 4 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 5 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 6 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 7 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 8 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 9 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 10 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 11 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 12 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 13 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 14 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 15 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 16 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 17 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 18 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 19 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 20 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 21 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 22 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 23 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 24 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 25 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 26 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 27 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 28 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 29 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 30 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 31 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 32 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 33 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 34 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 35 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 36 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 37 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 38 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 39 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 40 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 41 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 42 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 43 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 44 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 45 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 46 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 47 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 48 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 49 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 50 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 51 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 52 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 53 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 54 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 55 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 56 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 57 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 58 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 59 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 60 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 61 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 62 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 63 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 64 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 65 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 66 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 67 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 68 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 69 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 70 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 71 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 72 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 73 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 74 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 75 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 76 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 77 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 78 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 79 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 80 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 81 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 82 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 83 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 84 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 85 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 86 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 87 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 88 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 89 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 90 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 91 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 92 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 93 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 94 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 95 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 96 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 97 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 98 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > 99 2f: 6.23E-01 2x: 1.04E+06 2s: 1.16E+05 ir: 6.15E-01 iu: 1.16E+05 rsn: 0 > 100 2f: 6.18E-04 2x: 1.06E+06 2s: 1.16E+05 ir: 5.25E-04 iu: 1.16E+05 rsn: 0 > -> Cut time step: snes= -5 icut= 3[ 15] t= 3.05176E-05 dt= 3.81470E-06 > 1 2f: 8.15E-01 2x: 1.05E+06 2s: 6.36E+03 ir: 8.05E-01 iu: 6.33E+03 rsn: 0 > 2 2f: 1.34E-03 2x: 1.05E+06 2s: 4.20E+03 ir: 1.31E-03 iu: 4.19E+03 rsn: 0 > 3 2f: 4.65E-04 2x: 1.05E+06 2s: 2.99E+02 ir: 4.59E-04 iu: 2.50E+02 rsn: 0 > 4 2f: 1.43E-04 2x: 1.05E+06 2s: 1.18E+02 ir: 1.41E-04 iu: 1.18E+02 rsn: 0 > 5 2f: 4.35E-05 2x: 1.05E+06 2s: 8.37E+01 ir: 4.30E-05 iu: 8.37E+01 rsn: 0 > 6 2f: 1.27E-05 2x: 1.05E+06 2s: 5.80E+01 ir: 1.25E-05 iu: 5.80E+01 rsn: 0 > 7 2f: 3.26E-06 2x: 1.05E+06 2s: 3.69E+01 ir: 3.22E-06 iu: 3.69E+01 rsn: 0 > 8 2f: 5.57E-07 2x: 1.05E+06 2s: 1.83E+01 ir: 5.50E-07 iu: 1.83E+01 rsn: 0 > 9 2f: 2.98E-08 2x: 1.05E+06 2s: 4.65E+00 ir: 2.95E-08 iu: 4.65E+00 rsn: 0 > 10 2f: 9.92E-11 2x: 1.05E+06 2s: 2.79E-01 ir: 9.78E-11 iu: 2.79E-01 rsn: itol_res > > Step 2 Time= 3.43323E-05 Dt= 3.81470E-06 [y] snes_conv_reason: 10 > newton = 214 [ 261] linear = 214 [ 266] cuts = 3 [ 18] > --> SNES Linear/Non-Linear Iterations = 10 / 10 > --> SNES Residual: 9.921423E-11 9.921423E-13 9.778422E-11 > --> max chng: dpmx= 1.6597E+03 dtmpmx= 1.9780E-01 > > > -- > You received this message because you are subscribed to the Google Groups "pflotran-dev" group. > To view this discussion on the web visit https://groups.google.com/d/msgid/pflotran-dev/CAPz1TnfKQJov5Uxvghp_31uKpxz…<https://groups.google.com/d/msgid/pflotran-dev/CAPz1TnfKQJov5Uxvghp_31uKpxz…>. > > -- > You received this message because you are subscribed to the Google Groups "pflotran-dev" group. > To view this discussion on the web visit https://groups.google.com/d/msgid/pflotran-dev/92b8b510fceb47309cc0b8f7ada4….

1 0