Actual source code: petscdef.h

  1: !
  2: !
  3: !  Part of the base include file for Fortran use of PETSc.
  4: !  Note: This file should contain only define statements and
  5: !  not the declaration of variables.

  7: ! No spaces for #defines as some compilers (PGI) also adds
  8: ! those additional spaces during preprocessing - bad for fixed format
  9: !
 10: #if !defined (__PETSCDEF_H)
 12: !
 13: #include "petscconf.h"
 14: !
 15: #define MPI_Comm integer
 16: !
 17: #define PetscEnum integer
 18: !
 19: #define PetscTruth PetscEnum
 20: #define PetscDataType PetscEnum
 21: #define PetscFPTrap PetscEnum

 23: #define PetscErrorCode integer
 24: #define PetscCookie integer
 25: #define PetscEvent integer
 26: #define PetscMPIInt integer
 27: !
 28: !
 29: ! The real*8,complex*16 notatiton is used so that the
 30: ! PETSc double/complex variables are not affected by
 31: ! compiler options like -r4,-r8, sometimes invoked
 32: ! by the user. NAG compiler does not like integer*4,real*8

 34: #if defined(PETSC_USE_FORTRANKIND)
 35: #define integer8 integer(kind=selected_int_kind(10))
 36: #define integer4 integer(kind=selected_int_kind(5))
 37: #else
 38: #define integer8 integer*8
 39: #define integer4 integer*4
 40: #endif

 42: #if (PETSC_SIZEOF_VOID_P == 8)
 43: #define PetscFortranAddr integer8
 44: #define PetscOffset integer8
 45: #else
 46: #define PetscOffset integer4
 47: #define PetscFortranAddr integer4
 48: #endif

 50: #if defined(PETSC_USE_64BIT_INDICES)
 51: #define PetscInt integer8
 52: #else
 53: #define PetscInt integer4
 54: #endif

 56: #if defined (PETSC_USE_FORTRANKIND)
 57: #define PetscFortranFloat real(kind=selected_real_kind(5))
 58: #define PetscFortranDouble real(kind=selected_real_kind(10))
 59: #define PetscFortranComplex complex(kind=selected_real_kind(10))
 60: #define PetscChar(a) character(len = a) ::
 61: #else
 62: #define PetscFortranFloat real*4
 63: #define PetscFortranDouble real*8
 64: #define PetscFortranComplex complex*16
 65: #define PetscChar(a) character*(a)
 66: #endif

 68: #if defined(PETSC_USE_COMPLEX)
 69: #define PETSC_SCALAR PETSC_COMPLEX
 70: #else
 71: #if defined(PETSC_USE_SINGLE)
 72: #define PETSC_SCALAR PETSC_FLOAT
 73: #else
 74: #define PETSC_SCALAR PETSC_DOUBLE
 75: #endif     
 76: #endif
 77: !
 78: !     Macro for templating between real and complex
 79: !
 80: #if defined(PETSC_USE_COMPLEX)
 81: #define PetscScalar PetscFortranComplex
 82: !
 83: ! F90 uses real(), conjg() when KIND parameter is used.
 84: !
 85: #if defined (PETSC_MISSING_DREAL)
 86: #define PetscRealPart(a) real(a)
 87: #define PetscConj(a) conjg(a)
 88: #define PetscImaginaryPart(a) aimg(a)
 89: #else
 90: #define PetscRealPart(a) dreal(a)
 91: #define PetscConj(a) dconjg(a)
 92: #define PetscImaginaryPart(a) daimg(a)
 93: #endif
 94: #define MPIU_SCALAR MPI_DOUBLE_COMPLEX
 95: #else
 96: #if defined (PETSC_USE_SINGLE)
 97: #define PetscScalar PetscFortranFloat
 98: #else
 99: #define PetscScalar PetscFortranDouble
100: #endif
101: #define PetscRealPart(a) a
102: #define PetscConj(a) a
103: #define PetscImaginaryPart(a) a
104: #define MPIU_SCALAR MPI_DOUBLE_PRECISION
105: #endif

107: #if defined (PETSC_USE_SINGLE)
108: #define PetscReal PetscFortranFloat
109: #else
110: #define PetscReal PetscFortranDouble
111: #endif

113: !
114: !    Allows the matrix Fortran Kernels to work with single precision
115: !    matrix data structures
116: !
117: #if defined(PETSC_USE_COMPLEX)
118: #define MatScalar PetscScalar 
119: #elif defined(PETSC_USE_MAT_SINGLE)
120: #define MatScalar real*4
121: #else
122: #define MatScalar PetscScalar
123: #endif
124: !
125: !     PetscLogDouble variables are used to contain double precision numbers
126: !     that are not used in the numerical computations, but rather in logging,
127: !     timing etc.
128: !
129: #define PetscObject PetscFortranAddr
130: #define PetscLogDouble PetscFortranDouble
131: !
132: !     Macros for error checking
133: !
134: #if defined(PETSC_USE_DEBUG)
135: #define SETERRQ(n,s,ierr) call MPI_Abort(PETSC_COMM_WORLD,n,ierr)
136: #define CHKERRQ(n) if (n .ne. 0) call MPI_Abort(PETSC_COMM_WORLD,n,n)
137: #define CHKMEMQ call chkmemfortran(__LINE__,__FILE__,ierr)
138: #else
139: #define SETERRQ(n,s,ierr)
140: #define CHKERRQ(n)
141: #define CHKMEMQ
142: #endif

144: #define PetscMatlabEngine PetscFortranAddr

146: #if defined(PETSC_STDCALL)
147: #define PETSC_DEC_ATTRIBUTES(A,B) DEC$ ATTRIBUTESC, ALIAS:B ::A
148: #else
149: #define PETSC_DEC_ATTRIBUTES(A,B)
150: #endif

152: #endif