feat: add C implementation

---
type: pre_commit_static_analysis_report
description: Results of running static analysis checks when committing changes.
report:
  - task: lint_filenames
    status: passed
  - task: lint_editorconfig
    status: passed
  - task: lint_markdown
    status: na
  - task: lint_package_json
    status: na
  - task: lint_repl_help
    status: na
  - task: lint_javascript_src
    status: passed
  - task: lint_javascript_cli
    status: na
  - task: lint_javascript_examples
    status: na
  - task: lint_javascript_tests
    status: passed
  - task: lint_javascript_benchmarks
    status: na
  - task: lint_python
    status: na
  - task: lint_r
    status: na
  - task: lint_c_src
    status: passed
  - task: lint_c_examples
    status: na
  - task: lint_c_benchmarks
    status: na
  - task: lint_c_tests_fixtures
    status: na
  - task: lint_shell
    status: na
  - task: lint_typescript_declarations
    status: na
  - task: lint_typescript_tests
    status: na
  - task: lint_license_headers
    status: passed
---

Author: anandkaranubc

lib/node_modules/@stdlib/math/base/special/powf/include/stdlib/math/base/special/powf.h

@@ -0,0 +1,38 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2025 The Stdlib Authors.
+*
+* 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.
+*/
+
+#ifndef STDLIB_MATH_BASE_SPECIAL_POWF_H
+#define STDLIB_MATH_BASE_SPECIAL_POWF_H
+
+/*
+* If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
+*/
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+/**
+* Evaluates the exponential function for a single-precision floating-point number.
+*/
+float stdlib_base_powf( const float x, const float y );
+
+#ifdef __cplusplus
+}
+#endif
+
+#endif // !STDLIB_MATH_BASE_SPECIAL_POWF_H

lib/node_modules/@stdlib/math/base/special/powf/lib/log2axf.js

@@ -96,7 +96,7 @@ var THREE = f32( 3.0 );
 * @private
 * @param {Array} out - output array
 * @param {number} ax - absolute value of `x`
-* @param {number} ahx - high word of `ax`
+* @param {number} ahx - word of `ax`
 * @returns {Array} output array containing a tuple comprised of high and low parts
 *
 * @example

lib/node_modules/@stdlib/math/base/special/powf/lib/logxf.js

@@ -54,6 +54,8 @@ var INV_LN2_HI = f32( 1.4426879883e+00 ); // 0x3fb8aa00
 // Low: 1/LN2
 var INV_LN2_LO = f32( 7.0526075433e-06 ); // 0x36eca570
 
+var ONE = f32( 1.0 );
+
 
 // MAIN //
 
@@ -67,7 +69,7 @@ var INV_LN2_LO = f32( 7.0526075433e-06 ); // 0x36eca570
 *
 * @example
 * var t = logxf( [ 0.0, 0.0 ], 9.0 );
-* // Expected output: [ -1265.7236328125, -0.0008163940840404393 ]
+* // returns [ -1265.5, ~-0.224 ]
 */
 function logxf( out, ax ) {
 	var tmp;
@@ -78,7 +80,7 @@ function logxf( out, ax ) {
 	var u;
 	var v;
 
-	t = f32( ax - 1.0 ); // `t` has `20` trailing zeros
+	t = f32( ax - ONE ); // `t` has `20` trailing zeros
 	w = f32( f32( t * t ) * polyvalW( t ) );
 	u = f32( INV_LN2_HI * t ); // `INV_LN2_HI` has `16` significant bits
 	v = f32( f32( t*INV_LN2_LO ) - f32( w*INV_LN2 ) );

lib/node_modules/@stdlib/math/base/special/powf/lib/main.js

@@ -60,7 +60,6 @@ var HALF = f32( 0.5 );
 var ONE = f32( 1.0 );
 var TWO = f32( 2.0 );
 var THREE = f32( 3.0 );
-var FOUR = f32( 4.0 );
 var NEG_ZERO = f32( -0.0 );
 var NEG_HALF = f32( -0.5 );
 var NEG_ONE = f32( -1.0 );
@@ -141,9 +140,9 @@ function powf( x, y ) {
 	var ahx;  // absolute value word `x`
 	var ahy;  // absolute value word `y`
 	var ax;   // absolute value `x`
-	var hx;   // word representation of `x	`
+	var hx;   // word representation of `x`
 	var hy;   // word representation of `y`
-	var sn;
+	var sn;   // sign of the result
 	var y1;
 	var hp;
 	var lp;
@@ -183,10 +182,6 @@ function powf( x, y ) {
 	if ( y === THREE ) {
 		return f32( f32( x * x ) * x );
 	}
-	if ( y === FOUR ) {
-		x = f32( x * x );
-		return f32( x * x );
-	}
 	if ( isInfinitef( y ) ) { // y is +-inf
 		return yIsInfinitef( x, y );
 	}

lib/node_modules/@stdlib/math/base/special/powf/lib/pow2f.js

@@ -83,7 +83,6 @@ var TWO = f32( 2.0 );
 * // returns ~0.79
 */
 function pow2f( j, hp, lp ) {
-	var ahs;
 	var tmp;
 	var t1;
 	var t;
@@ -100,21 +99,21 @@ function pow2f( j, hp, lp ) {
 	k = ((i>>NUM_SIGNIFICAND_BITS) - BIAS)|0; // asm type annotation
 	n = 0;
 
-	// `|z| > 0.5`, set `n = [z+0.5]`
+	// `|z| > 0.5`, set `n = z+0.5`
 	if ( i > HALF_WORD ) {
 		n = (j + (MIN_NORM_WORD>>(k+1)))|0; // asm type annotation
 		k = (((n & ABS_MASK)>>NUM_SIGNIFICAND_BITS) - BIAS)|0; // new k for n
 		tmp = ((n & ~(SIGNIFICAND_MASK >> k)))|0; // asm type annotation
 		t = fromWordf( tmp );
-		n = (((n & SIGNIFICAND_MASK)|MIN_NORM_WORD) >> (NUM_SIGNIFICAND_BITS-k))>>>0; // eslint-disable-line max-len
+		n = (((n & SIGNIFICAND_MASK)|MIN_NORM_WORD) >> (NUM_SIGNIFICAND_BITS-k))|0; // eslint-disable-line max-len
 		if ( j < 0 ) {
 			n = -n;
 		}
 		hp = f32( hp - t );
 	}
 	t = f32( lp + hp );
-	ahs = toWordf( t );
-	t = fromWordf( ahs & TRUNC_MASK_15 );
+	tmp = toWordf( t )|0; // asm type annotation
+	t = fromWordf( tmp & TRUNC_MASK_15 );
 	u = f32( t * LN2_HI );
 	v = f32( f32( f32( lp - f32(t-hp) )*LN2 ) + f32( t*LN2_LO ) );
 	z = f32( u + v );
@@ -124,7 +123,7 @@ function pow2f( j, hp, lp ) {
 	r = f32( f32( f32(z*t1) / f32(t1-TWO) ) - f32( w + f32(z*w) ) );
 	z = f32( ONE - f32(r - z) );
 	j = toWordf( z ) |0; // asm type annotation
-	j += ( (n>>>0) << NUM_SIGNIFICAND_BITS) |0; // asm type annotation
+	j += ( (n>>>0) << NUM_SIGNIFICAND_BITS ) |0; // asm type annotation
 
 	// Check for subnormal output...
 	if ( (j>>NUM_SIGNIFICAND_BITS) <= 0 ) {

lib/node_modules/@stdlib/math/base/special/powf/src/Makefile

@@ -1,7 +1,7 @@
 #/
 # @license Apache-2.0
 #
-# Copyright (c) 2024 The Stdlib Authors.
+# Copyright (c) 2025 The Stdlib Authors.
 #
 # Licensed under the Apache License, Version 2.0 (the "License");
 # you may not use this file except in compliance with the License.

lib/node_modules/@stdlib/math/base/special/powf/src/addon.c

@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2024 The Stdlib Authors.
+* Copyright (c) 2025 The Stdlib Authors.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
@@ -16,7 +16,7 @@
 * limitations under the License.
 */
 
-#include "stdlib/math/base/special/pow.h"
+#include "stdlib/math/base/special/powf.h"
 #include "stdlib/math/base/napi/binary.h"
 
-STDLIB_MATH_BASE_NAPI_MODULE_DD_D( stdlib_base_pow )
+STDLIB_MATH_BASE_NAPI_MODULE_FF_F( stdlib_base_powf )