// fuzzy.cc -- template for fuzzy comparisons et al on floating point values
// $Header$
//
// jtutil::fuzzy::tolerance - comparison tolerance
// jtutil::fuzzy::EQ
// jtutil::fuzzy::is_integer
// jtutil::fuzzy::floor
// jtutil::fuzzy::ceiling
//
// ***** template instantiations *****
//

#include <algorithm>
#include "stdc.h"
#include "util.hh"

//******************************************************************************

namespace jtutil
	  {
template <typename fp>
bool fuzzy<fp>::EQ(fp x, fp y)
{
fp max_abs = std::max(jtutil::abs(x), jtutil::abs(y));
fp epsilon = std::max(tolerance, tolerance*max_abs);

return jtutil::abs(x-y) <= epsilon;
}
	  }	// namespace jtutil::

//******************************************************************************

namespace jtutil
	  {
template <typename fp>
bool fuzzy<fp>::is_integer(fp x)
{
int i = round<fp>::to_integer(x);
return EQ(x, fp(i));
}
	  }	// namespace jtutil::

//******************************************************************************

namespace jtutil
	  {
template <typename fp>
int fuzzy<fp>::floor(fp x)
{
return fuzzy<fp>::is_integer(x)
       ? round<fp>::to_integer(x)
       : round<fp>::floor(x);
}
	  }	// namespace jtutil::

//******************************************************************************

namespace jtutil
	  {
template <typename fp>
int fuzzy<fp>::ceiling(fp x)
{
return fuzzy<fp>::is_integer(x)
       ? round<fp>::to_integer(x)
       : round<fp>::ceiling(x);
}
	  }	// namespace jtutil::

//******************************************************************************
//******************************************************************************
//******************************************************************************

//
// ***** template instantiations *****
//

// instantiation for <float>
template class jtutil::fuzzy<float>;
float jtutil::fuzzy<float>::tolerance = 1.0e-5;	     // about 100 * FLT_EPSILON

// instantiations for <double>
template class jtutil::fuzzy<double>;
double jtutil::fuzzy<double>::tolerance = 1.0e-12;   // about 1e4 * DBL_EPSILON