@stdlib/math-base-special-ceil2 0.2.1 → 0.3.1

package/NOTICE

@@ -1 +1 @@
-Copyright (c) 2016-2024 The Stdlib Authors.
+Copyright (c) 2016-2026 The Stdlib Authors.

package/README.md

@@ -106,17 +106,96 @@ v = ceil2( NaN );
 <!-- eslint no-undef: "error" -->
 
 ```javascript
-var randu = require( '@stdlib/random-base-randu' );
+var uniform = require( '@stdlib/random-array-uniform' );
+var logEachMap = require( '@stdlib/console-log-each-map' );
 var ceil2 = require( '@stdlib/math-base-special-ceil2' );
 
-var x;
-var v;
-var i;
+var opts = {
+    'dtype': 'float64'
+};
+var x = uniform( 100, -50.0, 50.0, opts );
 
-for ( i = 0; i < 100; i++ ) {
-    x = (randu()*100.0) - 50.0;
-    v = ceil2( x );
-    console.log( 'Value: %d. Rounded: %d.', x, v );
+logEachMap( 'x: %0.4f. Rounded: %0.4f.', x, ceil2 );
+```
+
+</section>
+
+<!-- /.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/ceil2.h"
+```
+
+#### stdlib_base_ceil2( x )
+
+Rounds a `numeric` value to the nearest power of two toward positive infinity.
+
+```c
+double y = stdlib_base_ceil2( 3.14 );
+// returns 4.0
+```
+
+The function accepts the following arguments:
+
+-   **x**: `[in] double` input value.
+
+```c
+double stdlib_base_ceil2( 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/ceil2.h"
+#include <stdio.h>
+
+int main( void ) {
+    const double x[] = { 3.14, -3.14, 0.5, 0.0 / 0.0 };
+
+    double y;
+    int i;
+    for ( i = 0; i < 4; i++ ) {
+        y = stdlib_base_ceil2( x[ i ] );
+        printf( "ceil2(%lf) = %lf\n", x[ i ], y );
+    }
 }
 ```
 
@@ -124,6 +203,10 @@ for ( i = 0; i < 100; i++ ) {
 
 <!-- /.examples -->
 
+</section>
+
+<!-- /.c -->
+
 <!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. -->
 
 <section class="related">
@@ -167,7 +250,7 @@ See [LICENSE][stdlib-license].
 
 ## Copyright
 
-Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
+Copyright &copy; 2016-2026. The Stdlib [Authors][stdlib-authors].
 
 </section>
 
@@ -180,8 +263,8 @@ Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
 [npm-image]: http://img.shields.io/npm/v/@stdlib/math-base-special-ceil2.svg
 [npm-url]: https://npmjs.org/package/@stdlib/math-base-special-ceil2
 
-[test-image]: https://github.com/stdlib-js/math-base-special-ceil2/actions/workflows/test.yml/badge.svg?branch=v0.2.1
-[test-url]: https://github.com/stdlib-js/math-base-special-ceil2/actions/workflows/test.yml?query=branch:v0.2.1
+[test-image]: https://github.com/stdlib-js/math-base-special-ceil2/actions/workflows/test.yml/badge.svg?branch=v0.3.1
+[test-url]: https://github.com/stdlib-js/math-base-special-ceil2/actions/workflows/test.yml?query=branch:v0.3.1
 
 [coverage-image]: https://img.shields.io/codecov/c/github/stdlib-js/math-base-special-ceil2/main.svg
 [coverage-url]: https://codecov.io/github/stdlib-js/math-base-special-ceil2?branch=main
@@ -193,8 +276,8 @@ Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
 
 -->
 
-[chat-image]: https://img.shields.io/gitter/room/stdlib-js/stdlib.svg
-[chat-url]: https://app.gitter.im/#/room/#stdlib-js_stdlib:gitter.im
+[chat-image]: https://img.shields.io/badge/zulip-join_chat-brightgreen.svg
+[chat-url]: https://stdlib.zulipchat.com
 
 [stdlib]: https://github.com/stdlib-js/stdlib
 

package/include/stdlib/math/base/special/ceil2.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_CEIL2_H
+#define STDLIB_MATH_BASE_SPECIAL_CEIL2_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 toward positive infinity.
+*/
+double stdlib_base_ceil2( const double x );
+
+#ifdef __cplusplus
+}
+#endif
+
+#endif // !STDLIB_MATH_BASE_SPECIAL_CEIL2_H

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 toward positive infinity.
+*
+* @private
+* @param {number} x - input value
+* @returns {number} rounded value
+*
+* @example
+* var v = ceil2( 3.141592653589793 );
+* // returns 4.0
+*
+* @example
+* var v = ceil2( 9.0 );
+* // returns 16.0
+*
+* @example
+* var v = ceil2( -0.314 );
+* // returns -0.25
+*/
+function ceil2( x ) {
+	return addon( x );
+}
+
+
+// EXPORTS //
+
+module.exports = ceil2;

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-pow",
+        "@stdlib/math-base-special-floor",
+        "@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-pow",
+        "@stdlib/math-base-special-floor",
+        "@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-pow",
+        "@stdlib/math-base-special-floor",
+        "@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-ceil2",
-  "version": "0.2.1",
+  "version": "0.3.1",
   "description": "Round a numeric value to the nearest power of two toward positive infinity.",
   "license": "Apache-2.0",
   "author": {
@@ -14,9 +14,12 @@
     }
   ],
   "main": "./lib",
+  "gypfile": false,
   "directories": {
     "doc": "./docs",
+    "include": "./include",
     "lib": "./lib",
+    "src": "./src",
     "dist": "./dist"
   },
   "types": "./docs/types",
@@ -30,15 +33,17 @@
     "url": "https://github.com/stdlib-js/stdlib/issues"
   },
   "dependencies": {
-    "@stdlib/constants-float64-max-base2-exponent": "^0.2.1",
-    "@stdlib/constants-float64-min-base2-exponent-subnormal": "^0.2.0",
-    "@stdlib/constants-float64-pinf": "^0.2.1",
-    "@stdlib/math-base-assert-is-infinite": "^0.2.1",
-    "@stdlib/math-base-assert-is-nan": "^0.2.1",
-    "@stdlib/math-base-special-ceil": "^0.2.1",
-    "@stdlib/math-base-special-floor": "^0.2.1",
-    "@stdlib/math-base-special-log2": "^0.2.1",
-    "@stdlib/math-base-special-pow": "^0.2.1"
+    "@stdlib/constants-float64-max-base2-exponent": "^0.2.3",
+    "@stdlib/constants-float64-min-base2-exponent-subnormal": "^0.2.1",
+    "@stdlib/constants-float64-pinf": "^0.2.3",
+    "@stdlib/math-base-assert-is-infinite": "^0.2.3",
+    "@stdlib/math-base-assert-is-nan": "^0.2.3",
+    "@stdlib/math-base-napi-unary": "^0.2.7",
+    "@stdlib/math-base-special-ceil": "^0.2.3",
+    "@stdlib/math-base-special-floor": "^0.2.4",
+    "@stdlib/math-base-special-log2": "^0.3.1",
+    "@stdlib/math-base-special-pow": "^0.3.0",
+    "@stdlib/utils-library-manifest": "^0.2.4"
   },
   "devDependencies": {},
   "engines": {

package/src/addon.c

@@ -0,0 +1,22 @@
+/**
+* @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/ceil2.h"
+#include "stdlib/math/base/napi/unary.h"
+
+STDLIB_MATH_BASE_NAPI_MODULE_D_D( stdlib_base_ceil2 )

package/src/main.c

@@ -0,0 +1,76 @@
+/**
+* @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/ceil2.h"
+#include "stdlib/math/base/assert/is_nan.h"
+#include "stdlib/math/base/assert/is_infinite.h"
+#include "stdlib/math/base/special/pow.h"
+#include "stdlib/math/base/special/floor.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 toward positive infinity.
+*
+* @param x    input value
+* @return     rounded value
+*
+* @example
+* double y = stdlib_base_ceil2( 3.141592653589793 );
+* // returns 4.0
+*/
+double stdlib_base_ceil2( const double x ) {
+	double sign;
+	double xc;
+	double p;
+
+	xc = x;
+	if ( stdlib_base_is_nan( xc ) || stdlib_base_is_infinite( xc ) || xc == 0.0 ) {
+		return xc;
+	}
+	if ( xc < 0.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;
+	}
+
+	// Determine a power of two which rounds the input value toward positive infinity:
+	if ( sign == -1.0 ) {
+		p = stdlib_base_floor( p );
+	} else {
+		p = stdlib_base_ceil( p );
+	}
+
+	// Handle overflow:
+	if ( p > STDLIB_CONSTANT_FLOAT64_MAX_BASE2_EXPONENT ) {
+		return STDLIB_CONSTANT_FLOAT64_PINF;
+	}
+	return sign * stdlib_base_pow( 2.0, p );
+}