@stdlib/math-base-special-ceilsd 0.2.1 → 0.3.0

package/README.md

@@ -53,25 +53,21 @@ npm install @stdlib/math-base-special-ceilsd
 var ceilsd = require( '@stdlib/math-base-special-ceilsd' );
 ```
 
-#### ceilsd( x, n\[, b] )
+#### ceilsd( x, n, b )
 
 Rounds a `numeric` value to the nearest `number` toward positive infinity with `n` significant figures.
 
 ```javascript
-var v = ceilsd( 3.141592653589793, 5 );
+var v = ceilsd( 3.141592653589793, 5, 10 );
 // returns 3.1416
 
-v = ceilsd( 3.141592653589793, 1 );
+v = ceilsd( 3.141592653589793, 1, 10 );
 // returns 4.0
 
-v = ceilsd( 12368.0, 2 );
+v = ceilsd( 12368.0, 2, 10 );
 // returns 13000.0
-```
-
-The default base is `10` (decimal). To round using a different base, provide a third argument.
 
-```javascript
-var v = ceilsd( 0.0313, 2, 2 );
+v = ceilsd( 0.0313, 2, 2 );
 // returns 0.046875
 ```
 
@@ -101,7 +97,7 @@ var i;
 
 for ( i = 0; i < 100; i++ ) {
     x = (randu()*10000.0) - 5000.0;
-    y = ceilsd( x, 5 );
+    y = ceilsd( x, 5, 10 );
     console.log( 'x: %d. Rounded: %d.', x, y );
 }
 ```
@@ -110,6 +106,96 @@ 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/ceilsd.h"
+```
+
+#### stdlib_base_ceilsd( x, n, b )
+
+Rounds a `numeric` value to the nearest `number` toward negative infinity with `n` significant figures.
+
+```c
+double out = stdlib_base_ceilsd( 0.0313, 2, 2 );
+// returns 0.046875
+```
+
+The function accepts the following arguments:
+
+-   **x**: `[in] double` input value.
+-   **n**: `[in] int32_t` number of significant figures.
+-   **b**: `[in] int32_t` base.
+
+```c
+double stdlib_base_ceilsd( const double x, const int32_t n, const int32_t b );
+```
+
+</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/ceilsd.h"
+#include <stdio.h>
+#include <stdint.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 };
+    const int32_t n[] = { 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 };
+    const int32_t b[] = { 20, 19, 18, 17, 16, 15, 14, 13, 12, 11 };
+
+    double v;
+    int i;
+    for ( i = 0; i < 10; i++ ) {
+        v = stdlib_base_ceilsd( x[ i ], n[ i ], b[ i ] );
+        printf( "ceilsd(%lf, %d, %d) = %lf\n", x[ i ], n[ i ], b[ 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">
@@ -166,8 +252,8 @@ Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
 [npm-image]: http://img.shields.io/npm/v/@stdlib/math-base-special-ceilsd.svg
 [npm-url]: https://npmjs.org/package/@stdlib/math-base-special-ceilsd
 
-[test-image]: https://github.com/stdlib-js/math-base-special-ceilsd/actions/workflows/test.yml/badge.svg?branch=v0.2.1
-[test-url]: https://github.com/stdlib-js/math-base-special-ceilsd/actions/workflows/test.yml?query=branch:v0.2.1
+[test-image]: https://github.com/stdlib-js/math-base-special-ceilsd/actions/workflows/test.yml/badge.svg?branch=v0.3.0
+[test-url]: https://github.com/stdlib-js/math-base-special-ceilsd/actions/workflows/test.yml?query=branch:v0.3.0
 
 [coverage-image]: https://img.shields.io/codecov/c/github/stdlib-js/math-base-special-ceilsd/main.svg
 [coverage-url]: https://codecov.io/github/stdlib-js/math-base-special-ceilsd?branch=main

package/dist/index.js

@@ -1,5 +1,5 @@
-"use strict";var p=function(r,e){return function(){return e||r((e={exports:{}}).exports,e),e.exports}};var c=p(function(I,o){
-var f=require('@stdlib/math-base-assert-is-nan/dist'),s=require('@stdlib/math-base-assert-is-infinite/dist'),N=require('@stdlib/math-base-special-pow/dist'),g=require('@stdlib/math-base-special-log10/dist'),l=require('@stdlib/math-base-special-ln/dist'),v=require('@stdlib/math-base-special-abs/dist'),m=require('@stdlib/math-base-special-floor/dist'),d=require('@stdlib/number-float64-base-exponent/dist'),q=require('@stdlib/math-base-special-ceil/dist');function h(r,e,t){var a,i,u,n;if(f(r)||f(e)||e<1||s(e))return NaN;if(arguments.length>2){if(f(t)||t<=0||s(t))return NaN;a=t}else a=10;return s(r)||r===0||(a===10?i=g(v(r)):a===2?i=d(v(r)):i=l(v(r))/l(a),i=m(i-e+1),u=N(a,v(i)),s(u))||(i<0?n=q(r*u)/u:n=q(r/u)*u,s(n))?r:n}o.exports=h
-});var w=c();module.exports=w;
+"use strict";var c=function(r,e){return function(){return e||r((e={exports:{}}).exports,e),e.exports}};var o=c(function(I,l){
+var n=require('@stdlib/math-base-assert-is-nan/dist'),s=require('@stdlib/math-base-assert-is-infinite/dist'),p=require('@stdlib/math-base-special-pow/dist'),N=require('@stdlib/math-base-special-log10/dist'),f=require('@stdlib/math-base-special-ln/dist'),v=require('@stdlib/math-base-special-abs/dist'),d=require('@stdlib/math-base-special-floor/dist'),g=require('@stdlib/number-float64-base-exponent/dist'),q=require('@stdlib/math-base-special-ceil/dist');function m(r,e,i){var u,a,t;return n(r)||n(e)||e<1||s(e)||n(i)||i<=0||s(i)?NaN:s(r)||r===0||(i===10?u=N(v(r)):i===2?u=g(v(r)):u=f(v(r))/f(i),u=d(u-e+1),a=p(i,v(u)),s(a))||(u<0?t=q(r*a)/a:t=q(r/a)*a,s(t))?r:t}l.exports=m
+});var w=o();module.exports=w;
 /** @license Apache-2.0 */
 //# sourceMappingURL=index.js.map

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 log10 = require( '@stdlib/math-base-special-log10' );\nvar ln = require( '@stdlib/math-base-special-ln' );\nvar abs = require( '@stdlib/math-base-special-abs' );\nvar floor = require( '@stdlib/math-base-special-floor' );\nvar exponent = require( '@stdlib/number-float64-base-exponent' );\nvar ceil = require( '@stdlib/math-base-special-ceil' );\n\n\n// MAIN //\n\n/**\n* Rounds a numeric value to the nearest number toward positive infinity with \\\\(n\\\\) significant figures.\n*\n* @param {number} x - input value\n* @param {PositiveInteger} n - number of significant figures\n* @param {PositiveInteger} [b=10] - base\n* @returns {number} rounded value\n*\n* @example\n* var v = ceilsd( 3.141592653589793, 5 );\n* // returns 3.1416\n*\n* @example\n* var v = ceilsd( 3.141592653589793, 1 );\n* // returns 4.0\n*\n* @example\n* var v = ceilsd( 12368.0, 2 );\n* // returns 13000.0\n*\n* @example\n* var v = ceilsd( 0.0313, 2, 2 );\n* // returns 0.046875\n*/\nfunction ceilsd( x, n, b ) {\n\tvar base;\n\tvar exp;\n\tvar s;\n\tvar y;\n\tif (\n\t\tisnan( x ) ||\n\t\tisnan( n ) ||\n\t\tn < 1 ||\n\t\tisInfinite( n )\n\t) {\n\t\treturn NaN;\n\t}\n\tif ( arguments.length > 2 ) {\n\t\tif (\n\t\t\tisnan( b ) ||\n\t\t\tb <= 0 ||\n\t\t\tisInfinite( b )\n\t\t) {\n\t\t\treturn NaN;\n\t\t}\n\t\tbase = b;\n\t} else {\n\t\tbase = 10;\n\t}\n\tif ( isInfinite( x ) || x === 0.0 ) {\n\t\treturn x;\n\t}\n\tif ( base === 10 ) {\n\t\texp = log10( abs( x ) );\n\t}\n\telse if ( base === 2 ) {\n\t\texp = exponent( abs( x ) );\n\t}\n\telse {\n\t\texp = ln( abs(x) ) / ln( base );\n\t}\n\texp = floor( exp - n + 1.0 );\n\ts = pow( base, abs( exp ) );\n\n\t// Check for overflow:\n\tif ( isInfinite( s ) ) {\n\t\treturn x;\n\t}\n\t// To avoid numerical stability issues due to floating-point rounding error (e.g., 3.55/0.1-35.5 = -7.105427357601e-15 and 3.55*10-35.5 = 0), we must treat positive and negative exponents separately.\n\tif ( exp < 0 ) {\n\t\ty = ceil( x * s ) / s;\n\t} else {\n\t\ty = ceil( x / s ) * s;\n\t}\n\t// Check for overflow:\n\tif ( isInfinite( y ) ) {\n\t\treturn x;\n\t}\n\treturn y;\n}\n\n\n// EXPORTS //\n\nmodule.exports = ceilsd;\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 number toward positive infinity with `n` significant figures.\n*\n* @module @stdlib/math-base-special-ceilsd\n*\n* @example\n* var ceilsd = require( '@stdlib/math-base-special-ceilsd' );\n*\n* var v = ceilsd( 3.141592653589793, 5 );\n* // returns 3.1416\n*\n* v = ceilsd( 3.141592653589793, 1 );\n* // returns 4.0\n*\n* v = ceilsd( 12368.0, 2 );\n* // returns 13000.0\n*\n* v = ceilsd( 0.0313, 2, 2 );\n* // returns 0.046875\n*/\n\n// MODULES //\n\nvar main = require( './main.js' );\n\n\n// EXPORTS //\n\nmodule.exports = main;\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,EAAK,QAAS,8BAA+B,EAC7CC,EAAM,QAAS,+BAAgC,EAC/CC,EAAQ,QAAS,iCAAkC,EACnDC,EAAW,QAAS,sCAAuC,EAC3DC,EAAO,QAAS,gCAAiC,EA6BrD,SAASC,EAAQC,EAAGC,EAAGC,EAAI,CAC1B,IAAIC,EACAC,EACAC,EACAC,EACJ,GACChB,EAAOU,CAAE,GACTV,EAAOW,CAAE,GACTA,EAAI,GACJV,EAAYU,CAAE,EAEd,MAAO,KAER,GAAK,UAAU,OAAS,EAAI,CAC3B,GACCX,EAAOY,CAAE,GACTA,GAAK,GACLX,EAAYW,CAAE,EAEd,MAAO,KAERC,EAAOD,CACR,MACCC,EAAO,GA4BR,OA1BKZ,EAAYS,CAAE,GAAKA,IAAM,IAGzBG,IAAS,GACbC,EAAMX,EAAOE,EAAKK,CAAE,CAAE,EAEbG,IAAS,EAClBC,EAAMP,EAAUF,EAAKK,CAAE,CAAE,EAGzBI,EAAMV,EAAIC,EAAIK,CAAC,CAAE,EAAIN,EAAIS,CAAK,EAE/BC,EAAMR,EAAOQ,EAAMH,EAAI,CAAI,EAC3BI,EAAIb,EAAKW,EAAMR,EAAKS,CAAI,CAAE,EAGrBb,EAAYc,CAAE,KAIdD,EAAM,EACVE,EAAIR,EAAME,EAAIK,CAAE,EAAIA,EAEpBC,EAAIR,EAAME,EAAIK,CAAE,EAAIA,EAGhBd,EAAYe,CAAE,GACXN,EAEDM,CACR,CAKAjB,EAAO,QAAUU,IC5EjB,IAAIQ,EAAO,IAKX,OAAO,QAAUA",
-  "names": ["require_main", "__commonJSMin", "exports", "module", "isnan", "isInfinite", "pow", "log10", "ln", "abs", "floor", "exponent", "ceil", "ceilsd", "x", "n", "b", "base", "exp", "s", "y", "main"]
+  "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 log10 = require( '@stdlib/math-base-special-log10' );\nvar ln = require( '@stdlib/math-base-special-ln' );\nvar abs = require( '@stdlib/math-base-special-abs' );\nvar floor = require( '@stdlib/math-base-special-floor' );\nvar exponent = require( '@stdlib/number-float64-base-exponent' );\nvar ceil = require( '@stdlib/math-base-special-ceil' );\n\n\n// MAIN //\n\n/**\n* Rounds a numeric value to the nearest number toward positive infinity with \\\\(n\\\\) significant figures.\n*\n* @param {number} x - input value\n* @param {PositiveInteger} n - number of significant figures\n* @param {PositiveInteger} b - base\n* @returns {number} rounded value\n*\n* @example\n* var v = ceilsd( 3.141592653589793, 5, 10 );\n* // returns 3.1416\n*\n* @example\n* var v = ceilsd( 3.141592653589793, 1, 10 );\n* // returns 4.0\n*\n* @example\n* var v = ceilsd( 12368.0, 2, 10 );\n* // returns 13000.0\n*\n* @example\n* var v = ceilsd( 0.0313, 2, 2 );\n* // returns 0.046875\n*/\nfunction ceilsd( x, n, b ) {\n\tvar exp;\n\tvar s;\n\tvar y;\n\tif (\n\t\tisnan( x ) ||\n\t\tisnan( n ) ||\n\t\tn < 1 ||\n\t\tisInfinite( n ) ||\n\t\tisnan( b ) ||\n\t\tb <= 0 ||\n\t\tisInfinite( b )\n\t) {\n\t\treturn NaN;\n\t}\n\tif ( isInfinite( x ) || x === 0.0 ) {\n\t\treturn x;\n\t}\n\tif ( b === 10 ) {\n\t\texp = log10( abs( x ) );\n\t}\n\telse if ( b === 2 ) {\n\t\texp = exponent( abs( x ) );\n\t}\n\telse {\n\t\texp = ln( abs(x) ) / ln( b );\n\t}\n\texp = floor( exp - n + 1.0 );\n\ts = pow( b, abs( exp ) );\n\n\t// Check for overflow:\n\tif ( isInfinite( s ) ) {\n\t\treturn x;\n\t}\n\t// To avoid numerical stability issues due to floating-point rounding error (e.g., 3.55/0.1-35.5 = -7.105427357601e-15 and 3.55*10-35.5 = 0), we must treat positive and negative exponents separately.\n\tif ( exp < 0 ) {\n\t\ty = ceil( x * s ) / s;\n\t} else {\n\t\ty = ceil( x / s ) * s;\n\t}\n\t// Check for overflow:\n\tif ( isInfinite( y ) ) {\n\t\treturn x;\n\t}\n\treturn y;\n}\n\n\n// EXPORTS //\n\nmodule.exports = ceilsd;\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 number toward positive infinity with `n` significant figures.\n*\n* @module @stdlib/math-base-special-ceilsd\n*\n* @example\n* var ceilsd = require( '@stdlib/math-base-special-ceilsd' );\n*\n* var v = ceilsd( 3.141592653589793, 5, 10 );\n* // returns 3.1416\n*\n* v = ceilsd( 3.141592653589793, 1, 10 );\n* // returns 4.0\n*\n* v = ceilsd( 12368.0, 2, 10 );\n* // returns 13000.0\n*\n* v = ceilsd( 0.0313, 2, 2 );\n* // returns 0.046875\n*/\n\n// MODULES //\n\nvar main = require( './main.js' );\n\n\n// EXPORTS //\n\nmodule.exports = main;\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,EAAK,QAAS,8BAA+B,EAC7CC,EAAM,QAAS,+BAAgC,EAC/CC,EAAQ,QAAS,iCAAkC,EACnDC,EAAW,QAAS,sCAAuC,EAC3DC,EAAO,QAAS,gCAAiC,EA6BrD,SAASC,EAAQC,EAAGC,EAAGC,EAAI,CAC1B,IAAIC,EACAC,EACAC,EACJ,OACCf,EAAOU,CAAE,GACTV,EAAOW,CAAE,GACTA,EAAI,GACJV,EAAYU,CAAE,GACdX,EAAOY,CAAE,GACTA,GAAK,GACLX,EAAYW,CAAE,EAEP,IAEHX,EAAYS,CAAE,GAAKA,IAAM,IAGzBE,IAAM,GACVC,EAAMV,EAAOE,EAAKK,CAAE,CAAE,EAEbE,IAAM,EACfC,EAAMN,EAAUF,EAAKK,CAAE,CAAE,EAGzBG,EAAMT,EAAIC,EAAIK,CAAC,CAAE,EAAIN,EAAIQ,CAAE,EAE5BC,EAAMP,EAAOO,EAAMF,EAAI,CAAI,EAC3BG,EAAIZ,EAAKU,EAAGP,EAAKQ,CAAI,CAAE,EAGlBZ,EAAYa,CAAE,KAIdD,EAAM,EACVE,EAAIP,EAAME,EAAII,CAAE,EAAIA,EAEpBC,EAAIP,EAAME,EAAII,CAAE,EAAIA,EAGhBb,EAAYc,CAAE,GACXL,EAEDK,CACR,CAKAhB,EAAO,QAAUU,IClEjB,IAAIO,EAAO,IAKX,OAAO,QAAUA",
+  "names": ["require_main", "__commonJSMin", "exports", "module", "isnan", "isInfinite", "pow", "log10", "ln", "abs", "floor", "exponent", "ceil", "ceilsd", "x", "n", "b", "exp", "s", "y", "main"]
 }

package/docs/types/index.d.ts

@@ -23,26 +23,26 @@
 *
 * @param x - input value
 * @param n - number of significant figures
-* @param b - base (default: 10)
+* @param b - base
 * @returns rounded value
 *
 * @example
-* var v = ceilsd( 3.141592653589793, 5 );
+* var v = ceilsd( 3.141592653589793, 5, 10 );
 * // returns 3.1416
 *
 * @example
-* var v = ceilsd( 3.141592653589793, 1 );
+* var v = ceilsd( 3.141592653589793, 1, 10 );
 * // returns 4.0
 *
 * @example
-* var v = ceilsd( 12368.0, 2 );
+* var v = ceilsd( 12368.0, 2, 10 );
 * // returns 13000.0
 *
 * @example
 * var v = ceilsd( 0.0313, 2, 2 );
 * // returns 0.046875
 */
-declare function ceilsd( x: number, n: number, b?: number ): number;
+declare function ceilsd( x: number, n: number, b: number ): number;
 
 
 // EXPORTS //

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

@@ -0,0 +1,40 @@
+/**
+* @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_CEILSD_H
+#define STDLIB_MATH_BASE_SPECIAL_CEILSD_H
+
+#include <stdint.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 number toward positive infinity with \\(n\\) significant figures.
+*/
+double stdlib_base_ceilsd( const double x, const int32_t n, const int32_t b );
+
+#ifdef __cplusplus
+}
+#endif
+
+#endif // !STDLIB_MATH_BASE_SPECIAL_CEILSD_H

package/lib/index.js

@@ -26,13 +26,13 @@
 * @example
 * var ceilsd = require( '@stdlib/math-base-special-ceilsd' );
 *
-* var v = ceilsd( 3.141592653589793, 5 );
+* var v = ceilsd( 3.141592653589793, 5, 10 );
 * // returns 3.1416
 *
-* v = ceilsd( 3.141592653589793, 1 );
+* v = ceilsd( 3.141592653589793, 1, 10 );
 * // returns 4.0
 *
-* v = ceilsd( 12368.0, 2 );
+* v = ceilsd( 12368.0, 2, 10 );
 * // returns 13000.0
 *
 * v = ceilsd( 0.0313, 2, 2 );

package/lib/main.js

@@ -38,19 +38,19 @@ var ceil = require( '@stdlib/math-base-special-ceil' );
 *
 * @param {number} x - input value
 * @param {PositiveInteger} n - number of significant figures
-* @param {PositiveInteger} [b=10] - base
+* @param {PositiveInteger} b - base
 * @returns {number} rounded value
 *
 * @example
-* var v = ceilsd( 3.141592653589793, 5 );
+* var v = ceilsd( 3.141592653589793, 5, 10 );
 * // returns 3.1416
 *
 * @example
-* var v = ceilsd( 3.141592653589793, 1 );
+* var v = ceilsd( 3.141592653589793, 1, 10 );
 * // returns 4.0
 *
 * @example
-* var v = ceilsd( 12368.0, 2 );
+* var v = ceilsd( 12368.0, 2, 10 );
 * // returns 13000.0
 *
 * @example
@@ -58,7 +58,6 @@ var ceil = require( '@stdlib/math-base-special-ceil' );
 * // returns 0.046875
 */
 function ceilsd( x, n, b ) {
-	var base;
 	var exp;
 	var s;
 	var y;
@@ -66,36 +65,27 @@ function ceilsd( x, n, b ) {
 		isnan( x ) ||
 		isnan( n ) ||
 		n < 1 ||
-		isInfinite( n )
+		isInfinite( n ) ||
+		isnan( b ) ||
+		b <= 0 ||
+		isInfinite( b )
 	) {
 		return NaN;
 	}
-	if ( arguments.length > 2 ) {
-		if (
-			isnan( b ) ||
-			b <= 0 ||
-			isInfinite( b )
-		) {
-			return NaN;
-		}
-		base = b;
-	} else {
-		base = 10;
-	}
 	if ( isInfinite( x ) || x === 0.0 ) {
 		return x;
 	}
-	if ( base === 10 ) {
+	if ( b === 10 ) {
 		exp = log10( abs( x ) );
 	}
-	else if ( base === 2 ) {
+	else if ( b === 2 ) {
 		exp = exponent( abs( x ) );
 	}
 	else {
-		exp = ln( abs(x) ) / ln( base );
+		exp = ln( abs(x) ) / ln( b );
 	}
 	exp = floor( exp - n + 1.0 );
-	s = pow( base, abs( exp ) );
+	s = pow( b, abs( exp ) );
 
 	// Check for overflow:
 	if ( isInfinite( s ) ) {

package/lib/native.js

@@ -0,0 +1,60 @@
+/**
+* @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 number toward positive infinity with \\(n\\) significant figures.
+*
+* @private
+* @param {number} x - input value
+* @param {PositiveInteger} n - number of significant figures
+* @param {PositiveInteger} b - base
+* @returns {number} rounded value
+*
+* @example
+* var v = ceilsd( 3.141592653589793, 5, 10 );
+* // returns 3.1416
+*
+* @example
+* var v = ceilsd( 3.141592653589793, 1, 10 );
+* // returns 4.0
+*
+* @example
+* var v = ceilsd( 12368.0, 2, 10 );
+* // returns 13000.0
+*
+* @example
+* var v = ceilsd( 0.0313, 2, 2 );
+* // returns 0.046875
+*/
+function ceilsd( x, n, b ) {
+	return addon( x, n, b );
+}
+
+
+// EXPORTS //
+
+module.exports = ceilsd;

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-ternary",
+        "@stdlib/math-base-assert-is-nan",
+        "@stdlib/math-base-assert-is-infinite",
+        "@stdlib/math-base-special-pow",
+        "@stdlib/math-base-special-log10",
+        "@stdlib/math-base-special-ln",
+        "@stdlib/math-base-special-abs",
+        "@stdlib/math-base-special-floor",
+        "@stdlib/math-base-special-ceil",
+        "@stdlib/number-float64-base-exponent"
+      ]
+    },
+    {
+      "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-log10",
+        "@stdlib/math-base-special-ln",
+        "@stdlib/math-base-special-abs",
+        "@stdlib/math-base-special-floor",
+        "@stdlib/math-base-special-ceil",
+        "@stdlib/number-float64-base-exponent"
+      ]
+    },
+    {
+      "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-log10",
+        "@stdlib/math-base-special-ln",
+        "@stdlib/math-base-special-abs",
+        "@stdlib/math-base-special-floor",
+        "@stdlib/math-base-special-ceil",
+        "@stdlib/number-float64-base-exponent"
+      ]
+    }
+  ]
+}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@stdlib/math-base-special-ceilsd",
-  "version": "0.2.1",
+  "version": "0.3.0",
   "description": "Round a numeric value to the nearest number toward positive infinity with N significant figures.",
   "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/math-base-assert-is-infinite": "^0.2.1",
-    "@stdlib/math-base-assert-is-nan": "^0.2.1",
-    "@stdlib/math-base-special-abs": "^0.2.1",
-    "@stdlib/math-base-special-ceil": "^0.2.1",
-    "@stdlib/math-base-special-floor": "^0.2.1",
-    "@stdlib/math-base-special-ln": "^0.2.1",
-    "@stdlib/math-base-special-log10": "^0.2.1",
-    "@stdlib/math-base-special-pow": "^0.2.1",
-    "@stdlib/number-float64-base-exponent": "^0.2.1"
+    "@stdlib/math-base-assert-is-infinite": "^0.2.2",
+    "@stdlib/math-base-assert-is-nan": "^0.2.2",
+    "@stdlib/math-base-napi-ternary": "^0.3.0",
+    "@stdlib/math-base-special-abs": "^0.2.2",
+    "@stdlib/math-base-special-ceil": "^0.2.2",
+    "@stdlib/math-base-special-floor": "^0.2.3",
+    "@stdlib/math-base-special-ln": "^0.2.4",
+    "@stdlib/math-base-special-log10": "^0.3.0",
+    "@stdlib/math-base-special-pow": "^0.3.0",
+    "@stdlib/number-float64-base-exponent": "^0.2.2",
+    "@stdlib/utils-library-manifest": "^0.2.2"
   },
   "devDependencies": {},
   "engines": {

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/ceilsd.h"
+#include "stdlib/math/base/napi/ternary.h"
+
+// cppcheck-suppress shadowFunction
+STDLIB_MATH_BASE_NAPI_MODULE_DII_D( stdlib_base_ceilsd )

package/src/main.c

@@ -0,0 +1,81 @@
+/**
+* @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/ceilsd.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/log10.h"
+#include "stdlib/math/base/special/ln.h"
+#include "stdlib/math/base/special/abs.h"
+#include "stdlib/math/base/special/floor.h"
+#include "stdlib/math/base/special/ceil.h"
+#include "stdlib/number/float64/base/exponent.h"
+#include <stdint.h>
+
+/**
+* Rounds a numeric value to the nearest number toward positive infinity with \\(n\\) significant figures.
+*
+* @param x    input value
+* @param n    number of significant figures
+* @param b    base
+* @return     rounded value
+*
+* @example
+* double out = stdlib_base_ceilsd( 3.141592653589793, 5, 10 );
+* // returns 3.1416
+*/
+double stdlib_base_ceilsd( const double x, const int32_t n, const int32_t b ) {
+	double exp;
+	double s;
+	double y;
+
+	if ( stdlib_base_is_nan( x ) || n < 1 || b <= 0 ) {
+		return 0.0 / 0.0; // NaN
+	}
+	if ( stdlib_base_is_infinite( x ) || x == 0.0 ) {
+		return x;
+	}
+	if ( b == 10 ) {
+		exp = stdlib_base_log10( stdlib_base_abs( x ) );
+	} else if ( b == 2 ) {
+		exp = stdlib_base_float64_exponent( stdlib_base_abs( x ) );
+	} else {
+		exp = stdlib_base_ln( stdlib_base_abs( x ) ) / stdlib_base_ln( (double)b );
+	}
+	exp = stdlib_base_floor( exp - (double)n + 1.0 );
+	s = stdlib_base_pow( (double)b, stdlib_base_abs( exp ) );
+
+	// Check for overflow:
+	if ( stdlib_base_is_infinite( s ) ) {
+		return x;
+	}
+
+	// To avoid numerical stability issues due to floating-point rounding error (e.g., 3.55/0.1-35.5 = -7.105427357601e-15 and 3.55*10-35.5 = 0), we must treat positive and negative exponents separately.
+	if ( exp < 0 ) {
+		y = stdlib_base_ceil( x * s ) / s;
+	} else {
+		y = stdlib_base_ceil( x / s ) * s;
+	}
+
+	// Check for overflow:
+	if ( stdlib_base_is_infinite( y ) ) {
+		return x;
+	}
+	return y;
+}