<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<span style="font-family: STSongti-SC-Regular;" class="">Are there other approaches to solve this kind of systems in PETSc except for field-split methods? </span>
<div class="" style="font-family: STSongti-SC-Regular;"><br class="">
</div>
<div class="" style="font-family: STSongti-SC-Regular;">Thanks,</div>
<div class="" style="font-family: STSongti-SC-Regular;">Xiaofeng</div>
<div><br class="">
<blockquote type="cite" class="">
<div class="">On Sep 26, 2022, at 14:13, Jed Brown <<a href="mailto:jed@jedbrown.org" class="">jed@jedbrown.org</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class=""><span style="caret-color: rgb(0, 0, 0); font-family: STSongti-SC-Regular; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">This
is the joy of factorization field-split methods. The actual Schur complement is dense, so we represent it implicitly. A common strategy is to assemble the mass matrix and drop it in the 11 block of the Pmat. You can check out some examples in the repository
for incompressible flow (Stokes problems). The LSC (least squares commutator) is another option. You'll likely find that lumping diag(A00)^{-1} works poorly because the resulting operator behaves like a Laplacian rather than like a mass matrix.</span><br style="caret-color: rgb(0, 0, 0); font-family: STSongti-SC-Regular; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<br style="caret-color: rgb(0, 0, 0); font-family: STSongti-SC-Regular; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="caret-color: rgb(0, 0, 0); font-family: STSongti-SC-Regular; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">晓峰
何 <</span><a href="mailto:tlanyan@hotmail.com" style="font-family: STSongti-SC-Regular; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">tlanyan@hotmail.com</a><span style="caret-color: rgb(0, 0, 0); font-family: STSongti-SC-Regular; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">>
writes:</span><br style="caret-color: rgb(0, 0, 0); font-family: STSongti-SC-Regular; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<br style="caret-color: rgb(0, 0, 0); font-family: STSongti-SC-Regular; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<blockquote type="cite" style="font-family: STSongti-SC-Regular; font-size: 18px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
If assigned a preconditioner to A11 with this cmd options:<br class="">
<br class="">
-fieldsplit_0_ksp_type gmres -fieldsplit_0_pc_type ilu -fieldsplit_1_ksp_type gmres -fieldsplit_1_pc_type ilu<br class="">
<br class="">
Then I got this error:<br class="">
<br class="">
"Could not locate a solver type for factorization type ILU and matrix type schurcomplement"<br class="">
<br class="">
How could I specify a preconditioner for A11?<br class="">
<br class="">
BR,<br class="">
Xiaofeng<br class="">
<br class="">
<br class="">
On Sep 26, 2022, at 11:02, 晓峰 何 <<a href="mailto:tlanyan@hotmail.com" class="">tlanyan@hotmail.com</a><<a href="mailto:tlanyan@hotmail.com" class="">mailto:tlanyan@hotmail.com</a>>> wrote:<br class="">
<br class="">
-fieldsplit_0_ksp_type gmres -fieldsplit_0_pc_type ilu -fieldsplit_1_ksp_type gmres -fieldsplit_1_pc_type none</blockquote>
</div>
</blockquote>
</div>
<br class="">
</body>
</html>