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