<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Thanks a lot.<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 19-May-2018, at 9:16 PM, Smith, Barry F. <<a href="mailto:bsmith@mcs.anl.gov" class="">bsmith@mcs.anl.gov</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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=""> Praveen,</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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=""> Ahh, we didn't have support for passing PETSC_NULL_MAT from Fortran for this routine. I have added it in the branch barry/fix-null-matrix-set-jacobian/maint which after testing will go into the maint branch.<span class="Apple-converted-space"> </span></span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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=""> Could you please try this branch and let us know if it resolves the problem (or a new problem pops up)?</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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=""> Barry</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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: Helvetica; font-size: 13px; font-style: normal; font-variant-caps: normal; font-weight: normal; 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=""> Unfortunately handling null objects from Fortran requires us to manually tweak the Fortran interface functions and sometimes we forget or don't realize they need fixing until someone reports a problem.</span></div></blockquote></div><br class=""><div class="">I meanwhile passed a Mat object to TSSetRHSJacobian function. It seems that with -snes_mf and PCSHELL, my RHSJacobian is never invoked. Is this the correct behaviour ?</div><div class=""><br class=""></div><div class="">In RHSJacobian, I am storing the current solution into my usercontext so that I can use it in my shell preconditioner</div><div class=""><br class=""></div><div class=""><div class="">subroutine RHSJacobian(ts, time, u, J, P, ctx, ierr)</div><div class=""><br class=""></div><div class=""> call VecCopy(u, ctx%p%v_u, ierr); CHKERRQ(ierr)</div></div><div class=""><br class=""></div><div class="">end</div><div class=""><br class=""></div><div class="">However, since my RHSJacobian is never invoked, I dont have access to the current solution needed in my ApplyPC. How can I access the current solution in my ApplyPC function ?</div><div class=""><br class=""></div><div class="">Another issue I am facing is this. I attach my user context to the PC</div><div class=""><br class=""></div><div class=""><div class=""> call TSGetSNES(ts, snes, ierr); CHKERRQ(ierr)</div><div class=""> call SNESGetKSP(snes, ksp, ierr); CHKERRQ(ierr)</div><div class=""> call KSPGetPC(ksp, pc, ierr); CHKERRQ(ierr)</div><div class=""> call PCSetType(pc, PCSHELL, ierr); CHKERRQ(ierr)</div><div class=""> call PCShellSetContext(pc, ctx, ierr); CHKERRQ(ierr)</div><div class=""> call PCShellSetApply(pc, ApplyPC, ierr); CHKERRQ(ierr)</div></div><div class=""><br class=""></div><div class="">Inside my ApplyPC, I access ctx using </div><div class=""><br class=""></div><div class=""><div class="">subroutine ApplyPC(pc, x, y, ierr)</div></div><div class=""><div class=""> use petscpc</div><div class=""> use mtsdata</div><div class=""> implicit none</div><div class=""> PC :: pc</div><div class=""> Vec :: x, y</div><div class=""> PetscErrorCode :: ierr</div><div class=""> ! Local variables</div><div class=""> type(tsdata) :: ctx</div></div><div class=""><br class=""></div><div class=""><div class=""> call PCShellGetContext(pc, ctx, ierr); CHKERRQ(ierr)</div></div><div class=""><br class=""></div><div class="">end</div><div class=""><br class=""></div><div class="">But the ctx I get from above is not correct, since it does not have the data I stored in my original ctx.</div><div class=""><br class=""></div><div class="">I am rather clueless at this stage.</div><div class=""><br class=""></div><div class="">Thanks for your help.</div><div class="">Best regards</div><div class="">praveen</div></body></html>