--- /dev/null
+///////////////////////////////////////////////////////////////////////////
+//
+// Copyright (c) 2002, Industrial Light & Magic, a division of Lucas
+// Digital Ltd. LLC
+//
+// All rights reserved.
+//
+// Redistribution and use in source and binary forms, with or without
+// modification, are permitted provided that the following conditions are
+// met:
+// * Redistributions of source code must retain the above copyright
+// notice, this list of conditions and the following disclaimer.
+// * Redistributions in binary form must reproduce the above
+// copyright notice, this list of conditions and the following disclaimer
+// in the documentation and/or other materials provided with the
+// distribution.
+// * Neither the name of Industrial Light & Magic nor the names of
+// its contributors may be used to endorse or promote products derived
+// from this software without specific prior written permission.
+//
+// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+// "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+// LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+// A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
+// OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
+// SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
+// LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
+// DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
+// THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
+// (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+// OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+//
+///////////////////////////////////////////////////////////////////////////
+
+
+
+#ifndef INCLUDED_IMATHMATH_H
+#define INCLUDED_IMATHMATH_H
+
+//----------------------------------------------------------------------------
+//
+// ImathMath.h
+//
+// This file contains template functions which call the double-
+// precision math functions defined in math.h (sin(), sqrt(),
+// exp() etc.), with specializations that call the faster
+// single-precision versions (sinf(), sqrtf(), expf() etc.)
+// when appropriate.
+//
+// Example:
+//
+// double x = Math<double>::sqrt (3); // calls ::sqrt(double);
+// float y = Math<float>::sqrt (3); // calls ::sqrtf(float);
+//
+// When would I want to use this?
+//
+// You may be writing a template which needs to call some function
+// defined in math.h, for example to extract a square root, but you
+// don't know whether to call the single- or the double-precision
+// version of this function (sqrt() or sqrtf()):
+//
+// template <class T>
+// T
+// glorp (T x)
+// {
+// return sqrt (x + 1); // should call ::sqrtf(float)
+// } // if x is a float, but we
+// // don't know if it is
+//
+// Using the templates in this file, you can make sure that
+// the appropriate version of the math function is called:
+//
+// template <class T>
+// T
+// glorp (T x, T y)
+// {
+// return Math<T>::sqrt (x + 1); // calls ::sqrtf(float) if x
+// } // is a float, ::sqrt(double)
+// // otherwise
+//
+//----------------------------------------------------------------------------
+
+#include "ImathPlatform.h"
+#include <math.h>
+
+namespace Imath {
+
+
+template <class T>
+struct Math
+{
+ static T acos (T x) {return ::acos (double(x));}
+ static T asin (T x) {return ::asin (double(x));}
+ static T atan (T x) {return ::atan (double(x));}
+ static T atan2 (T x, T y) {return ::atan2 (double(x), double(y));}
+ static T cos (T x) {return ::cos (double(x));}
+ static T sin (T x) {return ::sin (double(x));}
+ static T tan (T x) {return ::tan (double(x));}
+ static T cosh (T x) {return ::cosh (double(x));}
+ static T sinh (T x) {return ::sinh (double(x));}
+ static T tanh (T x) {return ::tanh (double(x));}
+ static T exp (T x) {return ::exp (double(x));}
+ static T log (T x) {return ::log (double(x));}
+ static T log10 (T x) {return ::log10 (double(x));}
+ static T modf (T x, T *iptr)
+ {
+ double ival;
+ T rval( ::modf (double(x),&ival));
+ *iptr = ival;
+ return rval;
+ }
+ static T pow (T x, T y) {return ::pow (double(x), double(y));}
+ static T sqrt (T x) {return ::sqrt (double(x));}
+ static T ceil (T x) {return ::ceil (double(x));}
+ static T fabs (T x) {return ::fabs (double(x));}
+ static T floor (T x) {return ::floor (double(x));}
+ static T fmod (T x, T y) {return ::fmod (double(x), double(y));}
+ static T hypot (T x, T y) {return ::hypot (double(x), double(y));}
+};
+
+
+template <>
+struct Math<float>
+{
+ static float acos (float x) {return ::acosf (x);}
+ static float asin (float x) {return ::asinf (x);}
+ static float atan (float x) {return ::atanf (x);}
+ static float atan2 (float x, float y) {return ::atan2f (x, y);}
+ static float cos (float x) {return ::cosf (x);}
+ static float sin (float x) {return ::sinf (x);}
+ static float tan (float x) {return ::tanf (x);}
+ static float cosh (float x) {return ::coshf (x);}
+ static float sinh (float x) {return ::sinhf (x);}
+ static float tanh (float x) {return ::tanhf (x);}
+ static float exp (float x) {return ::expf (x);}
+ static float log (float x) {return ::logf (x);}
+ static float log10 (float x) {return ::log10f (x);}
+ static float modf (float x, float *y) {return ::modff (x, y);}
+ static float pow (float x, float y) {return ::powf (x, y);}
+ static float sqrt (float x) {return ::sqrtf (x);}
+ static float ceil (float x) {return ::ceilf (x);}
+ static float fabs (float x) {return ::fabsf (x);}
+ static float floor (float x) {return ::floorf (x);}
+ static float fmod (float x, float y) {return ::fmodf (x, y);}
+#if !defined(_MSC_VER)
+ static float hypot (float x, float y) {return ::hypotf (x, y);}
+#else
+ static float hypot (float x, float y) {return ::sqrtf(x*x + y*y);}
+#endif
+};
+
+
+//--------------------------------------------------------------------------
+// Compare two numbers and test if they are "approximately equal":
+//
+// equalWithAbsError (x1, x2, e)
+//
+// Returns true if x1 is the same as x2 with an absolute error of
+// no more than e,
+//
+// abs (x1 - x2) <= e
+//
+// equalWithRelError (x1, x2, e)
+//
+// Returns true if x1 is the same as x2 with an relative error of
+// no more than e,
+//
+// abs (x1 - x2) <= e * x1
+//
+//--------------------------------------------------------------------------
+
+template <class T>
+inline bool
+equalWithAbsError (T x1, T x2, T e)
+{
+ return ((x1 > x2)? x1 - x2: x2 - x1) <= e;
+}
+
+
+template <class T>
+inline bool
+equalWithRelError (T x1, T x2, T e)
+{
+ return ((x1 > x2)? x1 - x2: x2 - x1) <= e * ((x1 > 0)? x1: -x1);
+}
+
+
+
+} // namespace Imath
+
+#endif