Update to 2.0.0 tree from current Fremantle build

[opencv] / 3rdparty / include / OpenEXR / halfFunction.h

diff --git a/3rdparty/include/OpenEXR/halfFunction.h b/3rdparty/include/OpenEXR/halfFunction.h
new file mode 100644 (file)
index 0000000..1983880
--- /dev/null
+++ b/3rdparty/include/OpenEXR/halfFunction.h
@@ -0,0 +1,159 @@
+///////////////////////////////////////////////////////////////////////////
+//
+// 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.
+//
+///////////////////////////////////////////////////////////////////////////
+
+// Primary authors:
+//     Florian Kainz <kainz@ilm.com>
+//     Rod Bogart <rgb@ilm.com>
+
+
+//---------------------------------------------------------------------------
+//
+//     halfFunction<T> -- a class for fast evaluation
+//                        of half --> T functions
+//
+//     The constructor for a halfFunction object,
+//
+//         halfFunction (function,
+//                       domainMin, domainMax,
+//                       defaultValue,
+//                       posInfValue, negInfValue,
+//                       nanValue);
+//
+//     evaluates the function for all finite half values in the interval
+//     [domainMin, domainMax], and stores the results in a lookup table.
+//     For finite half values that are not in [domainMin, domainMax], the
+//     constructor stores defaultValue in the table.  For positive infinity,
+//     negative infinity and NANs, posInfValue, negInfValue and nanValue
+//     are stored in the table.
+//
+//     The tabulated function can then be evaluated quickly for arbitrary
+//     half values by calling the the halfFunction object's operator()
+//     method.
+//
+//     Example:
+//
+//         #include <math.h>
+//         #include <halfFunction.h>
+//
+//         halfFunction<half> hsin (sin);
+//
+//         halfFunction<half> hsqrt (sqrt,             // function
+//                                   0, HALF_MAX,      // domain
+//                                   half::qNan(),     // sqrt(x) for x < 0
+//                                   half::posInf(),   // sqrt(+inf)
+//                                   half::qNan(),     // sqrt(-inf)
+//                                   half::qNan());    // sqrt(nan)
+//
+//         half x = hsin (1);
+//         half y = hsqrt (3.5);
+//
+//---------------------------------------------------------------------------
+
+#ifndef _HALF_FUNCTION_H_
+#define _HALF_FUNCTION_H_
+
+#include <float.h>
+#include "half.h"
+
+
+template <class T>
+class halfFunction
+{
+  public:
+
+    //------------
+    // Constructor
+    //------------
+
+    template <class Function>
+    halfFunction (Function f,
+                 half domainMin = -HALF_MAX,
+                 half domainMax =  HALF_MAX,
+                 T defaultValue = 0,
+                 T posInfValue  = 0,
+                 T negInfValue  = 0,
+                 T nanValue     = 0);
+
+    //-----------
+    // Evaluation
+    //-----------
+
+    T          operator () (half x) const;
+
+  private:
+
+    T          _lut[1 << 16];
+};
+
+
+//---------------
+// Implementation
+//---------------
+
+template <class T>
+template <class Function>
+halfFunction<T>::halfFunction (Function f,
+                              half domainMin,
+                              half domainMax,
+                              T defaultValue,
+                              T posInfValue,
+                              T negInfValue,
+                              T nanValue)
+{
+    for (int i = 0; i < (1 << 16); i++)
+    {
+       half x;
+       x.setBits (i);
+
+       if (x.isNan())
+           _lut[i] = nanValue;
+       else if (x.isInfinity())
+           _lut[i] = x.isNegative()? negInfValue: posInfValue;
+       else if (x < domainMin || x > domainMax)
+           _lut[i] = defaultValue;
+       else
+           _lut[i] = f (x);
+    }
+}
+
+
+template <class T>
+inline T
+halfFunction<T>::operator () (half x) const
+{
+    return _lut[x.bits()];
+}
+
+
+#endif