hypot function (#118)

zpl-c · Aug 18, 2024 · 05c6cb0 · 05c6cb0
1 parent 5b69b50
commit 05c6cb0
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/code/header/math.h b/code/header/math.h
@@ -192,6 +192,7 @@ ZPL_DEF zpl_f32 zpl_rsqrt(zpl_f32 a);
 ZPL_DEF zpl_f32 zpl_quake_rsqrt(zpl_f32 a); /* NOTE: It's probably better to use 1.0f/zpl_sqrt(a)
                                      * And for simd, there is usually isqrt functions too!
                                      */
+ZPL_DEF zpl_f32 zpl_hypot(zpl_f32 x, zpl_f32 y);
 ZPL_DEF zpl_f32 zpl_sin(zpl_f32 radians);
 ZPL_DEF zpl_f32 zpl_cos(zpl_f32 radians);
 ZPL_DEF zpl_f32 zpl_tan(zpl_f32 radians);

diff --git a/code/source/math.c b/code/source/math.c
@@ -207,7 +207,6 @@ zpl_f32 zpl_quake_rsqrt(zpl_f32 a) {
             r *= f;
             return flipped ? 1.0f / r : r;
         }
-
 #    else
 
         zpl_f32 zpl_rsqrt(zpl_f32 a)              { return 1.0f / __builtin_sqrt(a); }
@@ -225,6 +224,7 @@ zpl_f32 zpl_quake_rsqrt(zpl_f32 a) {
 
         // TODO: Should this be zpl_exp(y * zpl_log(x)) ???
         zpl_f32 zpl_pow(zpl_f32 x, zpl_f32 y) { return __builtin_powf(x, y); }
+        zpl_f32 zpl_hypot(zpl_f32 x, zpl_f32 y){ return __builtin_sqrt(zpl_square(x) + zpl_square(y)); }
 
 #    endif
 #else
@@ -241,6 +241,7 @@ zpl_f32 zpl_quake_rsqrt(zpl_f32 a) {
     zpl_f32 zpl_exp(zpl_f32 x)            { return expf(x); }
     zpl_f32 zpl_log(zpl_f32 x)            { return logf(x); }
     zpl_f32 zpl_pow(zpl_f32 x, zpl_f32 y) { return powf(x, y); }
+    zpl_f32 zpl_hypot(zpl_f32 x, zpl_f32 y) { return sqrtf(zpl_square(x) + zpl_square(y)); }
 #endif
 
 zpl_f32 zpl_exp2(zpl_f32 x) { return zpl_exp(ZPL_LOG_TWO * x); }