Actual source code: zmfnf.c
slepc-3.23.1 2025-05-01
1: /*
2: - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
3: SLEPc - Scalable Library for Eigenvalue Problem Computations
4: Copyright (c) 2002-, Universitat Politecnica de Valencia, Spain
6: This file is part of SLEPc.
7: SLEPc is distributed under a 2-clause BSD license (see LICENSE).
8: - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
9: */
11: #include <petsc/private/ftnimpl.h>
12: #include <slepcmfn.h>
14: #if defined(PETSC_HAVE_FORTRAN_CAPS)
15: #define mfnmonitordefault_ MFNMONITORDEFAULT
16: #define mfnmonitorset_ MFNMONITORSET
17: #elif !defined(PETSC_HAVE_FORTRAN_UNDERSCORE)
18: #define mfnmonitordefault_ mfnmonitordefault
19: #define mfnmonitorset_ mfnmonitorset
20: #endif
22: /*
23: These cannot be called from Fortran but allow Fortran users
24: to transparently set these monitors from .F code
25: */
26: SLEPC_EXTERN void mfnmonitordefault_(MFN*,PetscInt*,PetscReal*,PetscViewerAndFormat*,PetscErrorCode*);
28: static struct {
29: PetscFortranCallbackId monitor;
30: PetscFortranCallbackId monitordestroy;
31: } _cb;
33: /* These are not extern C because they are passed into non-extern C user level functions */
34: static PetscErrorCode ourmonitor(MFN mfn,PetscInt i,PetscReal d,void* ctx)
35: {
36: PetscObjectUseFortranCallback(mfn,_cb.monitor,(MFN*,PetscInt*,PetscReal*,void*,PetscErrorCode*),(&mfn,&i,&d,_ctx,&ierr));
37: }
39: static PetscErrorCode ourdestroy(void** ctx)
40: {
41: MFN mfn = (MFN)*ctx;
42: PetscObjectUseFortranCallback(mfn,_cb.monitordestroy,(void*,PetscErrorCode*),(_ctx,&ierr));
43: }
45: SLEPC_EXTERN void mfnmonitorset_(MFN *mfn,void (*monitor)(MFN*,PetscInt*,PetscReal*,void*,PetscErrorCode*),void *mctx,void (*monitordestroy)(void *,PetscErrorCode*),PetscErrorCode *ierr)
46: {
47: CHKFORTRANNULLOBJECT(mctx);
48: CHKFORTRANNULLFUNCTION(monitordestroy);
49: if ((PetscVoidFunction)monitor == (PetscVoidFunction)mfnmonitordefault_) {
50: *ierr = MFNMonitorSet(*mfn,(PetscErrorCode (*)(MFN,PetscInt,PetscReal,void*))MFNMonitorDefault,*(PetscViewerAndFormat**)mctx,(PetscCtxDestroyFn*)PetscViewerAndFormatDestroy);
51: } else {
52: *ierr = PetscObjectSetFortranCallback((PetscObject)*mfn,PETSC_FORTRAN_CALLBACK_CLASS,&_cb.monitor,(PetscVoidFunction)monitor,mctx); if (*ierr) return;
53: *ierr = PetscObjectSetFortranCallback((PetscObject)*mfn,PETSC_FORTRAN_CALLBACK_CLASS,&_cb.monitordestroy,(PetscVoidFunction)monitordestroy,mctx); if (*ierr) return;
54: *ierr = MFNMonitorSet(*mfn,ourmonitor,*mfn,ourdestroy);
55: }
56: }