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