Dear All,
I am having some issues upgrading to PETSc 3.7 due to some changes in the
KSPMonitor routines.
I need to configure my KSP solver to output its convergence log to an existing
ASCII viewer through a custom monitor, printing only one every 10 iterations.
The calling code is Fortran, we wrote a small C code that just calls the
default monitor every 10 iterations, passing the viewer as last argument. If I
understood correctly it is now necessary to set up a PetscViewerAndFormat
structure and pass that as last argument to the monitor routine.
I tried to create one with PetscViewerAndFormatCreate but I'm getting a runtime
error (see below).
Could you please help me understand what I did wrong?
Thanks in advance,
Federico
Error:
[1]PETSC ERROR: #1 PetscViewerPushFormat() line 149 in
/PETSc/petsc-3.7.0/src/sys/classes/viewer/interface/viewa.c
Error Message --------------------------------------------------------------
[8]PETSC ERROR: Argument out of range
[8]PETSC ERROR: Too many PetscViewerPushFormat(), perhaps you forgot
PetscViewerPopFormat()?
Fortran code:
PetscViewer :: viewer
PetscViewerAndFormat :: vf
external ShellKSPMonitor
...
call PetscViewerAndFormatCreate(viewer, PETSC_VIEWER_DEFAULT, vf,ierr)
call KSPMonitorSet(ksp, ShellKSPMonitor, vf, PetscViewerAndFormatDestroy,
ierr)
C code:
PetscErrorCode shellkspmonitor_(KSP *ksp, PetscInt *n, PetscReal *rnorm, void
*ptr)
{
PetscErrorCode ierr=0;
if (*n % 10 == 0) {
ierr =
KSPMonitorTrueResidualNorm(*ksp,*n,*rnorm,(PetscViewerAndFormat*)ptr);CHKERRQ(ierr);
}
return ierr;
}
______ ______ ______
Federico Miorelli
Senior R&D Geophysicist
Subsurface Imaging - General Geophysics Italy
This email and any accompanying attachments are confidential. If you received
this email by mistake, please delete
it from your system. Any review, disclosure, copying, distribution, or use of
the email by others is strictly prohibited.
