<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Sep 11, 2014 at 10:34 AM, Eric Chamberland <span dir="ltr"><<a href="mailto:Eric.Chamberland@giref.ulaval.ca" target="_blank">Eric.Chamberland@giref.ulaval.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I was just curious to know why the prefix of the sub-matrix a_{11} in a matnest is forced to the ksp prefix by PCSetUp_FieldSplit?<br>
<br>
In my case it changes from "gcrSchur_fieldsplit_a_11_" to "gcrSchur_fieldsplit_schur_" after the PCSetUp_FieldSplit, by this (I think) lines of code in fieldsplit.c (petsc 3.5.2):<br>
<br>
const char *prefix;<br>
ierr = MatGetSubMatrix(pc->pmat,<u></u>ilink->is,ilink->is_col,MAT_<u></u>INITIAL_MATRIX,&jac->pmat[i]);<u></u>CHKERRQ(ierr);<br>
ierr = KSPGetOptionsPrefix(ilink-><u></u>ksp,&prefix);CHKERRQ(ierr);<br>
ierr = MatSetOptionsPrefix(jac->pmat[<u></u>i],prefix);CHKERRQ(ierr);<br>
ierr = MatViewFromOptions(jac->pmat[<u></u>i],NULL,"-mat_view");CHKERRQ(<u></u>ierr);<br></blockquote><div><br></div><div>I did this. My thinking was that if someone wanted to view the Schur complement matrix, then</div><div>this was the correct prefix.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
We wanted to pass options to the a_11 matrix, but using the PC_COMPOSITE_SCHUR, we have to give the a_{11} matrix a unique prefix different from the ksp prefix used to solve the shur complement (which we named "schur") and have MatSetFromOptions use this unique prefixe. It all worked, but we saw this "curiosity" in the kspview (view attached log) about the PCSetUp_FieldSplit which rename our matrix (after all the precautions we took to name it differently)!!!!<br></blockquote><div><br></div><div>I am trying to understand the problem. This matrix is created by by the MatGetSubMatrix() call.</div><div><br></div><div>  1) What options are you trying to pass to it?</div><div><br></div><div>  2) Why would they interfere with the KSP options?</div><div><br></div><div>  Thanks,</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">
The code is working, so maybe this is not an issue, but I can't tell from my knowledge if it can be harmful?<br>
<br>
thank you!<span class="HOEnZb"><font color="#888888"><br>
<br>
Eric<br>
<br>
<br>
</font></span></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>