From 9c9ed39870107e246e7bcd64e1238f42485514f5 Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 10:11:45 +0530 Subject: [PATCH 01/20] add c-api documentation in readme --- .../@stdlib/math/base/special/log10/README.md | 94 ++++++++++++++++++- 1 file changed, 92 insertions(+), 2 deletions(-) diff --git a/lib/node_modules/@stdlib/math/base/special/log10/README.md b/lib/node_modules/@stdlib/math/base/special/log10/README.md index f8f9f577e3e4..7f91e9c0f4a2 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/README.md +++ b/lib/node_modules/@stdlib/math/base/special/log10/README.md @@ -2,7 +2,7 @@ @license Apache-2.0 -Copyright (c) 2018 The Stdlib Authors. +Copyright (c) 2022 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. @@ -57,7 +57,7 @@ Evaluates the [common logarithm][common-logarithm]. ```javascript var v = log10( 100.0 ); -// returns 2.0 +// returns ~2.0 v = log10( 8.0 ); // returns ~0.903 @@ -107,6 +107,96 @@ for ( i = 0; i < 100; i++ ) { + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/math/base/special/log10.h" +``` + +#### stdlib_base_log10( x ) + +Evaluates the [common logarithm][common-logarithm]. + +```c +double out = stdlib_base_log10( 100.0 ); +// returns ~2.0 + +out = stdlib_base_log10( 8.0 ); +// returns ~0.903 +``` + +The function accepts the following arguments: + +- **x**: `[in] double` input value. + +```c +double stdlib_base_log10( const double x ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/math/base/special/log10.h" +#include +#include + +int main() { + double x; + double v; + int i; + + for ( i = 0; i < 100; i++ ) { + x = ( (double) rand() / (double) RAND_MAX ) * 100.0; + v = stdlib_base_log10( x ); + printf( "log10(%lf) = %lf\n", x, v ); + } +} +``` + +
+ + + +
+ + +
From 7c09bdb2f945259fbb0e0a09accb91b2acb663f7 Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 10:17:48 +0530 Subject: [PATCH 02/20] add include.gypi and binding.gyp --- .../math/base/special/log10/include.gypi | 53 +++++++++++++++++++ 1 file changed, 53 insertions(+) create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/include.gypi diff --git a/lib/node_modules/@stdlib/math/base/special/log10/include.gypi b/lib/node_modules/@stdlib/math/base/special/log10/include.gypi new file mode 100644 index 000000000000..3b437d524797 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/include.gypi @@ -0,0 +1,53 @@ +# @license Apache-2.0 +# +# Copyright (c) 2022 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. + +# A GYP include file for building a Node.js native add-on. +# +# Main documentation: +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Source directory: + 'src_dir': './src', + + # Include directories: + 'include_dirs': [ + ' Date: Tue, 20 Dec 2022 10:19:06 +0530 Subject: [PATCH 03/20] Add binding.gyp --- .../math/base/special/log10/binding.gyp | 170 ++++++++++++++++++ 1 file changed, 170 insertions(+) create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/binding.gyp diff --git a/lib/node_modules/@stdlib/math/base/special/log10/binding.gyp b/lib/node_modules/@stdlib/math/base/special/log10/binding.gyp new file mode 100644 index 000000000000..1058b57bab16 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/binding.gyp @@ -0,0 +1,170 @@ +# @license Apache-2.0 +# +# Copyright (c) 2022 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. + +# A `.gyp` file for building a Node.js native add-on. +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # List of files to include in this file: + 'includes': [ + './include.gypi', + ], + + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Target name should match the add-on export name: + 'addon_target_name%': 'addon', + + # Set variables based on the host OS: + 'conditions': [ + [ + 'OS=="win"', + { + # Define the object file suffix: + 'obj': 'obj', + }, + { + # Define the object file suffix: + 'obj': 'o', + } + ], # end condition (OS=="win") + ], # end conditions + }, # end variables + + # Define compile targets: + 'targets': [ + + # Target to generate an add-on: + { + # The target name should match the add-on export name: + 'target_name': '<(addon_target_name)', + + # Define dependencies: + 'dependencies': [], + + # Define directories which contain relevant include headers: + 'include_dirs': [ + # Local include directory: + '<@(include_dirs)', + ], + + # List of source files: + 'sources': [ + '<@(src_files)', + ], + + # Settings which should be applied when a target's object files are used as linker input: + 'link_settings': { + # Define libraries: + 'libraries': [ + '<@(libraries)', + ], + + # Define library directories: + 'library_dirs': [ + '<@(library_dirs)', + ], + }, + + # C/C++ compiler flags: + 'cflags': [ + # Enable commonly used warning options: + '-Wall', + + # Aggressive optimization: + '-O3', + ], + + # C specific compiler flags: + 'cflags_c': [ + # Specify the C standard to which a program is expected to conform: + '-std=c99', + ], + + # C++ specific compiler flags: + 'cflags_cpp': [ + # Specify the C++ standard to which a program is expected to conform: + '-std=c++11', + ], + + # Linker flags: + 'ldflags': [], + + # Apply conditions based on the host OS: + 'conditions': [ + [ + 'OS=="mac"', + { + # Linker flags: + 'ldflags': [ + '-undefined dynamic_lookup', + '-Wl,-no-pie', + '-Wl,-search_paths_first', + ], + }, + ], # end condition (OS=="mac") + [ + 'OS!="win"', + { + # C/C++ flags: + 'cflags': [ + # Generate platform-independent code: + '-fPIC', + ], + }, + ], # end condition (OS!="win") + ], # end conditions + }, # end target <(addon_target_name) + + # Target to copy a generated add-on to a standard location: + { + 'target_name': 'copy_addon', + + # Declare that the output of this target is not linked: + 'type': 'none', + + # Define dependencies: + 'dependencies': [ + # Require that the add-on be generated before building this target: + '<(addon_target_name)', + ], + + # Define a list of actions: + 'actions': [ + { + 'action_name': 'copy_addon', + 'message': 'Copying addon...', + + # Explicitly list the inputs in the command-line invocation below: + 'inputs': [], + + # Declare the expected outputs: + 'outputs': [ + '<(addon_output_dir)/<(addon_target_name).node', + ], + + # Define the command-line invocation: + 'action': [ + 'cp', + '<(PRODUCT_DIR)/<(addon_target_name).node', + '<(addon_output_dir)/<(addon_target_name).node', + ], + }, + ], # end actions + }, # end target copy_addon + ], # end targets +} From 0d05417b47a3a308d32f5e3a208e887deaf523db Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 10:20:35 +0530 Subject: [PATCH 04/20] Add log10.h --- .../include/stdlib/math/base/special/log10.h | 38 +++++++++++++++++++ 1 file changed, 38 insertions(+) create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h diff --git a/lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h b/lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h new file mode 100644 index 000000000000..6ab17d8dc34f --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h @@ -0,0 +1,38 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2022 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_LOG10_H +#define STDLIB_MATH_BASE_SPECIAL_LOG10_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 common logarithm. +*/ +double stdlib_base_log10( const double x ); + +#ifdef __cplusplus +} +#endif + +#endif // !STDLIB_MATH_BASE_SPECIAL_LOG10_H \ No newline at end of file From 9a0a4f2347753e60f4881717c753a059cd9ce2d1 Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 10:47:12 +0530 Subject: [PATCH 05/20] Add Makefile, addon.c and main.c in src/ --- .../math/base/special/log10/src/Makefile | 70 ++++++ .../math/base/special/log10/src/addon.c | 23 ++ .../math/base/special/log10/src/main.c | 209 ++++++++++++++++++ 3 files changed, 302 insertions(+) create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/src/Makefile create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/src/addon.c create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/src/main.c diff --git a/lib/node_modules/@stdlib/math/base/special/log10/src/Makefile b/lib/node_modules/@stdlib/math/base/special/log10/src/Makefile new file mode 100644 index 000000000000..f79b87238713 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/src/Makefile @@ -0,0 +1,70 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2022 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. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + + +# RULES # + +#/ +# Removes generated files for building an add-on. +# +# @example +# make clean-addon +#/ +clean-addon: + $(QUIET) -rm -f *.o *.node + +.PHONY: clean-addon + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: clean-addon + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/math/base/special/log10/src/addon.c b/lib/node_modules/@stdlib/math/base/special/log10/src/addon.c new file mode 100644 index 000000000000..068a8d89162a --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/src/addon.c @@ -0,0 +1,23 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2022 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/log10.h" +#include "stdlib/math/base/napi/unary.h" + +// cppcheck-suppress shadowFunction +STDLIB_MATH_BASE_NAPI_MODULE_D_D( stdlib_base_log10 ) diff --git a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c new file mode 100644 index 000000000000..ce201e77a09a --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c @@ -0,0 +1,209 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2022 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. +* +* +* ## Notice +* +* The following copyright and license were part of the original implementation available as part of [FreeBSD]{@link https://svnweb.freebsd.org/base/release/9.3.0/lib/msun/src/e_log10.c}. The implementation follows the original, but has been modified for JavaScript. +* +* ```text +* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. +* +* Developed at SunPro, a Sun Microsystems, Inc. business. +* Permission to use, copy, modify, and distribute this +* software is freely granted, provided that this notice +* is preserved. +* ``` +*/ + +#include "stdlib/math/base/special/log10.h" +#include "stdlib/math/base/assert/is_nan.h" +#include "stdlib/number/float64/base/get_high_word.h" +#include "stdlib/number/float64/base/set_high_word.h" +#include "stdlib/number/float64/base/set_low_word.h" +#include "stdlib/constants/float64/exponent_bias.h" +#include "stdlib/constants/float64/ninf.h" + +static const double TWO54 = 1.80143985094819840000e+16; // 0x43500000, 0x00000000 +static const double IVLN10HI = 4.34294481878168880939e-01; // 0x3fdbcb7b, 0x15200000 +static const double IVLN10LO = 2.50829467116452752298e-11; // 0x3dbb9438, 0xca9aadd5 +static const double LOG10_2HI = 3.01029995663611771306e-01; // 0x3FD34413, 0x509F6000 +static const double LOG10_2LO = 3.69423907715893078616e-13; // 0x3D59FEF3, 0x11F12B36 + +// 0x000fffff = 1048575 => 0 00000000000 11111111111111111111 +static const int64_t HIGH_SIGNIFICAND_MASK = 0x000fffff; + +// 0x7ff00000 = 2146435072 => 0 11111111111 00000000000000000000 => biased exponent: 2047 = 1023+1023 => 2^1023 +static const int64_t HIGH_MAX_NORMAL_EXP = 0x7ff00000; + +// 0x00100000 = 1048576 => 0 00000000001 00000000000000000000 => biased exponent: 1 = -1022+1023 => 2^-1022 +static const int64_t HIGH_MIN_NORMAL_EXP = 0x00100000; + +// 0x3ff00000 = 1072693248 => 0 01111111111 00000000000000000000 => biased exponent: 1023 = 0+1023 => 2^0 = 1 +static const int64_t HIGH_BIASED_EXP_0 = 0x3ff00000; + +// 1/3 +static const double ONE_THIRD = 0.33333333333333333; + +/* Begin auto-generated functions. The following functions are auto-generated. Do not edit directly. */ + +// BEGIN: polyval_p + +/** +* Evaluates a polynomial. +* +* ## Notes +* +* - The implementation uses [Horner's rule][horners-method] for efficient computation. +* +* [horners-method]: https://en.wikipedia.org/wiki/Horner%27s_method +* +* @param x value at which to evaluate the polynomial +* @return evaluated polynomial +*/ +static double polyval_p( const double x ) { + return 0.3999999999940942 + (x * (0.22222198432149784 + (x * 0.15313837699209373))); +} + +// END: polyval_p + +// BEGIN: polyval_q + +/** +* Evaluates a polynomial. +* +* ## Notes +* +* - The implementation uses [Horner's rule][horners-method] for efficient computation. +* +* [horners-method]: https://en.wikipedia.org/wiki/Horner%27s_method +* +* @param x value at which to evaluate the polynomial +* @return evaluated polynomial +*/ +static double polyval_q( const double x ) { + return 0.6666666666666735 + (x * (0.2857142874366239 + (x * (0.1818357216161805 + (x * 0.14798198605116586))))); +} + +// END: polyval_q + +/* End auto-generated functions. */ + +/** +* Return `log(x) - (x-1)` for `x` in `~[sqrt(2)/2, sqrt(2)]`. +* +* @param {number} x - input value +* @returns {number} function value +*/ +double klog( const double x ) { + double hfsq; + double t1; + double t2; + uint32_t hx; + double f; + double s; + double z; + double R; + double w; + int32_t i; + int32_t j; + + stdlib_base_float64_get_high_word( x, &hx ); + f = x - 1.0; + if ( ( HIGH_SIGNIFICAND_MASK & ( 2 + hx ) ) < 3 ) { + // Case: -2**-20 <= f < 2**-20 + if ( f == 0.0 ) { + return 0.0; + } + return f * f * ( ( ONE_THIRD * f ) - 0.5 ); + } + s = f / ( 2.0 + f ); + z = s * s; + hx &= HIGH_SIGNIFICAND_MASK; + i = (hx - 0x6147a); + w = z * z; + j = ( 0x6b851 - hx ); + t1 = w * polyval_p( w ); + t2 = z * polyval_q( w ); + i |= j; + R = t2 + t1; + if ( i > 0 ) { + hfsq = 0.5 * f * f; + return ( s * ( hfsq + R ) ) - hfsq; + } + return s * ( R - f ); +} + +/** +* Evaluates the common logarithm. +* +* @param x input value +* @return output value +* +* @example +* double out = stdlib_base_log10( 4.0 ); +* // returns 0.6020599913279624 +*/ +double stdlib_base_log10( const double x ) { + double xc = x; + double hi; + uint32_t hx; + double lo; + double f; + int32_t i; + int32_t k; + int32_t y; + double z; + + if ( stdlib_base_is_nan( xc ) || xc < 0.0 ) { + return 0.0 / 0.0; // NaN + } + if ( xc == 0.0 ) { + return STDLIB_CONSTANT_FLOAT64_NINF; + } + stdlib_base_float64_get_high_word( xc, &hx ); + k = 0; + + // Case: 0 < x < 2**-1022 + if ( hx < HIGH_MIN_NORMAL_EXP ) { + // Subnormal number, scale up `x`... + k -= 54; + xc *= TWO54; + stdlib_base_float64_get_high_word( xc, &hx ); + } + if ( hx >= HIGH_MAX_NORMAL_EXP ) { + return x + x; + } + k += ( ( hx>>20 ) - STDLIB_CONSTANT_FLOAT64_EXPONENT_BIAS ); + hx &= HIGH_SIGNIFICAND_MASK; + i = ( hx + 0x95f64 ) & 0x100000; + + // Normalize `x` or `x/2`... + stdlib_base_float64_set_high_word( xc, hx | ( i ^ HIGH_BIASED_EXP_0 ) ); + k += i >> 20; + y = k; + f = klog( xc ); + xc -= 1.0; + stdlib_base_float64_set_low_word( 0, &hi ); + lo = xc - hi; + z = ( y * LOG10_2LO ) + ( ( xc + f ) * IVLN10LO ); + z += ( ( lo + f ) * IVLN10HI ) + ( hi * IVLN10HI ); + return z + ( y * LOG10_2HI ); +} + + + From cba12c63d90cec50bab7ae62794bed67ecd1cc7e Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 10:50:16 +0530 Subject: [PATCH 06/20] Add manifest.json --- .../math/base/special/log10/manifest.json | 87 +++++++++++++++++++ 1 file changed, 87 insertions(+) create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/manifest.json diff --git a/lib/node_modules/@stdlib/math/base/special/log10/manifest.json b/lib/node_modules/@stdlib/math/base/special/log10/manifest.json new file mode 100644 index 000000000000..379083a0ccdd --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/manifest.json @@ -0,0 +1,87 @@ +{ + "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/constants/float64/ninf", + "@stdlib/constants/float64/exponent-bias", + "@stdlib/number/float64/base/get-high-word", + "@stdlib/number/float64/base/set-high-word", + "@stdlib/number/float64/base/set-low-word" + ] + }, + { + "task": "benchmark", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/assert/is-nan", + "@stdlib/constants/float64/ninf", + "@stdlib/constants/float64/exponent-bias", + "@stdlib/number/float64/base/get-high-word", + "@stdlib/number/float64/base/set-high-word", + "@stdlib/number/float64/base/set-low-word" + ] + }, + { + "task": "examples", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/assert/is-nan", + "@stdlib/constants/float64/ninf", + "@stdlib/constants/float64/exponent-bias", + "@stdlib/number/float64/base/get-high-word", + "@stdlib/number/float64/base/set-high-word", + "@stdlib/number/float64/base/set-low-word" + ] + } + ] + } \ No newline at end of file From f04eb499cff20b9e2f53d7fd8fb726e11da9740f Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 10:52:15 +0530 Subject: [PATCH 07/20] update manifest.json --- .../math/base/special/log10/manifest.json | 134 +++++++++--------- 1 file changed, 67 insertions(+), 67 deletions(-) diff --git a/lib/node_modules/@stdlib/math/base/special/log10/manifest.json b/lib/node_modules/@stdlib/math/base/special/log10/manifest.json index 379083a0ccdd..6a8b58c4b818 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/manifest.json +++ b/lib/node_modules/@stdlib/math/base/special/log10/manifest.json @@ -1,87 +1,87 @@ { - "options": { - "task": "build" + "options": { + "task": "build" + }, + "fields": [ + { + "field": "src", + "resolve": true, + "relative": true }, - "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", + { + "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/constants/float64/ninf", "@stdlib/constants/float64/exponent-bias", "@stdlib/number/float64/base/get-high-word", "@stdlib/number/float64/base/set-high-word", "@stdlib/number/float64/base/set-low-word" - ] - }, - { - "task": "benchmark", - "src": [ - "./src/main.c" - ], - "include": [ - "./include" - ], - "libraries": [], - "libpath": [], - "dependencies": [ + ] + }, + { + "task": "benchmark", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ "@stdlib/math/base/assert/is-nan", "@stdlib/constants/float64/ninf", "@stdlib/constants/float64/exponent-bias", "@stdlib/number/float64/base/get-high-word", "@stdlib/number/float64/base/set-high-word", "@stdlib/number/float64/base/set-low-word" - ] - }, - { - "task": "examples", - "src": [ - "./src/main.c" - ], - "include": [ - "./include" - ], - "libraries": [], - "libpath": [], - "dependencies": [ + ] + }, + { + "task": "examples", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ "@stdlib/math/base/assert/is-nan", "@stdlib/constants/float64/ninf", "@stdlib/constants/float64/exponent-bias", "@stdlib/number/float64/base/get-high-word", "@stdlib/number/float64/base/set-high-word", "@stdlib/number/float64/base/set-low-word" - ] - } - ] - } \ No newline at end of file + ] + } + ] +} From 3290ff70e5d98844c659e4099076d300284046a1 Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 10:53:16 +0530 Subject: [PATCH 08/20] fix indentation in manifest.json --- .../math/base/special/log10/manifest.json | 36 +++++++++---------- 1 file changed, 18 insertions(+), 18 deletions(-) diff --git a/lib/node_modules/@stdlib/math/base/special/log10/manifest.json b/lib/node_modules/@stdlib/math/base/special/log10/manifest.json index 6a8b58c4b818..90494d753517 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/manifest.json +++ b/lib/node_modules/@stdlib/math/base/special/log10/manifest.json @@ -37,12 +37,12 @@ "libpath": [], "dependencies": [ "@stdlib/math/base/napi/unary", - "@stdlib/math/base/assert/is-nan", - "@stdlib/constants/float64/ninf", - "@stdlib/constants/float64/exponent-bias", - "@stdlib/number/float64/base/get-high-word", - "@stdlib/number/float64/base/set-high-word", - "@stdlib/number/float64/base/set-low-word" + "@stdlib/math/base/assert/is-nan", + "@stdlib/constants/float64/ninf", + "@stdlib/constants/float64/exponent-bias", + "@stdlib/number/float64/base/get-high-word", + "@stdlib/number/float64/base/set-high-word", + "@stdlib/number/float64/base/set-low-word" ] }, { @@ -56,12 +56,12 @@ "libraries": [], "libpath": [], "dependencies": [ - "@stdlib/math/base/assert/is-nan", - "@stdlib/constants/float64/ninf", - "@stdlib/constants/float64/exponent-bias", - "@stdlib/number/float64/base/get-high-word", - "@stdlib/number/float64/base/set-high-word", - "@stdlib/number/float64/base/set-low-word" + "@stdlib/math/base/assert/is-nan", + "@stdlib/constants/float64/ninf", + "@stdlib/constants/float64/exponent-bias", + "@stdlib/number/float64/base/get-high-word", + "@stdlib/number/float64/base/set-high-word", + "@stdlib/number/float64/base/set-low-word" ] }, { @@ -75,12 +75,12 @@ "libraries": [], "libpath": [], "dependencies": [ - "@stdlib/math/base/assert/is-nan", - "@stdlib/constants/float64/ninf", - "@stdlib/constants/float64/exponent-bias", - "@stdlib/number/float64/base/get-high-word", - "@stdlib/number/float64/base/set-high-word", - "@stdlib/number/float64/base/set-low-word" + "@stdlib/math/base/assert/is-nan", + "@stdlib/constants/float64/ninf", + "@stdlib/constants/float64/exponent-bias", + "@stdlib/number/float64/base/get-high-word", + "@stdlib/number/float64/base/set-high-word", + "@stdlib/number/float64/base/set-low-word" ] } ] From 1a0186f27c4ee1845c2ad3dabd2bf03f4d4ab0d4 Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 10:54:35 +0530 Subject: [PATCH 09/20] add native.js --- .../math/base/special/log10/lib/native.js | 66 +++++++++++++++++++ 1 file changed, 66 insertions(+) create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/lib/native.js diff --git a/lib/node_modules/@stdlib/math/base/special/log10/lib/native.js b/lib/node_modules/@stdlib/math/base/special/log10/lib/native.js new file mode 100644 index 000000000000..a613782bbd8b --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/lib/native.js @@ -0,0 +1,66 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2022 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 // + +/** +* Evaluates the common logarithm (base ten). +* +* @private +* @param {NonNegativeNumber} x - input value +* @returns {number} function value +* +* @example +* var v = log10( 4.0 ); +* // returns ~0.602 +* +* @example +* var v = log10( 8.0 ); +* // returns ~0.903 +* +* @example +* var v = log10( 0.0 ); +* // returns -Infinity +* +* @example +* var v = log10( Infinity ); +* // returns Infinity +* +* @example +* var v = log10( NaN ); +* // returns NaN +* +* @example +* var v = log10( -4.0 ); +* // returns NaN +*/ +function log10( x ) { + return addon( x ); +} + + +// EXPORTS // + +module.exports = log10; From a630b5dda7aa5ac37983e77a94c14ba2c4692a0f Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 10:55:43 +0530 Subject: [PATCH 10/20] sync description with main.js --- .../base/special/log10/include/stdlib/math/base/special/log10.h | 2 +- lib/node_modules/@stdlib/math/base/special/log10/src/main.c | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h b/lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h index 6ab17d8dc34f..44362869af17 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h +++ b/lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h @@ -27,7 +27,7 @@ extern "C" { #endif /** -* Evaluates the common logarithm. +* Evaluates the common logarithm (base ten). */ double stdlib_base_log10( const double x ); diff --git a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c index ce201e77a09a..a868cca09c24 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c +++ b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c @@ -149,7 +149,7 @@ double klog( const double x ) { } /** -* Evaluates the common logarithm. +* Evaluates the common logarithm (base ten). * * @param x input value * @return output value From 949f1de6a2226a5281450d4cec98b89ea93d38ac Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 11:03:59 +0530 Subject: [PATCH 11/20] updates in main.c --- lib/node_modules/@stdlib/math/base/special/log10/src/main.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c index a868cca09c24..92aef1963665 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c +++ b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c @@ -193,7 +193,7 @@ double stdlib_base_log10( const double x ) { i = ( hx + 0x95f64 ) & 0x100000; // Normalize `x` or `x/2`... - stdlib_base_float64_set_high_word( xc, hx | ( i ^ HIGH_BIASED_EXP_0 ) ); + stdlib_base_float64_set_high_word( hx | ( i ^ HIGH_BIASED_EXP_0 ), &xc ); k += i >> 20; y = k; f = klog( xc ); From 64c111b1902803ae694f4092fe5b76267f66c3db Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 11:04:06 +0530 Subject: [PATCH 12/20] add c examples --- .../base/special/log10/examples/c/Makefile | 146 ++++++++++++++++++ .../base/special/log10/examples/c/example.c | 33 ++++ 2 files changed, 179 insertions(+) create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/examples/c/Makefile create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c diff --git a/lib/node_modules/@stdlib/math/base/special/log10/examples/c/Makefile b/lib/node_modules/@stdlib/math/base/special/log10/examples/c/Makefile new file mode 100644 index 000000000000..91d364d19fc3 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/examples/c/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2022 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. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + +# Define the program used for compiling C source files: +ifdef C_COMPILER + CC := $(C_COMPILER) +else + CC := gcc +endif + +# Define the command-line options when compiling C files: +CFLAGS ?= \ + -std=c99 \ + -O3 \ + -Wall \ + -pedantic + +# Determine whether to generate position independent code ([1][1], [2][2]). +# +# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options +# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option +ifeq ($(OS), WINNT) + fPIC ?= +else + fPIC ?= -fPIC +endif + +# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`): +INCLUDE ?= + +# List of source files: +SOURCE_FILES ?= + +# List of libraries (e.g., `-lopenblas -lpthread`): +LIBRARIES ?= + +# List of library paths (e.g., `-L /foo/bar -L /beep/boop`): +LIBPATH ?= + +# List of C targets: +c_targets := example.out + + +# RULES # + +#/ +# Compiles source files. +# +# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`) +# @param {string} [CFLAGS] - C compiler options +# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`) +# @param {string} [SOURCE_FILES] - list of source files +# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`) +# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`) +# +# @example +# make +# +# @example +# make all +#/ +all: $(c_targets) + +.PHONY: all + +#/ +# Compiles C source files. +# +# @private +# @param {string} CC - C compiler (e.g., `gcc`) +# @param {string} CFLAGS - C compiler options +# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`) +# @param {string} SOURCE_FILES - list of source files +# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`) +# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`) +#/ +$(c_targets): %.out: %.c + $(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES) + +#/ +# Runs compiled examples. +# +# @example +# make run +#/ +run: $(c_targets) + $(QUIET) ./$< + +.PHONY: run + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: + $(QUIET) -rm -f *.o *.out + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c b/lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c new file mode 100644 index 000000000000..7970a7971945 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c @@ -0,0 +1,33 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2022 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/log10.h" +#include +#include + +int main() { + double x; + double v; + int i; + + for ( i = 0; i < 100; i++ ) { + x = ( (double) rand() / (double) RAND_MAX ) * 100.0; + v = stdlib_base_log10( x ); + printf( "log10(%lf) = %lf\n", x, v ); + } +} \ No newline at end of file From 76eb082a6a521f6dab171496328c6db444f8945a Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 11:08:13 +0530 Subject: [PATCH 13/20] Add native c and js benchmarks --- .../log10/benchmark/benchmark.native.js | 61 ++++++++ .../special/log10/benchmark/c/native/Makefile | 146 ++++++++++++++++++ .../log10/benchmark/c/native/benchmark.c | 136 ++++++++++++++++ 3 files changed, 343 insertions(+) create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/benchmark/benchmark.native.js create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/benchmark/c/native/Makefile create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/benchmark/c/native/benchmark.c diff --git a/lib/node_modules/@stdlib/math/base/special/log10/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/math/base/special/log10/benchmark/benchmark.native.js new file mode 100644 index 000000000000..f69d7b157ce4 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/benchmark/benchmark.native.js @@ -0,0 +1,61 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2022 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 resolve = require( 'path' ).resolve; +var bench = require( '@stdlib/bench' ); +var randu = require( '@stdlib/random/base/randu' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var log10 = tryRequire( resolve( __dirname, './../lib/native.js' ) ); +var opts = { + 'skip': ( log10 instanceof Error ) +}; + + + +// MAIN // + +bench( pkg+'::native', opts, function benchmark( b ) { + var x; + var y; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + x = ( randu()*1000.0 ) - 0.0; + y = log10( x ); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); diff --git a/lib/node_modules/@stdlib/math/base/special/log10/benchmark/c/native/Makefile b/lib/node_modules/@stdlib/math/base/special/log10/benchmark/c/native/Makefile new file mode 100644 index 000000000000..d6b58c7f5d3e --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/benchmark/c/native/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2022 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. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + +# Define the program used for compiling C source files: +ifdef C_COMPILER + CC := $(C_COMPILER) +else + CC := gcc +endif + +# Define the command-line options when compiling C files: +CFLAGS ?= \ + -std=c99 \ + -O3 \ + -Wall \ + -pedantic + +# Determine whether to generate position independent code ([1][1], [2][2]). +# +# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options +# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option +ifeq ($(OS), WINNT) + fPIC ?= +else + fPIC ?= -fPIC +endif + +# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`): +INCLUDE ?= + +# List of source files: +SOURCE_FILES ?= + +# List of libraries (e.g., `-lopenblas -lpthread`): +LIBRARIES ?= + +# List of library paths (e.g., `-L /foo/bar -L /beep/boop`): +LIBPATH ?= + +# List of C targets: +c_targets := benchmark.out + + +# RULES # + +#/ +# Compiles source files. +# +# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`) +# @param {string} [CFLAGS] - C compiler options +# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`) +# @param {string} [SOURCE_FILES] - list of source files +# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`) +# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`) +# +# @example +# make +# +# @example +# make all +#/ +all: $(c_targets) + +.PHONY: all + +#/ +# Compiles C source files. +# +# @private +# @param {string} CC - C compiler (e.g., `gcc`) +# @param {string} CFLAGS - C compiler options +# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`) +# @param {string} SOURCE_FILES - list of source files +# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`) +# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`) +#/ +$(c_targets): %.out: %.c + $(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES) + +#/ +# Runs compiled benchmarks. +# +# @example +# make run +#/ +run: $(c_targets) + $(QUIET) ./$< + +.PHONY: run + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: + $(QUIET) -rm -f *.o *.out + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/math/base/special/log10/benchmark/c/native/benchmark.c b/lib/node_modules/@stdlib/math/base/special/log10/benchmark/c/native/benchmark.c new file mode 100644 index 000000000000..66f92c35aa10 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/benchmark/c/native/benchmark.c @@ -0,0 +1,136 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2022 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. +*/ + +/** +* Benchmark `log10`. +*/ +#include "stdlib/math/base/special/log10.h" +#include +#include +#include +#include +#include + +#define NAME "log10" +#define ITERATIONS 1000000 +#define REPEATS 3 + +/** +* Prints the TAP version. +*/ +void print_version() { + printf( "TAP version 13\n" ); +} + +/** +* Prints the TAP summary. +* +* @param total total number of tests +* @param passing total number of passing tests +*/ +void print_summary( int total, int passing ) { + printf( "#\n" ); + printf( "1..%d\n", total ); // TAP plan + printf( "# total %d\n", total ); + printf( "# pass %d\n", passing ); + printf( "#\n" ); + printf( "# ok\n" ); +} + +/** +* Prints benchmarks results. +* +* @param elapsed elapsed time in seconds +*/ +void print_results( double elapsed ) { + double rate = (double)ITERATIONS / elapsed; + printf( " ---\n" ); + printf( " iterations: %d\n", ITERATIONS ); + printf( " elapsed: %0.9f\n", elapsed ); + printf( " rate: %0.9f\n", rate ); + printf( " ...\n" ); +} + +/** +* Returns a clock time. +* +* @return clock time +*/ +double tic() { + struct timeval now; + gettimeofday( &now, NULL ); + return (double)now.tv_sec + (double)now.tv_usec/1.0e6; +} + +/** +* Generates a random number on the interval [0,1]. +* +* @return random number +*/ +double rand_double() { + int r = rand(); + return (double)r / ( (double)RAND_MAX + 1.0 ); +} + +/** +* Runs a benchmark. +* +* @return elapsed time in seconds +*/ +double benchmark() { + double elapsed; + double x; + double y; + double t; + int i; + + t = tic(); + for ( i = 0; i < ITERATIONS; i++ ) { + x = ( 1000.0*rand_double() ) - 0.0; + y = stdlib_base_log10( x ); + if ( y != y ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( y != y ) { + printf( "should not return NaN\n" ); + } + return elapsed; +} + +/** +* Main execution sequence. +*/ +int main( void ) { + double elapsed; + int i; + + // Use the current time to seed the random number generator: + srand( time( NULL ) ); + + print_version(); + for ( i = 0; i < REPEATS; i++ ) { + printf( "# c::native::%s\n", NAME ); + elapsed = benchmark(); + print_results( elapsed ); + printf( "ok %d benchmark finished\n", i+1 ); + } + print_summary( REPEATS, REPEATS ); +} From a67fe0f370fff1d1e327b3152a9a8fd2c314c264 Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 11:10:33 +0530 Subject: [PATCH 14/20] Add native tests --- .../base/special/log10/test/test.native.js | 235 ++++++++++++++++++ 1 file changed, 235 insertions(+) create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/test/test.native.js diff --git a/lib/node_modules/@stdlib/math/base/special/log10/test/test.native.js b/lib/node_modules/@stdlib/math/base/special/log10/test/test.native.js new file mode 100644 index 000000000000..cb3b02afd26d --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/test/test.native.js @@ -0,0 +1,235 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2022 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 resolve = require( 'path' ).resolve; +var tape = require( 'tape' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var PINF = require( '@stdlib/constants/float64/pinf' ); +var NINF = require( '@stdlib/constants/float64/ninf' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var EPS = require( '@stdlib/constants/float64/eps' ); +var abs = require( '@stdlib/math/base/special/abs' ); + + +// FIXTURES // + +var veryLargePositive = require( './fixtures/julia/very_large_positive.json' ); +var largePositive = require( './fixtures/julia/large_positive.json' ); +var mediumPositive = require( './fixtures/julia/medium_positive.json' ); +var smallPositive = require( './fixtures/julia/small_positive.json' ); +var smaller = require( './fixtures/julia/smaller.json' ); +var tinyPositive = require( './fixtures/julia/tiny_positive.json' ); +var subnormal = require( './fixtures/julia/subnormal.json' ); + + +// VARIABLES // + +var log10 = tryRequire( resolve( __dirname, './../lib/native.js' ) ); +var opts = { + 'skip': ( log10 instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.equal( typeof log10, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function evaluates the common logarithm of `x` (very large positive values)', opts, function test( t ) { + var expected; + var delta; + var tol; + var x; + var y; + var i; + + expected = veryLargePositive.expected; + x = veryLargePositive.x; + for ( i = 0; i < x.length; i++ ) { + y = log10( x[i] ); + if ( y === expected[i] ) { + t.equal( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); + } else { + delta = abs( y - expected[i] ); + tol = EPS * abs( expected[i] ); + t.equal( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); + } + } + t.end(); +}); + +tape( 'the function evaluates the common logarithm of `x` (large positive values)', opts, function test( t ) { + var expected; + var delta; + var tol; + var x; + var y; + var i; + + expected = largePositive.expected; + x = largePositive.x; + for ( i = 0; i < x.length; i++ ) { + y = log10( x[i] ); + if ( y === expected[i] ) { + t.equal( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); + } else { + delta = abs( y - expected[i] ); + tol = EPS * abs( expected[i] ); + t.equal( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); + } + } + t.end(); +}); + +tape( 'the function evaluates the common logarithm of `x` (medium positive values)', opts, function test( t ) { + var expected; + var delta; + var tol; + var x; + var y; + var i; + + expected = mediumPositive.expected; + x = mediumPositive.x; + for ( i = 0; i < x.length; i++ ) { + y = log10( x[i] ); + if ( y === expected[i] ) { + t.equal( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); + } else { + delta = abs( y - expected[i] ); + tol = EPS * abs( expected[i] ); + t.equal( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); + } + } + t.end(); +}); + +tape( 'the function evaluates the common logarithm of `x` (small positive values)', opts, function test( t ) { + var expected; + var delta; + var tol; + var x; + var y; + var i; + + expected = smallPositive.expected; + x = smallPositive.x; + for ( i = 0; i < x.length; i++ ) { + y = log10( x[i] ); + if ( y === expected[i] ) { + t.equal( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); + } else { + delta = abs( y - expected[i] ); + tol = EPS * abs( expected[i] ); + t.equal( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); + } + } + t.end(); +}); + +tape( 'the function evaluates the common logarithm of `x` (smaller values)', opts, function test( t ) { + var expected; + var delta; + var tol; + var x; + var y; + var i; + + expected = smaller.expected; + x = smaller.x; + for ( i = 0; i < x.length; i++ ) { + y = log10( x[i] ); + if ( y === expected[i] ) { + t.equal( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); + } else { + delta = abs( y - expected[i] ); + tol = EPS * abs( expected[i] ); + t.equal( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); + } + } + t.end(); +}); + +tape( 'the function evaluates the common logarithm of `x` (tiny positive values)', opts, function test( t ) { + var expected; + var delta; + var tol; + var x; + var y; + var i; + + expected = tinyPositive.expected; + x = tinyPositive.x; + for ( i = 0; i < x.length; i++ ) { + y = log10( x[i] ); + if ( y === expected[i] ) { + t.equal( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); + } else { + delta = abs( y - expected[i] ); + tol = EPS * abs( expected[i] ); + t.equal( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); + } + } + t.end(); +}); + +tape( 'the function evaluates the common logarithm of `x` (subnormal values)', opts, function test( t ) { + var expected; + var delta; + var tol; + var x; + var y; + var i; + + expected = subnormal.expected; + x = subnormal.x; + for ( i = 0; i < x.length; i++ ) { + y = log10( x[i] ); + if ( y === expected[i] ) { + t.equal( y, expected[i], 'x: '+x[i]+', y: '+y+', expected: '+expected[i] ); + } else { + delta = abs( y - expected[i] ); + tol = EPS * abs( expected[i] ); + t.equal( delta <= tol, true, 'within tolerance. x: '+x[i]+'. y: '+y+'. E: '+expected[i]+'. Δ: '+delta+'. Tolerance: '+tol+'.' ); + } + } + t.end(); +}); + +tape( 'the function returns `-infinity` if provided `0`', opts, function test( t ) { + t.equal( log10( 0.0 ), NINF, 'equals -infinity' ); + t.end(); +}); + +tape( 'the function returns `+infinity` if provided `+infinity`', opts, function test( t ) { + t.equal( log10( PINF ), PINF, 'equals +infinity' ); + t.end(); +}); + +tape( 'the function returns `NaN` if provided a negative number', opts, function test( t ) { + var v = log10( -1.0 ); + t.equal( isnan( v ), true, 'returns NaN' ); + t.end(); +}); From 2699dd153a923e2955e03b958e63ba8598225bd4 Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Tue, 20 Dec 2022 11:14:53 +0530 Subject: [PATCH 15/20] fix spacing issue in benchmark.native.js, example.c, log10.h and main.c --- .../math/base/special/log10/benchmark/benchmark.native.js | 1 - .../@stdlib/math/base/special/log10/examples/c/example.c | 2 +- .../special/log10/include/stdlib/math/base/special/log10.h | 2 +- lib/node_modules/@stdlib/math/base/special/log10/src/main.c | 3 --- 4 files changed, 2 insertions(+), 6 deletions(-) diff --git a/lib/node_modules/@stdlib/math/base/special/log10/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/math/base/special/log10/benchmark/benchmark.native.js index f69d7b157ce4..8818b8d85d54 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/math/base/special/log10/benchmark/benchmark.native.js @@ -36,7 +36,6 @@ var opts = { }; - // MAIN // bench( pkg+'::native', opts, function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c b/lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c index 7970a7971945..1aafa6a22f12 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c +++ b/lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c @@ -30,4 +30,4 @@ int main() { v = stdlib_base_log10( x ); printf( "log10(%lf) = %lf\n", x, v ); } -} \ No newline at end of file +} diff --git a/lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h b/lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h index 44362869af17..9a8bd513cecd 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h +++ b/lib/node_modules/@stdlib/math/base/special/log10/include/stdlib/math/base/special/log10.h @@ -35,4 +35,4 @@ double stdlib_base_log10( const double x ); } #endif -#endif // !STDLIB_MATH_BASE_SPECIAL_LOG10_H \ No newline at end of file +#endif // !STDLIB_MATH_BASE_SPECIAL_LOG10_H diff --git a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c index 92aef1963665..ce115959b425 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c +++ b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c @@ -204,6 +204,3 @@ double stdlib_base_log10( const double x ) { z += ( ( lo + f ) * IVLN10HI ) + ( hi * IVLN10HI ); return z + ( y * LOG10_2HI ); } - - - From 57c835123fb649df898dd2134fb2cb981c929901 Mon Sep 17 00:00:00 2001 From: Pranav <85227306+Pranavchiku@users.noreply.github.com> Date: Sun, 25 Dec 2022 23:35:38 +0530 Subject: [PATCH 16/20] Update lib/node_modules/@stdlib/math/base/special/log10/README.md --- lib/node_modules/@stdlib/math/base/special/log10/README.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/lib/node_modules/@stdlib/math/base/special/log10/README.md b/lib/node_modules/@stdlib/math/base/special/log10/README.md index 7f91e9c0f4a2..d4a960fe6898 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/README.md +++ b/lib/node_modules/@stdlib/math/base/special/log10/README.md @@ -57,7 +57,7 @@ Evaluates the [common logarithm][common-logarithm]. ```javascript var v = log10( 100.0 ); -// returns ~2.0 +// returns 2.0 v = log10( 8.0 ); // returns ~0.903 From 7f2872a8e4e4021e5acbbfaef8978775f49ed05c Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Fri, 26 Apr 2024 01:19:10 +0530 Subject: [PATCH 17/20] chore: apply code review --- .../@stdlib/math/base/special/log10/README.md | 6 +- .../base/special/log10/examples/c/example.c | 18 +- .../math/base/special/log10/lib/native.js | 2 +- .../math/base/special/log10/manifest.json | 12 +- .../math/base/special/log10/src/main.c | 185 +++++++++--------- 5 files changed, 119 insertions(+), 104 deletions(-) diff --git a/lib/node_modules/@stdlib/math/base/special/log10/README.md b/lib/node_modules/@stdlib/math/base/special/log10/README.md index d4a960fe6898..547a7a08c534 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/README.md +++ b/lib/node_modules/@stdlib/math/base/special/log10/README.md @@ -2,7 +2,7 @@ @license Apache-2.0 -Copyright (c) 2022 The Stdlib Authors. +Copyright (c) 2018 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. @@ -139,7 +139,7 @@ Evaluates the [common logarithm][common-logarithm]. ```c double out = stdlib_base_log10( 100.0 ); -// returns ~2.0 +// returns 2.0 out = stdlib_base_log10( 8.0 ); // returns ~0.903 @@ -182,7 +182,7 @@ int main() { int i; for ( i = 0; i < 100; i++ ) { - x = ( (double) rand() / (double) RAND_MAX ) * 100.0; + x = ( (double)rand() / (double)RAND_MAX ) * 100.0; v = stdlib_base_log10( x ); printf( "log10(%lf) = %lf\n", x, v ); } diff --git a/lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c b/lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c index 1aafa6a22f12..0a70bbbe0e1f 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c +++ b/lib/node_modules/@stdlib/math/base/special/log10/examples/c/example.c @@ -21,13 +21,13 @@ #include int main() { - double x; - double v; - int i; - - for ( i = 0; i < 100; i++ ) { - x = ( (double) rand() / (double) RAND_MAX ) * 100.0; - v = stdlib_base_log10( x ); - printf( "log10(%lf) = %lf\n", x, v ); - } + double x; + double v; + int i; + + for ( i = 0; i < 100; i++ ) { + x = ( (double)rand() / (double)RAND_MAX ) * 100.0; + v = stdlib_base_log10( x ); + printf( "log10(%lf) = %lf\n", x, v ); + } } diff --git a/lib/node_modules/@stdlib/math/base/special/log10/lib/native.js b/lib/node_modules/@stdlib/math/base/special/log10/lib/native.js index a613782bbd8b..c393a2df8127 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/lib/native.js +++ b/lib/node_modules/@stdlib/math/base/special/log10/lib/native.js @@ -57,7 +57,7 @@ var addon = require( './../src/addon.node' ); * // returns NaN */ function log10( x ) { - return addon( x ); + return addon( x ); } diff --git a/lib/node_modules/@stdlib/math/base/special/log10/manifest.json b/lib/node_modules/@stdlib/math/base/special/log10/manifest.json index 90494d753517..739ba3847e0d 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/manifest.json +++ b/lib/node_modules/@stdlib/math/base/special/log10/manifest.json @@ -37,11 +37,13 @@ "libpath": [], "dependencies": [ "@stdlib/math/base/napi/unary", - "@stdlib/math/base/assert/is-nan", "@stdlib/constants/float64/ninf", "@stdlib/constants/float64/exponent-bias", + "@stdlib/number/float64/base/to-words", "@stdlib/number/float64/base/get-high-word", "@stdlib/number/float64/base/set-high-word", + "@stdlib/constants/float64/high-word-abs-mask", + "@stdlib/constants/float64/high-word-significand-mask", "@stdlib/number/float64/base/set-low-word" ] }, @@ -56,11 +58,13 @@ "libraries": [], "libpath": [], "dependencies": [ - "@stdlib/math/base/assert/is-nan", "@stdlib/constants/float64/ninf", "@stdlib/constants/float64/exponent-bias", + "@stdlib/number/float64/base/to-words", "@stdlib/number/float64/base/get-high-word", "@stdlib/number/float64/base/set-high-word", + "@stdlib/constants/float64/high-word-abs-mask", + "@stdlib/constants/float64/high-word-significand-mask", "@stdlib/number/float64/base/set-low-word" ] }, @@ -75,11 +79,13 @@ "libraries": [], "libpath": [], "dependencies": [ - "@stdlib/math/base/assert/is-nan", "@stdlib/constants/float64/ninf", "@stdlib/constants/float64/exponent-bias", + "@stdlib/number/float64/base/to-words", "@stdlib/number/float64/base/get-high-word", "@stdlib/number/float64/base/set-high-word", + "@stdlib/constants/float64/high-word-abs-mask", + "@stdlib/constants/float64/high-word-significand-mask", "@stdlib/number/float64/base/set-low-word" ] } diff --git a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c index ce115959b425..e988638cd16e 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c +++ b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c @@ -31,10 +31,12 @@ */ #include "stdlib/math/base/special/log10.h" -#include "stdlib/math/base/assert/is_nan.h" +#include "stdlib/number/float64/base/to_words.h" #include "stdlib/number/float64/base/get_high_word.h" #include "stdlib/number/float64/base/set_high_word.h" #include "stdlib/number/float64/base/set_low_word.h" +#include "stdlib/constants/float64/high_word_abs_mask.h" +#include "stdlib/constants/float64/high_word_significand_mask.h" #include "stdlib/constants/float64/exponent_bias.h" #include "stdlib/constants/float64/ninf.h" @@ -44,17 +46,14 @@ static const double IVLN10LO = 2.50829467116452752298e-11; // 0x3dbb9438, 0xca9 static const double LOG10_2HI = 3.01029995663611771306e-01; // 0x3FD34413, 0x509F6000 static const double LOG10_2LO = 3.69423907715893078616e-13; // 0x3D59FEF3, 0x11F12B36 -// 0x000fffff = 1048575 => 0 00000000000 11111111111111111111 -static const int64_t HIGH_SIGNIFICAND_MASK = 0x000fffff; - // 0x7ff00000 = 2146435072 => 0 11111111111 00000000000000000000 => biased exponent: 2047 = 1023+1023 => 2^1023 -static const int64_t HIGH_MAX_NORMAL_EXP = 0x7ff00000; +static const int32_t HIGH_MAX_NORMAL_EXP = 0x7ff00000; // 0x00100000 = 1048576 => 0 00000000001 00000000000000000000 => biased exponent: 1 = -1022+1023 => 2^-1022 -static const int64_t HIGH_MIN_NORMAL_EXP = 0x00100000; +static const int32_t HIGH_MIN_NORMAL_EXP = 0x00100000; // 0x3ff00000 = 1072693248 => 0 01111111111 00000000000000000000 => biased exponent: 1023 = 0+1023 => 2^0 = 1 -static const int64_t HIGH_BIASED_EXP_0 = 0x3ff00000; +static const int32_t HIGH_BIASED_EXP_0 = 0x3ff00000; // 1/3 static const double ONE_THIRD = 0.33333333333333333; @@ -110,97 +109,107 @@ static double polyval_q( const double x ) { * @returns {number} function value */ double klog( const double x ) { - double hfsq; - double t1; - double t2; - uint32_t hx; - double f; - double s; - double z; - double R; - double w; - int32_t i; - int32_t j; - - stdlib_base_float64_get_high_word( x, &hx ); - f = x - 1.0; - if ( ( HIGH_SIGNIFICAND_MASK & ( 2 + hx ) ) < 3 ) { - // Case: -2**-20 <= f < 2**-20 - if ( f == 0.0 ) { - return 0.0; - } - return f * f * ( ( ONE_THIRD * f ) - 0.5 ); - } - s = f / ( 2.0 + f ); - z = s * s; - hx &= HIGH_SIGNIFICAND_MASK; - i = (hx - 0x6147a); - w = z * z; - j = ( 0x6b851 - hx ); - t1 = w * polyval_p( w ); - t2 = z * polyval_q( w ); - i |= j; - R = t2 + t1; - if ( i > 0 ) { - hfsq = 0.5 * f * f; - return ( s * ( hfsq + R ) ) - hfsq; - } - return s * ( R - f ); + double hfsq; + uint32_t hx; + int32_t ihx; + int32_t i; + int32_t j; + double t1; + double t2; + double f; + double s; + double z; + double R; + double w; + + stdlib_base_float64_get_high_word( x, &hx ); + ihx = (int32_t)hx; + f = x - 1.0; + if ( ( STDLIB_CONSTANT_FLOAT64_HIGH_WORD_SIGNIFICAND_MASK & ( 2 + ihx ) ) < 3 ) { + // Case: -2**-20 <= f < 2**-20 + if ( f == 0.0 ) { + return 0.0; + } + return f * f * ( ( ONE_THIRD * f ) - 0.5 ); + } + s = f / ( 2.0 + f ); + z = s * s; + ihx &= STDLIB_CONSTANT_FLOAT64_HIGH_WORD_SIGNIFICAND_MASK; + i = (ihx - 0x6147a); + w = z * z; + j = ( 0x6b851 - ihx ); + t1 = w * polyval_p( w ); + t2 = z * polyval_q( w ); + i |= j; + R = t2 + t1; + if ( i > 0 ) { + hfsq = 0.5 * f * f; + return ( s * ( hfsq + R ) ) - hfsq; + } + return s * ( R - f ); } /** * Evaluates the common logarithm (base ten). * * @param x input value -* @return output value +* @return output value * * @example * double out = stdlib_base_log10( 4.0 ); * // returns 0.6020599913279624 */ double stdlib_base_log10( const double x ) { - double xc = x; - double hi; - uint32_t hx; - double lo; - double f; - int32_t i; - int32_t k; - int32_t y; - double z; - - if ( stdlib_base_is_nan( xc ) || xc < 0.0 ) { - return 0.0 / 0.0; // NaN - } - if ( xc == 0.0 ) { - return STDLIB_CONSTANT_FLOAT64_NINF; - } - stdlib_base_float64_get_high_word( xc, &hx ); - k = 0; - - // Case: 0 < x < 2**-1022 - if ( hx < HIGH_MIN_NORMAL_EXP ) { - // Subnormal number, scale up `x`... - k -= 54; - xc *= TWO54; - stdlib_base_float64_get_high_word( xc, &hx ); - } - if ( hx >= HIGH_MAX_NORMAL_EXP ) { - return x + x; - } - k += ( ( hx>>20 ) - STDLIB_CONSTANT_FLOAT64_EXPONENT_BIAS ); - hx &= HIGH_SIGNIFICAND_MASK; - i = ( hx + 0x95f64 ) & 0x100000; - - // Normalize `x` or `x/2`... - stdlib_base_float64_set_high_word( hx | ( i ^ HIGH_BIASED_EXP_0 ), &xc ); - k += i >> 20; - y = k; - f = klog( xc ); - xc -= 1.0; - stdlib_base_float64_set_low_word( 0, &hi ); - lo = xc - hi; - z = ( y * LOG10_2LO ) + ( ( xc + f ) * IVLN10LO ); - z += ( ( lo + f ) * IVLN10HI ) + ( hi * IVLN10HI ); - return z + ( y * LOG10_2HI ); + int32_t ihx; + int32_t ilx; + uint32_t hx; + uint32_t lx; + int32_t i; + int32_t k; + double xc; + double hi; + double lo; + double f; + double y; + double z; + + xc = x; + stdlib_base_float64_to_words( x, &hx, &lx ); + ihx = (int32_t)hx; + ilx = (int32_t)lx; + k = 0; + + // Case: 0 < x < 2**-1022 + if ( ihx < HIGH_MIN_NORMAL_EXP ) { + if ( ( ( ihx & STDLIB_CONSTANT_FLOAT64_HIGH_WORD_ABS_MASK ) | ilx ) == 0 ) { + return STDLIB_CONSTANT_FLOAT64_NINF; + } + if ( ihx < 0 ) { + return 0.0 / 0.0; // NaN + } + // Subnormal number, scale up `x`... + k -= 54; + xc *= TWO54; + stdlib_base_float64_get_high_word( xc, &hx ); + ihx = (int32_t)hx; + } + if ( ihx >= HIGH_MAX_NORMAL_EXP ) { + return x + x; + } + k += ( ( ihx>>20 ) - STDLIB_CONSTANT_FLOAT64_EXPONENT_BIAS ); + ihx &= STDLIB_CONSTANT_FLOAT64_HIGH_WORD_SIGNIFICAND_MASK; + i = ( ihx + 0x95f64 ) & 0x100000; + + // Normalize `x` or `x/2`... + stdlib_base_float64_set_high_word( ihx | ( i ^ HIGH_BIASED_EXP_0 ), &xc ); + k += i >> 20; + y = (double)k; + f = klog( xc ); + xc -= 1.0; + hi = xc; + stdlib_base_float64_set_low_word( 0, &hi ); + lo = xc - hi; + z = ( y * LOG10_2LO ) + ( ( xc + f ) * IVLN10LO ); + z += ( ( lo + f ) * IVLN10HI ) + ( hi * IVLN10HI ); + return z + ( y * LOG10_2HI ); } From e52c3cfbb646ec75a5e7f53dcf107346361427dd Mon Sep 17 00:00:00 2001 From: Pranavchiku Date: Fri, 26 Apr 2024 10:11:28 +0530 Subject: [PATCH 18/20] chore: add evalpoly.js --- .../base/special/log10/scripts/evalpoly.js | 115 ++++++++++++++++++ 1 file changed, 115 insertions(+) create mode 100644 lib/node_modules/@stdlib/math/base/special/log10/scripts/evalpoly.js diff --git a/lib/node_modules/@stdlib/math/base/special/log10/scripts/evalpoly.js b/lib/node_modules/@stdlib/math/base/special/log10/scripts/evalpoly.js new file mode 100644 index 000000000000..bdb13f4573ef --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/log10/scripts/evalpoly.js @@ -0,0 +1,115 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2018 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. +*/ + +/* +* This script compiles modules for evaluating polynomial functions. If any polynomial coefficients change, this script should be rerun to update the compiled files. +*/ +'use strict'; + +// MODULES // + +var resolve = require( 'path' ).resolve; +var readFileSync = require( '@stdlib/fs/read-file' ).sync; +var writeFileSync = require( '@stdlib/fs/write-file' ).sync; +var substringBefore = require( '@stdlib/string/substring-before' ); +var substringAfter = require( '@stdlib/string/substring-after' ); +var format = require( '@stdlib/string/format' ); +var compileC = require( '@stdlib/math/base/tools/evalpoly-compile-c' ); + + +// VARIABLES // + +// Polynomial coefficients ordered in ascending degree... +var P = [ + 3.999999999940941908e-01, // 3FD99999 9997FA04 + 2.222219843214978396e-01, // 3FCC71C5 1D8E78AF + 1.531383769920937332e-01 // 3FC39A09 D078C69F +]; +var Q = [ + 6.666666666666735130e-01, // 3FE55555 55555593 + 2.857142874366239149e-01, // 3FD24924 94229359 + 1.818357216161805012e-01, // 3FC74664 96CB03DE + 1.479819860511658591e-01 // 3FC2F112 DF3E5244 +]; + + +// FUNCTIONS // + +/** +* Inserts a compiled function into file content. +* +* @private +* @param {string} text - source content +* @param {string} id - function identifier +* @param {string} str - function string +* @returns {string} updated content +*/ +function insert( text, id, str ) { + var before; + var after; + var begin; + var end; + + begin = '// BEGIN: '+id; + end = '// END: '+id; + + before = substringBefore( text, begin ); + after = substringAfter( text, end ); + + return format( '%s// BEGIN: %s\n\n%s\n%s%s', before, id, str, end, after ); +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var fpath; + var copts; + var opts; + var file; + var str; + + opts = { + 'encoding': 'utf8' + }; + + copts = { + 'dtype': 'double', + 'name': '' + }; + + fpath = resolve( __dirname, '..', 'src', 'main.c' ); + file = readFileSync( fpath, opts ); + + copts.name = 'polyval_p'; + str = compileC( P, copts ); + file = insert( file, copts.name, str ); + + copts.name = 'polyval_q'; + str = compileC( Q, copts ); + file = insert( file, copts.name, str ); + + writeFileSync( fpath, file, opts ); +} + +main(); From df079d559d96a87bbe521629d55f024b0f70eba3 Mon Sep 17 00:00:00 2001 From: Athan Date: Sat, 27 Apr 2024 14:48:57 -0700 Subject: [PATCH 19/20] Apply suggestions from code review Signed-off-by: Athan --- lib/node_modules/@stdlib/math/base/special/log10/src/main.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c index e988638cd16e..c826012f7c52 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c +++ b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c @@ -108,7 +108,7 @@ static double polyval_q( const double x ) { * @param {number} x - input value * @returns {number} function value */ -double klog( const double x ) { +static double klog( const double x ) { double hfsq; uint32_t hx; int32_t ihx; From 3c8add49dab207400cc1b3fa503b9c6c3454becd Mon Sep 17 00:00:00 2001 From: Athan Date: Sat, 27 Apr 2024 14:49:42 -0700 Subject: [PATCH 20/20] Apply suggestions from code review Signed-off-by: Athan --- lib/node_modules/@stdlib/math/base/special/log10/src/main.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c index c826012f7c52..db75c0fd5c03 100644 --- a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c +++ b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c @@ -157,7 +157,7 @@ static double klog( const double x ) { * * @example * double out = stdlib_base_log10( 4.0 ); -* // returns 0.6020599913279624 +* // returns ~0.602 */ double stdlib_base_log10( const double x ) { int32_t ihx;