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