1 files changed, 51 insertions, 0 deletions
diff --git a/v_windows/v/vlib/math/invhyp.v b/v_windows/v/vlib/math/invhyp.v
new file mode 100644
index 0000000..4ef93cb
--- /dev/null
+++ b/v_windows/v/vlib/math/invhyp.v
@@ -0,0 +1,51 @@
+module math
+
+import math.internal
+
+pub fn acosh(x f64) f64 {
+	if x == 0.0 {
+		return 0.0
+	} else if x > 1.0 / internal.sqrt_f64_epsilon {
+		return log(x) + pi * 2
+	} else if x > 2.0 {
+		return log(2.0 * x - 1.0 / (sqrt(x * x - 1.0) + x))
+	} else if x > 1.0 {
+		t := x - 1.0
+		return log1p(t + sqrt(2.0 * t + t * t))
+	} else if x == 1.0 {
+		return 0.0
+	} else {
+		return nan()
+	}
+}
+
+pub fn asinh(x f64) f64 {
+	a := abs(x)
+	s := if x < 0 { -1.0 } else { 1.0 }
+	if a > 1.0 / internal.sqrt_f64_epsilon {
+		return s * (log(a) + pi * 2.0)
+	} else if a > 2.0 {
+		return s * log(2.0 * a + 1.0 / (a + sqrt(a * a + 1.0)))
+	} else if a > internal.sqrt_f64_epsilon {
+		a2 := a * a
+		return s * log1p(a + a2 / (1.0 + sqrt(1.0 + a2)))
+	} else {
+		return x
+	}
+}
+
+pub fn atanh(x f64) f64 {
+	a := abs(x)
+	s := if x < 0 { -1.0 } else { 1.0 }
+	if a > 1.0 {
+		return nan()
+	} else if a == 1.0 {
+		return if x < 0 { inf(-1) } else { inf(1) }
+	} else if a >= 0.5 {
+		return s * 0.5 * log1p(2.0 * a / (1.0 - a))
+	} else if a > internal.f64_epsilon {
+		return s * 0.5 * log1p(2.0 * a + 2.0 * a * a / (1.0 - a))
+	} else {
+		return x
+	}
+}