Changeset 16905 for trunk/HeuristicLab.ExtLibs/HeuristicLab.NativeInterpreter

Timestamp:

05/07/19 14:20:33 (6 years ago)

Author:

gkronber

Message:

#2915: Corrected calculation of cuberoot function in all interpreters and updated all formatters.
I tested a run with all interpreters and got the same results with all of them.

Location:

trunk/HeuristicLab.ExtLibs/HeuristicLab.NativeInterpreter/0.1

Files:

• HeuristicLab.NativeInterpreter-0.1/hl-native-interpreter-msvc-x64.dll (modified) (previous)
• HeuristicLab.NativeInterpreter-0.1/hl-native-interpreter-msvc-x86.dll (modified) (previous)
• NativeInterpreter-0.1/src/interpreter.h (modified) (2 diffs)
• NativeInterpreter-0.1/src/vector_operations.h (modified) (4 diffs)

trunk/HeuristicLab.ExtLibs/HeuristicLab.NativeInterpreter/0.1/NativeInterpreter-0.1/src/interpreter.h

@@ -118 +118 @@
             case OpCodes::CubeRoot:
                 {
-                    in.value = hl_pow(code[in.childIndex].value, 1./3.);
+                    in.value = hl_cbrt(code[in.childIndex].value);
                     break;
                 }
@@ -272 +272 @@
             case OpCodes::CubeRoot:
                 {
-                    pow(in.buf, code[in.childIndex].buf, 1./3.);
+                    cbrt(in.buf, code[in.childIndex].buf);
                     break;
                 }

trunk/HeuristicLab.ExtLibs/HeuristicLab.NativeInterpreter/0.1/NativeInterpreter-0.1/src/vector_operations.h

@@ -17 +17 @@
 #define hl_sqrt vdt::fast_sqrt
 #define hl_pow vdt::fast_pow
+#define hl_cbrt std::cbrt
 #define hl_round vdt::fast_round
 #define hl_inv vdt::fast_inv
@@ -28 +29 @@
 #define hl_sqrt std::sqrt
 #define hl_pow std::pow
+#define hl_cbrt std::cbrt
 #define hl_round std::round
 #define hl_inv(x) 1. / x;
@@ -56 +58 @@
 inline void pow(double* __restrict a, double const * __restrict b) noexcept { FOR(i) a[i] = hl_pow(a[i], hl_round(b[i])); };
 inline void root(double* __restrict a, double const * __restrict b) noexcept { FOR(i) a[i] = hl_pow(a[i], 1. / hl_round(b[i])); };
+inline void cbrt(double* __restrict a, double const * __restrict b) noexcept { FOR(i) a[i] = hl_cbrt(b[i]); };
 inline void square(double* __restrict a, double const * __restrict b) noexcept { FOR(i) a[i] = hl_pow(b[i], 2.); };
 inline void inv(double* __restrict a, double const * __restrict b) noexcept { FOR(i) a[i] = hl_inv(b[i]); }
@@ -81 +84 @@
 inline void round(double* __restrict a) noexcept { FOR(i) a[i] = hl_round(a[i]); }
 inline void square(double* __restrict a) noexcept { FOR(i) a[i] = hl_pow(a[i], 2.); }
+inline void cbrt(double* __restrict a) noexcept { FOR(i) a[i] = hl_cbrt(a[i]); }
 
 #undef FOR