diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/README.md b/lib/node_modules/@stdlib/math/base/special/factorialf/README.md new file mode 100644 index 000000000000..b37035a2ddb9 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/README.md @@ -0,0 +1,227 @@ + + +# factorialf + +> Evaluate the [Factorial][factorial-function] function of a single-precision floating-point number. + +
+ +The [factorial][factorial-function] function may be defined as the product + +```math +n! = \prod_{k=1}^n k +``` + +or according to the recurrence relation + +```math +n! = \begin{cases}1 & \textrm{if } n = 0,\\(n-1)! \times n & \textrm{if } n > 1\end{cases} +``` + +Following the convention for an [empty product][empty-product], in all definitions, + +```math +0! = 1 +``` + +The [Gamma][@stdlib/math/base/special/gamma] function extends the [factorial][factorial-function] function for non-integer values. + +```math +n! = \Gamma(n+1) +``` + +The [factorial][factorial-function] of a **negative** integer is not defined. + +
+ + + +
+ +## Usage + +```javascript +var factorialf = require( '@stdlib/math/base/special/factorialf' ); +``` + +#### factorialf( x ) + +Evaluates the [factorial][factorial-function] function of a single-precision floating-point number. + +```javascript +var v = factorialf( 3.0 ); +// returns 6.0 + +v = factorialf( -1.5 ); +// returns ~-3.545 + +v = factorialf( -0.5 ); +// returns ~1.772 + +v = factorialf( 0.5 ); +// returns ~0.886 + +v = factorialf( -10.0 ); +// returns NaN + +v = factorialf( 34.0 ); +// returns ~2.952e38 + +v = factorialf( 35.0 ); +// returns Infinity + +v = factorialf( NaN ); +// returns NaN +``` + +
+ + + +
+ +## Examples + + + +```javascript +var incrspace = require( '@stdlib/array/base/incrspace' ); +var factorialf = require( '@stdlib/math/base/special/factorialf' ); + +var x = incrspace( -10.0, 50.0, 1.0 ); + +var i; +for ( i = 0; i < x.length; i++ ) { + console.log( 'x: %d, f(x): %d', x[ i ], factorialf( x[ i ] ) ); +} +``` + +
+ + + + + +* * * + +
+ +## C APIs + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/math/base/special/factorialf.h" +``` + +#### stdlib_base_factorialf( x ) + +Evaluates the [factorial][factorial-function] function of a single-precision floating-point number. + +```c +float out = stdlib_base_factorialf( 3.0f ); +// returns 6.0f + +out = stdlib_base_factorialf( -1.5f ); +// returns ~-3.545f +``` + +The function accepts the following arguments: + +- **x**: `[in] float` input value. + +```c +float stdlib_base_factorialf( const float x ); +``` + +
+ + + +
+ +
+ + + +
+ +### Examples + +```c +#include "stdlib/math/base/special/factorialf.h" +#include + +int main( void ) { + const float x[] = { 2.0f, 3.0f, 5.0f, 8.0f }; + + float y; + int i; + for ( i = 0; i < 4; i++ ) { + y = stdlib_base_factorialf( x[ i ] ); + printf( "factorialf(%f) = %f\n", x[ i ], y ); + } +} +``` + +
+ + + +
+ + + + + +
+ +
+ + + + + +
+ +[@stdlib/math/base/special/gamma]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/math/base/special/gamma + +[factorial-function]: https://en.wikipedia.org/wiki/Factorial + +[empty-product]: https://en.wikipedia.org/wiki/Empty_product + + + + + +
+ + diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/benchmark/benchmark.js b/lib/node_modules/@stdlib/math/base/special/factorialf/benchmark/benchmark.js new file mode 100644 index 000000000000..5768b7022e86 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/benchmark/benchmark.js @@ -0,0 +1,68 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 bench = require( '@stdlib/bench' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; +var factorialf = require( './../lib' ); + + +// MAIN // + +bench( format( '%s::integers', pkg ), function benchmark( b ) { + var y; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = factorialf( i%35 ); + if ( isnanf( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); + +bench( format( '%s::decimals', pkg ), function benchmark( b ) { + var y; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = factorialf( (i%35)*0.68163 ); + if ( isnanf( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/math/base/special/factorialf/benchmark/benchmark.native.js new file mode 100644 index 000000000000..60e8fff758cf --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/benchmark/benchmark.native.js @@ -0,0 +1,77 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var factorialf = tryRequire( resolve( __dirname, './../lib/native.js' ) ); +var opts = { + 'skip': ( factorialf instanceof Error ) +}; + + +// MAIN // + +bench( format( '%s::native,integers', pkg ), opts, function benchmark( b ) { + var y; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = factorialf( i % 171 ); + if ( isnanf( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); + +bench( format( '%s::native,decimals', pkg ), opts, function benchmark( b ) { + var y; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = factorialf( ( i % 171 ) * 0.68163 ); + if ( isnanf( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/benchmark/c/native/Makefile b/lib/node_modules/@stdlib/math/base/special/factorialf/benchmark/c/native/Makefile new file mode 100644 index 000000000000..979768abbcec --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/benchmark/c/native/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 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/factorialf/benchmark/c/native/benchmark.c b/lib/node_modules/@stdlib/math/base/special/factorialf/benchmark/c/native/benchmark.c new file mode 100644 index 000000000000..6e0d4095c74c --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/benchmark/c/native/benchmark.c @@ -0,0 +1,139 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/factorialf.h" +#include +#include +#include +#include +#include +#include + +#define NAME "factorialf" +#define ITERATIONS 1000000 +#define REPEATS 3 + +/** +* Prints the TAP version. +*/ +static void print_version( void ) { + printf( "TAP version 13\n" ); +} + +/** +* Prints the TAP summary. +* +* @param total total number of tests +* @param passing total number of passing tests +*/ +static 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 +*/ +static 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 +*/ +static double tic( void ) { + struct timeval now; + gettimeofday( &now, NULL ); + return (double)now.tv_sec + (double)now.tv_usec/1.0e6; +} + +/** +* Generates a random number on the interval [min,max). +* +* @param min minimum value (inclusive) +* @param max maximum value (exclusive) +* @return random number +*/ +static float random_uniform( const float min, const float max ) { + float v = (float)rand() / ( (float)RAND_MAX + 1.0f ); + return min + ( v*(max-min) ); +} + +/** +* Runs a benchmark. +* +* @return elapsed time in seconds +*/ +static double benchmark( void ) { + float x[ 100 ]; + double elapsed; + double t; + float y; + int i; + + for ( i = 0; i < 100; i++ ) { + x[ i ] = random_uniform( 0.0f, 34.0f ); + } + + t = tic(); + for ( i = 0; i < ITERATIONS; i++ ) { + y = stdlib_base_factorialf( x[ i%100 ] ); + 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 ); +} diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/binding.gyp b/lib/node_modules/@stdlib/math/base/special/factorialf/binding.gyp new file mode 100644 index 000000000000..0d6508a12e99 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/binding.gyp @@ -0,0 +1,170 @@ +# @license Apache-2.0 +# +# Copyright (c) 2026 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 +} diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_factorial_function.svg b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_factorial_function.svg new file mode 100644 index 000000000000..b1736ddb4b25 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_factorial_function.svg @@ -0,0 +1,26 @@ + +n factorial equals product Underscript k equals 1 Overscript n Endscripts k + + + \ No newline at end of file diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_factorial_function_and_gamma.svg b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_factorial_function_and_gamma.svg new file mode 100644 index 000000000000..bec85bef1e9e --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_factorial_function_and_gamma.svg @@ -0,0 +1,24 @@ + +n factorial equals normal upper Gamma left-parenthesis n plus 1 right-parenthesis + + + \ No newline at end of file diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_factorial_recurrence_relation.svg b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_factorial_recurrence_relation.svg new file mode 100644 index 000000000000..40b655b7c293 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_factorial_recurrence_relation.svg @@ -0,0 +1,60 @@ + +n factorial equals StartLayout Enlarged left-brace 1st Row 1st Column 1 2nd Column if n equals 0 comma 2nd Row 1st Column left-parenthesis n minus 1 right-parenthesis factorial times n 2nd Column if n greater-than 1 EndLayout + + + \ No newline at end of file diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_zero_factorial.svg b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_zero_factorial.svg new file mode 100644 index 000000000000..6bb43a008850 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/img/equation_zero_factorial.svg @@ -0,0 +1,15 @@ + +0 factorial equals 1 + + + \ No newline at end of file diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/docs/repl.txt b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/repl.txt new file mode 100644 index 000000000000..eaa9257b0428 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/repl.txt @@ -0,0 +1,45 @@ + +{{alias}}( x ) + Evaluates the factorial function of a single-precision floating-point + number. + + If `x` is greater than `34`, the function returns `+Infinity`, as larger + factorial values overflow single-precision floating-point format. + + For non-integer inputs, the function returns $\Gamma(x+1)$, where $\Gamma$ + is the Gamma function. + + The factorial of a negative integer is not defined and results in `NaN`. + + If provided `NaN`, the function returns `NaN`. + + Parameters + ---------- + x: number + Input value. + + Returns + ------- + y: number + Factorial. + + Examples + -------- + > var y = {{alias}}( 3.0 ) + 6.0 + > y = {{alias}}( -1.5 ) + ~-3.545 + > y = {{alias}}( -0.5 ) + ~1.772 + > y = {{alias}}( 0.5 ) + ~0.886 + > y = {{alias}}( -10.0 ) + NaN + > y = {{alias}}( 35.0 ) + Infinity + > y = {{alias}}( NaN ) + NaN + + See Also + -------- + diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/docs/types/index.d.ts b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/types/index.d.ts new file mode 100644 index 000000000000..78286f0d9497 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/types/index.d.ts @@ -0,0 +1,66 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2026 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. +*/ + +// TypeScript Version: 4.1 + +/** +* Evaluates the factorial of `x` as a single-precision floating-point number. +* +* ## Notes +* +* - If `x` is greater than `34` and an integer, the function returns `+Infinity`, as larger factorial values overflow single-precision floating-point format. +* - For non-integer inputs, the function returns $\Gamma(x+1)$, where $\Gamma$ is the Gamma function. +* - The factorial of a negative integer is not defined and results in `NaN`. +* +* @param x - input value +* @returns factorial +* +* @example +* var v = factorialf( 3.0 ); +* // returns 6.0 +* +* @example +* var v = factorialf( -1.5 ); +* // returns ~-3.545 +* +* @example +* var v = factorialf( -0.5 ); +* // returns ~1.772 +* +* @example +* var v = factorialf( 0.5 ); +* // returns ~0.886 +* +* @example +* var v = factorialf( -10.0 ); +* // returns NaN +* +* @example +* var v = factorialf( 35.0 ); +* // returns Infinity +* +* @example +* var v = factorialf( NaN ); +* // returns NaN +*/ +declare function factorialf( x: number ): number; + + +// EXPORTS // + +export = factorialf; diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/docs/types/test.ts b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/types/test.ts new file mode 100644 index 000000000000..b6bf86120794 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/docs/types/test.ts @@ -0,0 +1,44 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2026 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. +*/ + +import factorialf = require( './index' ); + + +// TESTS // + +// The function returns a number... +{ + factorialf( 0.5 ); // $ExpectType number +} + +// The compiler throws an error if the function is provided a value other than a number... +{ + factorialf( true ); // $ExpectError + factorialf( false ); // $ExpectError + factorialf( null ); // $ExpectError + factorialf( undefined ); // $ExpectError + factorialf( '5' ); // $ExpectError + factorialf( [] ); // $ExpectError + factorialf( {} ); // $ExpectError + factorialf( ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the function is provided insufficient arguments... +{ + factorialf(); // $ExpectError +} diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/examples/c/Makefile b/lib/node_modules/@stdlib/math/base/special/factorialf/examples/c/Makefile new file mode 100644 index 000000000000..c8f8e9a1517b --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/examples/c/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 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/factorialf/examples/c/example.c b/lib/node_modules/@stdlib/math/base/special/factorialf/examples/c/example.c new file mode 100644 index 000000000000..64cf831d266c --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/examples/c/example.c @@ -0,0 +1,31 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/factorialf.h" +#include + +int main( void ) { + const float x[] = { 2.0f, 3.0f, 5.0f, 8.0f }; + + float y; + int i; + for ( i = 0; i < 4; i++ ) { + y = stdlib_base_factorialf( x[ i ] ); + printf( "factorialf(%f) = %f\n", x[ i ], y ); + } +} diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/examples/index.js b/lib/node_modules/@stdlib/math/base/special/factorialf/examples/index.js new file mode 100644 index 000000000000..9f26b3e88376 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/examples/index.js @@ -0,0 +1,29 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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'; + +var incrspace = require( '@stdlib/array/base/incrspace' ); +var factorialf = require( './../lib' ); + +var x = incrspace( -10.0, 50.0, 1.0 ); + +var i; +for ( i = 0; i < x.length; i++ ) { + console.log( 'x: %d, f(x): %d', x[ i ], factorialf( x[ i ] ) ); +} diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/include.gypi b/lib/node_modules/@stdlib/math/base/special/factorialf/include.gypi new file mode 100644 index 000000000000..bee8d41a2caf --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/include.gypi @@ -0,0 +1,53 @@ +# @license Apache-2.0 +# +# Copyright (c) 2026 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': [ + '=0.10.0", + "npm": ">2.7.0" + }, + "os": [ + "aix", + "darwin", + "freebsd", + "linux", + "macos", + "openbsd", + "sunos", + "win32", + "windows" + ], + "keywords": [ + "stdlib", + "stdmath", + "mathematics", + "math", + "special functions", + "special", + "function", + "factorial", + "fact", + "combinatorics", + "gamma", + "number" + ], + "__stdlib__": { + "scaffold": { + "$schema": "math/base@v1.0", + "base_alias": "factorial", + "alias": "factorialf", + "pkg_desc": "compute the factorial of a single-precision floating-point number", + "desc": "computes the factorial of a single-precision floating-point number", + "short_desc": "factorial", + "parameters": [ + { + "name": "x", + "desc": "input value", + "type": { + "javascript": "number", + "jsdoc": "number", + "c": "float", + "dtype": "float32" + }, + "domain": [ + { + "min": 0, + "max": "infinity" + } + ], + "rand": { + "prng": "random/base/discrete-uniform", + "parameters": [ + 0, + 100 + ] + }, + "example_values": [ + 1, + 2, + 3, + 4, + 5, + 6, + 7, + 8, + 9, + 10, + 11, + 12, + 13, + 14, + 15, + 16, + 17, + 21, + 34 + ] + } + ], + "output_policy": "real_floating_point_and_generic", + "returns": { + "desc": "factorial", + "type": { + "javascript": "number", + "jsdoc": "number", + "c": "float", + "dtype": "float32" + } + }, + "keywords": [ + "factorial", + "fact" + ], + "extra_keywords": [] + } + } +} diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/src/Makefile b/lib/node_modules/@stdlib/math/base/special/factorialf/src/Makefile new file mode 100644 index 000000000000..2caf905cedbe --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/src/Makefile @@ -0,0 +1,70 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 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/factorialf/src/addon.c b/lib/node_modules/@stdlib/math/base/special/factorialf/src/addon.c new file mode 100644 index 000000000000..100a45e68fc9 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/src/addon.c @@ -0,0 +1,22 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/factorialf.h" +#include "stdlib/math/base/napi/unary.h" + +STDLIB_MATH_BASE_NAPI_MODULE_F_F( stdlib_base_factorialf ) diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/src/main.c b/lib/node_modules/@stdlib/math/base/special/factorialf/src/main.c new file mode 100644 index 000000000000..9ba35e8e4486 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/src/main.c @@ -0,0 +1,90 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/factorialf.h" +#include "stdlib/math/base/assert/is_nanf.h" +#include "stdlib/math/base/assert/is_integerf.h" +#include "stdlib/math/base/special/gamma.h" +#include "stdlib/constants/float32/pinf.h" +#include "stdlib/constants/float32/max_nth_factorial.h" +#include + +static const float FACTORIALS[ 35 ] = { + 1.0f, + 1.0f, + 2.0f, + 6.0f, + 24.0f, + 120.0f, + 720.0f, + 5040.0f, + 40320.0f, + 362880.0f, + 3628800.0f, + 39916800.0f, + 479001600.0f, + 6227020800.0f, + 87178291200.0f, + 1307674368000.0f, + 20922789888000.0f, + 355687428096000.0f, + 6402373705728000.0f, + 0.12164510e18f, + 0.24329020e19f, + 0.51090942e20f, + 0.11240007e22f, + 0.25852017e23f, + 0.62044840e24f, + 0.15511210e26f, + 0.40329146e27f, + 0.10888869e29f, + 0.30488834e30f, + 0.88417620e31f, + 0.26525286e33f, + 0.82228387e34f, + 0.26313084e36f, + 0.86833176e37f, + 0.29523280e39f +}; + + +/** +* Evaluates the factorial of `x` as a single-precision floating-point number. +* +* @param x input value +* @return factorial +* +* @example +* float out = stdlib_base_factorialf( 3.0f ); +* // returns 6.0f +*/ +float stdlib_base_factorialf( const float x ) { + if ( stdlib_base_is_nanf( x ) ) { + return 0.0f / 0.0f; // NaN + } + if ( stdlib_base_is_integerf( x ) ) { + if ( x < 0.0f ) { + return 0.0f / 0.0f; // NaN + } + if ( x <= STDLIB_CONSTANT_FLOAT32_MAX_NTH_FACTORIAL ) { + return FACTORIALS[ (int32_t)x ]; + } + return STDLIB_CONSTANT_FLOAT32_PINF; + } + return (float)stdlib_base_gamma( (double)x + 1.0 ); +} diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/REQUIRE b/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/REQUIRE new file mode 100644 index 000000000000..401d9d45732f --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/REQUIRE @@ -0,0 +1,2 @@ +JSON +SpecialFunctions diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/decimals.json b/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/decimals.json new file mode 100644 index 000000000000..b2c6ed102b5f --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/decimals.json @@ -0,0 +1 @@ +{"expected":[-8.379156073085479e-307,2.1722521399420905e-306,8.848780641324387e-306,2.0167215420957614e-304,-3.8511398573246966e-304,-1.708431205503738e-303,-5.226292007175996e-302,6.843897604264987e-302,3.2949164116299437e-301,1.5741451988986375e-299,-1.217934412484103e-299,-6.349254256498959e-299,-7.518114721948931e-297,2.1687646687495735e-297,1.22282225340586e-296,-2.464767095480898e-294,-3.861886596783345e-295,-2.3547056435463546e-294,1.2059550370539137e-292,6.873351037412174e-293,4.535889495400768e-292,-1.2759268240010285e-290,-1.2222069652075396e-290,-8.7462989575318e-290,1.6370857207694982e-288,2.1706388459179076e-288,1.689621539979808e-287,-2.293242208173914e-286,-3.84932637987738e-286,-3.2736973390445827e-285,3.374060773580803e-284,6.814725147940859e-284,6.371042776717613e-283,-5.117895951893993e-282,-1.2042370317630222e-281,-1.2478797377895167e-280,7.920330095678027e-280,2.123874345957367e-279,2.466809277694776e-278,-1.2425703822159714e-277,-3.738269139374552e-277,-4.941408956558281e-276,1.967813620650582e-275,6.566395831778869e-275,1.009203237775541e-273,-3.1365538385683894e-273,-1.1510828864926268e-272,-2.1224251837111807e-271,5.021137980058961e-271,2.0138956002367626e-270,4.678317716928544e-269,-8.060128880847373e-269,-3.516958057098194e-268,-1.1207886752741209e-266,1.2958191586763275e-266,6.131551984792092e-266,3.2126277006067587e-264,-2.0844796321628222e-264,-1.0674504857593635e-263,-1.7802508779502918e-261,3.3525525510742503e-262,1.856238578765711e-261,-2.2923931973858397e-259,-5.387866462032551e-260,-3.2255546848573236e-259,1.3753224509630425e-257,8.647892104667782e-258,5.603860560539522e-257,-1.3870643507296015e-255,-1.3857499586389032e-255,-9.740411705483564e-255,1.640570980596484e-253,2.2161660602014776e-253,1.695344750905834e-252,-2.094424879492834e-251,-3.5363199731075124e-251,-2.958244789042545e-250,2.7928852741675717e-249,5.629181318361077e-249,5.18298734588412e-248,-3.8268501771544535e-247,-8.937517379990193e-247,-9.137300285201595e-246,5.3377160701781425e-245,1.4152048913798343e-244,1.625713088698119e-243,-7.534282651980737e-243,-2.2347223955527005e-242,-2.93192507755061e-241,1.071996909236231e-240,3.5189963919625765e-240,5.395829069139462e-239,-1.5332453508070398e-238,-5.5260230804867554e-238,-1.0246225684386004e-236,2.1999971781571104e-236,8.65431175060507e-236,2.0484863458316603e-234,-3.16204437259928e-234,-1.3518464045609678e-233,-4.5018996475239544e-232,4.547192473860956e-232,2.1065452378685717e-231,1.2304487173156384e-229,-6.5365890454357945e-230,-3.2754078414695346e-229,-9.812920299921034e-227,9.38587983807058e-228,5.083331654464993e-227,-4.4032955114235166e-225,-1.345425193148173e-225,-7.877715133862347e-225,2.832756635462507e-223,1.9244094821288897e-223,1.2197073116348132e-222,-2.6619664671000147e-221,-2.7454884600099386e-221,-1.8880783558983376e-220,2.858453409740057e-219,3.905603295938481e-219,2.9247801023578337e-218,-3.280709064452583e-217,-5.538492196236871e-217,-4.539460633487126e-216,3.913018429031141e-215,7.827827762611232e-215,7.070739842168505e-214,-4.7803747803950456e-213,-1.1024757589362582e-212,-1.1077830175521957e-211,5.9310408716712e-211,1.5471254954491933e-210,1.7513127437798484e-209,-7.433124861898638e-209,-2.1631141140986784e-208,-2.8070132668447455e-207,9.375498385649275e-207,3.01312644994807e-206,4.595109233725879e-205,-1.187059604818343e-204,-4.1816194982480564e-204,-7.778390980583175e-203,1.5058320989288564e-202,5.782104673291267e-202,1.3934504384399645e-200,-1.9110763684814888e-200,-7.966888870365663e-200,-2.782254568176531e-198,2.4237590025205378e-198,1.0940200120324885e-197,7.288303342718714e-196,-3.0691991975699945e-196,-1.497605957139244e-195,-3.554752256581815e-192,3.877708781110762e-194,2.0442960833288521e-193,-1.3291128910208012e-191,-4.885263406846251e-192,-2.7838573802534393e-191,8.511358059713116e-190,6.134188004224469e-190,3.783966333279612e-189,-7.28241468438634e-188,-7.673836639664125e-188,-5.1375900204092e-187,6.973782497738859e-186,9.56124507299609e-186,6.974264835451697e-185,-7.076348122196084e-184,-1.1861734453008249e-183,-9.478011812688437e-183,7.42615495517674e-182,1.464944799527305e-181,1.2917117455318216e-180,-7.956500394709099e-180,-1.800781061612006e-179,-1.7696368655786212e-178,8.636380040339214e-178,2.202991622135916e-177,2.445476568719537e-176,-9.449865913561455e-176,-2.68188710340233e-175,-3.4263397489616223e-174,1.0387667519920923e-173,3.248815710516482e-173,4.9068899583450716e-172,-1.144309735546926e-171,-3.9162132292647796e-171,-7.28322083795426e-170,1.2609872104952429e-169,4.697685784075478e-169,1.150863736300878e-167,-1.3880778487737913e-167,-5.6081552110454455e-167,-2.0619740170344791e-165,1.5246803032487332e-165,6.664112519137411e-165,5.237399704831706e-163,-1.6696553404851894e-163,-7.884048181874147e-163,2.7912238994920994e-160,1.821589837541054e-161,9.289163845379027e-161,-4.710091139159297e-159,-1.9787798875690506e-159,-1.0904527473434245e-158,2.846088754682745e-157,2.1392322736393302e-157,1.2760930033667355e-156,-2.161029006088634e-155,-2.3006842760694162e-155,-1.489787000921675e-154,1.80521392917136e-153,2.460641731603809e-153,1.736840068630618e-152,-1.5853065459641284e-151,-2.616434751086319e-151,-2.0247171233598554e-150,1.4329993341408279e-149,2.7652897048498255e-149,2.3644210755774952e-148,-1.3179584616050458e-147,-2.904415139293546e-147,-2.7729662872002488e-146,1.2246654529077252e-145,3.0310983176876295e-145,3.278095446233646e-144,-1.144408651545471e-143,-3.1428082022939096e-143,-3.928021521502476e-142,1.0720000398470278e-141,3.237300027985969e-141,4.813792330905615e-140,-1.0042545601920785e-139,-3.3127162650823377e-139,-6.128664629859317e-138,9.392155261164062e-138,3.367678821525536e-137,8.363331517063275e-136,-8.757268274239039e-136,-3.401367863265234e-135,-1.32238620934838e-133,8.131781767970379e-134,3.4135836901383037e-133,3.3427096187653646e-131,-7.513416374272746e-132,-3.404789738328826e-131,5.2318109785985314e-129,6.9025865260091516e-130,3.376137054591683e-129,-1.3643824313041047e-127,-6.30156429886776e-128,-3.329473643129236e-127,7.417169913764835e-126,5.7138251054866895e-126,3.267346257055205e-125,-4.844267163506203e-124,-5.143516568657621e-124,-3.19300826020857e-123,3.429634830202082e-122,4.595020772989101e-122,3.1104568516792313e-121,-2.533360058326788e-120,-4.072595718360417e-120,-3.024540156373404e-119,1.9165976035969175e-118,3.5800897549245033e-118,2.9412087016112633e-117,-1.4696664608404235e-116,-3.120726129247681e-116,-2.8680597977551765e-115,1.1348583988780456e-114,2.6969553339150504e-114,2.81550036119958e-113,-8.786505605788642e-113,-2.310371907354975e-112,-2.7993248826675467e-111,6.800069095053091e-111,1.9616904984494525e-110,2.8469376182708037e-109,-5.248723390015473e-109,-1.650773987924886e-108,-3.014656128281357e-107,4.033567252053917e-107,1.376704674900718e-106,3.4481416417924408e-105,-3.0819923344325144e-105,-1.1378882806556247e-104,-4.698954694588333e-103,2.3388652153739226e-103,9.321799542078056e-103,1.2599802666745998e-100,-1.7612375527888305e-101,-7.570216337769496e-101,7.005600543240901e-99,1.3150422776776367e-99,6.095809765540968e-99,-1.9825940953056094e-97,-9.729448064148285e-98,-4.868834947008163e-97,9.20896479844182e-96,7.128834039787898e-96,3.859313593646435e-95,-4.962898798342298e-94,-5.170299571603466e-94,-3.0380436976642865e-93,2.8603705733224084e-92,3.7101331639275414e-92,2.3774065545038543e-91,-1.7066431885430767e-90,-2.6331175961725228e-90,-1.8519665334345327e-89,1.036989113512067e-88,1.8475929512353584e-88,1.438877272099772e-87,-6.356128263024029e-87,-1.2813330391349678e-86,-1.1181258860537838e-85,3.906378210961296e-85,8.780451422013893e-85,8.726693698066437e-84,-2.3973387050170598e-83,-5.943807993028969e-83,-6.885612974299635e-82,1.4647678230207493e-81,3.973879649554017e-81,5.5529244164443535e-80,-8.890404346158968e-80,-2.623562664895348e-79,-4.6704111563403055e-78,5.350917933604819e-78,1.7101626651413792e-77,4.2792741525396744e-76,-3.189156565092746e-76,-1.1005646134410267e-75,-4.845866089006215e-74,1.8799858489061437e-74,6.992189211637791e-74,1.68436571665312e-71,-1.0950410246781485e-72,-4.385789699155792e-72,2.733790370482116e-70,6.296876844438747e-71,2.71628325828793e-70,-7.095652637639809e-69,-3.5719329605734964e-69,-1.6614983657061648e-68,2.6335216963813735e-67,1.9973836693156822e-67,1.0041366087817824e-66,-1.1014370884612229e-65,-1.1003273651510507e-65,-5.999470268603561e-65,4.859933011241968e-64,5.9679710485725875e-64,3.546851826356016e-63,-2.199700361435585e-62,-3.1851811513884897e-62,-2.0775103128627424e-61,1.0063934824328733e-60,1.6719208823379374e-60,1.2079139045313878e-59,-4.613428021240374e-59,-8.62688947406943e-59,-6.991172475428158e-58,2.1068251117535515e-57,4.373633644084899e-57,4.0455222784641608e-56,-9.54600384884919e-56,-2.177620852112662e-55,-2.3569427747962757e-54,4.278500617982483e-54,1.06435173315725e-53,1.399279620571961e-52,-1.892399466067966e-52,-5.1047259461203104e-52,-8.662244027897601e-51,8.244205119111406e-51,2.401467885927561e-50,5.891104531057422e-49,-3.5317650481079795e-49,-1.1077577742706341e-48,-5.201504894752087e-47,1.4856737185045072e-47,5.0088796035107605e-47,2.4543358603456696e-43,-6.128973852936303e-46,-2.2194793781572853e-45,9.685414412111026e-44,2.476695010965836e-44,9.635937374690966e-44,-1.9838871066856808e-42,-9.79254302000384e-43,-4.098534179585286e-42,5.293858032449176e-41,3.7843756180551835e-41,1.7079391632577913e-40,-1.547913406527303e-39,-1.427967556339901e-39,-6.974695136599793e-39,4.696339102631557e-38,5.255609994461314e-38,2.7925803819881297e-37,-1.442343822968692e-36,-1.884797175385593e-36,-1.097257352826971e-35,4.423388423005951e-35,6.579536087115991e-35,4.237380769925742e-34,-1.3430856027594178e-33,-2.2333681338620563e-33,-1.6123669286501844e-32,4.0136932557417026e-32,7.363673625671829e-32,6.070341417350373e-31,-1.1753063161266218e-30,-2.3556966425334753e-30,-2.277339438048587e-29,3.360334321605555e-29,7.303655710872043e-29,8.623816454580651e-28,-9.352575614180496e-28,-2.1920198767191015e-27,-3.382861291501576e-26,2.5271658912903497e-26,6.36060954127026e-26,1.4637963906177559e-24,-6.613187914991939e-25,-1.7821644092222266e-24,-8.784191820452554e-23,1.671962484545739e-23,4.8151867051740776e-23,-9.351548277914864e-21,-4.0743403868556156e-22,-1.2528323839501552e-21,3.747713210922599e-20,9.547023950698036e-21,3.134463190330884e-20,-4.801269194027617e-19,-2.145807500144421e-19,-7.529754611988379e-19,7.405765982858303e-18,4.614297243855689e-18,1.7341689275997814e-17,-1.205197519674421e-16,-9.467188116666065e-17,-3.823334666728837e-16,1.9736942378590472e-15,1.847794517930657e-15,8.057629139919885e-15,-3.176358790329746e-14,-3.419835877092127e-14,-1.621174752009691e-13,4.950287350829181e-13,5.980477310879885e-13,3.111015058728449e-12,-7.392014346040976e-12,-9.84309458950673e-12,-5.692610439289876e-11,1.0484072469723386e-10,1.5179768175252612e-10,9.943217894592053e-10,-1.4011359153491299e-9,-2.1824488856823623e-9,-1.6641803889541555e-8,1.750721893137372e-8,2.908340936688993e-8,2.6950455059776027e-7,-2.028560131802823e-7,-3.568010046452482e-7,-4.327804942669424e-6,2.160195911808437e-6,3.997626713632721e-6,7.397499200754782e-5,-2.092576096848763e-5,-4.051130136344611e-5,-0.0018597657188592117,0.00018217537339921874,0.000366931050584885,-0.023523529713407765,-0.001404375678903569,-0.0029262558079613914,0.04882709978402207,0.009408154555202568,0.020149931567742918,-0.17679380858844643,-0.0534277865527404,-0.1166663307937594,0.6404074460440544,0.24844473646104276,0.5467269254661202,-2.0072247527824203,-0.8980862044971847,-1.9539561046501777,4.878793885482199,2.3128479576291805,4.790511473331294,-7.88162874195985,-3.5494925084096987,-6.3168350769542245,5.831212166494657,1.772453850905516,1.1216476257258134,0.9234267991754131,0.8871157602049371,0.9493352227643639,1.1033618352130108,1.370877739826718,1.8012496382870498,2.483415943332259,3.5715564047148405,5.333211782896224,8.238140091244544,13.123619888497574,21.505796217615796,36.17400900603283,62.34069828111632,109.89598235693734,197.88600958543833,363.52229563018955,680.5326929884206,1296.9958005869191,2514.2721575057717,4953.5825924444725,9911.53629415503,20127.307062127744,41455.90376904078,86556.02848364945,183101.34432626006,392246.29817187064,850558.3133107716,1.866145481817764e6,4.1410779363501086e6,9.29071372469531e6,2.1067048251811147e7,4.826556724607287e7,1.1169111872082755e8,2.6098963958351523e8,6.156502767605731e8,1.4656834621906223e9,3.520760920764776e9,8.53146168854503e9,2.085000766514558e10,5.138007975838699e10,1.276448544786957e11,3.196318316871095e11,8.065973890316306e11,2.0509215576902158e12,5.253580036778337e12,1.355524497592884e13,3.5224075168081086e13,9.21699190788961e13,2.4282638300860294e14,6.440233842130081e14,1.7192964819099988e15,4.619443472671949e15,1.249012089997775e16,3.398067210112702e16,9.301171130419837e16,2.561163757159423e17,7.093921725606889e17,1.9762485924120095e18,5.536829317896093e18,1.5599284947548684e19,4.4190931332854874e19,1.2586607241238518e20,3.604081525173044e20,1.0374218639369305e21,3.0016260977250927e21,8.729005434590555e21,2.551217365873386e22,7.493301045503555e22,2.211624088270904e23,6.558916599156119e23,1.9543693924235424e24,5.850686334976575e24,1.7595646564117852e25,5.315881089475517e25,1.61321020814564e26,4.9173042182366236e26,1.5054245348076276e27,4.6287391259705193e27,1.4292718141757418e28,4.431926479183303e28,1.3799825637540872e29,4.314557500130318e29,1.3544386162689752e30,4.268960046293447e30,1.3508412799014446e31,4.291274209279003e31,1.3685119900466135e32,4.3809920610910946e32,1.4077948947346063e33,4.540788907311517e33,1.4700448538659695e34,4.776611552139449e34,1.5576949609279884e35,5.0980195864147685e35,1.6744067502449167e36,5.518803563753132e36,1.8253151357475918e37,6.057933238743228e37,2.0173900950248048e38,6.740923821134393e38,2.2599495030007822e39,7.601753070289913e39,2.565373974894949e40,8.685522984856887e40,2.950097306002222e41,1.0052145101639081e42,3.435978221044128e42,1.1781449894247433e43,4.0522052147845434e43,1.39802959085657e44,4.837952998134798e44,1.67925090273155e45,5.846106493977773e45,2.041285550715801e46,7.148511551440194e46,2.510680288423983e47,8.84342342937994e47,3.123863967191474e48,1.1066139471440092e49,3.9311744230329836e49,1.4004274630551676e50,5.002662465307714e50,1.7919849854616184e51,6.436511858072379e51,2.3181441297887677e52,8.371334744521159e52,3.031128156887444e53,1.1004244987490978e54,4.0054724214573775e54,1.4617600954192625e55,5.348337033824471e55,1.9618839708724105e56,7.21491900170525e56,2.660020697214137e57,9.831639206300934e57,3.642891817426233e58,1.353125936909032e59,5.038416405082198e59,1.880642261467387e60,7.03666911096889e60,2.63918288608507e61,9.922154373357598e61,3.739120216699995e62,1.4123838741803545e63,5.347477434060994e63,2.0293229592137862e64,7.718849591948246e64,2.9427033487642237e65,1.124414060509204e66,4.306112803279511e66,1.6527891184861042e67,6.3579343104223e67,2.45117828900327e68,9.47082570333816e68,3.667319515214014e69,1.4231529576707267e70,5.534652691979889e70,2.1570451934796704e71,8.424659941760089e71,3.297342339182071e72,1.2932680627460147e73,5.083002469152159e73,2.0019526940100485e74,7.901027234505343e74,3.1246700547234493e75,1.238255587500414e76,4.9169600651223496e76,1.956402057723652e77,7.799907259390171e77,3.1159148299799025e78,1.2472155039444695e79,5.0020984743387334e79,2.010077835962817e80,8.093168029886788e80,3.2648594000254446e81,1.3196109356353663e82,5.34389764804692e82,2.1681859151334187e83,8.813663763062142e83,3.589495978584689e84,1.4646134161770216e85,5.987157034117714e85,2.452008677528756e86,1.0060573184302912e87,4.1354094623507865e87,1.7029626727030544e88,7.025518037114276e88,2.903581551625408e89,1.2021759528418253e90,4.98627587061348e90,2.071833228837072e91,8.623816770503222e91,3.5958911106210504e92,1.5020068729137283e93,6.284805108452709e93,2.634282742692048e94,1.1060629496860177e95,4.652007815602594e95,1.9599294735686944e96,8.27134604987033e96,3.4965845509479654e97,1.4806089428073397e98,6.280028596061646e98,2.668113238222093e99,1.1354415313761737e100,4.8399369806139184e100,2.066453405181873e101,8.837289956671872e101,3.7854437756015145e102,1.6241097216315224e103,6.9792940520216325e103,3.0040137527934785e104,1.2950405948659784e105,5.591813511689578e105,2.41828021719749e106,1.0474711232487806e107,4.544183507675853e107,1.9744445050330043e108,8.592236084403253e108,3.744868987726936e109,1.634681695277226e110,7.146492243546689e110,3.1290524986712003e111,1.3721134911517056e112,6.025895179791671e112,2.650358253570006e113,1.1674443975903327e114,5.150075667546943e114,2.275275489000857e115,1.0066872240912007e116,4.4605925704508815e116,1.97936210617205e117,8.796095373198749e117,3.9145685382858845e118,1.7446346765122577e119,7.786621823551856e119,3.480286737863696e120,1.5577570747984138e121,6.982334038127475e121,3.1341165176290984e122,1.4087721130956504e123,6.341250806665016e123,2.858348626042305e124,1.2902060478091571e125,5.831819908877259e125,2.639659179169261e126,1.196431712050708e127,5.430275797685077e127,2.4680127333940403e128,1.12321290542243e129,5.118746422566715e129,2.3358753790150335e130,1.0673772547522588e131,4.883893186679894e131,2.237649409756108e132,1.0265813127892774e133,4.71593554215361e133,2.169268692304077e134,9.991429356307452e134,4.607957897200337e135,2.1279127975134066e136,9.839237607623212e136,4.555428540569263e137,2.11181116901069e138,9.80250071680559e138,4.555882943501406e139,2.1201180263334454e140,9.878642667782403e140,4.608740272484814e141,2.1528456332672737e142,1.0069020838005772e143,4.715234432876256e143,2.210848913101871e144,1.0378930418787497e145,4.878451026699412e145,2.2958588162837134e146,1.0817835275931254e147,5.103470425946921e147,2.4105658490497695e148,1.139983789939527e149,5.397625704291866e149,2.558759185831499e150,1.2144420165198401e151,5.770893328981602e151,2.7455312288619505e152,1.3077508355039954e153,6.236445225835546e153,2.977562797735454e154,1.423294246145542e155,6.811404187080667e155,3.2635108749083713e156,1.565446398926099e157,7.517861980901653e157,3.6145297059536905e158,1.7398381855912544e159,8.384242795664792e159,4.044968010203822e160,1.9537137525193862e161,9.447126364562083e161,4.573301423554996e162,2.2164075070447032e163,1.075368885533529e164,5.223381940963026e164,2.539983921734467e165,1.2364980475368792e166,6.026117719061113e166,2.9400988524891653e167,1.436034404918899e168,7.021740973832351e168,3.437164184226607e169,1.6843399132899545e170,8.262884411267258e170,4.057930318419161e171,1.9950186872659177e172,9.818775730246538e172,4.8376475691731635e173,2.3860314180464948e174,1.1780987039271963e175,5.823034164622061e175,2.881228474870887e176,1.4271358884135376e177,7.076374426655318e177,3.5124707427480625e178,1.7452982633074276e179,8.681209395658523e179,4.322580141595531e180,2.154550829222353e181,1.0750283868747008e182,5.369467967950497e182,2.684660427817604e183,1.3436708722564649e184,6.731944743696817e184,3.376221656045714e185,1.694973087637221e186,8.517944303126845e186,4.284948556179227e187,2.1577143880309373e188,1.0876233958633914e189,5.487797531597675e189,2.7717324009611554e190,1.4013187372018448e191,7.09175478676187e191,3.5925298856866007e192,1.821695227091982e193,9.246523208367835e193,4.697939165570857e194,2.389248212680435e195,1.2162949110282099e196,6.197821120740016e196,3.16126213854506e197,1.6139947201092047e198,8.248265781564954e198,4.219305166973688e199,2.160408141875779e200,1.1072509972704594e201,5.68029117179991e201,2.9168104184727104e202,1.499193679089154e203,7.712906625495843e203,3.9718062817983747e204,2.047229659280739e205,1.0562153324131863e206,5.454370657556086e206,2.819304545960926e207,1.4586241265179083e208,7.55349159110206e208,3.915200148783325e209,2.0312386848274856e210,1.0547939704337874e211,5.482427909294707e211,2.85217232063906e212,1.4851662328156853e213,7.740513745298792e213,4.0379307422622245e214,2.108342951133213e215,1.1018333530740247e216,5.763440351610192e216,3.017434758685767e217,1.5811866745596913e218,8.293093312318467e218,4.353485393802315e219,2.2874086232060636e220,1.2029166422603673e221,6.331570245046247e221,3.335572980780946e222,1.7587803171864766e223,9.281836050114486e223,4.902711706798969e224,2.5918986278625282e225,1.3714434028766153e226,7.262982960163224e226,3.8497139291676e227,2.0422891694224533e228,1.0843772618900855e229,5.762580661166326e229,3.0649693703864494e230,1.631574204597715e231,8.692769445930831e231,4.635314485580186e232,2.47382536510874e233,1.3213770058647047e234,7.064011739356308e234,3.77956772289142e235,2.0239417021227102e236,1.0847211671904852e237,5.818373345763043e237,3.123543210892627e238,1.678243712039733e239,9.0245017643187e239,4.856811601508528e240,2.616002921075256e241,1.4102083019124004e242,7.608263052075697e242,4.108129467006156e243,2.2200267980595414e244,1.2006791034547101e245,6.499044035417996e245,3.520667040594475e246,1.908765324244226e247,1.0356938319646272e248,5.624199033258932e248,3.0566068569926948e249,1.6625214992227097e250,9.049882116615464e250,4.930213449741808e251,2.688034615758029e252,1.4667277811758835e253,8.009564430495456e253,4.37735933754999e254,2.3941900421606133e255,1.3105312085795264e256,7.179226428067843e256,3.9359419386552926e257,2.159532186558495e258,1.1857956841653292e259,6.516260821324504e259,3.5836425964495023e260,1.9723672793315063e261,1.0863930801579475e262,5.9885479671591226e262,3.303624331051728e263,1.8238688493383373e264,1.0076962219308059e265,5.571831583250964e265,3.0831731280867376e266,1.7073740748151603e267,9.46214401789977e267,5.247830395026179e268,2.9127194834713784e269,1.6178767576580238e270,8.993306412336206e270,5.002876733214417e271,2.785133084315428e272,1.551661999776744e273,8.651125374343582e273,4.826937631540426e274,2.6952165330682747e275,1.5060445766354067e276,8.421772012632872e276,4.712918463593597e277,2.6393468682231876e278,1.479184746515137e279,8.295968369224977e279,4.656180287474308e280,2.6152296915477904e281,1.4699635274565894e282,8.268358033672469e282,4.65422429629532e283,2.6217422514617174e284,1.4779054765514027e285,8.337131273239216e285,4.706513130830841e286,2.658850406520796e287,1.503141514592792e288,8.503868163392297e288,4.8144114080795705e289,2.7275910924515268e290,1.5464079181271777e291,8.773582416899671e291,4.981238188822435e292,2.830117543731086e293,1.6090806731922014e294,9.154965486089249e294,5.21243343707357e295,2.9698097123632627e296,1.6932472931130934e297,9.660846910150978e297,5.515849830404352e298,3.151457641731085e299,1.8018212572424988e300,1.0308904535961347e301,5.902191543208586e301,3.381531549964242e302,1.938707748135469e303,1.1122679002484526e304,6.385633918190735e304,3.668559683164094e305,2.109033721105841e306,1.213297293421517e307],"x":[-171.1,-170.75948103792416,-170.4189620758483,-170.07844311377247,-169.7379241516966,-169.39740518962077,-169.0568862275449,-168.71636726546907,-168.3758483033932,-168.03532934131738,-167.69481037924152,-167.35429141716565,-167.01377245508982,-166.67325349301396,-166.33273453093813,-165.99221556886226,-165.65169660678643,-165.31117764471057,-164.97065868263473,-164.63013972055887,-164.28962075848304,-163.94910179640718,-163.60858283433134,-163.26806387225548,-162.92754491017965,-162.5870259481038,-162.24650698602795,-161.9059880239521,-161.56546906187626,-161.2249500998004,-160.88443113772456,-160.5439121756487,-160.20339321357287,-159.862874251497,-159.52235528942117,-159.1818363273453,-158.84131736526948,-158.5007984031936,-158.16027944111775,-157.81976047904192,-157.47924151696606,-157.13872255489022,-156.79820359281436,-156.45768463073853,-156.11716566866266,-155.77664670658683,-155.43612774451097,-155.09560878243514,-154.75508982035927,-154.41457085828344,-154.07405189620758,-153.73353293413174,-153.39301397205588,-153.05249500998005,-152.7119760479042,-152.37145708582835,-152.0309381237525,-151.69041916167666,-151.3499001996008,-151.00938123752496,-150.6688622754491,-150.32834331337327,-149.9878243512974,-149.64730538922154,-149.3067864271457,-148.96626746506985,-148.625748502994,-148.28522954091815,-147.94471057884232,-147.60419161676646,-147.26367265469062,-146.92315369261476,-146.58263473053893,-146.24211576846307,-145.90159680638723,-145.56107784431137,-145.22055888223554,-144.88003992015967,-144.53952095808384,-144.19900199600798,-143.85848303393215,-143.51796407185628,-143.17744510978045,-142.8369261477046,-142.49640718562875,-142.1558882235529,-141.81536926147706,-141.4748502994012,-141.13433133732534,-140.7938123752495,-140.45329341317364,-140.1127744510978,-139.77225548902194,-139.4317365269461,-139.09121756487025,-138.75069860279442,-138.41017964071855,-138.06966067864272,-137.72914171656686,-137.38862275449102,-137.04810379241516,-136.70758483033933,-136.36706586826347,-136.02654690618763,-135.68602794411177,-135.34550898203594,-135.00499001996008,-134.66447105788424,-134.32395209580838,-133.98343313373255,-133.64291417165668,-133.30239520958085,-132.961876247505,-132.62135728542916,-132.2808383233533,-131.94031936127743,-131.5998003992016,-131.25928143712574,-130.9187624750499,-130.57824351297404,-130.2377245508982,-129.89720558882235,-129.5566866267465,-129.21616766467065,-128.87564870259482,-128.53512974051895,-128.19461077844312,-127.85409181636726,-127.51357285429141,-127.17305389221556,-126.83253493013972,-126.49201596806387,-126.15149700598802,-125.81097804391217,-125.47045908183632,-125.12994011976048,-124.78942115768463,-124.44890219560878,-124.10838323353293,-123.76786427145709,-123.42734530938124,-123.08682634730539,-122.74630738522954,-122.4057884231537,-122.06526946107785,-121.724750499002,-121.38423153692615,-121.0437125748503,-120.70319361277446,-120.36267465069861,-120.02215568862276,-119.68163672654691,-119.34111776447106,-119.0005988023952,-118.66007984031936,-118.31956087824351,-117.97904191616766,-117.63852295409181,-117.29800399201596,-116.95748502994012,-116.61696606786427,-116.27644710578842,-115.93592814371257,-115.59540918163673,-115.25489021956088,-114.91437125748503,-114.57385229540918,-114.23333333333333,-113.89281437125749,-113.55229540918164,-113.21177644710579,-112.87125748502994,-112.5307385229541,-112.19021956087825,-111.8497005988024,-111.50918163672655,-111.1686626746507,-110.82814371257486,-110.48762475049901,-110.14710578842315,-109.8065868263473,-109.46606786427145,-109.1255489021956,-108.78502994011976,-108.44451097804391,-108.10399201596806,-107.76347305389221,-107.42295409181637,-107.08243512974052,-106.74191616766467,-106.40139720558882,-106.06087824351297,-105.72035928143713,-105.37984031936128,-105.03932135728543,-104.69880239520958,-104.35828343313374,-104.01776447105789,-103.67724550898204,-103.33672654690619,-102.99620758483034,-102.6556886227545,-102.31516966067865,-101.9746506986028,-101.63413173652695,-101.29361277445109,-100.95309381237524,-100.6125748502994,-100.27205588822355,-99.9315369261477,-99.59101796407185,-99.250499001996,-98.90998003992016,-98.56946107784431,-98.22894211576846,-97.88842315369261,-97.54790419161677,-97.20738522954092,-96.86686626746507,-96.52634730538922,-96.18582834331337,-95.84530938123753,-95.50479041916168,-95.16427145708583,-94.82375249500998,-94.48323353293414,-94.14271457085829,-93.80219560878244,-93.4616766467066,-93.12115768463075,-92.78063872255488,-92.44011976047904,-92.09960079840319,-91.75908183632734,-91.41856287425149,-91.07804391217564,-90.7375249500998,-90.39700598802395,-90.0564870259481,-89.71596806387225,-89.3754491017964,-89.03493013972056,-88.69441117764471,-88.35389221556886,-88.01337325349301,-87.67285429141717,-87.33233532934132,-86.99181636726547,-86.65129740518962,-86.31077844311378,-85.97025948103793,-85.62974051896208,-85.28922155688623,-84.94870259481038,-84.60818363273454,-84.26766467065869,-83.92714570858283,-83.58662674650698,-83.24610778443113,-82.90558882235528,-82.56506986027944,-82.22455089820359,-81.88403193612774,-81.5435129740519,-81.20299401197605,-80.8624750499002,-80.52195608782435,-80.1814371257485,-79.84091816367265,-79.5003992015968,-79.15988023952096,-78.81936127744511,-78.47884231536926,-78.13832335329342,-77.79780439121757,-77.45728542914172,-77.11676646706587,-76.77624750499002,-76.43572854291418,-76.09520958083833,-75.75469061876248,-75.41417165668663,-75.07365269461077,-74.73313373253492,-74.39261477045908,-74.05209580838323,-73.71157684630738,-73.37105788423153,-73.03053892215569,-72.69001996007984,-72.34950099800399,-72.00898203592814,-71.6684630738523,-71.32794411177645,-70.9874251497006,-70.64690618762475,-70.3063872255489,-69.96586826347306,-69.62534930139721,-69.28483033932136,-68.94431137724551,-68.60379241516966,-68.26327345309382,-67.92275449101797,-67.58223552894212,-67.24171656686627,-66.90119760479043,-66.56067864271458,-66.22015968063872,-65.87964071856287,-65.53912175648702,-65.19860279441117,-64.85808383233532,-64.51756487025948,-64.17704590818363,-63.83652694610778,-63.496007984031934,-63.155489021956086,-62.81497005988024,-62.47445109780439,-62.13393213572854,-61.793413173652695,-61.45289421157685,-61.112375249501,-60.77185628742515,-60.431337325349304,-60.090818363273456,-59.7502994011976,-59.409780439121754,-59.069261477045906,-58.72874251497006,-58.38822355289421,-58.04770459081836,-57.707185628742515,-57.36666666666667,-57.02614770459082,-56.68562874251497,-56.345109780439124,-56.004590818363276,-55.66407185628743,-55.32355289421157,-54.983033932135726,-54.64251497005988,-54.30199600798403,-53.96147704590818,-53.620958083832335,-53.28043912175649,-52.93992015968064,-52.59940119760479,-52.258882235528944,-51.918363273453096,-51.57784431137725,-51.2373253493014,-50.896806387225546,-50.5562874251497,-50.21576846307385,-49.875249500998,-49.534730538922155,-49.19421157684631,-48.85369261477046,-48.51317365269461,-48.172654690618764,-47.832135728542916,-47.49161676646707,-47.15109780439122,-46.81057884231537,-46.47005988023952,-46.12954091816367,-45.78902195608782,-45.448502994011974,-45.10798403193613,-44.76746506986028,-44.42694610778443,-44.08642714570858,-43.745908183632736,-43.40538922155689,-43.06487025948104,-42.72435129740519,-42.383832335329345,-42.04331337325349,-41.70279441117764,-41.362275449101794,-41.02175648702595,-40.6812375249501,-40.34071856287425,-40.0001996007984,-39.659680638722556,-39.31916167664671,-38.97864271457086,-38.63812375249501,-38.297604790419165,-37.95708582834332,-37.61656686626746,-37.276047904191614,-36.935528942115766,-36.59500998003992,-36.25449101796407,-35.91397205588822,-35.573453093812375,-35.23293413173653,-34.89241516966068,-34.55189620758483,-34.211377245508984,-33.87085828343314,-33.53033932135729,-33.189820359281434,-32.849301397205586,-32.50878243512974,-32.16826347305389,-31.827744510978043,-31.487225548902195,-31.146706586826348,-30.8061876247505,-30.465668662674652,-30.1251497005988,-29.784630738522953,-29.444111776447105,-29.103592814371257,-28.76307385229541,-28.422554890219562,-28.082035928143714,-27.741516966067863,-27.400998003992015,-27.060479041916167,-26.71996007984032,-26.379441117764472,-26.038922155688624,-25.698403193612773,-25.357884231536925,-25.017365269461077,-24.67684630738523,-24.336327345309382,-23.995808383233534,-23.655289421157686,-23.314770459081835,-22.974251497005987,-22.63373253493014,-22.29321357285429,-21.952694610778444,-21.612175648702596,-21.271656686626745,-20.931137724550897,-20.59061876247505,-20.2500998003992,-19.909580838323354,-19.569061876247506,-19.22854291417166,-18.888023952095807,-18.54750499001996,-18.20698602794411,-17.866467065868264,-17.525948103792416,-17.18542914171657,-16.844910179640717,-16.50439121756487,-16.16387225548902,-15.823353293413174,-15.482834331337326,-15.142315369261476,-14.801796407185629,-14.461277445109781,-14.120758483033931,-13.780239520958084,-13.439720558882236,-13.099201596806386,-12.758682634730539,-12.418163672654691,-12.077644710578843,-11.737125748502994,-11.396606786427146,-11.056087824351298,-10.715568862275449,-10.3750499001996,-10.034530938123753,-9.694011976047904,-9.353493013972056,-9.012974051896208,-8.672455089820359,-8.33193612774451,-7.991417165668663,-7.650898203592814,-7.310379241516966,-6.969860279441118,-6.629341317365269,-6.288822355289422,-5.948303393213573,-5.607784431137724,-5.2672654690618765,-4.926746506986028,-4.586227544910179,-4.2457085828343315,-3.905189620758483,-3.5646706586826347,-3.2241516966067865,-2.8836327345309383,-2.5431137724550896,-2.2025948103792414,-1.8620758483033932,-1.5215568862275448,-1.1810379241516966,-0.8405189620758483,-0.5,-0.1594810379241517,0.18103792415169662,0.5215568862275449,0.8620758483033932,1.2025948103792414,1.5431137724550898,1.883632734530938,2.2241516966067865,2.5646706586826347,2.905189620758483,3.2457085828343315,3.5862275449101797,3.926746506986028,4.2672654690618765,4.607784431137724,4.948303393213573,5.288822355289422,5.629341317365269,5.969860279441118,6.310379241516966,6.650898203592814,6.991417165668663,7.331936127744511,7.672455089820359,8.012974051896208,8.353493013972056,8.694011976047904,9.034530938123753,9.3750499001996,9.715568862275449,10.056087824351298,10.396606786427146,10.737125748502994,11.077644710578843,11.418163672654691,11.758682634730539,12.099201596806386,12.439720558882236,12.780239520958084,13.120758483033931,13.461277445109781,13.801796407185629,14.142315369261476,14.482834331337326,14.823353293413174,15.163872255489022,15.504391217564871,15.844910179640719,16.18542914171657,16.525948103792416,16.866467065868264,17.20698602794411,17.54750499001996,17.888023952095807,18.22854291417166,18.569061876247506,18.909580838323354,19.2500998003992,19.59061876247505,19.931137724550897,20.271656686626745,20.612175648702596,20.952694610778444,21.29321357285429,21.63373253493014,21.974251497005987,22.314770459081835,22.655289421157686,22.995808383233534,23.336327345309382,23.67684630738523,24.017365269461077,24.357884231536925,24.698403193612773,25.038922155688624,25.379441117764472,25.71996007984032,26.060479041916167,26.400998003992015,26.741516966067863,27.082035928143714,27.422554890219562,27.76307385229541,28.103592814371257,28.444111776447105,28.784630738522953,29.1251497005988,29.465668662674652,29.8061876247505,30.146706586826348,30.487225548902195,30.827744510978043,31.16826347305389,31.508782435129742,31.84930139720559,32.189820359281434,32.53033932135729,32.87085828343314,33.211377245508984,33.55189620758483,33.89241516966068,34.23293413173653,34.573453093812375,34.91397205588822,35.25449101796407,35.59500998003992,35.935528942115766,36.276047904191614,36.61656686626746,36.95708582834332,37.297604790419165,37.63812375249501,37.97864271457086,38.31916167664671,38.659680638722556,39.0001996007984,39.34071856287425,39.6812375249501,40.02175648702595,40.362275449101794,40.70279441117764,41.04331337325349,41.383832335329345,41.72435129740519,42.06487025948104,42.40538922155689,42.745908183632736,43.08642714570858,43.42694610778443,43.76746506986028,44.10798403193613,44.448502994011974,44.78902195608782,45.12954091816367,45.47005988023952,45.81057884231537,46.15109780439122,46.49161676646707,46.832135728542916,47.172654690618764,47.51317365269461,47.85369261477046,48.19421157684631,48.534730538922155,48.875249500998,49.21576846307385,49.5562874251497,49.896806387225546,50.2373253493014,50.57784431137725,50.918363273453096,51.258882235528944,51.59940119760479,51.93992015968064,52.28043912175649,52.620958083832335,52.96147704590818,53.30199600798403,53.64251497005988,53.983033932135726,54.32355289421157,54.66407185628743,55.004590818363276,55.345109780439124,55.68562874251497,56.02614770459082,56.36666666666667,56.707185628742515,57.04770459081836,57.38822355289421,57.72874251497006,58.069261477045906,58.409780439121754,58.7502994011976,59.090818363273456,59.431337325349304,59.77185628742515,60.112375249501,60.45289421157685,60.793413173652695,61.13393213572854,61.47445109780439,61.81497005988024,62.155489021956086,62.496007984031934,62.83652694610778,63.17704590818363,63.517564870259484,63.85808383233533,64.19860279441117,64.53912175648702,64.87964071856287,65.22015968063872,65.56067864271458,65.90119760479043,66.24171656686627,66.58223552894212,66.92275449101797,67.26327345309382,67.60379241516966,67.94431137724551,68.28483033932136,68.62534930139721,68.96586826347306,69.3063872255489,69.64690618762475,69.9874251497006,70.32794411177645,70.6684630738523,71.00898203592814,71.34950099800399,71.69001996007984,72.03053892215569,72.37105788423153,72.71157684630738,73.05209580838323,73.39261477045908,73.73313373253492,74.07365269461077,74.41417165668663,74.75469061876248,75.09520958083833,75.43572854291418,75.77624750499002,76.11676646706587,76.45728542914172,76.79780439121757,77.13832335329342,77.47884231536926,77.81936127744511,78.15988023952096,78.5003992015968,78.84091816367265,79.1814371257485,79.52195608782435,79.8624750499002,80.20299401197605,80.5435129740519,80.88403193612774,81.22455089820359,81.56506986027944,81.90558882235528,82.24610778443113,82.58662674650698,82.92714570858283,83.26766467065869,83.60818363273454,83.94870259481038,84.28922155688623,84.62974051896208,84.97025948103793,85.31077844311378,85.65129740518962,85.99181636726547,86.33233532934132,86.67285429141717,87.01337325349301,87.35389221556886,87.69441117764471,88.03493013972056,88.3754491017964,88.71596806387225,89.0564870259481,89.39700598802395,89.7375249500998,90.07804391217564,90.41856287425149,90.75908183632734,91.09960079840319,91.44011976047904,91.78063872255488,92.12115768463075,92.4616766467066,92.80219560878244,93.14271457085829,93.48323353293414,93.82375249500998,94.16427145708583,94.50479041916168,94.84530938123753,95.18582834331337,95.52634730538922,95.86686626746507,96.20738522954092,96.54790419161677,96.88842315369261,97.22894211576846,97.56946107784431,97.90998003992016,98.250499001996,98.59101796407185,98.9315369261477,99.27205588822355,99.6125748502994,99.95309381237524,100.29361277445109,100.63413173652695,100.9746506986028,101.31516966067865,101.6556886227545,101.99620758483034,102.33672654690619,102.67724550898204,103.01776447105789,103.35828343313374,103.69880239520958,104.03932135728543,104.37984031936128,104.72035928143713,105.06087824351297,105.40139720558882,105.74191616766467,106.08243512974052,106.42295409181637,106.76347305389221,107.10399201596806,107.44451097804391,107.78502994011976,108.1255489021956,108.46606786427145,108.8065868263473,109.14710578842315,109.48762475049901,109.82814371257486,110.1686626746507,110.50918163672655,110.8497005988024,111.19021956087825,111.5307385229541,111.87125748502994,112.21177644710579,112.55229540918164,112.89281437125749,113.23333333333333,113.57385229540918,113.91437125748503,114.25489021956088,114.59540918163673,114.93592814371257,115.27644710578842,115.61696606786427,115.95748502994012,116.29800399201596,116.63852295409181,116.97904191616766,117.31956087824351,117.66007984031936,118.0005988023952,118.34111776447106,118.68163672654691,119.02215568862276,119.36267465069861,119.70319361277446,120.0437125748503,120.38423153692615,120.724750499002,121.06526946107785,121.4057884231537,121.74630738522954,122.08682634730539,122.42734530938124,122.76786427145709,123.10838323353293,123.44890219560878,123.78942115768463,124.12994011976048,124.47045908183632,124.81097804391217,125.15149700598802,125.49201596806387,125.83253493013972,126.17305389221556,126.51357285429141,126.85409181636726,127.19461077844312,127.53512974051897,127.87564870259482,128.21616766467065,128.5566866267465,128.89720558882235,129.2377245508982,129.57824351297404,129.9187624750499,130.25928143712574,130.5998003992016,130.94031936127743,131.2808383233533,131.62135728542916,131.961876247505,132.30239520958085,132.64291417165668,132.98343313373255,133.32395209580838,133.66447105788424,134.00499001996008,134.34550898203594,134.68602794411177,135.02654690618763,135.36706586826347,135.70758483033933,136.04810379241516,136.38862275449102,136.72914171656686,137.06966067864272,137.41017964071855,137.75069860279442,138.09121756487025,138.4317365269461,138.77225548902194,139.1127744510978,139.45329341317364,139.7938123752495,140.13433133732534,140.4748502994012,140.81536926147706,141.1558882235529,141.49640718562875,141.8369261477046,142.17744510978045,142.51796407185628,142.85848303393215,143.19900199600798,143.53952095808384,143.88003992015967,144.22055888223554,144.56107784431137,144.90159680638723,145.24211576846307,145.58263473053893,145.92315369261476,146.26367265469062,146.60419161676646,146.94471057884232,147.28522954091815,147.625748502994,147.96626746506985,148.3067864271457,148.64730538922154,148.9878243512974,149.32834331337327,149.6688622754491,150.00938123752496,150.3499001996008,150.69041916167666,151.0309381237525,151.37145708582835,151.7119760479042,152.05249500998005,152.39301397205588,152.73353293413174,153.07405189620758,153.41457085828344,153.75508982035927,154.09560878243514,154.43612774451097,154.77664670658683,155.11716566866266,155.45768463073853,155.79820359281436,156.13872255489022,156.47924151696606,156.81976047904192,157.16027944111775,157.5007984031936,157.84131736526948,158.1818363273453,158.52235528942117,158.862874251497,159.20339321357287,159.5439121756487,159.88443113772456,160.2249500998004,160.56546906187626,160.9059880239521,161.24650698602795,161.5870259481038,161.92754491017965,162.26806387225548,162.60858283433134,162.94910179640718,163.28962075848304,163.63013972055887,163.97065868263473,164.31117764471057,164.65169660678643,164.99221556886226,165.33273453093813,165.67325349301396,166.01377245508982,166.35429141716565,166.69481037924152,167.03532934131738,167.3758483033932,167.71636726546907,168.0568862275449,168.39740518962077,168.7379241516966,169.07844311377247,169.4189620758483,169.75948103792416,170.1]} diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/integers.json b/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/integers.json new file mode 100644 index 000000000000..b9adae4fb316 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/integers.json @@ -0,0 +1 @@ +{"expected":[1.0,1.0,2.0,6.0,24.0,120.0,720.0,5040.0,40320.0,362880.0,3.6288e6,3.99168e7,4.790016e8,6.2270208e9,8.71782912e10,1.307674368e12,2.0922789888e13,3.55687428096e14,6.402373705728e15,1.21645100408832e17,2.43290200817664e18,5.109094217170944e19,1.1240007277776077e21,2.585201673888498e22,6.204484017332394e23,1.5511210043330986e25,4.0329146112660565e26,1.0888869450418352e28,3.0488834461171387e29,8.841761993739702e30,2.6525285981219107e32,8.222838654177922e33,2.631308369336935e35,8.683317618811887e36,2.9523279903960416e38,1.0333147966386145e40,3.719933267899012e41,1.3763753091226346e43,5.230226174666011e44,2.0397882081197444e46,8.159152832478977e47,3.345252661316381e49,1.40500611775288e51,6.041526306337383e52,2.658271574788449e54,1.1962222086548019e56,5.502622159812089e57,2.5862324151116818e59,1.2413915592536073e61,6.082818640342676e62,3.0414093201713376e64,1.5511187532873822e66,8.065817517094388e67,4.2748832840600255e69,2.308436973392414e71,1.2696403353658276e73,7.109985878048635e74,4.0526919504877214e76,2.3505613312828785e78,1.3868311854568984e80,8.32098711274139e81,5.075802138772248e83,3.146997326038794e85,1.98260831540444e87,1.2688693218588417e89,8.247650592082472e90,5.443449390774431e92,3.647111091818868e94,2.4800355424368305e96,1.7112245242814133e98,1.1978571669969892e100,8.504785885678623e101,6.1234458376886085e103,4.4701154615126844e105,3.307885441519386e107,2.48091408113954e109,1.88549470166605e111,1.4518309202828587e113,1.13242811782063e115,8.946182130782976e116,7.156945704626381e118,5.797126020747368e120,4.753643337012842e122,3.945523969720659e124,3.314240134565353e126,2.8171041143805505e128,2.4227095383672734e130,2.107757298379528e132,1.8548264225739844e134,1.650795516090846e136,1.4857159644817615e138,1.352001527678403e140,1.2438414054641308e142,1.1567725070816416e144,1.087366156656743e146,1.032997848823906e148,9.916779348709496e149,9.619275968248212e151,9.426890448883248e153,9.332621544394415e155,9.332621544394415e157,9.42594775983836e159,9.614466715035127e161,9.90290071648618e163,1.0299016745145628e166,1.081396758240291e168,1.1462805637347084e170,1.226520203196138e172,1.324641819451829e174,1.4438595832024937e176,1.5882455415227428e178,1.7629525510902446e180,1.974506857221074e182,2.2311927486598138e184,2.5435597334721877e186,2.9250936934930154e188,3.393108684451898e190,3.969937160808721e192,4.684525849754291e194,5.574585761207606e196,6.689502913449127e198,8.094298525273444e200,9.875044200833601e202,1.214630436702533e205,1.506141741511141e207,1.882677176888926e209,2.372173242880047e211,3.0126600184576594e213,3.856204823625804e215,4.974504222477287e217,6.466855489220474e219,8.47158069087882e221,1.1182486511960043e224,1.4872707060906857e226,1.9929427461615188e228,2.6904727073180504e230,3.6590428819525483e232,5.012888748274991e234,6.917786472619489e236,9.615723196941089e238,1.3462012475717526e241,1.898143759076171e243,2.695364137888163e245,3.854370717180073e247,5.5502938327393044e249,8.047926057471992e251,1.1749972043909107e254,1.727245890454639e256,2.556323917872866e258,3.80892263763057e260,5.713383956445855e262,8.62720977423324e264,1.3113358856834527e267,2.0063439050956823e269,3.0897696138473508e271,4.789142901463394e273,7.471062926282894e275,1.1729568794264145e278,1.853271869493735e280,2.9467022724950384e282,4.7147236359920616e284,7.590705053947219e286,1.2296942187394494e289,2.0044015765453026e291,3.287218585534296e293,5.423910666131589e295,9.003691705778438e297,1.503616514864999e300,2.5260757449731984e302,4.269068009004705e304,7.257415615307999e306],"x":[0.0,1.0,2.0,3.0,4.0,5.0,6.0,7.0,8.0,9.0,10.0,11.0,12.0,13.0,14.0,15.0,16.0,17.0,18.0,19.0,20.0,21.0,22.0,23.0,24.0,25.0,26.0,27.0,28.0,29.0,30.0,31.0,32.0,33.0,34.0,35.0,36.0,37.0,38.0,39.0,40.0,41.0,42.0,43.0,44.0,45.0,46.0,47.0,48.0,49.0,50.0,51.0,52.0,53.0,54.0,55.0,56.0,57.0,58.0,59.0,60.0,61.0,62.0,63.0,64.0,65.0,66.0,67.0,68.0,69.0,70.0,71.0,72.0,73.0,74.0,75.0,76.0,77.0,78.0,79.0,80.0,81.0,82.0,83.0,84.0,85.0,86.0,87.0,88.0,89.0,90.0,91.0,92.0,93.0,94.0,95.0,96.0,97.0,98.0,99.0,100.0,101.0,102.0,103.0,104.0,105.0,106.0,107.0,108.0,109.0,110.0,111.0,112.0,113.0,114.0,115.0,116.0,117.0,118.0,119.0,120.0,121.0,122.0,123.0,124.0,125.0,126.0,127.0,128.0,129.0,130.0,131.0,132.0,133.0,134.0,135.0,136.0,137.0,138.0,139.0,140.0,141.0,142.0,143.0,144.0,145.0,146.0,147.0,148.0,149.0,150.0,151.0,152.0,153.0,154.0,155.0,156.0,157.0,158.0,159.0,160.0,161.0,162.0,163.0,164.0,165.0,166.0,167.0,168.0,169.0,170.0]} diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/runner.jl b/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/runner.jl new file mode 100644 index 000000000000..01a81ae81627 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/test/fixtures/julia/runner.jl @@ -0,0 +1,78 @@ +#!/usr/bin/env julia +# +# @license Apache-2.0 +# +# Copyright (c) 2026 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. + +import JSON +using SpecialFunctions + +""" + gen( x, name ) + +Generate fixture data and write to file. + +# Arguments + +* `x`: domain +* `name::AbstractString`: output filename + +# Examples + +``` julia +julia> x = range( 0, stop = 34, length = 35 ); +julia> gen( x, \"data.json\" ); +``` +""" +function gen( x, name ) + y = Array{Float64}( undef, length( x ) ); + for i in eachindex(x) + if isinteger( x[i] ) && x[i] >= 0 + # Use big float for accuracy before rounding to float64/float32: + y[i] = Float64( factorial( big( x[i] ) ) ); + else + y[i] = Float64( gamma( Float32( x[i] ) + 1.0f0 ) ); + end + end + + # Store data to be written to file as a collection: + data = Dict([ + ("x", x), + ("expected", y) + ]); + + # Based on the script directory, create an output filepath: + filepath = joinpath( dir, name ); + + # Write the data to the output filepath as JSON: + outfile = open( filepath, "w" ); + write( outfile, JSON.json(data) ); + write( outfile, "\n" ); + close( outfile ); +end + +# Get the filename: +file = @__FILE__; + +# Extract the directory in which this file resides: +dir = dirname( file ); + +# Integer values (up to the limit for float32): +x = range( 0, stop = 34, length = 35 ); +gen( x, "integers.json" ); + +# Decimal values: +x = range( -35.0, stop = 34.0, length = 1001 ); +gen( x, "decimals.json" ); diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/test/test.js b/lib/node_modules/@stdlib/math/base/special/factorialf/test/test.js new file mode 100644 index 000000000000..f0285da9b152 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/test/test.js @@ -0,0 +1,144 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 tape = require( 'tape' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var incrspace = require( '@stdlib/array/base/incrspace' ); +var f32 = require( '@stdlib/number/float64/base/to-float32' ); +var isAlmostSameValue = require( '@stdlib/number/float32/base/assert/is-almost-same-value' ); +var PINF = require( '@stdlib/constants/float32/pinf' ); +var NINF = require( '@stdlib/constants/float32/ninf' ); +var factorialf = require( './../lib' ); + + +// FIXTURES // + +var integers = require( './fixtures/julia/integers.json' ); +var decimals = require( './fixtures/julia/decimals.json' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof factorialf, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'if provided a negative integer, the function returns `NaN`', function test( t ) { + var values; + var v; + var i; + + values = incrspace( -1.0, -100.0, -1.0 ); + for ( i = 0; i < values.length; i++ ) { + v = factorialf( values[ i ] ); + t.strictEqual( isnanf( v ), true, 'returns expected value when provided ' + values[ i ] ); + } + t.end(); +}); + +tape( 'if provided negative infinity, the function returns `NaN`', function test( t ) { + var v = factorialf( NINF ); + t.strictEqual( isnanf( v ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'if provided positive infinity, the function returns `+infinity`', function test( t ) { + var v = factorialf( PINF ); + t.strictEqual( v, PINF, 'returns expected value' ); + t.end(); +}); + +tape( 'if provided `NaN`, the function returns `NaN`', function test( t ) { + var v = factorialf( NaN ); + t.strictEqual( isnanf( v ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'if `x > 34.038...`, the function returns positive infinity', function test( t ) { + var values; + var v; + var i; + + values = incrspace( 35.0, 1000.0, 10.1234 ); + for ( i = 0; i < values.length; i++ ) { + v = factorialf( values[ i ] ); + t.strictEqual( v, PINF, 'returns expected value when provided ' + values[ i ] ); + } + t.end(); +}); + +tape( 'the function evaluates the factorial function (positive integers < 35)', function test( t ) { + var expected; + var x; + var v; + var i; + + x = integers.x; + expected = integers.expected; + for ( i = 0; i < x.length; i++ ) { + if ( x[ i ] > 34 ) { + continue; + } + v = factorialf( x[ i ] ); + t.strictEqual( isAlmostSameValue( v, f32( expected[ i ] ), 1 ), true, 'returns expected value when provided ' + x[ i ] + '. Value: ' + v + '. Expected: ' + expected[ i ] + '.' ); + } + t.end(); +}); + +tape( 'if provided positive integers greater than `34`, the function returns positive infinity', function test( t ) { + var i; + var v; + for ( i = 35; i < 500; i++ ) { + v = factorialf( i ); + t.strictEqual( v, PINF, 'returns expected value when provided ' + i ); + } + t.end(); +}); + +tape( 'the function evaluates the factorial function (decimal values)', function test( t ) { + var expected; + var x; + var v; + var i; + + x = decimals.x; + expected = decimals.expected; + for ( i = 0; i < x.length; i++ ) { + v = factorialf( x[ i ] ); + if ( isnanf( v ) ) { + t.strictEqual( isnanf( v ), isnanf( f32( expected[ i ] ) ), 'returns NaN for x='+x[i] ); + continue; + } + if ( v === PINF ) { + t.strictEqual( v, PINF, 'returns Infinity for x='+x[i] ); + continue; + } + if ( v === 0.0 && f32( expected[ i ] ) === 0.0 ) { + t.strictEqual( v, 0.0, 'returns 0.0 for x=' + x[ i ] ); + continue; + } + t.strictEqual( isAlmostSameValue( v, f32( expected[ i ] ), 3000 ), true, 'within tolerance. x: ' + x[ i ] + '. Value: ' + v + '. Expected: ' + expected[ i ] + '.' ); + } + t.end(); +}); diff --git a/lib/node_modules/@stdlib/math/base/special/factorialf/test/test.native.js b/lib/node_modules/@stdlib/math/base/special/factorialf/test/test.native.js new file mode 100644 index 000000000000..b8ada3985901 --- /dev/null +++ b/lib/node_modules/@stdlib/math/base/special/factorialf/test/test.native.js @@ -0,0 +1,155 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var incrspace = require( '@stdlib/array/base/incrspace' ); +var f32 = require( '@stdlib/number/float64/base/to-float32' ); +var isAlmostSameValue = require( '@stdlib/number/float32/base/assert/is-almost-same-value' ); +var PINF = require( '@stdlib/constants/float32/pinf' ); +var NINF = require( '@stdlib/constants/float32/ninf' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// FIXTURES // + +var integers = require( './fixtures/julia/integers.json' ); +var decimals = require( './fixtures/julia/decimals.json' ); + + +// VARIABLES // + +var factorialf = tryRequire( resolve( __dirname, './../lib/native.js' ) ); +var opts = { + 'skip': ( factorialf instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof factorialf, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'if provided a negative integer, the function returns `NaN`', opts, function test( t ) { + var values; + var v; + var i; + + values = incrspace( -1.0, -1000.0, -1.0 ); + for ( i = 0; i < values.length; i++ ) { + v = factorialf( values[ i ] ); + t.strictEqual( isnanf( v ), true, 'returns expected value when provided ' + values[ i ] ); + } + t.end(); +}); + +tape( 'if provided negative infinity, the function returns `NaN`', opts, function test( t ) { + var v = factorialf( NINF ); + t.strictEqual( isnanf( v ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'if provided positive infinity, the function returns `+infinity`', opts, function test( t ) { + var v = factorialf( PINF ); + t.strictEqual( v, PINF, 'returns expected value' ); + t.end(); +}); + +tape( 'if provided `NaN`, the function returns `NaN`', opts, function test( t ) { + var v = factorialf( NaN ); + t.strictEqual( isnanf( v ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'if `x > 34.038...`, the function returns positive infinity', opts, function test( t ) { + var values; + var v; + var i; + + values = incrspace( 35.0, 1000.0, 10.1234 ); + for ( i = 0; i < values.length; i++ ) { + v = factorialf( values[ i ] ); + t.strictEqual( v, PINF, 'returns expected value when provided ' + values[ i ] ); + } + t.end(); +}); + +// (removed negative range test as it is not applicable or different for float32/gamma) + +tape( 'the function evaluates the factorial function (positive integers < 35)', opts, function test( t ) { + var expected; + var x; + var v; + var i; + + x = integers.x; + expected = integers.expected; + for ( i = 0; i < x.length; i++ ) { + if ( x[ i ] > 34 ) { + continue; + } + v = factorialf( x[ i ] ); + t.strictEqual( isAlmostSameValue( v, f32( expected[ i ] ), 1 ), true, 'returns expected value when provided ' + x[ i ] + '. Value: ' + v + '. Expected: ' + expected[ i ] + '.' ); + } + t.end(); +}); + +tape( 'if provided positive integers greater than `34`, the function returns positive infinity', opts, function test( t ) { + var i; + var v; + for ( i = 35; i < 500; i++ ) { + v = factorialf( i ); + t.strictEqual( v, PINF, 'returns expected value when provided ' + i ); + } + t.end(); +}); + +tape( 'the function evaluates the factorial function (decimal values)', opts, function test( t ) { + var expected; + var x; + var v; + var i; + + x = decimals.x; + expected = decimals.expected; + for ( i = 0; i < x.length; i++ ) { + v = factorialf( x[ i ] ); + if ( isnanf( v ) ) { + t.strictEqual( isnanf( v ), isnanf( f32( expected[ i ] ) ), 'returns NaN for x='+x[i] ); + continue; + } + if ( v === PINF ) { + t.strictEqual( v, PINF, 'returns Infinity for x='+x[i] ); + continue; + } + if ( v === 0.0 && f32( expected[ i ] ) === 0.0 ) { + t.strictEqual( v, 0.0, 'returns 0.0 for x=' + x[ i ] ); + continue; + } + t.strictEqual( isAlmostSameValue( v, f32( expected[ i ] ), 3000 ), true, 'within tolerance. x: ' + x[ i ] + '. Value: ' + v + '. Expected: ' + expected[ i ] + '.' ); + } + t.end(); +});