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