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 PetscFortranLongDouble real(kind=selected_real_kind(16))
     60: #define PetscFortranComplex complex(kind=selected_real_kind(10))
     61: #define PetscChar(a) character(len = a) ::
     62: #else
     63: #define PetscFortranFloat real*4
     64: #define PetscFortranDouble real*8
     65: #define PetscFortranLongDouble real*16
     66: #define PetscFortranComplex complex*16
     67: #define PetscChar(a) character*(a)
     68: #endif

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

    115: #if defined (PETSC_USE_SINGLE)
    116: #define PetscReal PetscFortranFloat
    117: #elif defined(PETSC_USE_LONG_DOUBLE)
    118: #define PetscReal PetscFortranLongDouble
    119: #else
    120: #define PetscReal PetscFortranDouble
    121: #endif

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

    154: #define PetscMatlabEngine PetscFortranAddr

    156: #if defined(PETSC_STDCALL)
    157: #define PETSC_DEC_ATTRIBUTES(A,B) DEC$ ATTRIBUTESC, ALIAS:B ::A
    158: #else
    159: #define PETSC_DEC_ATTRIBUTES(A,B)
    160: #endif

    162: #endif