diff --git a/src/MNH/gamma.f90 b/src/MNH/gamma.f90 index 763042054b9d5526e09a1ef81b9a10723a165cfc..626aab326fb079165592123ace4e2c7c263f6f8d 100644 --- a/src/MNH/gamma.f90 +++ b/src/MNH/gamma.f90 @@ -72,7 +72,7 @@ END MODULE MODI_GAMMA !* 0. DECLARATIONS ! ------------ ! -#ifdef MNH_BITREP +#if defined(MNH_BITREP) || defined(MNH_BITREP_OMP) USE MODI_BITREP #endif ! @@ -119,7 +119,7 @@ ELSE END IF ZY = ZX ZTMP = ZX + 5.5 -#ifndef MNH_BITREP +#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP) ZTMP = (ZX + 0.5) * ALOG(ZTMP) - ZTMP #else ZTMP = (ZX + 0.5) * BR_LOG(ZTMP) - ZTMP @@ -132,13 +132,13 @@ DO JJ = 1, 6 END DO ! IF (PX .LT. 0.) THEN -#ifndef MNH_BITREP +#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP) PGAMMA = ZPI / SIN(ZPI*PX) / EXP(ZTMP + ALOG(ZSTP*ZSER/ZX)) #else PGAMMA = ZPI / SIN(ZPI*PX) / BR_EXP(ZTMP + BR_LOG(ZSTP*ZSER/ZX)) #endif ELSE -#ifndef MNH_BITREP +#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP) PGAMMA = EXP(ZTMP + ALOG(ZSTP*ZSER/ZX)) #else PGAMMA = BR_EXP(ZTMP + BR_LOG(ZSTP*ZSER/ZX)) @@ -197,7 +197,7 @@ END FUNCTION GAMMA_X0D !* 0. DECLARATIONS ! ------------ ! -#ifdef MNH_BITREP +#if defined(MNH_BITREP) || defined(MNH_BITREP_OMP) USE MODI_BITREP #endif ! @@ -235,7 +235,7 @@ WHERE ( PX(:)<0.0 ) END WHERE ZY(:) = ZX(:) ZTMP(:) = ZX(:) + 5.5 -#ifndef MNH_BITREP +#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP) ZTMP(:) = (ZX(:) + 0.5)*ALOG(ZTMP(:)) - ZTMP(:) #else ZTMP(:) = (ZX(:) + 0.5)*BR_LOG(ZTMP(:)) - ZTMP(:) @@ -247,7 +247,7 @@ DO JJ = 1 , 6 ZSER(:) = ZSER(:) + ZCOEF(JJ)/ZY(:) END DO ! -#ifndef MNH_BITREP +#if !defined(MNH_BITREP) && !defined(MNH_BITREP_OMP) PGAMMA(:) = EXP( ZTMP(:) + ALOG( ZSTP*ZSER(:)/ZX(:) ) ) #else PGAMMA(:) = BR_EXP( ZTMP(:) + BR_LOG( ZSTP*ZSER(:)/ZX(:) ) )