<div dir="ltr"><div><div><div><div>Barry:<br><br></div>Thanks for the quick response!<br><br></div>I assumed that the Mat object created by MatCreateSchurComplement can somehow be solved by a KSP option (like plain CG). <br><br>I am aware of the fieldsplit pc. My situation is that I have to solve a problem with S = L - D K^-1 D^T, wherein L, D, K  are all sparse. Of course, the last option for me is to write a simple CG solver for S myself. Do you have any recommendations?<br><br></div>Thanks,<br><br></div>Luis<br></div><div class="gmail_extra"><br><div class="gmail_quote">2018-01-23 11:25 GMT-08:00 Smith, Barry F. <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
  Luis,<br>
<br>
   This doesn't really have any thing to do with the Schur complement per se, you are asking do an ILU on this abstract Schur complement matrix which is not possible.<br>
<br>
   Almost for sure you want to use PCFIELDSPLIT and not directly create the MatCreateSchurComplement() yourself. PCFIELDSPLIT has many possible variants for solving systems with segregated Schur complement methods.<br>
<span class="HOEnZb"><font color="#888888"><br>
   Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
> On Jan 23, 2018, at 1:13 PM, saturday luis <<a href="mailto:luis.saturday@gmail.com">luis.saturday@gmail.com</a>> wrote:<br>
><br>
> Hi PETSc team:<br>
><br>
> I am trying to solve a Stokes problem in a segregated manner, in which the pressure field is obtained by solving the pressure Schur complement. This is what I did.<br>
><br>
> I created a Schur complement by calling<br>
><br>
> Mat S;<br>
> MatCreateSchurComplement(<wbr>gloAssem_ptr->K_11, gloAssem_ptr->K_11, gloAssem_ptr->K_10, gloAssem_ptr->K_01, gloAssem_ptr->K_00, &S);<br>
><br>
> Then I use the Mat S as an input for KSP:<br>
> KSPSetOperators(ksp, S, S);<br>
><br>
> I got error message as follows.<br>
> [0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/linearsolvertable.html" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>documentation/<wbr>linearsolvertable.html</a> for possible LU and Cholesky solvers<br>
> [0]PETSC ERROR: Could not locate a solver package. Perhaps you must ./configure with --download-<package><br>
> [0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>documentation/faq.html</a> for trouble shooting.<br>
> [0]PETSC ERROR: Petsc Release Version 3.8.3, Dec, 09, 2017<br>
> [0]PETSC ERROR: ./bipn3d on a arch-linux2-cxx-debug named bacon by luis Tue Jan 23 11:02:34 2018<br>
> [0]PETSC ERROR: Configure options --prefix=~/lib/petsc-3.8.3-<wbr>debug --with-x=0 -with-pic --with-mpi-compilers=1 --with-mpi-dir=/home/lib/<wbr>mpich-3.2/ --with-clanguage=C++ --with-scalar-type=real --with-shared-libraries=1 --with-precision=double --with-hypre=1 --download-hypre --with-spai=1 --download-spai --with-superlu_dist=1 --download-superlu_dist --with-superlu=1 --download-superlu --with-parmetis=1 --download-parmetis --with-metis=1 --download-metis --with-mumps=1 --download-mumps --with-parmetis=1 --download-parmetis --with-metis=1 --download-metis --with-scalapack=1 --download-scalapack --with-blacs=1 --download-blacs --with-spooles=1 --download-spooles --with-debugging=yes --download-fblaslapack<br>
> [0]PETSC ERROR: #1 MatGetFactor() line 4346 in /home/lib/petsc-3.8.3/src/mat/<wbr>interface/matrix.c<br>
> [0]PETSC ERROR: #2 PCSetUp_ILU() line 142 in /home/lib/petsc-3.8.3/src/ksp/<wbr>pc/impls/factor/ilu/ilu.c<br>
> [0]PETSC ERROR: #3 PCSetUp() line 924 in /home/lib/petsc-3.8.3/src/ksp/<wbr>pc/interface/precon.c<br>
> [0]PETSC ERROR: #4 KSPSetUp() line 381 in /home/lib/petsc-3.8.3/src/ksp/<wbr>ksp/interface/itfunc.c<br>
> [0]PETSC ERROR: #5 KSPSolve() line 612 in /home/lib/petsc-3.8.3/src/ksp/<wbr>ksp/interface/itfunc.c<br>
><br>
> Let me know if you can see any issue in the usage of Schur complement creation.<br>
><br>
> Thanks,<br>
><br>
> Luis<br>
<br>
</div></div></blockquote></div><br></div>