From 885f5974cdf65b59415837ae97f5a14ef1350670 Mon Sep 17 00:00:00 2001
From: Kacper <kacper@mail.openlinux.dev>
Date: Tue, 9 Dec 2025 19:20:15 +0100
Subject: feat: add gzip and new headers

---
 lib/libm/csqrtf.c | 55 +++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 55 insertions(+)
 create mode 100644 lib/libm/csqrtf.c

(limited to 'lib/libm/csqrtf.c')

diff --git a/lib/libm/csqrtf.c b/lib/libm/csqrtf.c
new file mode 100644
index 00000000..4ea894ab
--- /dev/null
+++ b/lib/libm/csqrtf.c
@@ -0,0 +1,55 @@
+#include "__complex.h"
+
+float complex csqrtf(float complex z)
+{
+	float complex w;
+	float x, y, r, t, scale;
+
+	x = crealf(z);
+	y = cimagf(z);
+	if (y == 0.0f) {
+		if (x < 0.0f) {
+			return 0.0f + sqrtf(-x) * I;
+		} else if (x == 0.0f) {
+			return (0.0f + y * I);
+		} else {
+			return sqrtf(x) + y * I;
+		}
+	}
+
+	if (x == 0.0f) {
+		r = fabsf(y);
+		r = sqrtf(0.5f * r);
+		if (y > 0)
+			return r + r * I;
+		else
+			return r - r * I;
+	}
+
+	if ((fabsf(x) > 4.0f) || (fabsf(y) > 4.0f)) {
+		x *= 0.25f;
+		y *= 0.25f;
+		scale = 2.0f;
+	} else {
+		x *= 6.7108864e7f; /* 2^26 */
+		y *= 6.7108864e7f;
+		scale = 1.220703125e-4f; /* 2^-13 */
+	}
+
+	w = x + y * I;
+	r = cabsf(w);
+	if (x > 0) {
+		t = sqrtf(0.5f * r + 0.5f * x);
+		r = scale * fabsf((0.5f * y) / t);
+		t *= scale;
+	} else {
+		r = sqrtf(0.5f * r - 0.5f * x);
+		t = scale * fabsf((0.5f * y) / r);
+		r *= scale;
+	}
+
+	if (y < 0)
+		return t - r * I;
+	else
+		return t + r * I;
+}
-- 
cgit v1.2.3