1 2 #include <petsc-private/fortranimpl.h> 3 4 #if defined(PETSC_HAVE_FORTRAN_CAPS) 5 #define petscinitializefortran_ PETSCINITIALIZEFORTRAN 6 #define petscsetcommonblock_ PETSCSETCOMMONBLOCK 7 #define petscsetfortranbasepointers_ PETSCSETFORTRANBASEPOINTERS 8 #define petsc_null_function_ PETSC_NULL_FUNCTION 9 #elif !defined(PETSC_HAVE_FORTRAN_UNDERSCORE) 10 #define petscinitializefortran_ petscinitializefortran 11 #define petscsetcommonblock_ petscsetcommonblock 12 #define petscsetfortranbasepointers_ petscsetfortranbasepointers 13 #define petsc_null_function_ petsc_null_function 14 #endif 15 16 #if defined(PETSC_HAVE_FORTRAN_UNDERSCORE_UNDERSCORE) 17 #define petsc_null_function_ petsc_null_function__ 18 #endif 19 20 PETSC_EXTERN void PETSC_STDCALL petscsetcommonblock_(MPI_Fint*,MPI_Fint*); 21 PETSC_EXTERN void PETSC_STDCALL petscsetcommonblockmpi_(MPI_Fint*,MPI_Fint*,MPI_Fint*); 22 23 /*@C 24 PetscInitializeFortran - Routine that should be called soon AFTER 25 the call to PetscInitialize() if one is using a C main program 26 that calls Fortran routines that in turn call PETSc routines. 27 28 Collective on PETSC_COMM_WORLD 29 30 Level: beginner 31 32 Notes: 33 PetscInitializeFortran() initializes some of the default viewers, 34 communicators, etc. for use in the Fortran if a user's main program is 35 written in C. PetscInitializeFortran() is NOT needed if a user's main 36 program is written in Fortran; in this case, just calling 37 PetscInitialize() in the main (Fortran) program is sufficient. 38 39 .seealso: PetscInitialize() 40 41 .keywords: Mixing C and Fortran, passing PETSc objects to Fortran 42 @*/ 43 PetscErrorCode PetscInitializeFortran(void) 44 { 45 MPI_Fint c1=0,c2=0; 46 47 if (PETSC_COMM_WORLD) c1 = MPI_Comm_c2f(PETSC_COMM_WORLD); 48 c2 = MPI_Comm_c2f(PETSC_COMM_SELF); 49 petscsetcommonblock_(&c1,&c2); 50 51 #if defined(PETSC_USE_REAL___FLOAT128) 52 { 53 MPI_Fint freal,fscalar,fsum; 54 freal = MPI_Type_c2f(MPIU_REAL); 55 fscalar = MPI_Type_c2f(MPIU_SCALAR); 56 fsum = MPI_Op_c2f(MPIU_SUM); 57 petscsetcommonblockmpi_(&freal,&fscalar,&fsum); 58 } 59 #endif 60 return 0; 61 } 62 63 PETSC_EXTERN void PETSC_STDCALL petscinitializefortran_(int *ierr) 64 { 65 *ierr = PetscInitializeFortran(); 66 } 67 68 PETSC_EXTERN void PETSC_STDCALL petscsetfortranbasepointers_(char *fnull_character PETSC_MIXED_LEN(len), 69 void *fnull_integer,void *fnull_scalar,void * fnull_double, 70 void *fnull_real,void *fnull_object, 71 void* fnull_truth,void (*fnull_function)(void) PETSC_END_LEN(len)) 72 { 73 PETSC_NULL_CHARACTER_Fortran = fnull_character; 74 PETSC_NULL_INTEGER_Fortran = fnull_integer; 75 PETSC_NULL_SCALAR_Fortran = fnull_scalar; 76 PETSC_NULL_DOUBLE_Fortran = fnull_double; 77 PETSC_NULL_REAL_Fortran = fnull_real; 78 PETSC_NULL_OBJECT_Fortran = fnull_object; 79 PETSC_NULL_BOOL_Fortran = fnull_truth; 80 PETSC_NULL_FUNCTION_Fortran = fnull_function; 81 } 82 83 /* 84 A valid address for the fortran variable PETSC_NULL_FUNCTION 85 */ 86 PETSC_EXTERN void petsc_null_function_(void) 87 { 88 return; 89 } 90 91 92