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: #include "petscconf.h"
 13:  #include finclude/petscviewerdef.h
 14:  #include finclude/petscerrordef.h
 15:  #include finclude/petsclogdef.h
 16:  #include finclude/petscsysdef.h
 17:  #include finclude/petscdrawdef.h

 19: !
 20: ! The real*8,complex*16 notatiton is used so that the
 21: ! PETSc double/complex variables are not affected by
 22: ! compiler options like -r4,-r8, sometimes invoked
 23: ! by the user. NAG compiler does not like integer*4,real*8

 25: #if defined(PETSC_USE_FORTRANKIND)
 26: #define integer8 integer(kind=selected_int_kind(10))
 27: #define integer4 integer(kind=selected_int_kind(5))
 28: #define PetscTruth logical(kind=4)
 29: #else
 30: #define integer8 integer*8
 31: #define integer4 integer*4
 32: #define PetscTruth logical*4
 33: #endif

 35: #if (PETSC_SIZEOF_VOID_P == 8)
 36: #define PetscFortranAddr integer8
 37: #define PetscOffset integer8
 38: #else
 39: #define PetscOffset integer4
 40: #define PetscFortranAddr integer4
 41: #endif

 43: #if defined(PETSC_USE_64BIT_INDICES)
 44: #define PetscInt integer8
 45: #else
 46: #define PetscInt integer4
 47: #endif

 49: #if (PETSC_SIZEOF_INT == 4)
 50: #define PetscFortranInt integer4
 51: #elif (PETSC_SIZEOF_INT == 8)
 52: #define PetscFortranInt integer8
 53: #endif
 54: !
 55: #if (PETSC_SIZEOF_SIZE_T == 8)
 56: #define PetscSizeT integer8
 57: #else
 58: #define PetscSizeT integer4
 59: #endif
 60: !
 61: #if defined(PETSC_HAVE_MPIUNI)
 62: #define MPI_Comm PetscFortranInt
 63: #define PetscMPIInt PetscFortranInt
 64: #else
 65: #define MPI_Comm integer
 66: #define PetscMPIInt integer
 67: #endif
 68: !
 69: #define PetscEnum PetscFortranInt
 70: #define PetscErrorCode PetscFortranInt
 71: #define PetscCookie PetscFortranInt
 72: #define PetscLogEvent PetscFortranInt
 73: #define PetscLogStage PetscFortranInt
 74: #define PetscVoid PetscFortranAddr
 75: !
 76: #if defined(PETSC_FORTRAN_PETSCTRUTH_INT)
 77: #undef PetscTruth
 78: #define PetscTruth PetscEnum
 79: #endif
 80: !
 81: #define PetscDataType PetscEnum
 82: #define PetscFPTrap PetscEnum
 83: !
 84: #if defined (PETSC_USE_FORTRANKIND)
 85: #define PetscFortranFloat real(kind=selected_real_kind(5))
 86: #define PetscFortranDouble real(kind=selected_real_kind(10))
 87: #define PetscFortranLongDouble real(kind=selected_real_kind(16))
 88: #define PetscFortranComplex complex(kind=selected_real_kind(10))
 89: #define PetscChar(a) character(len = a) ::
 90: #else
 91: #define PetscFortranFloat real*4
 92: #define PetscFortranDouble real*8
 93: #define PetscFortranLongDouble real*16
 94: #define PetscFortranComplex complex*16
 95: #define PetscChar(a) character*(a)
 96: #endif

 98: #if defined(PETSC_USE_COMPLEX)
 99: #define PETSC_SCALAR PETSC_COMPLEX
100: #else
101: #if defined(PETSC_USE_SINGLE)
102: #define PETSC_SCALAR PETSC_FLOAT
103: #elif defined(PETSC_USE_LONG_DOUBLE)
104: #define PETSC_SCALAR PETSC_LONG_DOUBLE
105: #else
106: #define PETSC_SCALAR PETSC_DOUBLE
107: #endif     
108: #endif
109: !
110: !     Macro for templating between real and complex
111: !
112: #if defined(PETSC_USE_COMPLEX)
113: #define PetscScalar PetscFortranComplex
114: !
115: ! F90 uses real(), conjg() when KIND parameter is used.
116: !
117: #if defined (PETSC_MISSING_DREAL)
118: #define PetscRealPart(a) real(a)
119: #define PetscConj(a) conjg(a)
120: #define PetscImaginaryPart(a) aimg(a)
121: #else
122: #define PetscRealPart(a) dreal(a)
123: #define PetscConj(a) dconjg(a)
124: #define PetscImaginaryPart(a) daimg(a)
125: #endif
126: #else
127: #if defined (PETSC_USE_SINGLE)
128: #define PetscScalar PetscFortranFloat
129: #elif defined(PETSC_USE_LONG_DOUBLE)
130: #define PetscScalar PetscFortranLongDouble
131: #else
132: #define PetscScalar PetscFortranDouble
133: #endif
134: #define PetscRealPart(a) a
135: #define PetscConj(a) a
136: #define PetscImaginaryPart(a) a
137: #endif

139: #if defined (PETSC_USE_SINGLE)
140: #define PetscReal PetscFortranFloat
141: #elif defined(PETSC_USE_LONG_DOUBLE)
142: #define PetscReal PetscFortranLongDouble
143: #else
144: #define PetscReal PetscFortranDouble
145: #endif

147: !
148: !    Allows the matrix Fortran Kernels to work with single precision
149: !    matrix data structures
150: !
151: #if defined(PETSC_USE_COMPLEX)
152: #define MatScalar PetscScalar 
153: #elif defined(PETSC_USE_MAT_SINGLE)
154: #define MatScalar real*4
155: #else
156: #define MatScalar PetscScalar
157: #endif
158: !
159: !     PetscLogDouble variables are used to contain double precision numbers
160: !     that are not used in the numerical computations, but rather in logging,
161: !     timing etc.
162: !
163: #define PetscObject PetscFortranAddr
164: #define PetscLogDouble PetscFortranDouble
165: !
166: !     Macros for error checking
167: !
168: #if defined(PETSC_USE_ERRORCHECKING)
169: #define SETERRQ(n,s,ierr) call MPI_Abort(PETSC_COMM_WORLD,n,ierr)
170: #define CHKERRQ(n) if (n .ne. 0) call MPI_Abort(PETSC_COMM_WORLD,n,n)
171: #define CHKMEMQ call chkmemfortran(__LINE__,__FILE__,ierr)
172: #else
173: #define SETERRQ(n,s,ierr)
174: #define CHKERRQ(n)
175: #define CHKMEMQ
176: #endif

178: #define PetscMatlabEngine PetscFortranAddr

180: #if defined(PETSC_STDCALL)
181: #define PETSC_DEC_ATTRIBUTES(A,B) DEC$ ATTRIBUTESC, ALIAS:B ::A
182: #else
183: #define PETSC_DEC_ATTRIBUTES(A,B)
184: #endif

186: #if !defined(PetscFlush)
187: #if defined(PETSC_HAVE_FLUSH)
188: #define PetscFlush(a)    call flush(a)
189: #elif defined(PETSC_HAVE_FLUSH_)
190: #define PetscFlush(a)    call flush_(a)
191: #else
192: #define PetscFlush(a)
193: #endif
194: #endif


197: #endif