<div dir="ltr"><div dir="ltr">On Fri, Jan 13, 2023 at 11:22 AM Alexis Marboeuf <<a href="mailto:alexis.marboeuf@hotmail.fr">alexis.marboeuf@hotmail.fr</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="msg5485660062342652653">
<div dir="ltr">
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">Hi all,</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)"><br>
</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)"><span style="background-color:rgb(255,255,255);display:inline">In
a variational approach of brittle fracture setting, </span>I try to solve a bound constraint minimization problem using TAO. I checkout on the main branch of Petsc. Minimization with respect to the bounded variable (damage) is achieved through the Bounded
Newton Trust Region (TAOBNTR). All other TAO parameters are set by default. On a Linux machine, I get the following error with a 4 processors run:</span></div></div></div></blockquote><div><br></div><div>Can you view the solver?</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </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="msg5485660062342652653"><div dir="ltr">
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)"><br>
</span></div>
<div><span style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
<div style="font-weight:normal;color:rgb(212,212,212);background-color:rgb(30,30,30);font-family:Menlo,Monaco,"Courier New",monospace;font-size:12px;line-height:18px">
<span><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: ---------------------
</span><span style="color:rgb(206,145,120);font-weight:bold">Error</span><span> Message --------------------------------------------------------------</span></span>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Nonconforming object sizes</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Preconditioner number of local rows
</span><span style="color:rgb(86,156,214)">1122</span><span> does not equal input vector size
</span><span style="color:rgb(86,156,214)">1161</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: See
</span><span style="color:rgb(86,156,214)"><a href="https://petsc.org/release/faq/" target="_blank">https://petsc.org/release/faq/</a></span><span> for trouble shooting.</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Petsc Development GIT revision:
</span><span style="color:rgb(86,156,214)">v3.18.3-342-gdab44c92d91</span><span> GIT Date:
</span><span style="color:rgb(106,153,85)">2023-01-04</span><span> </span><span style="color:rgb(106,153,85)">13:37:04</span><span> +</span><span style="color:rgb(86,156,214)">0000</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: /home/marboeua/Developpement/mef90/arch-darwin-c/bin/vDefTAO on a arch-darwin-c named bb01 by marboeua Thu Jan
</span><span style="color:rgb(86,156,214)">12</span><span> </span><span style="color:rgb(106,153,85)">16:55:18</span><span> </span><span style="color:rgb(86,156,214)">2023</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: ---------------------
</span><span style="color:rgb(206,145,120);font-weight:bold">Error</span><span> Message --------------------------------------------------------------</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Configure options --FFLAGS=-ffree-line-length-none --COPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3
-march=znver3 -g"</span><span> --CXXOPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3 -march=znver3 -g"</span><span> --FOPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3
-march=znver3 -g"</span><span> --download-fblaslapack=</span><span style="color:rgb(86,156,214)">1</span><span> --download-mumps=</span><span style="color:rgb(86,156,214)">1</span><span> --download-chaco=</span><span style="color:rgb(86,156,214)">1</span><span> --download-exodusii=</span><span style="color:rgb(86,156,214)">1</span><span> --download-hypre=</span><span style="color:rgb(86,156,214)">1</span><span> --download-ml=</span><span style="color:rgb(86,156,214)">1</span><span> --download-triangle
--download-scalapack=</span><span style="color:rgb(86,156,214)">1</span><span> --download-superlu=</span><span style="color:rgb(86,156,214)">1</span><span> --download-sowing=</span><span style="color:rgb(86,156,214)">1</span><span> --download-sowing-cc=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/gcc
--download-sowing-cxx=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/g++ --download-sowing-cpp=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/cpp
--download-sowing-cxxcpp=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/cpp --download-yaml=</span><span style="color:rgb(86,156,214)">1</span><span> --download-bison=</span><span style="color:rgb(86,156,214)">1</span><span> --download-hdf5=</span><span style="color:rgb(86,156,214)">1</span><span> --download-metis=</span><span style="color:rgb(86,156,214)">1</span><span> --download-parmetis=</span><span style="color:rgb(86,156,214)">1</span><span> --download-netcdf=</span><span style="color:rgb(86,156,214)">1</span><span> --download-pnetcdf=</span><span style="color:rgb(86,156,214)">1</span><span> --download-zlib=</span><span style="color:rgb(86,156,214)">1</span><span> --with-cmake=</span><span style="color:rgb(86,156,214)">1</span><span> --with-debugging=</span><span style="color:rgb(86,156,214)">0</span><span> --with-mpi-dir=/opt/HPC/mvapich2/</span><span style="color:rgb(86,156,214)">2</span><span>.</span><span style="color:rgb(86,156,214)">3</span><span>.</span><span style="color:rgb(86,156,214)">7</span><span>-</span><span style="color:rgb(86,156,214)">gcc11.2.1</span><span> --with-ranlib=ranlib
--with-shared-libraries=</span><span style="color:rgb(86,156,214)">1</span><span> --with-sieve=</span><span style="color:rgb(86,156,214)">1</span><span> --download-p4est=</span><span style="color:rgb(86,156,214)">1</span><span> --with-pic
--with-mpiexec=srun --with-x11=</span><span style="color:rgb(86,156,214)">0</span><span> PETSC_ARCH=arch-darwin-c</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">1</span><span> PCApply() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/pc/interface/precon.c:</span><span style="color:rgb(86,156,214)">434</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">2</span><span> KSP_PCApply() at /home/marboeua/Developpement/petsc/include/petsc/private/kspimpl.h:</span><span style="color:rgb(86,156,214)">380</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">3</span><span> KSPCGSolve_STCG() at
/</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/impls/cg/stcg/stcg.c:</span><span style="color:rgb(86,156,214)">76</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">4</span><span> KSPSolve_Private() at
/</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:</span><span style="color:rgb(86,156,214)">898</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">5</span><span> KSPSolve() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:</span><span style="color:rgb(86,156,214)">1070</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">6</span><span> TaoBNKComputeStep() at
/</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bnk.c:</span><span style="color:rgb(86,156,214)">459</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">7</span><span> TaoSolve_BNTR() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bntr.c:</span><span style="color:rgb(86,156,214)">138</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">8</span><span> TaoSolve() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/interface/taosolver.c:</span><span style="color:rgb(86,156,214)">177</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Nonconforming object sizes</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Preconditioner number of local rows
</span><span style="color:rgb(86,156,214)">1229</span><span> does not equal input vector size
</span><span style="color:rgb(86,156,214)">1254</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: See
</span><span style="color:rgb(86,156,214)"><a href="https://petsc.org/release/faq/" target="_blank">https://petsc.org/release/faq/</a></span><span> for trouble shooting.</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Petsc Development GIT revision:
</span><span style="color:rgb(86,156,214)">v3.18.3-342-gdab44c92d91</span><span> GIT Date:
</span><span style="color:rgb(106,153,85)">2023-01-04</span><span> </span><span style="color:rgb(106,153,85)">13:37:04</span><span> +</span><span style="color:rgb(86,156,214)">0000</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: /home/marboeua/Developpement/mef90/arch-darwin-c/bin/vDefTAO on a arch-darwin-c named bb01 by marboeua Thu Jan
</span><span style="color:rgb(86,156,214)">12</span><span> </span><span style="color:rgb(106,153,85)">16:55:18</span><span> </span><span style="color:rgb(86,156,214)">2023</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: Configure options --FFLAGS=-ffree-line-length-none --COPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3
-march=znver3 -g"</span><span> --CXXOPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3 -march=znver3 -g"</span><span> --FOPTFLAGS=</span><span style="color:rgb(206,145,120)">"-O3
-march=znver3 -g"</span><span> --download-fblaslapack=</span><span style="color:rgb(86,156,214)">1</span><span> --download-mumps=</span><span style="color:rgb(86,156,214)">1</span><span> --download-chaco=</span><span style="color:rgb(86,156,214)">1</span><span> --download-exodusii=</span><span style="color:rgb(86,156,214)">1</span><span> --download-hypre=</span><span style="color:rgb(86,156,214)">1</span><span> --download-ml=</span><span style="color:rgb(86,156,214)">1</span><span> --download-triangle
--download-scalapack=</span><span style="color:rgb(86,156,214)">1</span><span> --download-superlu=</span><span style="color:rgb(86,156,214)">1</span><span> --download-sowing=</span><span style="color:rgb(86,156,214)">1</span><span> --download-sowing-cc=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/gcc
--download-sowing-cxx=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/g++ --download-sowing-cpp=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/cpp
--download-sowing-cxxcpp=/opt/rh/devtoolset-</span><span style="color:rgb(86,156,214)">9</span><span>/root/usr/bin/cpp --download-yaml=</span><span style="color:rgb(86,156,214)">1</span><span> --download-bison=</span><span style="color:rgb(86,156,214)">1</span><span> --download-hdf5=</span><span style="color:rgb(86,156,214)">1</span><span> --download-metis=</span><span style="color:rgb(86,156,214)">1</span><span> --download-parmetis=</span><span style="color:rgb(86,156,214)">1</span><span> --download-netcdf=</span><span style="color:rgb(86,156,214)">1</span><span> --download-pnetcdf=</span><span style="color:rgb(86,156,214)">1</span><span> --download-zlib=</span><span style="color:rgb(86,156,214)">1</span><span> --with-cmake=</span><span style="color:rgb(86,156,214)">1</span><span> --with-debugging=</span><span style="color:rgb(86,156,214)">0</span><span> --with-mpi-dir=/opt/HPC/mvapich2/</span><span style="color:rgb(86,156,214)">2</span><span>.</span><span style="color:rgb(86,156,214)">3</span><span>.</span><span style="color:rgb(86,156,214)">7</span><span>-</span><span style="color:rgb(86,156,214)">gcc11.2.1</span><span> --with-ranlib=ranlib
--with-shared-libraries=</span><span style="color:rgb(86,156,214)">1</span><span> --with-sieve=</span><span style="color:rgb(86,156,214)">1</span><span> --download-p4est=</span><span style="color:rgb(86,156,214)">1</span><span> --with-pic
--with-mpiexec=srun --with-x11=</span><span style="color:rgb(86,156,214)">0</span><span> PETSC_ARCH=arch-darwin-c</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">1</span><span> PCApply() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/pc/interface/precon.c:</span><span style="color:rgb(86,156,214)">434</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">2</span><span> KSP_PCApply() at /home/marboeua/Developpement/petsc/include/petsc/private/kspimpl.h:</span><span style="color:rgb(86,156,214)">380</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">3</span><span> KSPCGSolve_STCG() at
/</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/impls/cg/stcg/stcg.c:</span><span style="color:rgb(86,156,214)">76</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">4</span><span> KSPSolve_Private() at
/</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:</span><span style="color:rgb(86,156,214)">898</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">5</span><span> KSPSolve() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/ksp/ksp/interface/itfunc.c:</span><span style="color:rgb(86,156,214)">1070</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">6</span><span> TaoBNKComputeStep() at
/</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bnk.c:</span><span style="color:rgb(86,156,214)">459</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">7</span><span> TaoSolve_BNTR() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/bound/impls/bnk/bntr.c:</span><span style="color:rgb(86,156,214)">138</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">8</span><span> TaoSolve() at /</span><span style="color:rgb(86,156,214)">1</span><span>/home/marboeua/Developpement/petsc/src/tao/interface/taosolver.c:</span><span style="color:rgb(86,156,214)">177</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">3</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">9</span><span> /home/marboeua/Developpement/mef90/vDef/vDefTAO.</span><span style="color:rgb(86,156,214)">F90:370</span></div>
<div><span>application called MPI_Abort(MPI_COMM_SELF, </span>
<span style="color:rgb(86,156,214)">60</span><span>) - process
</span><span style="color:rgb(86,156,214)">0</span></div>
<div><span>[</span><span style="color:rgb(86,156,214)">2</span><span>]PETSC
</span><span style="color:rgb(206,145,120);font-weight:bold">ERROR</span><span>: #</span><span style="color:rgb(86,156,214)">9</span><span> /home/marboeua/Developpement/mef90/vDef/vDefTAO.</span><span style="color:rgb(86,156,214)">F90:370</span></div>
<div><span>application called MPI_Abort(MPI_COMM_SELF, </span>
<span style="color:rgb(86,156,214)">60</span><span>) - process
</span><span style="color:rgb(86,156,214)">0</span></div>
<div><span>slurmstepd: </span><span style="color:rgb(206,145,120);font-weight:bold">error:</span><span> *** STEP
</span><span style="color:rgb(86,156,214)">5034</span><span>.</span><span style="color:rgb(86,156,214)">0</span><span> ON bb01 CANCELLED AT
</span><span style="color:rgb(106,153,85)">2023-01-12T17:21:07</span><span> ***</span></div>
<div><span>srun: Job step aborted: Waiting up to </span><span style="color:rgb(86,156,214)">32</span><span> seconds for job step to finish.</span></div>
<div><span>srun: </span><span style="color:rgb(206,145,120);font-weight:bold">error:</span><span> bb01: tasks
</span><span style="color:rgb(86,156,214)">0</span><span>-</span><span style="color:rgb(86,156,214)">1</span><span>: Killed</span></div>
<span><span>srun: </span><span style="color:rgb(206,145,120);font-weight:bold">error:</span><span> bb01: tasks
</span><span style="color:rgb(86,156,214)">2</span><span>-</span><span style="color:rgb(86,156,214)">3</span><span>: Exited with exit code
</span><span style="color:rgb(86,156,214)">1</span></span></div>
<br>
</span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
The error is raised in the middle of the computation after many successful calls of TAOSolve and TAO iterations. My guess is that TAO computes the preconditioner during its first iteration with all variables in the active set. But the preconditioner is never
updated when some variables are moved to the inactive set during the next TAO iterations. Am I right? Can you help me with that?</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
Thanks a lot for your help and your time.</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
Regards,</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
Alexis</div>
</div>
</div></blockquote></div><br clear="all"><div><br></div>-- <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="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>