[petsc-users] Access to 'step' in 'toil_ims.f90' and 'pm_toil_ims.f90'

[Jed Brown]

The files and functions you reference are not part of PETSc.  Perhaps
you can look at how they call PETSc.

Robert Annewandter <robert.annewandter at opengosim.com> writes:

> Hi,
>
> Is there a way to access 'step' as used in the output log file :
>
>  Step     35 Time=  1.00000E-03 Dt=  2.73903E-05 [d] snes_conv_reason:   10,
>
> from within the subroutines ' TOilImsResidual(snes,xx,r,realization,ierr)' and
> 'TOilImsJacobian(snes,xx,A,B,realization,ierr)' in 'toil_ims.f90', as well as
> from within 'PMTOilImsCheckUpdatePre(this,line_search,X,dX,changed,ierr)' in
> 'pm_toil_ims.f90'?
>
> I like to make sure that J, b and dx are taken from the same step. I can do
> that if the filenames of the exports reflect step and NR iteration. The problem
> is that, e.g. 'TOilImsResidual()' is called more often than is 'TOilImsJacobian
> ()', so using a simple counter which gets increased every time when the
> subroutines are called doesn't help unfortunately.
>
> Thank you for any pointers!
>
> Robert
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 800 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20161121/ca655dd3/attachment.pgp>