On Mon, Jun 24, 2024 at 2:50 PM MIGUEL MOLINOS PEREZ <mmoli...@us.es> wrote:
> 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) > ZjQcmQRYFpfptBannerStart > This Message Is From an External Sender > This message came from outside your organization. > > ZjQcmQRYFpfptBannerEnd > 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! > There is something inconsistent among the processes. First, I would try running with a constant A. Then execute your nonlinear solve, but return the constant A. If that passes, then likely you are returning inconsistent results across processes with your solve. Thanks, Matt > 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!YlgDM4WJ1LZd1CuLZGNlMThj7OBlwkCu03P9vWxFt-dOODaCkEboruflECm3ntB9Grn1TGgk5Ic8lfRVlzr2$ > > <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!YlgDM4WJ1LZd1CuLZGNlMThj7OBlwkCu03P9vWxFt-dOODaCkEboruflECm3ntB9Grn1TGgk5Ic8lfRVlzr2$ > > <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!YlgDM4WJ1LZd1CuLZGNlMThj7OBlwkCu03P9vWxFt-dOODaCkEboruflECm3ntB9Grn1TGgk5Ic8lfRVlzr2$ > > <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!YlgDM4WJ1LZd1CuLZGNlMThj7OBlwkCu03P9vWxFt-dOODaCkEboruflECm3ntB9Grn1TGgk5Ic8lfRVlzr2$ > > <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!YlgDM4WJ1LZd1CuLZGNlMThj7OBlwkCu03P9vWxFt-dOODaCkEboruflECm3ntB9Grn1TGgk5Ic8lfRVlzr2$ > > <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!YlgDM4WJ1LZd1CuLZGNlMThj7OBlwkCu03P9vWxFt-dOODaCkEboruflECm3ntB9Grn1TGgk5Ic8lfRVlzr2$ > > <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!YlgDM4WJ1LZd1CuLZGNlMThj7OBlwkCu03P9vWxFt-dOODaCkEboruflECm3ntB9Grn1TGgk5Ic8lfRVlzr2$ > > <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!YlgDM4WJ1LZd1CuLZGNlMThj7OBlwkCu03P9vWxFt-dOODaCkEboruflECm3ntB9Grn1TGgk5Ic8lfRVlzr2$ > > <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 > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener https://urldefense.us/v3/__https://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!YlgDM4WJ1LZd1CuLZGNlMThj7OBlwkCu03P9vWxFt-dOODaCkEboruflECm3ntB9Grn1TGgk5Ic8lUlphb0M$ <https://urldefense.us/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!YlgDM4WJ1LZd1CuLZGNlMThj7OBlwkCu03P9vWxFt-dOODaCkEboruflECm3ntB9Grn1TGgk5Ic8lXZDmbEm$ >