G4Log.hh File Reference

#include <limits>
#include <stdint.h>
#include "G4Types.hh"

Include dependency graph for G4Log.hh:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
union	G4LogConsts::ieee754

Namespaces
	G4LogConsts

Functions
G4double	G4LogConsts::get_log_px (const G4double x)
G4double	G4LogConsts::get_log_qx (const G4double x)
uint64_t	G4LogConsts::dp2uint64 (G4double x)
G4double	G4LogConsts::uint642dp (uint64_t ll)
G4float	G4LogConsts::uint322sp (G4int x)
uint32_t	G4LogConsts::sp2uint32 (G4float x)
G4double	G4LogConsts::getMantExponent (const G4double x, G4double &fe)
	Like frexp but vectorising and the exponent is a double. More...
G4float	G4LogConsts::getMantExponentf (const G4float x, G4float &fe)
	Like frexp but vectorising and the exponent is a float. More...
G4double	G4Log (G4double x)
G4float	G4LogConsts::get_log_poly (const G4float x)
G4float	G4Logf (G4float x)
void	logv (const uint32_t size, G4double const *__restrict__ iarray, G4double *__restrict__ oarray)
void	G4Logv (const uint32_t size, G4double const *__restrict__ iarray, G4double *__restrict__ oarray)
void	logfv (const uint32_t size, G4float const *__restrict__ iarray, G4float *__restrict__ oarray)
void	G4Logfv (const uint32_t size, G4float const *__restrict__ iarray, G4float *__restrict__ oarray)

Variables
const G4double	G4LogConsts::LOG_UPPER_LIMIT = 1e307
const G4double	G4LogConsts::LOG_LOWER_LIMIT = 0
const G4double	G4LogConsts::SQRTH = 0.70710678118654752440
const G4float	G4LogConsts::MAXNUMF = 3.4028234663852885981170418348451692544e38f
const G4float	G4LogConsts::LOGF_UPPER_LIMIT = MAXNUMF
const G4float	G4LogConsts::LOGF_LOWER_LIMIT = 0
const G4float	G4LogConsts::PX1logf = 7.0376836292E-2f
const G4float	G4LogConsts::PX2logf = -1.1514610310E-1f
const G4float	G4LogConsts::PX3logf = 1.1676998740E-1f
const G4float	G4LogConsts::PX4logf = -1.2420140846E-1f
const G4float	G4LogConsts::PX5logf = 1.4249322787E-1f
const G4float	G4LogConsts::PX6logf = -1.6668057665E-1f
const G4float	G4LogConsts::PX7logf = 2.0000714765E-1f
const G4float	G4LogConsts::PX8logf = -2.4999993993E-1f
const G4float	G4LogConsts::PX9logf = 3.3333331174E-1f
const G4float	G4LogConsts::SQRTHF = 0.707106781186547524f

Function Documentation

G4double G4Log ( G4double  x )

inline

Definition at line 230 of file G4Log.hh.

{
        const G4double original_x = x;

        /* separate mantissa from exponent */
        G4double fe;
        x = G4LogConsts::getMantExponent(x,fe);

        // blending
        x > G4LogConsts::SQRTH? fe+=1. : x+=x ;
        x -= 1.0;

        /* rational form */
        G4double px =  G4LogConsts::get_log_px(x);

        //for the final formula
        const G4double x2 = x*x;
        px *= x;
        px *= x2;

        const G4double qx = G4LogConsts::get_log_qx(x);

        G4double res = px / qx ;

        res -= fe * 2.121944400546905827679e-4;
        res -= 0.5 * x2  ;

        res = x + res;
        res += fe * 0.693359375;

        if (original_x > G4LogConsts::LOG_UPPER_LIMIT)
                res = std::numeric_limits<G4double>::infinity();
        if (original_x < G4LogConsts::LOG_LOWER_LIMIT) // THIS IS NAN!
                res =  - std::numeric_limits<G4double>::quiet_NaN();

        return res;
}

Here is the call graph for this function:

G4float G4Logf ( G4float  x )

inline

Definition at line 311 of file G4Log.hh.

{
        const G4float original_x = x;

        G4float fe;
        x = G4LogConsts::getMantExponentf( x, fe);

        x > G4LogConsts::SQRTHF? fe+=1.f : x+=x ;
        x -= 1.0f;

        const G4float x2 = x*x;

        G4float res = G4LogConsts::get_log_poly(x);
        res *= x2*x;

        res += -2.12194440e-4f * fe;
        res +=  -0.5f * x2;

        res= x + res;

        res += 0.693359375f * fe;

        if (original_x > G4LogConsts::LOGF_UPPER_LIMIT)
                res = std::numeric_limits<G4float>::infinity();
        if (original_x < G4LogConsts::LOGF_LOWER_LIMIT)
                res = -std::numeric_limits<G4float>::quiet_NaN();

        return res;
}

Here is the call graph for this function:

void G4Logfv	(	const uint32_t	size,
		G4float const *__restrict__	iarray,
		G4float *__restrict__	oarray
	)

void G4Logv	(	const uint32_t	size,
		G4double const *__restrict__	iarray,
		G4double *__restrict__	oarray
	)

void logfv	(	const uint32_t	size,
		G4float const *__restrict__	iarray,
		G4float *__restrict__	oarray
	)

void logv	(	const uint32_t	size,
		G4double const *__restrict__	iarray,
		G4double *__restrict__	oarray
	)