aboutsummaryrefslogtreecommitdiff
path: root/v_windows/v/vlib/math/sqrt.v
diff options
context:
space:
mode:
authorIndrajith K L2022-12-03 17:00:20 +0530
committerIndrajith K L2022-12-03 17:00:20 +0530
commitf5c4671bfbad96bf346bd7e9a21fc4317b4959df (patch)
tree2764fc62da58f2ba8da7ed341643fc359873142f /v_windows/v/vlib/math/sqrt.v
downloadcli-tools-windows-master.tar.gz
cli-tools-windows-master.tar.bz2
cli-tools-windows-master.zip
Adds most of the toolsHEADmaster
Diffstat (limited to 'v_windows/v/vlib/math/sqrt.v')
-rw-r--r--v_windows/v/vlib/math/sqrt.v37
1 files changed, 37 insertions, 0 deletions
diff --git a/v_windows/v/vlib/math/sqrt.v b/v_windows/v/vlib/math/sqrt.v
new file mode 100644
index 0000000..6513b79
--- /dev/null
+++ b/v_windows/v/vlib/math/sqrt.v
@@ -0,0 +1,37 @@
+module math
+
+// special cases are:
+// sqrt(+inf) = +inf
+// sqrt(±0) = ±0
+// sqrt(x < 0) = nan
+// sqrt(nan) = nan
+[inline]
+pub fn sqrt(a f64) f64 {
+ mut x := a
+ if x == 0.0 || is_nan(x) || is_inf(x, 1) {
+ return x
+ }
+ if x < 0.0 {
+ return nan()
+ }
+ z, ex := frexp(x)
+ w := x
+ // approximate square root of number between 0.5 and 1
+ // relative error of approximation = 7.47e-3
+ x = 4.173075996388649989089e-1 + 5.9016206709064458299663e-1 * z // adjust for odd powers of 2
+ if (ex & 1) != 0 {
+ x *= sqrt2
+ }
+ x = ldexp(x, ex >> 1)
+ // newton iterations
+ x = 0.5 * (x + w / x)
+ x = 0.5 * (x + w / x)
+ x = 0.5 * (x + w / x)
+ return x
+}
+
+// sqrtf calculates square-root of the provided value. (float32)
+[inline]
+pub fn sqrtf(a f32) f32 {
+ return f32(sqrt(a))
+}