<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Dave, and other PETSc users, sorry for bothering you all again. However continuing this topic, I have another question. Within this ex70.c, there is a user defined schur complement StokesSetupApproxSchur,
defined as
<div><span style="color: rgb(0, 132, 0); font-family: Menlo; font-size: 11px;"><br>
</span></div>
<div>
<p style="margin-right: 0px; margin-left: 0px; font-size: 11px; font-family: Menlo; color: rgb(0, 132, 0);">
/* Schur complement approximation: myS = A11 - A10 diag(A00)^(-1) A01 */</p>
<div><br>
</div>
<div>The preconditioner here can be written out explicitly as a matrix, and is not changed between iterations. Each time it is applied to the residue is simply matrix times vector. However, if I want to define a preconditioner which is allowed to vary between
iterations, for example, if I need to solve a linear system each time I apply the preconditioner to the residue, then is there a way I can do this? Any hint will be greatly appreciated. Thank you very much! </div>
<div><br>
</div>
<div>Best,</div>
<div>Hui</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF670824" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> Dave May [dave.mayhem23@gmail.com]<br>
<b>Sent:</b> Wednesday, July 16, 2014 3:40 PM<br>
<b>To:</b> Sun, Hui<br>
<b>Cc:</b> Matthew Knepley; petsc-users@mcs.anl.gov<br>
<b>Subject:</b> Re: [petsc-users] why a certain option cannot be used<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">You need to specify the fieldsplit preconditioning type. The choices are
<br>
<br>
-pc_fieldsplit_type <MULTIPLICATIVE> (choose one of) ADDITIVE MULTIPLICATIVE SYMMETRIC_MULTIPLICATIVE SPECIAL SCHUR<br>
<br>
Maybe what you want to run is the following<br>
<br>
./ex70 -ksp_view -pc_type fieldsplit -fieldsplit_0_pc_type jacobi -fieldsplit_1_pc_type jacobi -pc_fieldsplit_type SCHUR -pc_fieldsplit_schur_precondition user<br>
<div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">The option " -fieldsplit_1_user_pc" seems meaningless to me.<br>
I don't see anything in ex70.c which would cause this option to have any effect. The lines<br>
<br>
ierr = PetscOptionsHasName(NULL, "-user_pc", &s.userPC);CHKERRQ(ierr);<br>
ierr = PetscOptionsHasName(NULL, "-user_ksp", &s.userKSP);CHKERRQ(ierr);<br>
<br>
</div>
<div class="gmail_extra">imply that the option -user_pc and -user_ksp are meaningful and modify how the schur complement preconditioner will be configured. The first option triggers this setup to occur;<br>
<br>
if (s->userPC) {<br>
ierr = PCFieldSplitSchurPrecondition(pc, PC_FIELDSPLIT_SCHUR_PRE_USER, s->myS);CHKERRQ(ierr);<br>
}<br>
<br>
the second option will force this setup to occur<br>
<br>
if (s->userKSP) {<br>
ierr = PCSetUp(pc);CHKERRQ(ierr);<br>
ierr = PCFieldSplitGetSubKSP(pc, &n, &subksp);CHKERRQ(ierr);<br>
ierr = KSPSetOperators(subksp[1], s->myS, s->myS, SAME_PRECONDITIONER);CHKERRQ(ierr);<br>
ierr = PetscFree(subksp);CHKERRQ(ierr);<br>
}<br>
</div>
<div class="gmail_extra"><br>
<br>
<br>
<br>
<br>
<div class="gmail_quote">On 17 July 2014 00:24, Sun, Hui <span dir="ltr"><<a href="mailto:hus003@ucsd.edu" target="_blank">hus003@ucsd.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
Thank you Dave. It does give me the other options. However I want to be able to use the user defined fieldsplit_1_pc_type. In example ex70.c, there is a user defined Schur complement preconditioner. I want to use that. Then how should I set the values for the
options?
<div><br>
</div>
<div>I tried to see which values could be set by typing </div>
<div><span style="font-family:Menlo; font-size:11px"><br>
</span></div>
<div><span style="font-family:Menlo; font-size:11px">./ex70 -pc_type fieldsplit -help 2>/dev/null | grep -A5 fieldsplit_1_pc</span></div>
<div><span style="font-size:10pt"><br>
</span></div>
<div><span style="font-size:10pt"> </span><font>And</font> I get the following output:</div>
<div><font><br>
</font></div>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> -fieldsplit_1_pc_type <bjacobi>: Preconditioner (one of) none jacobi pbjacobi bjacobi sor lu shell mg</p>
<div class="">
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> eisenstat ilu icc cholesky asm gasm ksp composite redundant nn mat fieldsplit galerkin exotic hmpi asa cp lsc redistribute svd gamg tfs (PCSetType)</p>
</div>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> -fieldsplit_1_pc_use_amat: <FALSE> use Amat (instead of Pmat) to define preconditioner in nested inner solves (PCSetUseAmat)</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> Block Jacobi options</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> -fieldsplit_1_pc_bjacobi_blocks <-1>: Total number of blocks (PCBJacobiSetTotalBlocks)</p>
<div class="">
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">Krylov Method (KSP) options -------------------------------------------------</p>
</div>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> -fieldsplit_1_ksp_type <preonly>: Krylov method (one of) cg groppcg pipecg cgne nash stcg gltr richardson</p>
<div class="">
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> chebyshev gmres tcqmr bcgs ibcgs fbcgs fbcgsr bcgsl cgs tfqmr cr pipecr lsqr preonly qcg bicg fgmres minres symmlq lgmres lcd gcr pgmres specest dgmres (KSPSetType)</p>
</div>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> -fieldsplit_1_ksp_max_it <10000>: Maximum number of iterations (KSPSetTolerances)</p>
<div><span style="font-family:Menlo; font-size:11px"> </span><span style="font-family:Menlo; font-size:11px">-fieldsplit_1_ksp_rtol <1e-05>: Relative decrease in residual norm (KSPSetTolerances)</span></div>
<div><font><br>
</font></div>
<div><font>I cannot find anything here about the user defined pc. Do you have any idea what might be the choice? </font></div>
<div>
<div><br>
</div>
<div><br>
</div>
<div><br>
<div style="font-family:Times New Roman; color:rgb(0,0,0); font-size:16px">
<hr>
<div style="direction:ltr"><font color="#000000" face="Tahoma"><b>From:</b> Dave May [<a href="mailto:dave.mayhem23@gmail.com" target="_blank">dave.mayhem23@gmail.com</a>]<br>
<b>Sent:</b> Wednesday, July 16, 2014 3:11 PM
<div>
<div class="h5"><br>
<b>To:</b> Sun, Hui<br>
<b>Cc:</b> Matthew Knepley; <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">
petsc-users@mcs.anl.gov</a><br>
<b>Subject:</b> Re: [petsc-users] why a certain option cannot be used<br>
</div>
</div>
</font><br>
</div>
<div>
<div class="h5">
<div></div>
<div>
<div dir="ltr">
<div>
<div>To see all the fieldsplit options, run with<br>
</div>
./ex70 -pc_type fieldsplit -help<br>
<br>
</div>
-help will only display the options relevant to the current configuration of your KSP and PC.<br>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On 17 July 2014 00:10, Dave May <span dir="ltr"><<a href="mailto:dave.mayhem23@gmail.com" target="_blank">dave.mayhem23@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div dir="ltr">
<div>Run with <br>
</div>
-help<br>
</div>
<div>
<div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On 17 July 2014 00:07, Sun, Hui <span dir="ltr"><<a href="mailto:hus003@ucsd.edu" target="_blank">hus003@ucsd.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
Thank you Dave, it does sol<font>ve the problem. Another question is: since some options like</font>
<div><font><font face="Times New Roman"><br>
</font></font></div>
<div><font><font face="Times New Roman">-</font></font><font face="Times New Roman" style="font-size:small">fieldsplit_0_pc_type</font><span style="font-size:small"> </span></div>
<div><span style="font-size:small"><br>
</span></div>
<div><span style="font-size:small">does not show up if I type</span></div>
<div><font face="Times New Roman" style="font-size:small"><br>
</font></div>
<div><span style="font-size:small"></span><font face="Times New Roman" style="font-size:small">./ex70 -help</font></div>
<div><span style="font-size:small"><br>
</span></div>
<div><span style="font-size:small">how then can I have a complete list of all the options? </span></div>
<div>
<div><br>
</div>
<div>Best,</div>
<div>Hui</div>
<div><br>
</div>
<div><br>
<div style="font-family:Times New Roman; color:rgb(0,0,0); font-size:16px">
<hr>
<div style="direction:ltr"><font color="#000000" face="Tahoma"><b>From:</b> Dave May [<a href="mailto:dave.mayhem23@gmail.com" target="_blank">dave.mayhem23@gmail.com</a>]<br>
<b>Sent:</b> Wednesday, July 16, 2014 2:57 PM<br>
<b>To:</b> Sun, Hui<br>
<b>Cc:</b> Matthew Knepley; <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">
petsc-users@mcs.anl.gov</a>
<div>
<div><br>
<b>Subject:</b> Re: [petsc-users] why a certain option cannot be used<br>
</div>
</div>
</font><br>
</div>
<div>
<div>
<div></div>
<div>You need to specify on the command line
<div>-pc_type fieldsplit -fieldsplit_0_pc_type xxx -fieldsplit_1_pc_type yyy</div>
<div><br>
</div>
<div>Where xxx and yyy are the desired preconditioners for the (1,1) block and (2,2) block. Note that the default PC (ilu(0)) applied to both splits will fail for this example.</div>
<div><br>
</div>
<div>The comments in the code are misleading and incorrect. Without the option <font><span>-pc_type fieldsplit, all fieldsplit options have no effect.</span></font></div>
<div><font><span><br>
</span></font></div>
<div><font><span></span><span>When debugging solvers, always run with - ksp_view (as Matt recommends) </span></font></div>
<div><font><span><br>
</span></font></div>
<div><font><span><br>
</span></font><br>
On Wednesday, 16 July 2014, Sun, Hui <<a href="mailto:hus003@ucsd.edu" target="_blank">hus003@ucsd.edu</a>> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
No I don't think I'm using fieldsplitting although I want to use it. The output of
<div><span style="font-family:Menlo; font-size:11px"><br>
</span></div>
<div><span style="font-family:Menlo; font-size:11px">./ex70 -nx 32 -ny 48 -fieldsplit_1_pc_type jacobi -ksp_view</span></div>
<div><br>
</div>
<div>is the following: <span style="font-family:'Segoe UI',Helvetica,Arial,sans-serif; font-size:medium"> </span>
<div><br>
</div>
<div>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">KSP Object: 1 MPI processes</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> type: gmres</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> GMRES: happy breakdown tolerance 1e-30</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> maximum iterations=10000, initial guess is zero</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> tolerances: relative=1e-05, absolute=1e-50, divergence=10000</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> left preconditioning</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> using PRECONDITIONED norm type for convergence test</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">PC Object: 1 MPI processes</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> type: none</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> linear system matrix = precond matrix:</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> Matrix Object: 1 MPI processes</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> type: nest</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> rows=4608, cols=4608</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> Matrix object: </p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> type=nest, rows=2, cols=2 </p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> MatNest structure: </p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> (0,0) : prefix="a00_", type=mpiaij, rows=3072, cols=3072 </p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> (0,1) : prefix="a01_", type=mpiaij, rows=3072, cols=1536 </p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> (1,0) : prefix="a10_", type=mpiaij, rows=1536, cols=3072 </p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> (1,1) : prefix="a11_", type=mpiaij, rows=1536, cols=1536 </p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> residual u = 2.6315e-05</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> residual p = 0.000229872</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> residual [u,p] = 0.000231373</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error u = 0.00633503</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error p = 0.121534</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error [u,p] = 0.121699</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">WARNING! There are options you set that were not used!</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">WARNING! could be spelling mistake, etc!</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">Option left: name:-fieldsplit_1_pc_type value: jacobi</p>
<div><br>
</div>
<div><br>
<div style="font-family:Times New Roman; color:rgb(0,0,0); font-size:16px">
<hr>
<div style="direction:ltr"><font color="#000000" face="Tahoma"><b>From:</b> Sun, Hui<br>
<b>Sent:</b> Wednesday, July 16, 2014 2:40 PM<br>
<b>To:</b> Matthew Knepley<br>
<b>Cc:</b> <a href="http://UrlBlockedError.aspx" target="_blank">petsc-users@mcs.anl.gov</a><br>
<b>Subject:</b> RE: [petsc-users] why a certain option cannot be used<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
Thank you Matt. I've corrected it to
<div><br>
</div>
<div>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">./ex70 -nx 32 -ny 48 -fieldsplit_1_pc_type jacobi</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><br>
</p>
<p style="margin-right:0px; margin-left:0px">However the output is still:</p>
<p style="margin-right:0px; margin-left:0px"><br>
</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">residual u = 2.6315e-05</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> residual p = 0.000229872</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> residual [u,p] = 0.000231373</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error u = 0.00633503</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error p = 0.121534</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error [u,p] = 0.121699</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">WARNING! There are options you set that were not used!</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">WARNING! could be spelling mistake, etc!</p>
<p style="margin-right:0px; margin-left:0px"></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">Option left: name:-fieldsplit_1_pc_type value: jacobi</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><br>
</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px">If I type</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px"><br>
</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">./ex70 -help | grep fieldsplit_1_</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px"><br>
</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px">the output is that nothing coming out. </span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px"><br>
</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px"><br>
</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><br>
</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><br>
</p>
<div style="font-family:Times New Roman; color:rgb(0,0,0); font-size:16px">
<hr>
<div style="direction:ltr"><font color="#000000" face="Tahoma"><b>From:</b> Matthew Knepley [<a href="http://UrlBlockedError.aspx" target="_blank">knepley@gmail.com</a>]<br>
<b>Sent:</b> Wednesday, July 16, 2014 2:32 PM<br>
<b>To:</b> Sun, Hui<br>
<b>Cc:</b> <a href="http://UrlBlockedError.aspx" target="_blank">petsc-users@mcs.anl.gov</a><br>
<b>Subject:</b> Re: [petsc-users] why a certain option cannot be used<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">On Wed, Jul 16, 2014 at 4:25 PM, Matthew Knepley <span dir="ltr">
<<a href="http://UrlBlockedError.aspx" target="_blank">knepley@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">
<div>On Wed, Jul 16, 2014 at 4:22 PM, Sun, Hui <span dir="ltr"><<a href="http://UrlBlockedError.aspx" target="_blank">hus003@ucsd.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
Thank you Matt for answering, however even with that I don't think it works. For example, the command
<div><br>
</div>
<div>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">./ex70 -nx 32 -ny 48 -fieldsplit_1_type jacobi</p>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
</div>
<div>You have mistyped. Look carefully at my last message:</div>
<div><br>
</div>
<div> -fieldsplit_1_pc_type jacobi</div>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>I also want to note that you can see all the available options with -help. It does produce</div>
<div>a lot of output, but you can segregate it by prefix ("fieldsplit_1_"). You can also see the</div>
<div>prefix of each solver component using -ksp_view.</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 dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">
<div> Matt</div>
<div>
<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>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
<div>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px">gives me the following output</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px"><br>
</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">residual u = 2.6315e-05</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> residual p = 0.000229872</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> residual [u,p] = 0.000231373</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error u = 0.00633503</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error p = 0.121534</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error [u,p] = 0.121699</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">WARNING! There are options you set that were not used!</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">WARNING! could be spelling mistake, etc!</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">Option left: name:-fieldsplit_1_type value: jacobi</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><br>
</p>
<p style="margin-right:0px; margin-left:0px">Furthermore, if I look into the possible options with keyword fieldsplit involved by the command:</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><br>
</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">./ex70 -help | grep -A5 fieldsplit</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px"><br>
</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px">I get this output:</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px"><br>
</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">eisenstat ilu icc cholesky asm gasm ksp composite redundant nn mat fieldsplit galerkin exotic hmpi asa cp lsc redistribute svd gamg tfs (PCSetType)</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> -pc_use_amat: <FALSE> use Amat (instead of Pmat) to define preconditioner in nested inner solves (PCSetUseAmat)</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">Krylov Method (KSP) options -------------------------------------------------</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> -ksp_type <gmres>: Krylov method (one of) cg groppcg pipecg cgne nash stcg gltr richardson</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> chebyshev gmres tcqmr bcgs ibcgs fbcgs fbcgsr bcgsl cgs tfqmr cr pipecr lsqr preonly qcg bicg fgmres minres symmlq lgmres lcd gcr pgmres specest dgmres (KSPSetType)</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> -ksp_max_it <10000>: Maximum number of iterations (KSPSetTolerances)</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><br>
</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px">So I don't see any option that's similar to -fieldsplit_1_type. </span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px"><br>
</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px">Hui</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><span style="font-family:Tahoma; font-size:13px"><br>
</span></p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><br>
</p>
<div style="font-family:Times New Roman; color:rgb(0,0,0); font-size:16px">
<hr>
<div style="direction:ltr"><font color="#000000" face="Tahoma"><b>From:</b> Matthew Knepley [<a href="http://UrlBlockedError.aspx" target="_blank">knepley@gmail.com</a>]<br>
<b>Sent:</b> Wednesday, July 16, 2014 2:16 PM<br>
<b>To:</b> Sun, Hui<br>
<b>Cc:</b> <a href="http://UrlBlockedError.aspx" target="_blank">petsc-users@mcs.anl.gov</a><br>
<b>Subject:</b> Re: [petsc-users] why a certain option cannot be used<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">On Wed, Jul 16, 2014 at 4:09 PM, Sun, Hui <span dir="ltr">
<<a href="http://UrlBlockedError.aspx" target="_blank">hus003@ucsd.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
<div>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
I want to solve Stokes equation using user defined schur complement preconditioner. So I'm reading and testing snes/examples/tutorial/ex55.c, and ex70.c. In these examples, there are comments about the usage, for example:
<div><span style="color:rgb(215,57,30); font-family:Menlo; font-size:11px"><br>
</span></div>
<div><span style="color:rgb(215,57,30); font-family:Menlo; font-size:11px">mpiexec -n 2 ./stokes -nx 32 -ny 48 -fieldsplit_1_user_pc. </span>
<div><span style="font-size:10pt"><br>
</span></div>
<div><span style="font-size:10pt">However the option </span><span style="color:rgb(215,57,30); font-family:Menlo; font-size:11px">-fieldsplit_1_user_pc</span><span style="font-size:10pt"> is not recognized by the executable. The output of the above command
is:</span></div>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div>This ("user_pc") is just a mnemonic. What it means is that you can use any PC in this slot. What he should have written is</div>
<div><br>
</div>
<div> -fieldsplit_1_pc_type <user pc></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>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
<div>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> residual u = 2.6315e-05</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> residual p = 0.000229872</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> residual [u,p] = 0.000231373</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error u = 0.00633503</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error p = 0.121534</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"> discretization error [u,p] = 0.121699</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">WARNING! There are options you set that were not used!</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">WARNING! could be spelling mistake, etc!</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo">Option left: name:-fieldsplit_1_user_pc (no value)</p>
<p style="margin-right:0px; margin-left:0px; font-size:11px; font-family:Menlo"><br>
</p>
<p style="margin-right:0px; margin-left:0px">I'm wondering what's going on. Is that because I need some specific options during my configuring of the package? By the way, the version of PETSc I'm using is 3.4.4. </p>
<p style="margin-right:0px; margin-left:0px"><br>
</p>
<p style="margin-right:0px; margin-left:0px">Best,</p>
<p style="margin-right:0px; margin-left:0px">Hui</p>
<span><font color="#888888"></font></span></div>
<span><font color="#888888"></font></span></div>
<span><font color="#888888"></font></span></div>
<span><font color="#888888"></font></span></blockquote>
<span><font color="#888888"></font></span></div>
<span><font color="#888888"><br>
<br clear="all">
<div><br>
</div>
-- <br>
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 </font></span></div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
<div>
<div><br>
<br clear="all">
<div><br>
</div>
-- <br>
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>
</div>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
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>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>