safeLogRatio.cxx File Reference

#include "TrigBjetHypo/safeLogRatio.h"
#include <cmath>
#include <limits>

Include dependency graph for safeLogRatio.cxx:

Go to the source code of this file.

Functions
float	safeLogRatio (float num, float denom)

Function Documentation

safeLogRatio()

float safeLogRatio	(	float	num,
		float	denom
	)

Definition at line 10 of file safeLogRatio.cxx.

                                           {
  // ep(silon) is the smallest non-subnormal number
  // the recyprocal of this should not overflow
  float ep = std::numeric_limits<float>::min();
  float ratio = (std::abs(denom) < ep ? INFINITY : num / denom);
  return std::abs(ratio) < ep ? -INFINITY : std::log( ratio );
}