stdlib: Support kotlin.math

licenses/third_party/boost_LICENSE.txt

@@ -0,0 +1,23 @@
+Boost Software License - Version 1.0 - August 17th, 2003
+
+Permission is hereby granted, free of charge, to any person or organization
+obtaining a copy of the software and accompanying documentation covered by
+this license (the "Software") to use, reproduce, display, distribute,
+execute, and transmit the Software, and to prepare derivative works of the
+Software, and to permit third-parties to whom the Software is furnished to
+do so, all subject to the following:
+
+The copyright notices in the Software and this entire statement, including
+the above license grant, this restriction and the following disclaimer,
+must be included in all copies of the Software, in whole or in part, and
+all derivative works of the Software, unless such copies or derivative
+works are solely in the form of machine-executable object code generated by
+a source language processor.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE, TITLE AND NON-INFRINGEMENT. IN NO EVENT
+SHALL THE COPYRIGHT HOLDERS OR ANYONE DISTRIBUTING THE SOFTWARE BE LIABLE
+FOR ANY DAMAGES OR OTHER LIABILITY, WHETHER IN CONTRACT, TORT OR OTHERWISE,
+ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
+DEALINGS IN THE SOFTWARE.

runtime/src/main/cpp/Exceptions.h

@@ -46,6 +46,8 @@ void ThrowArithmeticException();
 void ThrowNumberFormatException();
 // Throws out of memory error.
 void ThrowOutOfMemoryError();
+// Throws not implemented error.
+void ThrowNotImplementedError();
 // Prints out mesage of Throwable.
 void PrintThrowable(KRef);
 

runtime/src/main/cpp/Math.cpp

@@ -0,0 +1,358 @@
+/*
+ * Copyright 2010-2017 JetBrains s.r.o.
+ *
+ * Licensed under the Apache License, Version 2.0 (the "License");
+ * you may not use this file except in compliance with the License.
+ * You may obtain a copy of the License at
+ *
+ * http://www.apache.org/licenses/LICENSE-2.0
+ *
+ * Unless required by applicable law or agreed to in writing, software
+ * distributed under the License is distributed on an "AS IS" BASIS,
+ * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+ * See the License for the specific language governing permissions and
+ * limitations under the License.
+ */
+
+#include <float.h>
+#include <math.h>
+#include <stdlib.h>
+
+#include "Exceptions.h"
+#include "Types.h"
+
+#if (__MINGW32__ || __MINGW64__)
+#define KONAN_NEED_ASINH_ACOSH 1
+#else
+#define KONAN_NEED_ASINH_ACOSH 0
+#endif
+
+#ifdef KONAN_WASM
+#define KONAN_NO_MATH 1
+#else
+#define KONAN_NO_MATH 0
+#endif
+
+
+#ifdef KONAN_NEED_ASINH_ACOSH
+namespace {
+
+    // MinGW's implmenetation if asinh/acosh function returns NaN for large arguments so we use another implementation.
+    // Both implementations derived from boost special math functions and are also used by Kotlin/JVM.
+    // Copyright Eric Ford & Hubert Holin 2001.
+
+    constexpr KDouble LN2 = 0.69314718055994530942;
+    constexpr KDouble SQRT2 = 1.41421356237309504880;
+
+    KDouble taylor_2_bound = sqrt(DBL_EPSILON);
+    KDouble taylor_n_bound = sqrt(taylor_2_bound);
+    KDouble upper_taylor_2_bound = 1.0 / taylor_2_bound;
+    KDouble upper_taylor_n_bound = 1.0 / taylor_n_bound;
+
+    KDouble custom_asinh(KDouble x) {
+        if (x >= +taylor_n_bound) {
+            if (x > upper_taylor_n_bound) {
+                if (x > upper_taylor_2_bound) {
+                    // approximation by laurent series in 1/x at 0+ order from -1 to 0
+                    return log(x) + LN2;
+                } else {
+                    // approximation by laurent series in 1/x at 0+ order from -1 to 1
+                    return log(x * 2 + (1.0 / (x * 2)));
+                }
+            } else {
+                return log(x + sqrt(x * x + 1));
+            }
+        } else if (x <= -taylor_n_bound) {
+            return -custom_asinh(-x);
+        } else {
+            // approximation by taylor series in x at 0 up to order 2
+            KDouble result = x;
+            if (fabs(x) >= taylor_2_bound) {
+                // approximation by taylor series in x at 0 up to order 4
+                result -= (x * x * x) / 6;
+            }
+            return result;
+        }
+    }
+
+    KDouble custom_acosh(KDouble x) {
+        if (x < 1) {
+            return NAN;
+        } else if (x > upper_taylor_2_bound) {
+            // approximation by laurent series in 1/x at 0+ order from -1 to 0
+            return log(x) + LN2;
+        } else if (x - 1 >= taylor_n_bound) {
+            return log(x + sqrt(x * x - 1));
+        } else {
+            KDouble y = sqrt(x - 1);
+            // approximation by taylor series in y at 0 up to order 2
+            KDouble result = y;
+            if (y >= taylor_2_bound) {
+                // approximation by taylor series in y at 0 up to order 4
+                result -= (y * y * y) / 12;
+            }
+            return SQRT2 * result;
+        }
+    }
+}
+#endif
+
+extern "C" {
+
+#if !KONAN_NO_MATH // There is a platform math library.
+
+// region Double math.
+
+KDouble Kotlin_math_sin(KDouble x) { return sin(x); }
+KDouble Kotlin_math_cos(KDouble x) { return cos(x); }
+KDouble Kotlin_math_tan(KDouble x) { return tan(x); }
+KDouble Kotlin_math_asin(KDouble x) { return asin(x); }
+KDouble Kotlin_math_acos(KDouble x) { return acos(x); }
+KDouble Kotlin_math_atan(KDouble x) { return atan(x); }
+KDouble Kotlin_math_atan2(KDouble y, KDouble x) { return atan2(y, x); }
+
+KDouble Kotlin_math_sinh(KDouble x) { return sinh(x); }
+KDouble Kotlin_math_cosh(KDouble x) { return cosh(x); }
+KDouble Kotlin_math_tanh(KDouble x) { return tanh(x); }
+
+KDouble Kotlin_math_asinh(KDouble x) {
+#if (KONAN_NEED_ASINH_ACOSH)
+    return custom_asinh(x);
+#else
+    return asinh(x);
+#endif
+}
+
+KDouble Kotlin_math_acosh(KDouble x) {
+#if (KONAN_NEED_ASINH_ACOSH)
+    return custom_acosh(x);
+#else
+    return acosh(x);
+#endif
+}
+
+KDouble Kotlin_math_atanh(KDouble x) { return atanh(x); }
+
+KDouble Kotlin_math_hypot(KDouble x, KDouble y) { return hypot(x, y); }
+KDouble Kotlin_math_sqrt(KDouble x) { return sqrt(x); }
+KDouble Kotlin_math_exp(KDouble x) { return exp(x); }
+KDouble Kotlin_math_expm1(KDouble x) { return expm1(x); }
+
+KDouble Kotlin_math_ln(KDouble x) { return log(x); }
+KDouble Kotlin_math_log10(KDouble x) { return log10(x); }
+KDouble Kotlin_math_log2(KDouble x) { return log2(x); }
+KDouble Kotlin_math_ln1p(KDouble x) { return log1p(x); }
+
+KDouble Kotlin_math_ceil(KDouble x) { return ceil(x); }
+KDouble Kotlin_math_floor(KDouble x) { return floor(x); }
+KDouble Kotlin_math_round(KDouble x) { return rint(x); }
+
+KDouble Kotlin_math_abs(KDouble x) { return fabs(x); }
+
+// extensions
+
+KDouble Kotlin_math_Double_pow(KDouble thiz, KDouble x) {
+    if (isinf(x) && (thiz == 1.0 || thiz == -1.0)) return NAN; // Kotlin corner case
+    return pow(thiz, x);
+}
+
+KDouble Kotlin_math_Double_IEEErem(KDouble thiz, KDouble divisor) { return remainder(thiz, divisor); }
+KDouble Kotlin_math_Double_withSign(KDouble thiz, KDouble sign) { return copysign(thiz, sign); }
+
+KDouble Kotlin_math_Double_nextUp(KDouble thiz) { return nextafter(thiz, HUGE_VAL); }
+KDouble Kotlin_math_Double_nextDown(KDouble thiz) { return nextafter(thiz, -HUGE_VAL); }
+KDouble Kotlin_math_Double_nextTowards(KDouble thiz, KDouble to) {
+    return (thiz == to) ? to : nextafter(thiz, to);
+}
+
+KBoolean Kotlin_math_Double_signBit(KDouble thiz) { return signbit(thiz) != 0; }
+
+// endregion
+
+// region Float math.
+
+KFloat Kotlin_math_sinf(KFloat x) { return sinf(x); }
+KFloat Kotlin_math_cosf(KFloat x) { return cosf(x); }
+KFloat Kotlin_math_tanf(KFloat x) { return tanf(x); }
+KFloat Kotlin_math_asinf(KFloat x) { return asinf(x); }
+KFloat Kotlin_math_acosf(KFloat x) { return acosf(x); }
+KFloat Kotlin_math_atanf(KFloat x) { return atanf(x); }
+KFloat Kotlin_math_atan2f(KFloat y, KFloat x) { return atan2f(y, x); }
+
+KFloat Kotlin_math_sinhf(KFloat x) { return sinhf(x); }
+KFloat Kotlin_math_coshf(KFloat x) { return coshf(x); }
+KFloat Kotlin_math_tanhf(KFloat x) { return tanhf(x); }
+
+KFloat Kotlin_math_asinhf(KFloat x) {
+#if (KONAN_NEED_ASINH_ACOSH)
+    return (KFloat)custom_asinh((KDouble)x);
+#else
+    return asinhf(x);
+#endif
+}
+
+KFloat Kotlin_math_acoshf(KFloat x) {
+#if (KONAN_NEED_ASINH_ACOSH)
+    return (KFloat)custom_acosh((KDouble)x);
+#else
+    return acoshf(x);
+#endif
+}
+
+KFloat Kotlin_math_atanhf(KFloat x) { return atanhf(x); }
+
+KFloat Kotlin_math_hypotf(KFloat x, KFloat y) { return hypotf(x, y); }
+KFloat Kotlin_math_sqrtf(KFloat x) { return sqrtf(x); }
+KFloat Kotlin_math_expf(KFloat x) { return expf(x); }
+KFloat Kotlin_math_expm1f(KFloat x) { return expm1f(x); }
+
+KFloat Kotlin_math_lnf(KFloat x) { return logf(x); }
+KFloat Kotlin_math_log10f(KFloat x) { return log10f(x); }
+KFloat Kotlin_math_log2f(KFloat x) { return log2f(x); }
+KFloat Kotlin_math_ln1pf(KFloat x) { return log1pf(x); }
+
+KFloat Kotlin_math_ceilf(KFloat x) { return ceilf(x); }
+KFloat Kotlin_math_floorf(KFloat x) { return floorf(x); }
+KFloat Kotlin_math_roundf(KFloat x) { return rintf(x); }
+
+KFloat Kotlin_math_absf(KFloat x) { return fabsf(x); }
+
+// extensions
+
+KFloat Kotlin_math_Float_pow(KFloat thiz, KFloat x) {
+    if (isinf(x) && (thiz == 1.0 || thiz == -1.0)) return NAN; // Kotlin corner case
+    return powf(thiz, x);
+}
+
+KFloat Kotlin_math_Float_IEEErem(KFloat thiz, KFloat divisor) { return remainderf(thiz, divisor); }
+KFloat Kotlin_math_Float_withSign(KFloat thiz, KFloat sign) { return copysignf(thiz, sign); }
+
+KFloat Kotlin_math_Float_nextUp(KFloat thiz) { return nextafterf(thiz, HUGE_VALF); }
+KFloat Kotlin_math_Float_nextDown(KFloat thiz) { return nextafterf(thiz, -HUGE_VALF); }
+KFloat Kotlin_math_Float_nextTowards(KFloat thiz, KFloat to) {
+    return (thiz == to) ? to : nextafterf(thiz, to);
+}
+
+KBoolean Kotlin_math_Float_signBit(KFloat thiz) { return signbit(thiz) != 0; }
+
+// endregion
+
+// region Integer math.
+
+KInt Kotlin_math_absi(KInt x) { return labs(x); }
+KLong Kotlin_math_absl(KLong x) { return llabs(x); }
+
+// endregion
+
+#else // if there is no platform math library.
+
+// kotlin.math isn't supported for WASM.
+namespace {
+    void NotImplemented() { ThrowNotImplementedError(); }
+}
+
+KDouble Kotlin_math_sin(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_cos(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_tan(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_asin(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_acos(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_atan(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_atan2(KDouble y, KDouble x) { NotImplemented(); }
+
+KDouble Kotlin_math_sinh(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_cosh(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_tanh(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_asinh(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_acosh(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_atanh(KDouble x) { NotImplemented(); }
+
+KDouble Kotlin_math_hypot(KDouble x, KDouble y) { NotImplemented(); }
+KDouble Kotlin_math_sqrt(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_exp(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_expm1(KDouble x) { NotImplemented(); }
+
+KDouble Kotlin_math_ln(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_log10(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_log2(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_ln1p(KDouble x) { NotImplemented(); }
+
+KDouble Kotlin_math_ceil(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_floor(KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_round(KDouble x) { NotImplemented(); }
+
+KDouble Kotlin_math_abs(KDouble x) { NotImplemented(); }
+
+// extensions
+
+KDouble Kotlin_math_Double_pow(KDouble thiz, KDouble x) { NotImplemented(); }
+KDouble Kotlin_math_Double_IEEErem(KDouble thiz, KDouble divisor) { NotImplemented(); }
+KDouble Kotlin_math_Double_withSign(KDouble thiz, KDouble sign) { NotImplemented(); }
+
+KDouble Kotlin_math_Double_nextUp(KDouble thiz) { NotImplemented(); }
+KDouble Kotlin_math_Double_nextDown(KDouble thiz) { NotImplemented(); }
+KDouble Kotlin_math_Double_nextTowards(KDouble thiz, KDouble to) { NotImplemented(); }
+
+KBoolean Kotlin_math_Double_signBit(KDouble thiz) { NotImplemented(); }
+
+// endregion
+
+// region Float math.
+
+KFloat Kotlin_math_sinf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_cosf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_tanf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_asinf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_acosf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_atanf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_atan2f(KFloat y, KFloat x) { NotImplemented(); }
+
+KFloat Kotlin_math_sinhf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_coshf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_tanhf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_asinhf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_acoshf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_atanhf(KFloat x) { NotImplemented(); }
+
+KFloat Kotlin_math_hypotf(KFloat x, KFloat y) { NotImplemented(); }
+KFloat Kotlin_math_sqrtf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_expf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_expm1f(KFloat x) { NotImplemented(); }
+
+KFloat Kotlin_math_lnf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_log10f(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_log2f(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_ln1pf(KFloat x) { NotImplemented(); }
+
+KFloat Kotlin_math_ceilf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_floorf(KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_roundf(KFloat x) { NotImplemented(); }
+
+KFloat Kotlin_math_absf(KFloat x) { NotImplemented(); }
+
+// extensions
+
+KFloat Kotlin_math_Float_pow(KFloat thiz, KFloat x) { NotImplemented(); }
+KFloat Kotlin_math_Float_IEEErem(KFloat thiz, KFloat divisor) { NotImplemented(); }
+KFloat Kotlin_math_Float_withSign(KFloat thiz, KFloat sign) { NotImplemented(); }
+
+KFloat Kotlin_math_Float_nextUp(KFloat thiz) { NotImplemented(); }
+KFloat Kotlin_math_Float_nextDown(KFloat thiz) {  NotImplemented(); }
+KFloat Kotlin_math_Float_nextTowards(KFloat thiz, KFloat to) { NotImplemented(); }
+
+KBoolean Kotlin_math_Float_signBit(KFloat thiz) { NotImplemented(); }
+
+// endregion
+
+// region Integer math
+
+KInt Kotlin_math_absi(KInt x) { NotImplemented(); }
+KInt Kotlin_math_mini(KInt a, KInt b) { NotImplemented(); }
+KInt Kotlin_math_maxi(KInt a, KInt b) { NotImplemented(); }
+
+KLong Kotlin_math_absl(KLong x) { NotImplemented(); }
+KLong Kotlin_math_minl(KLong a, KLong b) { NotImplemented(); }
+KLong Kotlin_math_maxl(KLong a, KLong b) { NotImplemented(); }
+
+#endif // #if !KONAN_NO_MATH
+
+} // extern "C"

runtime/src/main/kotlin/konan/internal/RuntimeUtils.kt

@@ -61,6 +61,11 @@ fun ThrowUninitializedPropertyAccessException(): Nothing {
     throw UninitializedPropertyAccessException()
 }
 
+@ExportForCppRuntime
+internal fun ThrowNotImplementedError(): Nothing {
+    throw NotImplementedError("An operation is not implemented.")
+}
+
 @ExportForCppRuntime
 fun PrintThrowable(throwable: Throwable) {
     println(throwable)

runtime/src/main/kotlin/kotlin/Numbers.kt

@@ -68,7 +68,8 @@ public inline fun Double.toBits(): Long = if (isNaN()) Double.NaN.toRawBits() el
  * preserving `NaN` values exact layout.
  */
 @SinceKotlin("1.2")
-public fun Double.toRawBits(): Long = bits()
+@kotlin.internal.InlineOnly
+public inline fun Double.toRawBits(): Long = bits()
 
 /**
  * Returns the [Double] value corresponding to a given bit representation.
@@ -95,7 +96,8 @@ public inline fun Float.toBits(): Int = if (isNaN()) Float.NaN.toRawBits() else
  * preserving `NaN` values exact layout.
  */
 @SinceKotlin("1.2")
-public fun Float.toRawBits(): Int = bits()
+@kotlin.internal.InlineOnly
+public inline fun Float.toRawBits(): Int = bits()
 
 /**
  * Returns the [Float] value corresponding to a given bit representation.