1) Always use options. It makes you code generic and process more flexible<br><br>2) If after thinking for a few hours, you still really want to do it in the code:<br><br>  <a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/PC/PCBJacobiGetSubKSP.html">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/PC/PCBJacobiGetSubKSP.html</a><br>
<br>   Matt<br><br><div class="gmail_quote">On Mon, Dec 7, 2009 at 11:52 AM,  <span dir="ltr">&lt;<a href="mailto:hxie@umn.edu">hxie@umn.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi,<br>
<br>
I try to customize the subksp for block jacobi method. Does anybody have an fortran example for this?<br>
For c language, we can declare &quot;ksp *subksp&quot;. Then how for fortran? Thanks.<br>
<br>
Bests,<br>
Hui<br>
</blockquote></div><br><br clear="all"><br>-- <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<br>