@stdlib/math-base-special-round2 0.2.0 → 0.3.0

package/README.md

@@ -124,6 +124,91 @@ for ( i = 0; i < 100; i++ ) {
 
 <!-- /.examples -->
 
+<!-- C interface documentation. -->
+
+* * *
+
+<section class="c">
+
+## C APIs
+
+<!-- Section to include introductory text. Make sure to keep an empty line after the intro `section` element and another before the `/section` close. -->
+
+<section class="intro">
+
+</section>
+
+<!-- /.intro -->
+
+<!-- C usage documentation. -->
+
+<section class="usage">
+
+### Usage
+
+```c
+#include "stdlib/math/base/special/round2.h"
+```
+
+#### stdlib_base_round2( x )
+
+Rounds a `numeric` value to the nearest power of two on a linear scale.
+
+```c
+double out = stdlib_base_round2( -4.2 );
+// returns -4.0
+```
+
+The function accepts the following arguments:
+
+-   **x**: `[in] double` input value.
+
+```c
+double stdlib_base_round2( const double x );
+```
+
+</section>
+
+<!-- /.usage -->
+
+<!-- C API usage notes. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
+
+<section class="notes">
+
+</section>
+
+<!-- /.notes -->
+
+<!-- C API usage examples. -->
+
+<section class="examples">
+
+### Examples
+
+```c
+#include "stdlib/math/base/special/round2.h"
+#include <stdio.h>
+
+int main( void ) {
+    const double x[] = { -5.0, -3.89, -2.78, -1.67, -0.56, 0.56, 1.67, 2.78, 3.89, 5.0 };
+    
+    double v;
+    int i;
+    for ( i = 0; i < 10; i++ ) {
+        v = stdlib_base_round2( x[ i ] );
+        printf( "round2(%lf) = %lf\n", x[ i ], v );
+    }
+}
+```
+
+</section>
+
+<!-- /.examples -->
+
+</section>
+
+<!-- /.c -->
+
 <!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. -->
 
 <section class="related">
@@ -180,8 +265,8 @@ Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
 [npm-image]: http://img.shields.io/npm/v/@stdlib/math-base-special-round2.svg
 [npm-url]: https://npmjs.org/package/@stdlib/math-base-special-round2
 
-[test-image]: https://github.com/stdlib-js/math-base-special-round2/actions/workflows/test.yml/badge.svg?branch=v0.2.0
-[test-url]: https://github.com/stdlib-js/math-base-special-round2/actions/workflows/test.yml?query=branch:v0.2.0
+[test-image]: https://github.com/stdlib-js/math-base-special-round2/actions/workflows/test.yml/badge.svg?branch=v0.3.0
+[test-url]: https://github.com/stdlib-js/math-base-special-round2/actions/workflows/test.yml?query=branch:v0.3.0
 
 [coverage-image]: https://img.shields.io/codecov/c/github/stdlib-js/math-base-special-round2/main.svg
 [coverage-url]: https://codecov.io/github/stdlib-js/math-base-special-round2?branch=main

package/dist/index.js.map

@@ -1,7 +1,7 @@
 {
   "version": 3,
   "sources": ["../lib/main.js", "../lib/index.js"],
-  "sourcesContent": ["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar isnan = require( '@stdlib/math-base-assert-is-nan' );\nvar isInfinite = require( '@stdlib/math-base-assert-is-infinite' );\nvar pow = require( '@stdlib/math-base-special-pow' );\nvar floor = require( '@stdlib/math-base-special-floor' );\nvar ceil = require( '@stdlib/math-base-special-ceil' );\nvar log2 = require( '@stdlib/math-base-special-log2' );\nvar MAX_EXP = require( '@stdlib/constants-float64-max-base2-exponent' );\nvar MIN_EXP_SUBNORMAL = require( '@stdlib/constants-float64-min-base2-exponent-subnormal' );\nvar PINF = require( '@stdlib/constants-float64-pinf' );\n\n\n// VARIABLES //\n\n// 2^1023:\nvar HUGE = pow( 2.0, MAX_EXP );\nvar HALF_HUGE = HUGE / 2.0;\n\n\n// MAIN //\n\n/**\n* Rounds a numeric value to the nearest power of two on a linear scale.\n*\n* @param {number} x - input value\n* @returns {number} rounded value\n*\n* @example\n* var v = round2( 3.141592653589793 );\n* // returns 4.0\n*\n* @example\n* var v = round2( 13.0 );\n* // returns 16.0\n*\n* @example\n* var v = round2( -0.314 );\n* // returns -0.25\n*/\nfunction round2( x ) {\n\tvar sign;\n\tvar half;\n\tvar p1;\n\tvar p2;\n\tvar y1;\n\tvar y2;\n\tvar p;\n\tif (\n\t\tisnan( x ) ||\n\t\tisInfinite( x ) ||\n\t\tx === 0.0\n\t) {\n\t\treturn x;\n\t}\n\tif ( x < 0 ) {\n\t\tx = -x;\n\t\tsign = -1.0;\n\t} else {\n\t\tsign = 1.0;\n\t}\n\t// Solve the equation `2^p = x` for `p`:\n\tp = log2( x );\n\n\t// If provided the smallest subnormal, no rounding possible:\n\tif ( p === MIN_EXP_SUBNORMAL ) {\n\t\treturn x;\n\t}\n\t// Find the previous and next integer powers:\n\tp1 = floor( p );\n\tp2 = ceil( p );\n\n\t// Handle overflow:\n\tif ( p1 === MAX_EXP ) {\n\t\tif ( x - HUGE >= HALF_HUGE ) {\n\t\t\treturn sign * PINF; // sign-preserving\n\t\t}\n\t\treturn sign * HUGE;\n\t}\n\t// Compute previous and next powers of two:\n\ty1 = pow( 2.0, p1 );\n\ty2 = pow( 2.0, p2 );\n\n\t// Find the closest power of two:\n\thalf = ( y2 - y1 ) / 2.0;\n\tif ( y1+half > x ) {\n\t\treturn sign * y1;\n\t}\n\treturn sign * y2;\n}\n\n\n// EXPORTS //\n\nmodule.exports = round2;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n/**\n* Round a numeric value to the nearest power of two on a linear scale.\n*\n* @module @stdlib/math-base-special-round2\n*\n* @example\n* var round2 = require( '@stdlib/math-base-special-round2' );\n*\n* var v = round2( 3.141592653589793 );\n* // returns 4.0\n*\n* v = round2( 13.0 );\n* // returns 16.0\n*\n* v = round2( -0.314 );\n* // returns -0.25\n*/\n\n// MODULES //\n\nvar round2 = require( './main.js' );\n\n\n// EXPORTS //\n\nmodule.exports = round2;\n"],
-  "mappings": "uGAAA,IAAAA,EAAAC,EAAA,SAAAC,EAAAC,EAAA,cAsBA,IAAIC,EAAQ,QAAS,iCAAkC,EACnDC,EAAa,QAAS,sCAAuC,EAC7DC,EAAM,QAAS,+BAAgC,EAC/CC,EAAQ,QAAS,iCAAkC,EACnDC,EAAO,QAAS,gCAAiC,EACjDC,EAAO,QAAS,gCAAiC,EACjDC,EAAU,QAAS,8CAA+C,EAClEC,EAAoB,QAAS,wDAAyD,EACtFC,EAAO,QAAS,gCAAiC,EAMjDC,EAAOP,EAAK,EAAKI,CAAQ,EACzBI,EAAYD,EAAO,EAuBvB,SAASE,EAAQC,EAAI,CACpB,IAAIC,EACAC,EACAC,EACAC,EACAC,EACAC,EACAC,EAkBJ,OAhBCnB,EAAOY,CAAE,GACTX,EAAYW,CAAE,GACdA,IAAM,IAIFA,EAAI,GACRA,EAAI,CAACA,EACLC,EAAO,IAEPA,EAAO,EAGRM,EAAId,EAAMO,CAAE,EAGPO,IAAMZ,GACHK,GAGRG,EAAKZ,EAAOgB,CAAE,EACdH,EAAKZ,EAAMe,CAAE,EAGRJ,IAAOT,EACNM,EAAIH,GAAQC,EACTG,EAAOL,EAERK,EAAOJ,GAGfQ,EAAKf,EAAK,EAAKa,CAAG,EAClBG,EAAKhB,EAAK,EAAKc,CAAG,EAGlBF,GAASI,EAAKD,GAAO,EAChBA,EAAGH,EAAOF,EACPC,EAAOI,EAERJ,EAAOK,GACf,CAKAnB,EAAO,QAAUY,IC1EjB,IAAIS,EAAS,IAKb,OAAO,QAAUA",
+  "sourcesContent": ["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar isnan = require( '@stdlib/math-base-assert-is-nan' );\nvar isInfinite = require( '@stdlib/math-base-assert-is-infinite' );\nvar pow = require( '@stdlib/math-base-special-pow' );\nvar floor = require( '@stdlib/math-base-special-floor' );\nvar ceil = require( '@stdlib/math-base-special-ceil' );\nvar log2 = require( '@stdlib/math-base-special-log2' );\nvar MAX_EXP = require( '@stdlib/constants-float64-max-base2-exponent' );\nvar MIN_EXP_SUBNORMAL = require( '@stdlib/constants-float64-min-base2-exponent-subnormal' );\nvar PINF = require( '@stdlib/constants-float64-pinf' );\n\n\n// VARIABLES //\n\n// 2^1023:\nvar HUGE = pow( 2.0, MAX_EXP );\nvar HALF_HUGE = HUGE / 2.0;\n\n\n// MAIN //\n\n/**\n* Rounds a numeric value to the nearest power of two on a linear scale.\n*\n* @param {number} x - input value\n* @returns {number} rounded value\n*\n* @example\n* var v = round2( 3.141592653589793 );\n* // returns 4.0\n*\n* @example\n* var v = round2( 13.0 );\n* // returns 16.0\n*\n* @example\n* var v = round2( -0.314 );\n* // returns -0.25\n*/\nfunction round2( x ) {\n\tvar sign;\n\tvar half;\n\tvar p1;\n\tvar p2;\n\tvar y1;\n\tvar y2;\n\tvar p;\n\tif (\n\t\tisnan( x ) ||\n\t\tisInfinite( x ) ||\n\t\tx === 0.0\n\t) {\n\t\treturn x;\n\t}\n\tif ( x < 0 ) {\n\t\tx = -x;\n\t\tsign = -1.0;\n\t} else {\n\t\tsign = 1.0;\n\t}\n\n\t// Solve the equation `2^p = x` for `p`:\n\tp = log2( x );\n\n\t// If provided the smallest subnormal, no rounding possible:\n\tif ( p === MIN_EXP_SUBNORMAL ) {\n\t\treturn x;\n\t}\n\n\t// Find the previous and next integer powers:\n\tp1 = floor( p );\n\tp2 = ceil( p );\n\n\t// Handle overflow:\n\tif ( p1 === MAX_EXP ) {\n\t\tif ( x - HUGE >= HALF_HUGE ) {\n\t\t\treturn sign * PINF; // sign-preserving\n\t\t}\n\t\treturn sign * HUGE;\n\t}\n\n\t// Compute previous and next powers of two:\n\ty1 = pow( 2.0, p1 );\n\ty2 = pow( 2.0, p2 );\n\n\t// Find the closest power of two:\n\thalf = ( y2 - y1 ) / 2.0;\n\tif ( y1 + half > x ) {\n\t\treturn sign * y1;\n\t}\n\treturn sign * y2;\n}\n\n\n// EXPORTS //\n\nmodule.exports = round2;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n/**\n* Round a numeric value to the nearest power of two on a linear scale.\n*\n* @module @stdlib/math-base-special-round2\n*\n* @example\n* var round2 = require( '@stdlib/math-base-special-round2' );\n*\n* var v = round2( 3.141592653589793 );\n* // returns 4.0\n*\n* v = round2( 13.0 );\n* // returns 16.0\n*\n* v = round2( -0.314 );\n* // returns -0.25\n*/\n\n// MODULES //\n\nvar round2 = require( './main.js' );\n\n\n// EXPORTS //\n\nmodule.exports = round2;\n"],
+  "mappings": "uGAAA,IAAAA,EAAAC,EAAA,SAAAC,EAAAC,EAAA,cAsBA,IAAIC,EAAQ,QAAS,iCAAkC,EACnDC,EAAa,QAAS,sCAAuC,EAC7DC,EAAM,QAAS,+BAAgC,EAC/CC,EAAQ,QAAS,iCAAkC,EACnDC,EAAO,QAAS,gCAAiC,EACjDC,EAAO,QAAS,gCAAiC,EACjDC,EAAU,QAAS,8CAA+C,EAClEC,EAAoB,QAAS,wDAAyD,EACtFC,EAAO,QAAS,gCAAiC,EAMjDC,EAAOP,EAAK,EAAKI,CAAQ,EACzBI,EAAYD,EAAO,EAuBvB,SAASE,EAAQC,EAAI,CACpB,IAAIC,EACAC,EACAC,EACAC,EACAC,EACAC,EACAC,EAmBJ,OAjBCnB,EAAOY,CAAE,GACTX,EAAYW,CAAE,GACdA,IAAM,IAIFA,EAAI,GACRA,EAAI,CAACA,EACLC,EAAO,IAEPA,EAAO,EAIRM,EAAId,EAAMO,CAAE,EAGPO,IAAMZ,GACHK,GAIRG,EAAKZ,EAAOgB,CAAE,EACdH,EAAKZ,EAAMe,CAAE,EAGRJ,IAAOT,EACNM,EAAIH,GAAQC,EACTG,EAAOL,EAERK,EAAOJ,GAIfQ,EAAKf,EAAK,EAAKa,CAAG,EAClBG,EAAKhB,EAAK,EAAKc,CAAG,EAGlBF,GAASI,EAAKD,GAAO,EAChBA,EAAKH,EAAOF,EACTC,EAAOI,EAERJ,EAAOK,GACf,CAKAnB,EAAO,QAAUY,IC7EjB,IAAIS,EAAS,IAKb,OAAO,QAAUA",
   "names": ["require_main", "__commonJSMin", "exports", "module", "isnan", "isInfinite", "pow", "floor", "ceil", "log2", "MAX_EXP", "MIN_EXP_SUBNORMAL", "PINF", "HUGE", "HALF_HUGE", "round2", "x", "sign", "half", "p1", "p2", "y1", "y2", "p", "round2"]
 }

package/include/stdlib/math/base/special/round2.h

@@ -0,0 +1,38 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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_ROUND2_H
+#define STDLIB_MATH_BASE_SPECIAL_ROUND2_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
+
+/**
+* Rounds a numeric value to the nearest power of two on a linear scale.
+*/
+double stdlib_base_round2( const double x );
+
+#ifdef __cplusplus
+}
+#endif
+
+#endif // !STDLIB_MATH_BASE_SPECIAL_ROUND2_H

package/lib/main.js

@@ -79,6 +79,7 @@ function round2( x ) {
 	} else {
 		sign = 1.0;
 	}
+
 	// Solve the equation `2^p = x` for `p`:
 	p = log2( x );
 
@@ -86,6 +87,7 @@ function round2( x ) {
 	if ( p === MIN_EXP_SUBNORMAL ) {
 		return x;
 	}
+
 	// Find the previous and next integer powers:
 	p1 = floor( p );
 	p2 = ceil( p );
@@ -97,13 +99,14 @@ function round2( x ) {
 		}
 		return sign * HUGE;
 	}
+
 	// Compute previous and next powers of two:
 	y1 = pow( 2.0, p1 );
 	y2 = pow( 2.0, p2 );
 
 	// Find the closest power of two:
 	half = ( y2 - y1 ) / 2.0;
-	if ( y1+half > x ) {
+	if ( y1 + half > x ) {
 		return sign * y1;
 	}
 	return sign * y2;

package/lib/native.js

@@ -0,0 +1,54 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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.
+*/
+
+'use strict';
+
+// MODULES //
+
+var addon = require( './../src/addon.node' );
+
+
+// MAIN //
+
+/**
+* Rounds a numeric value to the nearest power of two on a linear scale.
+*
+* @private
+* @param {number} x - input value
+* @returns {number} rounded value
+*
+* @example
+* var v = round2( 3.141592653589793 );
+* // returns 4.0
+*
+* @example
+* var v = round2( 13.0 );
+* // returns 16.0
+*
+* @example
+* var v = round2( -0.314 );
+* // returns -0.25
+*/
+function round2( x ) {
+	return addon( x );
+}
+
+
+// EXPORTS //
+
+module.exports = round2;

package/manifest.json

@@ -0,0 +1,96 @@
+{
+  "options": {
+    "task": "build"
+  },
+  "fields": [
+    {
+      "field": "src",
+      "resolve": true,
+      "relative": true
+    },
+    {
+      "field": "include",
+      "resolve": true,
+      "relative": true
+    },
+    {
+      "field": "libraries",
+      "resolve": false,
+      "relative": false
+    },
+    {
+      "field": "libpath",
+      "resolve": true,
+      "relative": false
+    }
+  ],
+  "confs": [
+    {
+      "task": "build",
+      "src": [
+        "./src/main.c"
+      ],
+      "include": [
+        "./include"
+      ],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/math-base-napi-unary",
+        "@stdlib/math-base-assert-is-nan",
+        "@stdlib/math-base-assert-is-infinite",
+        "@stdlib/math-base-special-floor",
+        "@stdlib/math-base-special-pow",
+        "@stdlib/math-base-special-ceil",
+        "@stdlib/math-base-special-log2",
+        "@stdlib/constants-float64-max-base2-exponent",
+        "@stdlib/constants-float64-min-base2-exponent-subnormal",
+        "@stdlib/constants-float64-pinf"
+      ]
+    },
+    {
+      "task": "benchmark",
+      "src": [
+        "./src/main.c"
+      ],
+      "include": [
+        "./include"
+      ],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/math-base-assert-is-nan",
+        "@stdlib/math-base-assert-is-infinite",
+        "@stdlib/math-base-special-floor",
+        "@stdlib/math-base-special-pow",
+        "@stdlib/math-base-special-ceil",
+        "@stdlib/math-base-special-log2",
+        "@stdlib/constants-float64-max-base2-exponent",
+        "@stdlib/constants-float64-min-base2-exponent-subnormal",
+        "@stdlib/constants-float64-pinf"
+      ]
+    },
+    {
+      "task": "examples",
+      "src": [
+        "./src/main.c"
+      ],
+      "include": [
+        "./include"
+      ],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/math-base-assert-is-nan",
+        "@stdlib/math-base-assert-is-infinite",
+        "@stdlib/math-base-special-floor",
+        "@stdlib/math-base-special-pow",
+        "@stdlib/math-base-special-ceil",
+        "@stdlib/math-base-special-log2",
+        "@stdlib/constants-float64-max-base2-exponent",
+        "@stdlib/constants-float64-min-base2-exponent-subnormal",
+        "@stdlib/constants-float64-pinf"
+      ]
+    }
+  ]
+}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@stdlib/math-base-special-round2",
-  "version": "0.2.0",
+  "version": "0.3.0",
   "description": "Round a numeric value to the nearest power of two on a linear scale.",
   "license": "Apache-2.0",
   "author": {
@@ -14,20 +14,16 @@
     }
   ],
   "main": "./lib",
+  "gypfile": false,
   "directories": {
-    "benchmark": "./benchmark",
     "doc": "./docs",
-    "example": "./examples",
+    "include": "./include",
     "lib": "./lib",
-    "test": "./test"
+    "src": "./src",
+    "dist": "./dist"
   },
   "types": "./docs/types",
-  "scripts": {
-    "test": "make test",
-    "test-cov": "make test-cov",
-    "examples": "make examples",
-    "benchmark": "make benchmark"
-  },
+  "scripts": {},
   "homepage": "https://stdlib.io",
   "repository": {
     "type": "git",
@@ -37,29 +33,19 @@
     "url": "https://github.com/stdlib-js/stdlib/issues"
   },
   "dependencies": {
-    "@stdlib/constants-float64-max-base2-exponent": "^0.2.0",
-    "@stdlib/constants-float64-min-base2-exponent-subnormal": "^0.2.0",
-    "@stdlib/constants-float64-pinf": "^0.2.0",
-    "@stdlib/math-base-assert-is-infinite": "^0.2.0",
-    "@stdlib/math-base-assert-is-nan": "^0.2.0",
-    "@stdlib/math-base-special-ceil": "^0.2.0",
-    "@stdlib/math-base-special-floor": "^0.2.0",
-    "@stdlib/math-base-special-log2": "^0.2.0",
-    "@stdlib/math-base-special-pow": "^0.2.0"
-  },
-  "devDependencies": {
-    "@stdlib/constants-float64-ninf": "^0.1.1",
-    "@stdlib/constants-float64-smallest-subnormal": "^0.2.0",
-    "@stdlib/math-base-assert-is-negative-zero": "^0.1.1",
-    "@stdlib/math-base-assert-is-positive-zero": "^0.1.1",
-    "@stdlib/math-base-special-exp2": "^0.1.0",
-    "@stdlib/random-base-randu": "^0.1.0",
-    "tape": "git+https://github.com/kgryte/tape.git#fix/globby",
-    "istanbul": "^0.4.1",
-    "tap-min": "git+https://github.com/Planeshifter/tap-min.git",
-    "@stdlib/bench-harness": "^0.2.0",
-    "@stdlib/bench": "^0.3.1"
+    "@stdlib/constants-float64-max-base2-exponent": "^0.2.2",
+    "@stdlib/constants-float64-min-base2-exponent-subnormal": "^0.2.1",
+    "@stdlib/constants-float64-pinf": "^0.2.2",
+    "@stdlib/math-base-assert-is-infinite": "^0.2.2",
+    "@stdlib/math-base-assert-is-nan": "^0.2.2",
+    "@stdlib/math-base-napi-unary": "^0.2.3",
+    "@stdlib/math-base-special-ceil": "^0.2.2",
+    "@stdlib/math-base-special-floor": "^0.2.3",
+    "@stdlib/math-base-special-log2": "^0.3.0",
+    "@stdlib/math-base-special-pow": "^0.3.0",
+    "@stdlib/utils-library-manifest": "^0.2.2"
   },
+  "devDependencies": {},
   "engines": {
     "node": ">=0.10.0",
     "npm": ">2.7.0"

package/src/addon.c

@@ -0,0 +1,23 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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.
+*/
+
+#include "stdlib/math/base/special/round2.h"
+#include "stdlib/math/base/napi/unary.h"
+
+// cppcheck-suppress shadowFunction
+STDLIB_MATH_BASE_NAPI_MODULE_D_D( stdlib_base_round2 )

package/src/main.c

@@ -0,0 +1,97 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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.
+*/
+
+#include "stdlib/math/base/special/round2.h"
+#include "stdlib/math/base/assert/is_nan.h"
+#include "stdlib/math/base/assert/is_infinite.h"
+#include "stdlib/math/base/special/floor.h"
+#include "stdlib/math/base/special/pow.h"
+#include "stdlib/math/base/special/ceil.h"
+#include "stdlib/math/base/special/log2.h"
+#include "stdlib/constants/float64/max_base2_exponent.h"
+#include "stdlib/constants/float64/min_base2_exponent_subnormal.h"
+#include "stdlib/constants/float64/pinf.h"
+
+/**
+* Rounds a numeric value to the nearest power of two on a linear scale.
+*
+* @param x    input value
+* @return     rounded value
+*
+* @example
+* double out = stdlib_base_round2( 3.141592653589793 );
+* // returns 4.0
+*/
+double stdlib_base_round2( const double x ) {
+	double HALF_HUGE_VALUE;
+	double HUGE_VALUE;
+	double sign;
+	double half;
+	double p1;
+	double p2;
+	double y1;
+	double y2;
+	double xc;
+	double p;
+
+	if ( stdlib_base_is_nan( x ) || stdlib_base_is_infinite( x ) || x == 0.0 ) {
+		return x;
+	}
+	xc = x;
+	if ( xc < 0 ) {
+		xc = -xc;
+		sign = -1.0;
+	} else {
+		sign = 1.0;
+	}
+
+	// Solve the equation `2^p = x` for `p`:
+	p = stdlib_base_log2( xc );
+
+	// If provided the smallest subnormal, no rounding possible:
+	if ( p == STDLIB_CONSTANT_FLOAT64_MIN_BASE2_EXPONENT_SUBNORMAL ) {
+		return xc;
+	}
+
+	// Find the previous and next integer powers:
+	p1 = stdlib_base_floor( p );
+	p2 = stdlib_base_ceil( p );
+
+	// 2^1023:
+	HUGE_VALUE = stdlib_base_pow( 2.0, STDLIB_CONSTANT_FLOAT64_MAX_BASE2_EXPONENT );
+	HALF_HUGE_VALUE = HUGE_VALUE / 2.0;
+
+	// Handle overflow:
+	if ( p1 == STDLIB_CONSTANT_FLOAT64_MAX_BASE2_EXPONENT ) {
+		if ( xc - HUGE_VALUE >= HALF_HUGE_VALUE ) {
+			return sign * STDLIB_CONSTANT_FLOAT64_PINF; // sign-preserving
+		}
+		return sign * HUGE_VALUE;
+	}
+
+	// Compute previous and next powers of two:
+	y1 = stdlib_base_pow( 2.0, p1 );
+	y2 = stdlib_base_pow( 2.0, p2 );
+
+	// Find the closest power of two:
+	half = ( y2 - y1 ) / 2.0;
+	if ( y1 + half > xc ) {
+		return sign * y1;
+	}
+	return sign * y2;
+}