[petsc-users] Help with KSPSetConvergenceTest

[Barry Smith]

   The code will not work as written because KSPConvergedDefault() requires a context created with KSPConvergedDefaultCreate(). 

   Here is a starting point for what you need, in main

   integer*8 defaultctx
   extern MyKSPConverged, KSPConvergedDefaultDestroy

KSPDefaultConvergedCreate(defaultctx,ierr)
KSPSetConvergenceTest <https://petsc.org/main/manualpages/KSP/KSPSetConvergenceTest/>(ksp,MyKSPConverged,defaultctx,KSPConvergedDefaultDestroy, ierr)


subroutine MyKSPConverged(ksp,n,rnorm,flag,defaultctx,ierr)

       KSP              ksp
       PetscErrorCode ierr
       PetscInt n
       integer*8 defaultctx
       KSPConvergedReason flag
       PetscReal rnorm

       if (n>1) then
         call KSPConvergedDefault(ksp, n, rnorm, flag,defaultctx, ierr)
       else
         flag = 0
       endif
       ierr = 0

    end subroutine MyKSPConverged



> On May 7, 2023, at 10:09 AM, Matthew Knepley <knepley at gmail.com> wrote:
> 
> On Sun, May 7, 2023 at 10:02 AM Edoardo alinovi <edoardo.alinovi at gmail.com <mailto:edoardo.alinovi at gmail.com>> wrote:
>> Thanks, 
>> 
>> Is this a reasonable thing to do if I want to replicate what KSP is doing by default?
> 
> Yes. The other option is to pass along 'dummy'
> 
>   Thanks,
> 
>     Matt
> 
> -- 
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
> -- Norbert Wiener
> 
> https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230507/3dd0b172/attachment.html>