1 2 #include <petsc-private/fortranimpl.h> 3 #include <petscdmcomposite.h> 4 5 #if defined(PETSC_HAVE_FORTRAN_CAPS) 6 #define dmcompositegetentries1_ DMCOMPOSITEGETENTRIES1 7 #define dmcompositegetentries2_ DMCOMPOSITEGETENTRIES2 8 #define dmcompositegetentries3_ DMCOMPOSITEGETENTRIES3 9 #define dmcompositegetentries4_ DMCOMPOSITEGETENTRIES4 10 #define dmcompositegetentries5_ DMCOMPOSITEGETENTRIES5 11 #define dmcompositecreate_ DMCOMPOSITECREATE 12 #define dmcompositeadddm_ DMCOMPOSITEADDDM 13 #define dmcompositedestroy_ DMCOMPOSITEDESTROY 14 #define dmcompositegetaccess4_ DMCOMPOSITEGETACCESS4 15 #define dmcompositescatter4_ DMCOMPOSITESCATTER4 16 #define dmcompositerestoreaccess4_ DMCOMPOSITERESTOREACCESS4 17 #define dmcompositegetlocalvectors4_ DMCOMPOSITEGETLOCALVECTORS4 18 #define dmcompositerestorelocalvectors4_ DMCOMPOSITERESTORELOCALVECTORS4 19 #define dmcompositegetglobaliss_ DMCOMPOSITEGETGLOBALISS 20 #define dmcompositegetlocaliss_ DMCOMPOSITEGETLOCALISS 21 #elif !defined(PETSC_HAVE_FORTRAN_UNDERSCORE) 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 #define dmcompositegetglobaliss_ dmcompositegetglobaliss 36 #define dmcompositegetlocaliss_ dmcompositegetlocaliss 37 #endif 38 39 EXTERN_C_BEGIN 40 41 void PETSC_STDCALL dmcompositegetentries1_(DM *dm,DM *da1,PetscErrorCode *ierr) 42 { 43 *ierr = DMCompositeGetEntries(*dm,da1); 44 } 45 46 void PETSC_STDCALL dmcompositegetentries2_(DM *dm,DM *da1,DM *da2,PetscErrorCode *ierr) 47 { 48 *ierr = DMCompositeGetEntries(*dm,da1,da2); 49 } 50 51 void PETSC_STDCALL dmcompositegetentries3_(DM *dm,DM *da1,DM *da2,DM *da3,PetscErrorCode *ierr) 52 { 53 *ierr = DMCompositeGetEntries(*dm,da1,da2,da3); 54 } 55 56 void PETSC_STDCALL dmcompositegetentries4_(DM *dm,DM *da1,DM *da2,DM *da3,DM *da4,PetscErrorCode *ierr) 57 { 58 *ierr = DMCompositeGetEntries(*dm,da1,da2,da3,da4); 59 } 60 61 void PETSC_STDCALL dmcompositegetentries5_(DM *dm,DM *da1,DM *da2,DM *da3,DM *da4,DM *da5,PetscErrorCode *ierr) 62 { 63 *ierr = DMCompositeGetEntries(*dm,da1,da2,da3,da4,da5); 64 } 65 66 void PETSC_STDCALL dmcompositecreate_(MPI_Fint * comm,DM *A, int *ierr) 67 { 68 *ierr = DMCompositeCreate(MPI_Comm_f2c(*(comm)),A); 69 } 70 71 void PETSC_STDCALL dmcompositeadddm_(DM *dm,DM *da,PetscErrorCode *ierr) 72 { 73 *ierr = DMCompositeAddDM(*dm,*da); 74 } 75 76 void PETSC_STDCALL dmcompositedestroy_(DM *dm,PetscErrorCode *ierr) 77 { 78 *ierr = DMDestroy(dm); 79 } 80 81 void PETSC_STDCALL dmcompositegetaccess4_(DM *dm,Vec *v,void **v1,void **p1,void **v2,void **p2,PetscErrorCode *ierr) 82 { 83 Vec *vv1 = (Vec*)v1,*vv2 = (Vec*)v2; 84 *ierr = DMCompositeGetAccess(*dm,*v,vv1,(PetscScalar*)p1,vv2,(PetscScalar*)p2); 85 } 86 87 void PETSC_STDCALL dmcompositescatter4_(DM *dm,Vec *v,void *v1,void *p1,void *v2,void *p2,PetscErrorCode *ierr) 88 { 89 Vec *vv1 = (Vec*)v1,*vv2 = (Vec*)v2; 90 *ierr = DMCompositeScatter(*dm,*v,*vv1,(PetscScalar*)p1,*vv2,(PetscScalar*)p2); 91 } 92 93 void PETSC_STDCALL dmcompositerestoreaccess4_(DM *dm,Vec *v,void **v1,void **p1,void **v2,void **p2,PetscErrorCode *ierr) 94 { 95 *ierr = DMCompositeRestoreAccess(*dm,*v,(Vec*)v1,0,(Vec*)v2,0); 96 } 97 98 void PETSC_STDCALL dmcompositegetlocalvectors4_(DM *dm,void **v1,void **p1,void **v2,void **p2,PetscErrorCode *ierr) 99 { 100 Vec *vv1 = (Vec*)v1,*vv2 = (Vec*)v2; 101 *ierr = DMCompositeGetLocalVectors(*dm,vv1,(PetscScalar*)p1,vv2,(PetscScalar*)p2); 102 } 103 104 void PETSC_STDCALL dmcompositerestorelocalvectors4_(DM *dm,void **v1,void **p1,void **v2,void **p2,PetscErrorCode *ierr) 105 { 106 Vec *vv1 = (Vec*)v1,*vv2 = (Vec*)v2; 107 *ierr = DMCompositeRestoreLocalVectors(*dm,vv1,(PetscScalar*)p1,vv2,(PetscScalar*)p2); 108 } 109 110 EXTERN_C_END 111 112 PETSC_EXTERN_C void PETSC_STDCALL dmcompositegetglobaliss_(DM *dm,IS *iss,PetscErrorCode *ierr) 113 { 114 IS *ais; 115 PetscInt i,ndm; 116 *ierr = DMCompositeGetGlobalISs(*dm,&ais); if (*ierr) return; 117 *ierr = DMCompositeGetNumberDM(*dm,&ndm); if (*ierr) return; 118 for (i=0; i<ndm; i++) iss[i] = ais[i]; 119 *ierr = PetscFree(ais); 120 } 121 122 PETSC_EXTERN_C void PETSC_STDCALL dmcompositegetlocaliss_(DM *dm,IS *iss,PetscErrorCode *ierr) 123 { 124 IS *ais; 125 PetscInt i,ndm; 126 *ierr = DMCompositeGetLocalISs(*dm,&ais); if (*ierr) return; 127 *ierr = DMCompositeGetNumberDM(*dm,&ndm); if (*ierr) return; 128 for (i=0; i<ndm; i++) iss[i] = ais[i]; 129 *ierr = PetscFree(ais); 130 } 131