Note you can uses SNESSetOptionsPrefix(snes,"myprefix_") to attach a prefix for the one dimensional solves and then set options for the one dimensional problems with -myprefix_snes_type etc Now you can set different options for the different types of snes in your code
> On Jun 24, 2024, at 4:18 PM, MIGUEL MOLINOS PEREZ <mmoli...@us.es> wrote: > > Thanks Barry, that was the problem! > > Thank you, > Miguel > >> On Jun 24, 2024, at 12:56 PM, Barry Smith <bsm...@petsc.dev> wrote: >> >> >> The error is coming from the parallel SNES used by TS TSTheta_SNESSolve() >> not your small SNES solver but TSTheta_SNESSolve() >> >> Process 7 uses [7]PETSC ERROR: #4 SNESLineSearchApply_CP() while the rest >> uses _BT(). I think this causes the problem since they are different >> algorithms using different MPI. >> >> Are you perhaps setting somewhere the line search to use? Options >> database or someway that different processes will get a different value? >> >> >> >>> On Jun 24, 2024, at 2:50 PM, MIGUEL MOLINOS PEREZ <mmoli...@us.es> wrote: >>> >>> This Message Is From an External Sender >>> This message came from outside your organization. >>> Dear all, >>> >>> I am trying to assemble a matrix A with coefficients which I need to >>> assemble the RHS (F) and its Jacobian (J) in a TS type of problem. >>> >>> Determining each coefficient of A involves the resolution of a small >>> non-linear problem (1 dof) using the serial version of SNES. By the way, >>> the matrix A is of the type “MATMPIAIJ”. >>> >>> The weird part is, if I pass the matrix A to the TS routine inside of a >>> user-context structure, without even accessing to the values inside of A, I >>> got the following error message: >>> > [1]PETSC ERROR: Petsc has generated inconsistent data >>> > [1]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) >>> > on different processors >>> But if I comment out the line which calls the SNES routine used to evaluate >>> the coefficients inside of A, I don’t get the error message. >>> >>> Some additional context: >>> - The SNES routine is called once at a time inside of each rank. >>> - I use PetscCall(SNESCreate(PETSC_COMM_SELF, >>> &snes)); >>> - The vectors inside of the SNES function are defined as follows: >>> VecCreateSeq(PETSC_COMM_SELF, >>> 1, >>> &Y) >>> - All the input fields for SNES are also sequential. >>> >>> Any feedback is greatly appreciated! >>> >>> Thanks, >>> Miguel >>> >>> >>> >>> [test-Mass-Transport-Master-Equation-PETSc-Backward-Euler][MgHx-hcp-x5x5x5-cell] >>> t=0.0000e+00 dt=1.0000e-07 it=( 0, 0) >>> 0 KSP Residual norm 4.776631889125e-07 >>> 1 KSP Residual norm 6.807505564283e-17 >>> [0]PETSC ERROR: --------------------- Error Message >>> -------------------------------------------------------------- >>> [0]PETSC ERROR: Petsc has generated inconsistent data >>> [5]PETSC ERROR: --------------------- Error Message >>> -------------------------------------------------------------- >>> [5]PETSC ERROR: Petsc has generated inconsistent data >>> [7]PETSC ERROR: --------------------- Error Message >>> -------------------------------------------------------------- >>> [7]PETSC ERROR: Petsc has generated inconsistent data >>> [0]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) >>> on different processors >>> [1]PETSC ERROR: --------------------- Error Message >>> -------------------------------------------------------------- >>> [1]PETSC ERROR: Petsc has generated inconsistent data >>> [1]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) >>> on different processors >>> [2]PETSC ERROR: --------------------- Error Message >>> -------------------------------------------------------------- >>> [2]PETSC ERROR: Petsc has generated inconsistent data >>> [2]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) >>> on different processors >>> [3]PETSC ERROR: --------------------- Error Message >>> -------------------------------------------------------------- >>> [3]PETSC ERROR: Petsc has generated inconsistent data >>> [3]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) >>> on different processors >>> [4]PETSC ERROR: --------------------- Error Message >>> -------------------------------------------------------------- >>> [4]PETSC ERROR: Petsc has generated inconsistent data >>> [4]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) >>> on different processors >>> [5]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) >>> on different processors >>> [6]PETSC ERROR: --------------------- Error Message >>> -------------------------------------------------------------- >>> [6]PETSC ERROR: Petsc has generated inconsistent data >>> [6]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) >>> on different processors >>> [7]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) >>> on different processors >>> [0]PETSC ERROR: WARNING! There are unused option(s) set! Could be the >>> program crashed before usage or a spelling mistake, etc! >>> [1]PETSC ERROR: WARNING! There are unused option(s) set! Could be the >>> program crashed before usage or a spelling mistake, etc! >>> [2]PETSC ERROR: WARNING! There are unused option(s) set! Could be the >>> program crashed before usage or a spelling mistake, etc! >>> [2]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code >>> [2]PETSC ERROR: [3]PETSC ERROR: WARNING! There are unused option(s) set! >>> Could be the program crashed before usage or a spelling mistake, etc! >>> [3]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code >>> [3]PETSC ERROR: See >>> https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!bC0DolaKYXB2ZhjL209T3VwDP2ypey_YUEI5guqAfAa7gn_JtXRJQ3v5QhgLFf432jn25U0OHqIzoJXRmLdKKFA$ >>> >>> <https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$> >>> for trouble shooting. >>> [4]PETSC ERROR: WARNING! There are unused option(s) set! Could be the >>> program crashed before usage or a spelling mistake, etc! >>> [4]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code >>> [4]PETSC ERROR: See >>> https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!bC0DolaKYXB2ZhjL209T3VwDP2ypey_YUEI5guqAfAa7gn_JtXRJQ3v5QhgLFf432jn25U0OHqIzoJXRmLdKKFA$ >>> >>> <https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$> >>> for trouble shooting. >>> [5]PETSC ERROR: WARNING! There are unused option(s) set! Could be the >>> program crashed before usage or a spelling mistake, etc! >>> [5]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code >>> [5]PETSC ERROR: See >>> https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!bC0DolaKYXB2ZhjL209T3VwDP2ypey_YUEI5guqAfAa7gn_JtXRJQ3v5QhgLFf432jn25U0OHqIzoJXRmLdKKFA$ >>> >>> <https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$> >>> for trouble shooting. >>> [5]PETSC ERROR: Petsc Release Version 3.21.0, unknown >>> [6]PETSC ERROR: WARNING! There are unused option(s) set! Could be the >>> program crashed before usage or a spelling mistake, etc! >>> [6]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code >>> [6]PETSC ERROR: See >>> https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!bC0DolaKYXB2ZhjL209T3VwDP2ypey_YUEI5guqAfAa7gn_JtXRJQ3v5QhgLFf432jn25U0OHqIzoJXRmLdKKFA$ >>> >>> <https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$> >>> for trouble shooting. >>> [6]PETSC ERROR: Petsc Release Version 3.21.0, unknown >>> [6]PETSC ERROR: [7]PETSC ERROR: WARNING! There are unused option(s) set! >>> Could be the program crashed before usage or a spelling mistake, etc! >>> [7]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code >>> [0]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code >>> [0]PETSC ERROR: See >>> https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!bC0DolaKYXB2ZhjL209T3VwDP2ypey_YUEI5guqAfAa7gn_JtXRJQ3v5QhgLFf432jn25U0OHqIzoJXRmLdKKFA$ >>> >>> <https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$> >>> for trouble shooting. >>> [0]PETSC ERROR: Petsc Release Version 3.21.0, unknown >>> [0]PETSC ERROR: ./exe-tasting-SOLERA on a arch-darwin-c-debug named >>> mmp-laptop.local by migmolper Mon Jun 24 11:37:38 2024 >>> [1]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code >>> [1]PETSC ERROR: See >>> https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!bC0DolaKYXB2ZhjL209T3VwDP2ypey_YUEI5guqAfAa7gn_JtXRJQ3v5QhgLFf432jn25U0OHqIzoJXRmLdKKFA$ >>> >>> <https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$> >>> for trouble shooting. >>> See >>> https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!bC0DolaKYXB2ZhjL209T3VwDP2ypey_YUEI5guqAfAa7gn_JtXRJQ3v5QhgLFf432jn25U0OHqIzoJXRmLdKKFA$ >>> >>> <https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$> >>> for trouble shooting. >>> [2]PETSC ERROR: Petsc Release Version 3.21.0, unknown >>> [2]PETSC ERROR: ./exe-tasting-SOLERA on a arch-darwin-c-debug named >>> mmp-laptop.local by migmolper Mon Jun 24 11:37:38 2024 >>> [2]PETSC ERROR: Configure options --download-hdf5=1 --download-mpich=1 >>> --with-debugging=1 CC=gcc CXX=c++ PETSC_ARCH=arch-darwin-c-debug --with-x >>> [3]PETSC ERROR: Petsc Release Version 3.21.0, unknown >>> [3]PETSC ERROR: ./exe-tasting-SOLERA on a arch-darwin-c-debug named >>> mmp-laptop.local by migmolper Mon Jun 24 11:37:38 2024 >>> [3]PETSC ERROR: Configure options --download-hdf5=1 --download-mpich=1 >>> --with-debugging=1 CC=gcc CXX=c++ PETSC_ARCH=arch-darwin-c-debug --with-x >>> [4]PETSC ERROR: Petsc Release Version 3.21.0, unknown >>> [4]PETSC ERROR: ./exe-tasting-SOLERA on a arch-darwin-c-debug named >>> mmp-laptop.local by migmolper Mon Jun 24 11:37:38 2024 >>> [4]PETSC ERROR: Configure options --download-hdf5=1 --download-mpich=1 >>> --with-debugging=1 CC=gcc CXX=c++ PETSC_ARCH=arch-darwin-c-debug --with-x >>> [5]PETSC ERROR: ./exe-tasting-SOLERA on a arch-darwin-c-debug named >>> mmp-laptop.local by migmolper Mon Jun 24 11:37:38 2024 >>> [5]PETSC ERROR: Configure options --download-hdf5=1 --download-mpich=1 >>> --with-debugging=1 CC=gcc CXX=c++ PETSC_ARCH=arch-darwin-c-debug --with-x >>> ./exe-tasting-SOLERA on a arch-darwin-c-debug named mmp-laptop.local by >>> migmolper Mon Jun 24 11:37:38 2024 >>> [6]PETSC ERROR: Configure options --download-hdf5=1 --download-mpich=1 >>> --with-debugging=1 CC=gcc CXX=c++ PETSC_ARCH=arch-darwin-c-debug --with-x >>> See >>> https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!bC0DolaKYXB2ZhjL209T3VwDP2ypey_YUEI5guqAfAa7gn_JtXRJQ3v5QhgLFf432jn25U0OHqIzoJXRmLdKKFA$ >>> >>> <https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$> >>> for trouble shooting. >>> [7]PETSC ERROR: Petsc Release Version 3.21.0, unknown >>> [7]PETSC ERROR: ./exe-tasting-SOLERA on a arch-darwin-c-debug named >>> mmp-laptop.local by migmolper Mon Jun 24 11:37:38 2024 >>> [7]PETSC ERROR: Configure options --download-hdf5=1 --download-mpich=1 >>> --with-debugging=1 CC=gcc CXX=c++ PETSC_ARCH=arch-darwin-c-debug --with-x >>> Configure options --download-hdf5=1 --download-mpich=1 --with-debugging=1 >>> CC=gcc CXX=c++ PETSC_ARCH=arch-darwin-c-debug --with-x >>> [1]PETSC ERROR: Petsc Release Version 3.21.0, unknown >>> [1]PETSC ERROR: ./exe-tasting-SOLERA on a arch-darwin-c-debug named >>> mmp-laptop.local by migmolper Mon Jun 24 11:37:38 2024 >>> [1]PETSC ERROR: Configure options --download-hdf5=1 --download-mpich=1 >>> --with-debugging=1 CC=gcc CXX=c++ PETSC_ARCH=arch-darwin-c-debug --with-x >>> [0]PETSC ERROR: #1 PetscSplitReductionApply() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230 >>> [2]PETSC ERROR: #1 PetscSplitReductionApply() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230 >>> [2]PETSC ERROR: [3]PETSC ERROR: #1 PetscSplitReductionApply() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230 >>> [3]PETSC ERROR: #2 PetscSplitReductionEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172 >>> [4]PETSC ERROR: #1 PetscSplitReductionApply() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230 >>> [4]PETSC ERROR: #2 PetscSplitReductionEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172 >>> [5]PETSC ERROR: #1 PetscSplitReductionApply() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230 >>> [5]PETSC ERROR: #2 PetscSplitReductionEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172 >>> [0]PETSC ERROR: #2 PetscSplitReductionEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172 >>> [0]PETSC ERROR: #3 VecNormEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553 >>> [0]PETSC ERROR: #4 SNESLineSearchApply_BT() at >>> /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88 >>> [1]PETSC ERROR: #1 PetscSplitReductionApply() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230 >>> [1]PETSC ERROR: #2 PetscSplitReductionEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172 >>> [1]PETSC ERROR: #3 VecNormEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553 >>> #2 PetscSplitReductionEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172 >>> [2]PETSC ERROR: #3 VecNormEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553 >>> [2]PETSC ERROR: #4 SNESLineSearchApply_BT() at >>> /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88 >>> [2]PETSC ERROR: [3]PETSC ERROR: #3 VecNormEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553 >>> [3]PETSC ERROR: #4 SNESLineSearchApply_BT() at >>> /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88 >>> [3]PETSC ERROR: #5 SNESLineSearchApply() at >>> /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645 >>> [4]PETSC ERROR: #3 VecNormEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553 >>> [4]PETSC ERROR: #4 SNESLineSearchApply_BT() at >>> /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88 >>> [4]PETSC ERROR: #5 SNESLineSearchApply() at >>> /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645 >>> [5]PETSC ERROR: #3 VecNormEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553 >>> [5]PETSC ERROR: #4 SNESLineSearchApply_BT() at >>> /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88 >>> [5]PETSC ERROR: #5 SNESLineSearchApply() at >>> /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645 >>> [5]PETSC ERROR: [6]PETSC ERROR: #1 PetscSplitReductionApply() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230 >>> [6]PETSC ERROR: #2 PetscSplitReductionEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172 >>> [6]PETSC ERROR: #3 VecNormEnd() at >>> /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553 >>> [6]PETSC ERROR: #4 SNESLineSearchApply_BT() at >>> /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88 >>> [6]PETSC ERROR: #5 SNESLineSearchApply() at >>> /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645 >>> [7]PETSC ERROR: #1 VecXDot_MPI_Default() at >>> /Users/migmolper/petsc/include/../src/vec/vec/impls/mpi/pvecimpl.h:107 >>> [7]PETSC ERROR: #2 VecDot_MPI() at >>> /Users/migmolper/petsc/src/vec/vec/impls/mpi/pvec2.c:10 >>> [7]PETSC ERROR: #3 VecDot() at >>> /Users/migmolper/petsc/src/vec/vec/interface/rvector.c:120 >>> [7]PETSC ERROR: #4 SNESLineSearchApply_CP() at >>> /Users/migmolper/petsc/src/snes/linesearch/impls/cp/linesearchcp.c:28 >>> [0]PETSC ERROR: #5 SNESLineSearchApply() at >>> /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645 >>> [0]PETSC ERROR: #6 SNESSolve_NEWTONLS() at >>> /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234 >>> [0]PETSC ERROR: #7 SNESSolve() at >>> /Users/migmolper/petsc/src/snes/interface/snes.c:4738 >>> [0]PETSC ERROR: #8 TSTheta_SNESSolve() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174 >>> [1]PETSC ERROR: #4 SNESLineSearchApply_BT() at >>> /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88 >>> [1]PETSC ERROR: #5 SNESLineSearchApply() at >>> /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645 >>> [1]PETSC ERROR: #6 SNESSolve_NEWTONLS() at >>> /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234 >>> [1]PETSC ERROR: #7 SNESSolve() at >>> /Users/migmolper/petsc/src/snes/interface/snes.c:4738 >>> [1]PETSC ERROR: #8 TSTheta_SNESSolve() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174 >>> [1]PETSC ERROR: #9 TSStep_Theta() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225 >>> [1]PETSC ERROR: #10 TSStep() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:3391 >>> [1]PETSC ERROR: #11 TSSolve() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:4037 >>> #5 SNESLineSearchApply() at >>> /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645 >>> [2]PETSC ERROR: #6 SNESSolve_NEWTONLS() at >>> /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234 >>> [2]PETSC ERROR: #7 SNESSolve() at >>> /Users/migmolper/petsc/src/snes/interface/snes.c:4738 >>> [2]PETSC ERROR: #8 TSTheta_SNESSolve() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174 >>> [2]PETSC ERROR: #9 TSStep_Theta() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225 >>> [2]PETSC ERROR: #10 TSStep() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:3391 >>> [2]PETSC ERROR: #11 TSSolve() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:4037 >>> [3]PETSC ERROR: #6 SNESSolve_NEWTONLS() at >>> /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234 >>> [3]PETSC ERROR: #7 SNESSolve() at >>> /Users/migmolper/petsc/src/snes/interface/snes.c:4738 >>> [3]PETSC ERROR: #8 TSTheta_SNESSolve() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174 >>> [3]PETSC ERROR: #9 TSStep_Theta() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225 >>> [3]PETSC ERROR: #10 TSStep() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:3391 >>> [3]PETSC ERROR: #11 TSSolve() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:4037 >>> [4]PETSC ERROR: #6 SNESSolve_NEWTONLS() at >>> /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234 >>> [4]PETSC ERROR: #7 SNESSolve() at >>> /Users/migmolper/petsc/src/snes/interface/snes.c:4738 >>> [4]PETSC ERROR: #8 TSTheta_SNESSolve() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174 >>> [4]PETSC ERROR: #9 TSStep_Theta() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225 >>> [4]PETSC ERROR: #10 TSStep() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:3391 >>> [4]PETSC ERROR: #11 TSSolve() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:4037 >>> #6 SNESSolve_NEWTONLS() at /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234 >>> [5]PETSC ERROR: #7 SNESSolve() at >>> /Users/migmolper/petsc/src/snes/interface/snes.c:4738 >>> [5]PETSC ERROR: #8 TSTheta_SNESSolve() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174 >>> [5]PETSC ERROR: #9 TSStep_Theta() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225 >>> [5]PETSC ERROR: #10 TSStep() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:3391 >>> [5]PETSC ERROR: #11 TSSolve() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:4037 >>> [6]PETSC ERROR: #6 SNESSolve_NEWTONLS() at >>> /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234 >>> [6]PETSC ERROR: #7 SNESSolve() at >>> /Users/migmolper/petsc/src/snes/interface/snes.c:4738 >>> [6]PETSC ERROR: #8 TSTheta_SNESSolve() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174 >>> [6]PETSC ERROR: #9 TSStep_Theta() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225 >>> [6]PETSC ERROR: #10 TSStep() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:3391 >>> [6]PETSC ERROR: #11 TSSolve() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:4037 >>> [7]PETSC ERROR: #5 SNESLineSearchApply() at >>> /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645 >>> [7]PETSC ERROR: #6 SNESSolve_NEWTONLS() at >>> /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234 >>> [7]PETSC ERROR: #7 SNESSolve() at >>> /Users/migmolper/petsc/src/snes/interface/snes.c:4738 >>> [7]PETSC ERROR: #8 TSTheta_SNESSolve() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174 >>> [7]PETSC ERROR: #9 TSStep_Theta() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225 >>> [7]PETSC ERROR: #10 TSStep() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:3391 >>> [7]PETSC ERROR: #11 TSSolve() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:4037 >>> [0]PETSC ERROR: #9 TSStep_Theta() at >>> /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225 >>> [0]PETSC ERROR: #10 TSStep() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:3391 >>> [0]PETSC ERROR: #11 TSSolve() at >>> /Users/migmolper/petsc/src/ts/interface/ts.c:4037 >>> [0]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at >>> /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251 >>> [2]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at >>> /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251 >>> [3]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at >>> /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251 >>> [1]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at >>> /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251 >>> [4]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at >>> /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251 >>> [5]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at >>> /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251 >>> [6]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at >>> /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251 >>> [7]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at >>> /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251 >> >
