<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sun, Jul 31, 2016 at 4:56 PM, Stefano Zampini <span dir="ltr"><<a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div style="margin:0px">It’s the first time I’m trying to setup field split so maybe I’m missing something</div><div style="margin:0px"><br></div><div style="margin:0px">I have a MatNest with NULL (1,1) block. I setup the field split using the row ISs from the MatNest (a note aside: why they cannot be inferred by FieldSplit if pc->pmat is of type MatNEST?)</div><div style="margin:0px">Options used are -pc_fieldsplit_type additive and -pc_fieldsplit_detect_saddle_point</div></div></blockquote><div><br></div><div>You don't need ISes is you use detect_saddle_point, so use one or the other.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div style="margin:0px">I then get this error, because MatGetSubMatrix_Nest returns a NULL pointer.</div><div style="margin:0px">What is the solution? having MatGetSubMatrix_Nest return a zero-valued matrix of correct size? any command line option I’m not using? Or instead having Fieldsplit handle the case?</div></div></blockquote><div><br></div><div>I think have Nest return a real matrix is the right solution, since this is what will happen with other types.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div style="margin:0px;font-size:14px;font-family:Menlo">[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div><div style="margin:0px;font-size:14px;font-family:Menlo">[0]PETSC ERROR: Null argument, when expecting valid pointer</div><div style="margin:0px;font-size:14px;font-family:Menlo">[0]PETSC ERROR: Null Object: Parameter # 1</div><div style="margin:0px;font-size:14px;font-family:Menlo">[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</div><div style="margin:0px;font-size:14px;font-family:Menlo">[0]PETSC ERROR: Petsc Development GIT revision: v3.7.2-744-gb54d608  GIT Date: 2016-06-28 18:20:20 +0300</div><div style="margin:0px;font-size:14px;font-family:Menlo">[0]PETSC ERROR: Unknown Name on a debug named kl-13916.local by szampini Sun Jul 31 14:47:53 2016</div><div style="margin:0px;font-size:14px;font-family:Menlo">[0]PETSC ERROR: Configure options --download-chaco --download-ctetgen --download-exodusii --download-hdf5 --download-hypre --download-metis --download-mpich --download-mumps --download-netcdf --download-p4est --download-parmetis --download-ptscotch --download-scalapack --download-suitesparse --download-triangle --with-valgrind-dir=/usr/local/ CC=/usr/local/bin/gcc-5 CXX=/usr/local/bin/g++-5 FC=/usr/local/bin/gfortran-5 PETSC_ARCH=debug</div><div style="margin:0px;font-size:14px;font-family:Menlo">[0]PETSC ERROR: #1 MatSetOptionsPrefix() line 654 in /Users/szampini/software/petsc/src/mat/interface/matrix.c</div><div style="margin:0px;font-size:14px;font-family:Menlo">[0]PETSC ERROR: #2 PCSetUp_FieldSplit() line 558 in /Users/szampini/software/petsc/src/ksp/pc/impls/fieldsplit/fieldsplit.c</div><div style="margin:0px;font-size:14px;font-family:Menlo">[0]PETSC ERROR: #3 PCSetUp() line 968 in /Users/szampini/software/petsc/src/ksp/pc/interface/precon.c</div><div><br></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">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></div>