1 2 #include <private/fortranimpl.h> 3 #include <petscdmcomposite.h> 4 5 #if defined(PETSC_HAVE_FORTRAN_CAPS) 6 #define dmgetmatrix_ DMGETMATRIX 7 #define dmcompositegetentries1_ DMCOMPOSITEGETENTRIES1 8 #define dmcompositegetentries2_ DMCOMPOSITEGETENTRIES2 9 #define dmcompositegetentries3_ DMCOMPOSITEGETENTRIES3 10 #define dmcompositegetentries4_ DMCOMPOSITEGETENTRIES4 11 #define dmcompositegetentries5_ DMCOMPOSITEGETENTRIES5 12 #define dmcompositecreate_ DMCOMPOSITECREATE 13 #define dmcompositeadddm_ DMCOMPOSITEADDDM 14 #define dmcompositedestroy_ DMCOMPOSITEDESTROY 15 #define dmcompositegetaccess4_ DMCOMPOSITEGETACCESS4 16 #define dmcompositescatter4_ DMCOMPOSITESCATTER4 17 #define dmcompositerestoreaccess4_ DMCOMPOSITERESTOREACCESS4 18 #define dmcompositegetlocalvectors4_ DMCOMPOSITEGETLOCALVECTORS4 19 #define dmcompositerestorelocalvectors4_ DMCOMPOSITERESTORELOCALVECTORS4 20 #elif !defined(PETSC_HAVE_FORTRAN_UNDERSCORE) 21 #define dmgetmatrix_ dmgetmatrix 22 #define dmcompositegetentries1_ dmcompositegetentries1 23 #define dmcompositegetentries2_ dmcompositegetentries2 24 #define dmcompositegetentries3_ dmcompositegetentries3 25 #define dmcompositegetentries4_ dmcompositegetentries4 26 #define dmcompositegetentries5_ dmcompositegetentries5 27 #define dmcompositecreate_ dmcompositecreate 28 #define dmcompositeadddm_ dmcompositeadddm 29 #define dmcompositedestroy_ dmcompositedestroy 30 #define dmcompositegetaccess4_ dmcompositegetaccess4 31 #define dmcompositescatter4_ dmcompositescatter4 32 #define dmcompositerestoreaccess4_ dmcompositerestoreaccess4 33 #define dmcompositegetlocalvectors4_ dmcompositegetlocalvectors4 34 #define dmcompositerestorelocalvectors4_ dmcompositerestorelocalvectors4 35 #endif 36 37 EXTERN_C_BEGIN 38 39 void PETSC_STDCALL dmgetmatrix_(DM *dm,CHAR mat_type PETSC_MIXED_LEN(len),Mat *J,PetscErrorCode *ierr PETSC_END_LEN(len)) 40 { 41 char *t; 42 FIXCHAR(mat_type,len,t); 43 *ierr = DMGetMatrix(*dm,t,J); 44 FREECHAR(mat_type,t); 45 } 46 47 void PETSC_STDCALL dmcompositegetentries1_(DM *dm,DM *da1,PetscErrorCode *ierr) 48 { 49 *ierr = DMCompositeGetEntries(*dm,da1); 50 } 51 52 void PETSC_STDCALL dmcompositegetentries2_(DM *dm,DM *da1,DM *da2,PetscErrorCode *ierr) 53 { 54 *ierr = DMCompositeGetEntries(*dm,da1,da2); 55 } 56 57 void PETSC_STDCALL dmcompositegetentries3_(DM *dm,DM *da1,DM *da2,DM *da3,PetscErrorCode *ierr) 58 { 59 *ierr = DMCompositeGetEntries(*dm,da1,da2,da3); 60 } 61 62 void PETSC_STDCALL dmcompositegetentries4_(DM *dm,DM *da1,DM *da2,DM *da3,DM *da4,PetscErrorCode *ierr) 63 { 64 *ierr = DMCompositeGetEntries(*dm,da1,da2,da3,da4); 65 } 66 67 void PETSC_STDCALL dmcompositegetentries5_(DM *dm,DM *da1,DM *da2,DM *da3,DM *da4,DM *da5,PetscErrorCode *ierr) 68 { 69 *ierr = DMCompositeGetEntries(*dm,da1,da2,da3,da4,da5); 70 } 71 72 void PETSC_STDCALL dmcompositecreate_(MPI_Fint * comm,DM *A, int *ierr ){ 73 *ierr = DMCompositeCreate(MPI_Comm_f2c( *(comm) ),A); 74 } 75 76 void PETSC_STDCALL dmcompositeadddm_(DM *dm,DM *da,PetscErrorCode *ierr) 77 { 78 *ierr = DMCompositeAddDM(*dm,*da); 79 } 80 81 void PETSC_STDCALL dmcompositedestroy_(DM *dm,PetscErrorCode *ierr) 82 { 83 *ierr = DMDestroy(dm); 84 } 85 86 void PETSC_STDCALL dmcompositegetaccess4_(DM *dm,Vec *v,void **v1,void **p1,void **v2,void **p2,PetscErrorCode *ierr) 87 { 88 Vec *vv1 = (Vec*)v1,*vv2 = (Vec*)v2; 89 *ierr = DMCompositeGetAccess(*dm,*v,vv1,(PetscScalar*)p1,vv2,(PetscScalar*)p2); 90 } 91 92 void PETSC_STDCALL dmcompositescatter4_(DM *dm,Vec *v,void *v1,void *p1,void *v2,void *p2,PetscErrorCode *ierr) 93 { 94 Vec *vv1 = (Vec*)v1,*vv2 = (Vec*)v2; 95 *ierr = DMCompositeScatter(*dm,*v,*vv1,(PetscScalar*)p1,*vv2,(PetscScalar*)p2); 96 } 97 98 void PETSC_STDCALL dmcompositerestoreaccess4_(DM *dm,Vec *v,void **v1,void **p1,void **v2,void **p2,PetscErrorCode *ierr) 99 { 100 *ierr = DMCompositeRestoreAccess(*dm,*v,(Vec*)v1,0,(Vec*)v2,0); 101 } 102 103 void PETSC_STDCALL dmcompositegetlocalvectors4_(DM *dm,void **v1,void **p1,void **v2,void **p2,PetscErrorCode *ierr) 104 { 105 Vec *vv1 = (Vec*)v1,*vv2 = (Vec*)v2; 106 *ierr = DMCompositeGetLocalVectors(*dm,vv1,(PetscScalar*)p1,vv2,(PetscScalar*)p2); 107 } 108 109 void PETSC_STDCALL dmcompositerestorelocalvectors4_(DM *dm,void **v1,void **p1,void **v2,void **p2,PetscErrorCode *ierr) 110 { 111 Vec *vv1 = (Vec*)v1,*vv2 = (Vec*)v2; 112 *ierr = DMCompositeRestoreLocalVectors(*dm,vv1,(PetscScalar*)p1,vv2,(PetscScalar*)p2); 113 } 114 115 EXTERN_C_END 116 117