<head><!-- BaNnErBlUrFlE-HeAdEr-start -->
<style>
#pfptBanner7183y32 { all: revert !important; display: block !important;
visibility: visible !important; opacity: 1 !important;
background-color: #D0D8DC !important;
max-width: none !important; max-height: none !important }
.pfptPrimaryButton7183y32:hover, .pfptPrimaryButton7183y32:focus {
background-color: #b4c1c7 !important; }
.pfptPrimaryButton7183y32:active {
background-color: #90a4ae !important; }
</style>
<!-- BaNnErBlUrFlE-HeAdEr-end -->
</head><!-- BaNnErBlUrFlE-BoDy-start -->
<!-- Preheader Text : BEGIN -->
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">
On Mon, Jun 24, 2024 at 2: 50 PM MIGUEL MOLINOS PEREZ <mmolinos@ 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
</div>
<!-- Preheader Text : END -->
<!-- Email Banner : BEGIN -->
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">ZjQcmQRYFpfptBannerStart</div>
<!--[if ((ie)|(mso))]>
<table border="0" cellspacing="0" cellpadding="0" width="100%" style="padding: 16px 0px 16px 0px; direction: ltr" ><tr><td>
<table border="0" cellspacing="0" cellpadding="0" style="padding: 0px 10px 5px 6px; width: 100%; border-radius:4px; border-top:4px solid #90a4ae;background-color:#D0D8DC;"><tr><td valign="top">
<table align="left" border="0" cellspacing="0" cellpadding="0" style="padding: 4px 8px 4px 8px">
<tr><td style="color:#000000; font-family: 'Arial', sans-serif; font-weight:bold; font-size:14px; direction: ltr">
This Message Is From an External Sender
</td></tr>
<tr><td style="color:#000000; font-weight:normal; font-family: 'Arial', sans-serif; font-size:12px; direction: ltr">
This message came from outside your organization.
</td></tr>
</table>
</td></tr></table>
</td></tr></table>
<![endif]-->
<![if !((ie)|(mso))]>
<div dir="ltr" id="pfptBanner7183y32" style="all: revert !important; display:block !important; text-align: left !important; margin:16px 0px 16px 0px !important; padding:8px 16px 8px 16px !important; border-radius: 4px !important; min-width: 200px !important; background-color: #D0D8DC !important; background-color: #D0D8DC; border-top: 4px solid #90a4ae !important; border-top: 4px solid #90a4ae;">
<div id="pfptBanner7183y32" style="all: unset !important; float:left !important; display:block !important; margin: 0px 0px 1px 0px !important; max-width: 600px !important;">
<div id="pfptBanner7183y32" style="all: unset !important; display:block !important; visibility: visible !important; background-color: #D0D8DC !important; color:#000000 !important; color:#000000; font-family: 'Arial', sans-serif !important; font-family: 'Arial', sans-serif; font-weight:bold !important; font-weight:bold; font-size:14px !important; line-height:18px !important; line-height:18px">
This Message Is From an External Sender
</div>
<div id="pfptBanner7183y32" style="all: unset !important; display:block !important; visibility: visible !important; background-color: #D0D8DC !important; color:#000000 !important; color:#000000; font-weight:normal; font-family: 'Arial', sans-serif !important; font-family: 'Arial', sans-serif; font-size:12px !important; line-height:18px !important; line-height:18px; margin-top:2px !important;">
This message came from outside your organization.
</div>
</div>
<div style="clear: both !important; display: block !important; visibility: hidden !important; line-height: 0 !important; font-size: 0.01px !important; height: 0px"> </div>
</div>
<![endif]>
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">ZjQcmQRYFpfptBannerEnd</div>
<!-- Email Banner : END -->
<!-- BaNnErBlUrFlE-BoDy-end -->
<div dir="ltr"><div dir="ltr">On Mon, Jun 24, 2024 at 2:50 PM MIGUEL MOLINOS PEREZ <<a href="mailto:mmolinos@us.es">mmolinos@us.es</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg770091079445759359">
<div style="font-size:1px;color:rgb(255,255,255);line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;display:none">
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)
</div>
<div style="font-size:1px;color:rgb(255,255,255);line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;display:none">ZjQcmQRYFpfptBannerStart</div>
<u></u>
<div dir="ltr" id="m_770091079445759359pfptBanner8lct7dh" style="display:block;text-align:left;margin:16px 0px;padding:8px 16px;border-radius:4px;min-width:200px;background-color:rgb(208,216,220);border-top:4px solid rgb(144,164,174)">
<div id="m_770091079445759359pfptBanner8lct7dh" style="float:left;display:block;margin:0px 0px 1px;max-width:600px">
<div id="m_770091079445759359pfptBanner8lct7dh" style="display:block;background-color:rgb(208,216,220);color:rgb(0,0,0);font-family:Arial,sans-serif;font-weight:bold;font-size:14px;line-height:18px">
This Message Is From an External Sender
</div>
<div id="m_770091079445759359pfptBanner8lct7dh" style="font-weight:normal;display:block;background-color:rgb(208,216,220);color:rgb(0,0,0);font-family:Arial,sans-serif;font-size:12px;line-height:18px;margin-top:2px">
This message came from outside your organization.
</div>
</div>
<div style="height:0px;clear:both;display:block;line-height:0;font-size:0.01px"> </div>
</div>
<u></u>
<div style="font-size:1px;color:rgb(255,255,255);line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;display:none">ZjQcmQRYFpfptBannerEnd</div>
<div>
Dear all,
<div><br>
</div>
<div>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.</div>
<div><br>
</div>
<div>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, t<span style="color:rgb(0,0,0)">he matrix A is of the type</span><span style="color:rgb(0,0,0)"> </span> “<span style="color:rgb(46,46,46);font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255)">MATMPIAIJ</span>”.</div>
<div><br>
</div>
<div>The weird part is, if I pass the matrix A <span style="color:rgb(0,0,0)"> </span><span style="color:rgb(0,0,0)">to the TS routine</span> inside of a user-context structure, without even accessing
to the values inside of A, I got the following error message:</div>
<div>
<pre style="color:rgb(0,0,0);background-color:rgb(255,255,255)">><i> [1]PETSC ERROR: Petsc has generated inconsistent data
</i>><i> [1]PETSC ERROR: MPI_Allreduce() called in different locations (code lines)
</i>><i> on different processors</i></pre>
</div>
<div>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. </div>
<div><br>
</div>
<div>Some additional context:</div>
<div>- The SNES routine is called once at a time inside of each rank.</div>
<div>- I use <span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(99,99,36)">PetscCall</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(94,94,94)">(</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(99,99,36)">SNESCreate</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(94,94,94)">(</span><span style="color:rgb(46,46,46);font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255)">PETSC_COMM_SELF</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(94,94,94)">,</span><span style="color:rgb(46,46,46);font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255)">
</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(94,94,94)">&</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(9,89,132)">snes</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(94,94,94)">));</span></div>
<div><font color="#000000"><span style="background-color:rgb(255,255,255)"><span style="white-space:pre-wrap">- The vectors inside of the SNES function are defined as follows:
</span></span></font><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(99,99,36)">VecCreateSeq</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(94,94,94)">(</span><span style="color:rgb(46,46,46);font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255)">PETSC_COMM_SELF</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(94,94,94)">,</span><span style="color:rgb(46,46,46);font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255)">
</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(73,104,57)">1</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(94,94,94)">,</span><span style="color:rgb(46,46,46);font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255)">
</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(94,94,94)">&</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(9,89,132)">Y</span><span style="font-family:Menlo,Monaco,"Courier New",monospace;white-space:pre-wrap;background-color:rgb(255,255,255);color:rgb(94,94,94)">)</span></div>
<div><font color="#000000"><span style="background-color:rgb(255,255,255)"><span style="white-space:pre-wrap">- All the input fields for SNES are also sequential</span>.</span></font></div>
<div><font color="#000000"><span style="background-color:rgb(255,255,255)"><br>
</span></font></div>
<div><font color="#000000"><span style="background-color:rgb(255,255,255)">Any feedback is greatly appreciated!</span></font></div></div></div></blockquote><div><br></div><div>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.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg770091079445759359"><div>
<div>Thanks,</div>
<div>Miguel</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div>[test-Mass-Transport-Master-Equation-PETSc-Backward-Euler][MgHx-hcp-x5x5x5-cell]</div>
<div>t=0.0000e+00 dt=1.0000e-07 it=( 0, 0)</div>
<div> 0 KSP Residual norm 4.776631889125e-07</div>
<div> 1 KSP Residual norm 6.807505564283e-17</div>
<div>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div>
<div>[0]PETSC ERROR: Petsc has generated inconsistent data</div>
<div>[5]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div>
<div>[5]PETSC ERROR: Petsc has generated inconsistent data</div>
<div>[7]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div>
<div>[7]PETSC ERROR: Petsc has generated inconsistent data</div>
<div>[0]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) on different processors</div>
<div>[1]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div>
<div>[1]PETSC ERROR: Petsc has generated inconsistent data</div>
<div>[1]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) on different processors</div>
<div>[2]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div>
<div>[2]PETSC ERROR: Petsc has generated inconsistent data</div>
<div>[2]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) on different processors</div>
<div>[3]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div>
<div>[3]PETSC ERROR: Petsc has generated inconsistent data</div>
<div>[3]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) on different processors</div>
<div>[4]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div>
<div>[4]PETSC ERROR: Petsc has generated inconsistent data</div>
<div>[4]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) on different processors</div>
<div>[5]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) on different processors</div>
<div>[6]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div>
<div>[6]PETSC ERROR: Petsc has generated inconsistent data</div>
<div>[6]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) on different processors</div>
<div>[7]PETSC ERROR: MPI_Allreduce() called in different locations (code lines) on different processors</div>
<div>[0]PETSC ERROR: WARNING! There are unused option(s) set! Could be the program crashed before usage or a spelling mistake, etc!</div>
<div>[1]PETSC ERROR: WARNING! There are unused option(s) set! Could be the program crashed before usage or a spelling mistake, etc!</div>
<div>[2]PETSC ERROR: WARNING! There are unused option(s) set! Could be the program crashed before usage or a spelling mistake, etc!</div>
<div>[2]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code</div>
<div>[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!</div>
<div>[3]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code</div>
<div>[3]PETSC ERROR: See <a href="https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$" target="_blank">https://petsc.org/release/faq/</a> for trouble shooting.</div>
<div>[4]PETSC ERROR: WARNING! There are unused option(s) set! Could be the program crashed before usage or a spelling mistake, etc!</div>
<div>[4]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code</div>
<div>[4]PETSC ERROR: See <a href="https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$" target="_blank">https://petsc.org/release/faq/</a> for trouble shooting.</div>
<div>[5]PETSC ERROR: WARNING! There are unused option(s) set! Could be the program crashed before usage or a spelling mistake, etc!</div>
<div>[5]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code</div>
<div>[5]PETSC ERROR: See <a href="https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$" target="_blank">https://petsc.org/release/faq/</a> for trouble shooting.</div>
<div>[5]PETSC ERROR: Petsc Release Version 3.21.0, unknown </div>
<div>[6]PETSC ERROR: WARNING! There are unused option(s) set! Could be the program crashed before usage or a spelling mistake, etc!</div>
<div>[6]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code</div>
<div>[6]PETSC ERROR: See <a href="https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$" target="_blank">https://petsc.org/release/faq/</a> for trouble shooting.</div>
<div>[6]PETSC ERROR: Petsc Release Version 3.21.0, unknown </div>
<div>[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!</div>
<div>[7]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code</div>
<div>[0]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code</div>
<div>[0]PETSC ERROR: See <a href="https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$" target="_blank">https://petsc.org/release/faq/</a> for trouble shooting.</div>
<div>[0]PETSC ERROR: Petsc Release Version 3.21.0, unknown </div>
<div>[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</div>
<div>[1]PETSC ERROR: Option left: name:-sns_monitor (no value) source: code</div>
<div>[1]PETSC ERROR: See <a href="https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$" target="_blank">https://petsc.org/release/faq/</a> for trouble shooting.</div>
<div>See <a href="https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$" target="_blank">https://petsc.org/release/faq/</a> for trouble shooting.</div>
<div>[2]PETSC ERROR: Petsc Release Version 3.21.0, unknown </div>
<div>[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</div>
<div>[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</div>
<div>[3]PETSC ERROR: Petsc Release Version 3.21.0, unknown </div>
<div>[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</div>
<div>[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</div>
<div>[4]PETSC ERROR: Petsc Release Version 3.21.0, unknown </div>
<div>[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</div>
<div>[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</div>
<div>[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</div>
<div>[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</div>
<div>./exe-tasting-SOLERA on a arch-darwin-c-debug named mmp-laptop.local by migmolper Mon Jun 24 11:37:38 2024</div>
<div>[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</div>
<div>See <a href="https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!YCo3JjgVyVqwScfRf05FevGAOdJG2APEIkhxmpYQmmJFmnrOBYrCKeun20x6gytf7m2IpXZPhyib8zJmPUCYCw$" target="_blank">https://petsc.org/release/faq/</a> for trouble shooting.</div>
<div>[7]PETSC ERROR: Petsc Release Version 3.21.0, unknown </div>
<div>[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</div>
<div>[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</div>
<div>Configure options --download-hdf5=1 --download-mpich=1 --with-debugging=1 CC=gcc CXX=c++ PETSC_ARCH=arch-darwin-c-debug --with-x</div>
<div>[1]PETSC ERROR: Petsc Release Version 3.21.0, unknown </div>
<div>[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</div>
<div>[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</div>
<div>[0]PETSC ERROR: #1 PetscSplitReductionApply() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230</div>
<div>[2]PETSC ERROR: #1 PetscSplitReductionApply() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230</div>
<div>[2]PETSC ERROR: [3]PETSC ERROR: #1 PetscSplitReductionApply() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230</div>
<div>[3]PETSC ERROR: #2 PetscSplitReductionEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172</div>
<div>[4]PETSC ERROR: #1 PetscSplitReductionApply() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230</div>
<div>[4]PETSC ERROR: #2 PetscSplitReductionEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172</div>
<div>[5]PETSC ERROR: #1 PetscSplitReductionApply() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230</div>
<div>[5]PETSC ERROR: #2 PetscSplitReductionEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172</div>
<div>[0]PETSC ERROR: #2 PetscSplitReductionEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172</div>
<div>[0]PETSC ERROR: #3 VecNormEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553</div>
<div>[0]PETSC ERROR: #4 SNESLineSearchApply_BT() at /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88</div>
<div>[1]PETSC ERROR: #1 PetscSplitReductionApply() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230</div>
<div>[1]PETSC ERROR: #2 PetscSplitReductionEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172</div>
<div>[1]PETSC ERROR: #3 VecNormEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553</div>
<div>#2 PetscSplitReductionEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172</div>
<div>[2]PETSC ERROR: #3 VecNormEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553</div>
<div>[2]PETSC ERROR: #4 SNESLineSearchApply_BT() at /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88</div>
<div>[2]PETSC ERROR: [3]PETSC ERROR: #3 VecNormEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553</div>
<div>[3]PETSC ERROR: #4 SNESLineSearchApply_BT() at /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88</div>
<div>[3]PETSC ERROR: #5 SNESLineSearchApply() at /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645</div>
<div>[4]PETSC ERROR: #3 VecNormEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553</div>
<div>[4]PETSC ERROR: #4 SNESLineSearchApply_BT() at /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88</div>
<div>[4]PETSC ERROR: #5 SNESLineSearchApply() at /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645</div>
<div>[5]PETSC ERROR: #3 VecNormEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553</div>
<div>[5]PETSC ERROR: #4 SNESLineSearchApply_BT() at /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88</div>
<div>[5]PETSC ERROR: #5 SNESLineSearchApply() at /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645</div>
<div>[5]PETSC ERROR: [6]PETSC ERROR: #1 PetscSplitReductionApply() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:230</div>
<div>[6]PETSC ERROR: #2 PetscSplitReductionEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:172</div>
<div>[6]PETSC ERROR: #3 VecNormEnd() at /Users/migmolper/petsc/src/vec/vec/utils/comb.c:553</div>
<div>[6]PETSC ERROR: #4 SNESLineSearchApply_BT() at /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88</div>
<div>[6]PETSC ERROR: #5 SNESLineSearchApply() at /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645</div>
<div>[7]PETSC ERROR: #1 VecXDot_MPI_Default() at /Users/migmolper/petsc/include/../src/vec/vec/impls/mpi/pvecimpl.h:107</div>
<div>[7]PETSC ERROR: #2 VecDot_MPI() at /Users/migmolper/petsc/src/vec/vec/impls/mpi/pvec2.c:10</div>
<div>[7]PETSC ERROR: #3 VecDot() at /Users/migmolper/petsc/src/vec/vec/interface/rvector.c:120</div>
<div>[7]PETSC ERROR: #4 SNESLineSearchApply_CP() at /Users/migmolper/petsc/src/snes/linesearch/impls/cp/linesearchcp.c:28</div>
<div>[0]PETSC ERROR: #5 SNESLineSearchApply() at /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645</div>
<div>[0]PETSC ERROR: #6 SNESSolve_NEWTONLS() at /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234</div>
<div>[0]PETSC ERROR: #7 SNESSolve() at /Users/migmolper/petsc/src/snes/interface/snes.c:4738</div>
<div>[0]PETSC ERROR: #8 TSTheta_SNESSolve() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174</div>
<div>[1]PETSC ERROR: #4 SNESLineSearchApply_BT() at /Users/migmolper/petsc/src/snes/linesearch/impls/bt/linesearchbt.c:88</div>
<div>[1]PETSC ERROR: #5 SNESLineSearchApply() at /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645</div>
<div>[1]PETSC ERROR: #6 SNESSolve_NEWTONLS() at /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234</div>
<div>[1]PETSC ERROR: #7 SNESSolve() at /Users/migmolper/petsc/src/snes/interface/snes.c:4738</div>
<div>[1]PETSC ERROR: #8 TSTheta_SNESSolve() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174</div>
<div>[1]PETSC ERROR: #9 TSStep_Theta() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225</div>
<div>[1]PETSC ERROR: #10 TSStep() at /Users/migmolper/petsc/src/ts/interface/ts.c:3391</div>
<div>[1]PETSC ERROR: #11 TSSolve() at /Users/migmolper/petsc/src/ts/interface/ts.c:4037</div>
<div>#5 SNESLineSearchApply() at /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645</div>
<div>[2]PETSC ERROR: #6 SNESSolve_NEWTONLS() at /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234</div>
<div>[2]PETSC ERROR: #7 SNESSolve() at /Users/migmolper/petsc/src/snes/interface/snes.c:4738</div>
<div>[2]PETSC ERROR: #8 TSTheta_SNESSolve() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174</div>
<div>[2]PETSC ERROR: #9 TSStep_Theta() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225</div>
<div>[2]PETSC ERROR: #10 TSStep() at /Users/migmolper/petsc/src/ts/interface/ts.c:3391</div>
<div>[2]PETSC ERROR: #11 TSSolve() at /Users/migmolper/petsc/src/ts/interface/ts.c:4037</div>
<div>[3]PETSC ERROR: #6 SNESSolve_NEWTONLS() at /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234</div>
<div>[3]PETSC ERROR: #7 SNESSolve() at /Users/migmolper/petsc/src/snes/interface/snes.c:4738</div>
<div>[3]PETSC ERROR: #8 TSTheta_SNESSolve() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174</div>
<div>[3]PETSC ERROR: #9 TSStep_Theta() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225</div>
<div>[3]PETSC ERROR: #10 TSStep() at /Users/migmolper/petsc/src/ts/interface/ts.c:3391</div>
<div>[3]PETSC ERROR: #11 TSSolve() at /Users/migmolper/petsc/src/ts/interface/ts.c:4037</div>
<div>[4]PETSC ERROR: #6 SNESSolve_NEWTONLS() at /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234</div>
<div>[4]PETSC ERROR: #7 SNESSolve() at /Users/migmolper/petsc/src/snes/interface/snes.c:4738</div>
<div>[4]PETSC ERROR: #8 TSTheta_SNESSolve() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174</div>
<div>[4]PETSC ERROR: #9 TSStep_Theta() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225</div>
<div>[4]PETSC ERROR: #10 TSStep() at /Users/migmolper/petsc/src/ts/interface/ts.c:3391</div>
<div>[4]PETSC ERROR: #11 TSSolve() at /Users/migmolper/petsc/src/ts/interface/ts.c:4037</div>
<div>#6 SNESSolve_NEWTONLS() at /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234</div>
<div>[5]PETSC ERROR: #7 SNESSolve() at /Users/migmolper/petsc/src/snes/interface/snes.c:4738</div>
<div>[5]PETSC ERROR: #8 TSTheta_SNESSolve() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174</div>
<div>[5]PETSC ERROR: #9 TSStep_Theta() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225</div>
<div>[5]PETSC ERROR: #10 TSStep() at /Users/migmolper/petsc/src/ts/interface/ts.c:3391</div>
<div>[5]PETSC ERROR: #11 TSSolve() at /Users/migmolper/petsc/src/ts/interface/ts.c:4037</div>
<div>[6]PETSC ERROR: #6 SNESSolve_NEWTONLS() at /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234</div>
<div>[6]PETSC ERROR: #7 SNESSolve() at /Users/migmolper/petsc/src/snes/interface/snes.c:4738</div>
<div>[6]PETSC ERROR: #8 TSTheta_SNESSolve() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174</div>
<div>[6]PETSC ERROR: #9 TSStep_Theta() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225</div>
<div>[6]PETSC ERROR: #10 TSStep() at /Users/migmolper/petsc/src/ts/interface/ts.c:3391</div>
<div>[6]PETSC ERROR: #11 TSSolve() at /Users/migmolper/petsc/src/ts/interface/ts.c:4037</div>
<div>[7]PETSC ERROR: #5 SNESLineSearchApply() at /Users/migmolper/petsc/src/snes/linesearch/interface/linesearch.c:645</div>
<div>[7]PETSC ERROR: #6 SNESSolve_NEWTONLS() at /Users/migmolper/petsc/src/snes/impls/ls/ls.c:234</div>
<div>[7]PETSC ERROR: #7 SNESSolve() at /Users/migmolper/petsc/src/snes/interface/snes.c:4738</div>
<div>[7]PETSC ERROR: #8 TSTheta_SNESSolve() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:174</div>
<div>[7]PETSC ERROR: #9 TSStep_Theta() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225</div>
<div>[7]PETSC ERROR: #10 TSStep() at /Users/migmolper/petsc/src/ts/interface/ts.c:3391</div>
<div>[7]PETSC ERROR: #11 TSSolve() at /Users/migmolper/petsc/src/ts/interface/ts.c:4037</div>
<div>[0]PETSC ERROR: #9 TSStep_Theta() at /Users/migmolper/petsc/src/ts/impls/implicit/theta/theta.c:225</div>
<div>[0]PETSC ERROR: #10 TSStep() at /Users/migmolper/petsc/src/ts/interface/ts.c:3391</div>
<div>[0]PETSC ERROR: #11 TSSolve() at /Users/migmolper/petsc/src/ts/interface/ts.c:4037</div>
<div>[0]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251</div>
<div>[2]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251</div>
<div>[3]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251</div>
<div>[1]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251</div>
<div>[4]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251</div>
<div>[5]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251</div>
<div>[6]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251</div>
<div>[7]PETSC ERROR: #12 Mass_Transport_Master_Equation_PETSc() at /Users/migmolper/DMD/SOLERA/Chemical-eqs/Mass-Transport-PETSc.cpp:251</div>
</div>
</div>
</div></blockquote></div><br clear="all"><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="https://urldefense.us/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!YlgDM4WJ1LZd1CuLZGNlMThj7OBlwkCu03P9vWxFt-dOODaCkEboruflECm3ntB9Grn1TGgk5Ic8lXZDmbEm$" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>