1c6db04a5SJed Brown #include <private/fortranimpl.h> 2c6db04a5SJed Brown #include <petscdm.h> 39a42bb27SBarry Smith 49a42bb27SBarry Smith #if defined(PETSC_HAVE_FORTRAN_CAPS) 524e7ded0SBarry Smith #define dmview_ DMVIEW 6e727c939SJed Brown #define dmcreatecoloring_ DMCREATECOLORING 7*7e833e3aSBarry Smith #define dmsetinitialguess_ DMSETINITIALGUESS 8*7e833e3aSBarry Smith #define dmsetfunction_ DMSETFUNCTION 9*7e833e3aSBarry Smith #define dmsetjacobian_ DMSETJACOBIAN 10950540a4SJed Brown #define dmcreatematrix_ DMCREATEMATRIX 119a42bb27SBarry Smith #elif !defined(PETSC_HAVE_FORTRAN_UNDERSCORE) 1224e7ded0SBarry Smith #define dmview_ dmview 13e727c939SJed Brown #define dmcreatecoloring_ dmcreatecoloring 14*7e833e3aSBarry Smith #define dmsetinitialguess_ dmsetinitialguess 15*7e833e3aSBarry Smith #define dmsetfunction_ dmsetfunction 16*7e833e3aSBarry Smith #define dmsetjacobian_ dmsetjacobian 17950540a4SJed Brown #define dmcreatematrix_ dmcreatematrix 189a42bb27SBarry Smith #endif 199a42bb27SBarry Smith 20*7e833e3aSBarry Smith static PetscErrorCode ourdminitialguess(DM dm,Vec x) 21*7e833e3aSBarry Smith { 22*7e833e3aSBarry Smith PetscErrorCode ierr = 0; 23*7e833e3aSBarry Smith (*(void (PETSC_STDCALL *)(DM*,Vec*,PetscErrorCode*))(((PetscObject)dm)->fortran_func_pointers[0]))(&dm,&x,&ierr);CHKERRQ(ierr); 24*7e833e3aSBarry Smith return 0; 25*7e833e3aSBarry Smith } 26*7e833e3aSBarry Smith 27*7e833e3aSBarry Smith static PetscErrorCode ourdmfunction(DM dm,Vec x,Vec b) 28*7e833e3aSBarry Smith { 29*7e833e3aSBarry Smith PetscErrorCode ierr = 0; 30*7e833e3aSBarry Smith (*(void (PETSC_STDCALL *)(DM*,Vec*,Vec*,PetscErrorCode*))(((PetscObject)dm)->fortran_func_pointers[1]))(&dm,&x,&b,&ierr);CHKERRQ(ierr); 31*7e833e3aSBarry Smith return 0; 32*7e833e3aSBarry Smith } 33*7e833e3aSBarry Smith 34*7e833e3aSBarry Smith static PetscErrorCode ourdmjacobian(DM dm,Vec x,Mat A,Mat B,MatStructure *str) 35*7e833e3aSBarry Smith { 36*7e833e3aSBarry Smith PetscErrorCode ierr = 0; 37*7e833e3aSBarry Smith (*(void (PETSC_STDCALL *)(DM*,Vec*,Mat*,Mat*,MatStructure*,PetscErrorCode*))(((PetscObject)dm)->fortran_func_pointers[2]))(&dm,&x,&A,&B,str,&ierr);CHKERRQ(ierr); 38*7e833e3aSBarry Smith return 0; 39*7e833e3aSBarry Smith } 40*7e833e3aSBarry Smith 41*7e833e3aSBarry Smith EXTERN_C_BEGIN 42*7e833e3aSBarry Smith void PETSC_STDCALL dmsetinitialguess_(DM *dm,PetscErrorCode (*f)(DM*,Vec*,PetscErrorCode*), int *ierr ) 43*7e833e3aSBarry Smith { 44*7e833e3aSBarry Smith PetscObjectAllocateFortranPointers(*dm,12); 45*7e833e3aSBarry Smith ((PetscObject)*dm)->fortran_func_pointers[0] = (PetscVoidFunction)f; 46*7e833e3aSBarry Smith *ierr = DMSetInitialGuess(*dm,ourdminitialguess); 47*7e833e3aSBarry Smith } 48*7e833e3aSBarry Smith EXTERN_C_END 49*7e833e3aSBarry Smith 50*7e833e3aSBarry Smith EXTERN_C_BEGIN 51*7e833e3aSBarry Smith void PETSC_STDCALL dmsetfunction_(DM *dm,PetscErrorCode (*f)(DM*,Vec*,Vec*,PetscErrorCode*), int *ierr ) 52*7e833e3aSBarry Smith { 53*7e833e3aSBarry Smith PetscObjectAllocateFortranPointers(*dm,12); 54*7e833e3aSBarry Smith ((PetscObject)*dm)->fortran_func_pointers[1] = (PetscVoidFunction)f; 55*7e833e3aSBarry Smith *ierr = DMSetFunction(*dm,ourdmfunction); 56*7e833e3aSBarry Smith } 57*7e833e3aSBarry Smith EXTERN_C_END 58*7e833e3aSBarry Smith 59*7e833e3aSBarry Smith EXTERN_C_BEGIN 60*7e833e3aSBarry Smith void PETSC_STDCALL dmsetjacobian_(DM *dm,PetscErrorCode (*f)(DM*,Vec*,Mat*,Mat*,MatStructure*,PetscErrorCode*), int *ierr ) 61*7e833e3aSBarry Smith { 62*7e833e3aSBarry Smith PetscObjectAllocateFortranPointers(*dm,12); 63*7e833e3aSBarry Smith ((PetscObject)*dm)->fortran_func_pointers[2] = (PetscVoidFunction)f; 64*7e833e3aSBarry Smith *ierr = DMSetJacobian(*dm,ourdmjacobian); 65*7e833e3aSBarry Smith } 66*7e833e3aSBarry Smith EXTERN_C_END 67*7e833e3aSBarry Smith 689a42bb27SBarry Smith EXTERN_C_BEGIN 69e727c939SJed Brown void PETSC_STDCALL dmcreatecoloring_(DM *dm,ISColoringType *ctype, CHAR mtype PETSC_MIXED_LEN(len),ISColoring *coloring, int *ierr PETSC_END_LEN(len)) 70d1e2c406SBarry Smith { 71d1e2c406SBarry Smith char *t; 72d1e2c406SBarry Smith 73d1e2c406SBarry Smith FIXCHAR(mtype,len,t); 74e727c939SJed Brown *ierr = DMCreateColoring(*dm,*ctype,t,coloring); 75d1e2c406SBarry Smith FREECHAR(mtype,t); 76d1e2c406SBarry Smith } 77d1e2c406SBarry Smith EXTERN_C_END 78d1e2c406SBarry Smith 79d1e2c406SBarry Smith EXTERN_C_BEGIN 8024e7ded0SBarry Smith void PETSC_STDCALL dmview_(DM *da,PetscViewer *vin,PetscErrorCode *ierr) 819a42bb27SBarry Smith { 829a42bb27SBarry Smith PetscViewer v; 839a42bb27SBarry Smith PetscPatchDefaultViewers_Fortran(vin,v); 849a42bb27SBarry Smith *ierr = DMView(*da,v); 859a42bb27SBarry Smith } 86564755cdSBarry Smith EXTERN_C_END 87950540a4SJed Brown 88950540a4SJed Brown EXTERN_C_BEGIN 89950540a4SJed Brown void PETSC_STDCALL dmcreatematrix_(DM *dm,CHAR mat_type PETSC_MIXED_LEN(len),Mat *J,PetscErrorCode *ierr PETSC_END_LEN(len)) 90950540a4SJed Brown { 91950540a4SJed Brown char *t; 92950540a4SJed Brown FIXCHAR(mat_type,len,t); 93950540a4SJed Brown *ierr = DMCreateMatrix(*dm,t,J); 94950540a4SJed Brown FREECHAR(mat_type,t); 95950540a4SJed Brown } 96950540a4SJed Brown EXTERN_C_END 97