1 #include <petsc/private/fortranimpl.h> 2 #include <petscds.h> 3 #include <petscviewer.h> 4 5 #if defined(PETSC_HAVE_FORTRAN_CAPS) 6 #define petscdsviewfromoptions_ PETSCDSVIEWFROMOPTIONS 7 #define petscdsview_ PETSCDSVIEW 8 #define petscdssetcontext_ PETSCDSSETCONTEXT 9 #define petscdssetriemannsolver_ PETSCDSSETRIEMANNSOLVER 10 #elif !defined(PETSC_HAVE_FORTRAN_UNDERSCORE) 11 #define petscdsviewfromoptions_ petscdsviewfromoptions 12 #define petscdsview_ petscdsview 13 #define petscdssetcontext_ petscdssetcontext 14 #define petscdssetriemannsolver_ petscdssetriemannsolver 15 #endif 16 17 static PetscFortranCallbackId riemannsolver; 18 19 // We can't use PetscObjectUseFortranCallback() because this function returns void 20 static void ourriemannsolver(PetscInt dim,PetscInt Nf,const PetscReal x[],const PetscReal n[],const PetscScalar uL[],const PetscScalar uR[],PetscInt numConstants,const PetscScalar constants[],PetscScalar flux[],void *ctx) 21 { 22 void (*func)(PetscInt *dim,PetscInt *Nf,const PetscReal x[],const PetscReal n[],const PetscScalar uL[],const PetscScalar uR[],const PetscInt *numConstants,const PetscScalar constants[],PetscScalar flux[],void *ctx); 23 void *_ctx; 24 PetscCallAbort(PETSC_COMM_SELF,PetscObjectGetFortranCallback((PetscObject)ctx,PETSC_FORTRAN_CALLBACK_CLASS,riemannsolver,(PetscVoidFunction*)&func,&_ctx)); 25 if (func) { 26 (*func)(&dim,&Nf,x,n,uL,uR,&numConstants,constants,flux,_ctx); 27 } 28 } 29 30 PETSC_EXTERN void petscdsviewfromoptions_(PetscDS *ao,PetscObject obj,char* type,PetscErrorCode *ierr,PETSC_FORTRAN_CHARLEN_T len) 31 { 32 char *t; 33 34 FIXCHAR(type,len,t); 35 CHKFORTRANNULLOBJECT(obj); 36 *ierr = PetscDSViewFromOptions(*ao,obj,t);if (*ierr) return; 37 FREECHAR(type,t); 38 } 39 40 PETSC_EXTERN void petscdsview_(PetscDS *prob,PetscViewer *vin,PetscErrorCode *ierr) 41 { 42 PetscViewer v; 43 PetscPatchDefaultViewers_Fortran(vin,v); 44 *ierr = PetscDSView(*prob,v);if (*ierr) return; 45 } 46 47 PETSC_EXTERN void petscdssetcontext_(PetscDS *prob,PetscInt *f,void *ctx,PetscErrorCode *ierr) 48 { 49 *ierr = PetscDSSetContext(*prob,*f,*prob);if (*ierr) return; 50 } 51 52 PETSC_EXTERN void petscdssetriemannsolver_(PetscDS *prob,PetscInt *f,void (*rs)(PetscInt*,PetscInt*,PetscReal*,PetscReal*,PetscScalar*,PetscScalar*,PetscInt*,PetscScalar*,PetscScalar*,void*,PetscErrorCode*),PetscErrorCode *ierr) 53 { 54 *ierr = PetscObjectSetFortranCallback((PetscObject)*prob,PETSC_FORTRAN_CALLBACK_CLASS,&riemannsolver,(PetscVoidFunction)rs,NULL);if (*ierr) return; 55 *ierr = PetscDSSetRiemannSolver(*prob,*f,ourriemannsolver);if (*ierr) return; 56 } 57