diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/README.md b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/README.md new file mode 100644 index 000000000000..7fda60686cba --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/README.md @@ -0,0 +1,260 @@ + + +# Probability Density Function + +> [Anglit][anglit-distribution] distribution [probability density function][pdf] (PDF). + +
+ +The [probability density function][pdf] (PDF) for an [anglit][anglit-distribution] random variable is + + + +```math +f(x;\mu,\sigma)=\begin{cases} \frac{1}{\sigma} \cos\left(2 \frac{x-\mu}{\sigma}\right) & \text{for } -\frac{\pi}{4} \le \frac{x-\mu}{\sigma} \le \frac{\pi}{4} \\ 0 & \text{otherwise} \end{cases} +``` + + + + + +where `μ` is the location parameter and `σ > 0` is the scale parameter. + +
+ + + +
+ +## Usage + +```javascript +var pdf = require( '@stdlib/stats/base/dists/anglit/pdf' ); +``` + +#### pdf( x, mu, sigma ) + +Evaluates the [probability density function][pdf] (PDF) for an [anglit][anglit-distribution] distribution with location parameter `μ` and scale parameter `σ > 0`. + +```javascript +var y = pdf( 0.5, 0.0, 1.0 ); +// returns ~0.540 + +y = pdf( 2.0, 0.0, 1.0 ); +// returns 0.0 + +y = pdf( 0.0, 0.0, 1.0 ); +// returns 1.0 +``` + +If provided `NaN` as any argument, the function returns `NaN`. + +```javascript +var y = pdf( NaN, 0.0, 1.0 ); +// returns NaN + +y = pdf( 0.0, NaN, 1.0 ); +// returns NaN + +y = pdf( 0.0, 0.0, NaN ); +// returns NaN +``` + +If provided `σ <= 0`, the function returns `NaN`. + +```javascript +var y = pdf( 0.0, 0.0, -1.0 ); +// returns NaN + +y = pdf( 0.0, 0.0, 0.0 ); +// returns NaN +``` + +#### pdf.factory( mu, sigma ) + +Returns a `function` for evaluating the [PDF][pdf] of an [anglit][anglit-distribution] distribution with location parameter `μ` and scale parameter `σ`. + +```javascript +var myPDF = pdf.factory( 0.0, 1.0 ); +var y = myPDF( 0.0 ); +// returns 1.0 + +y = myPDF( 1.0 ); +// returns 0.0 +``` + +
+ + + +
+ +## Examples + + + +```javascript +var uniform = require( '@stdlib/random/array/uniform' ); +var logEachMap = require( '@stdlib/console/log-each-map' ); +var pdf = require( '@stdlib/stats/base/dists/anglit/pdf' ); + +var opts = { + 'dtype': 'float64' +}; +var x = uniform( 10, -5.0, 5.0, opts ); +var mu = uniform( x.length, -5.0, 5.0, opts ); +var sigma = uniform( x.length, 0.5, 5.0, opts ); + +logEachMap( 'x: %lf, mu: %lf, σ: %lf, f(x;mu,σ): %lf', x, mu, sigma, pdf ); +``` + +
+ + + + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/stats/base/dists/anglit/pdf.h" +``` + +#### stdlib_base_dists_anglit_pdf( x, mu, sigma ) + +Evaluates the probability density function (PDF) for Anglit distribution. + +```c +double out = stdlib_base_dists_anglit_pdf( 0.5, 0.0, 1.0 ); +// returns ~0.540 +``` + +The function accepts the following arguments: + +- **x**: `[in] double` input value. +- **mu**: `[in] double` location parameter. +- **sigma**: `[in] double` scale parameter. + +```c +double stdlib_base_dists_anglit_pdf( const double x, const double mu, const double sigma ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/stats/base/dists/anglit/pdf.h" +#include +#include + +static double random_uniform( const double min, const double max ) { + double v = (double)rand() / ( (double)RAND_MAX + 1.0 ); + return min + ( v*(max-min) ); +} + +int main( void ) { + double sigma; + double mu; + double x; + double y; + int i; + + for ( i = 0; i < 25; i++ ) { + x = random_uniform( -10.0, 10.0 ); + mu = random_uniform( -20.0, 0.0 ); + sigma = random_uniform( mu, mu+40.0 ); + y = stdlib_base_dists_anglit_pdf( x, mu, sigma ); + printf( "x: %lf, mu: %lf, σ: %lf, f(x;mu,σ): %lf\n", x, mu, sigma, y ); + } +} +``` + +
+ + + +
+ + + + + +
+ +
+ + + + + +
+ +
+ + + + + +
+ +[pdf]: https://en.wikipedia.org/wiki/Probability_density_function + +[anglit-distribution]: https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.anglit.html + +
+ + diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/benchmark/benchmark.js new file mode 100644 index 000000000000..9af120245a47 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/benchmark/benchmark.js @@ -0,0 +1,89 @@ +/** +* @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 uniform = require( '@stdlib/random/array/uniform' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; +var pdf = require( './../lib' ); + + +// MAIN // + +bench( pkg, function benchmark( b ) { + var sigma; + var len; + var mu; + var x; + var y; + var i; + + len = 100; + x = uniform( len, -10.0, 10.0 ); + mu = uniform( len, -5.0, 5.0 ); + sigma = uniform( len, 0.1, 10.0 ); + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = pdf( x[ i % len ], mu[ i % len ], sigma[ i % len ] ); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); + +bench( format( '%s::factory', pkg ), function benchmark( b ) { + var mypdf; + var sigma; + var len; + var mu; + var x; + var y; + var i; + + mu = 0.0; + sigma = 2.0; + mypdf = pdf.factory( mu, sigma ); + len = 100; + x = uniform( len, -5.0, 5.0 ); + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = mypdf( x[ i % len ] ); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/benchmark/benchmark.native.js new file mode 100644 index 000000000000..095999d63890 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/benchmark/benchmark.native.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 resolve = require( 'path' ).resolve; +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var pdf = tryRequire( resolve( __dirname, './../lib/native.js' ) ); +var opts = { + 'skip': ( pdf instanceof Error ) +}; + + +// MAIN // + +bench( format( '%s::native', pkg ), opts, function benchmark( b ) { + var sigma; + var len; + var mu; + var x; + var y; + var i; + + len = 100; + x = uniform( len, -10.0, 10.0, opts ); + mu = uniform( len, -5.0, 5.0, opts ); + sigma = uniform( len, 0.1, 10.0, opts ); + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + y = pdf( x[ i % len ], mu[ i % len ], sigma[ i % len ] ); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/benchmark/c/Makefile b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/benchmark/c/Makefile new file mode 100644 index 000000000000..979768abbcec --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/benchmark/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 := 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/stats/base/dists/anglit/pdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/benchmark/c/benchmark.c new file mode 100644 index 000000000000..cddc5b19a7eb --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/benchmark/c/benchmark.c @@ -0,0 +1,142 @@ +/** +* @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/stats/base/dists/anglit/pdf.h" +#include +#include +#include +#include +#include + +#define NAME "anglit-pdf" +#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 double random_uniform( const double min, const double max ) { + double v = (double)rand() / ( (double)RAND_MAX + 1.0 ); + return min + ( v*(max-min) ); +} + +/** +* Runs a benchmark. +* +* @return elapsed time in seconds +*/ +static double benchmark( void ) { + double elapsed; + double sigma[ 100 ]; + double mu[ 100 ]; + double x[ 100 ]; + double y; + double t; + int i; + + for ( i = 0; i < 100; i++ ) { + x[ i ] = random_uniform( -10.0, 10.0 ); + mu[ i ] = random_uniform( -5.0, 5.0 ); + sigma[ i ] = random_uniform( 0.1, 10.0 ); // Scale must be positive + } + + t = tic(); + for ( i = 0; i < ITERATIONS; i++ ) { + y = stdlib_base_dists_anglit_pdf( x[ i%100 ], mu[ i%100 ], sigma[ 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::%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/stats/base/dists/anglit/pdf/binding.gyp b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/binding.gyp new file mode 100644 index 000000000000..0d6508a12e99 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/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/stats/base/dists/anglit/pdf/docs/repl.txt b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/docs/repl.txt new file mode 100644 index 000000000000..307ad187f1e9 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/docs/repl.txt @@ -0,0 +1,72 @@ +{{alias}}( x, mu, sigma ) + Evaluates the probability density function (PDF) for anglit + distribution with location parameter `mu` and scale parameter + `sigma` at a value `x`. + + If provided `NaN` as any argument, the function returns `NaN`. + + If provided `sigma <= 0`, the function returns `NaN`. + + Parameters + ---------- + x: number + Input value. + + mu: number + Location parameter. + + sigma: number + Scale parameter. + + Returns + ------- + out: number + Evaluated PDF. + + Examples + -------- + > var y = {{alias}}( 0.0, 0.0, 1.0 ) + 1.0 + > y = {{alias}}( 0.5, 0.0, 1.0 ) + ~0.540 + > y = {{alias}}( 2.0, 0.0, 1.0 ) + 0.0 + > y = {{alias}}( NaN, 0.0, 1.0 ) + NaN + > y = {{alias}}( 0.0, NaN, 1.0 ) + NaN + > y = {{alias}}( 0.0, 0.0, NaN ) + NaN + > y = {{alias}}( 0.0, 0.0, -1.0 ) + NaN + + +{{alias}}.factory( mu, sigma ) + Returns a function for evaluating the probability density function + (PDF) of anglit distribution with location parameter `mu` and scale + parameter `sigma`. + + Parameters + ---------- + mu: number + Location parameter. + + sigma: number + Scale parameter. + + Returns + ------- + pdf: Function + Probability density function (PDF). + + Examples + -------- + > var myPDF = {{alias}}.factory( 0.0, 1.0 ); + > var y = myPDF( 0.0 ) + 1.0 + > y = myPDF( 1.0 ) + 0.0 + + See Also + -------- + diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/docs/types/index.d.ts b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/docs/types/index.d.ts new file mode 100644 index 000000000000..f7eaeb31e081 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/docs/types/index.d.ts @@ -0,0 +1,117 @@ +/** +* @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 probability density function (PDF) for Anglit distribution. +* +* @param x - input value +* @returns evaluated PDF +*/ +type Unary = ( x: number ) => number; + +/** +* Interface for the probability density function (PDF) of Anglit distribution. +*/ +interface PDF { + /** + * Evaluates the probability density function (PDF) for Anglit distribution with location parameter `mu` and scale `sigma` at a value `x`. + * + * ## Notes + * + * - If provided `sigma <= 0`, the function returns `NaN`. + * + * @param x - input value + * @param mu - location parameter + * @param sigma - scale parameter + * @returns evaluated PDF + * + * @example + * var y = pdf( 0.0, 0.0, 1.0 ); + * // returns 1.0 + * + * @example + * var y = pdf( 0.5, 0.0, 1.0 ); + * // returns ~0.540 + * + * @example + * var y = pdf( 2.0, 0.0, 1.0 ); + * // returns 0.0 + * + * @example + * var y = pdf( NaN, 0.0, 1.0 ); + * // returns NaN + * + * @example + * var y = pdf( 0.0, NaN, 1.0 ); + * // returns NaN + * + * @example + * var y = pdf( 0.0, 0.0, NaN ); + * // returns NaN + * + * @example + * var y = pdf( 0.0, 0.0, -1.0 ); + * // returns NaN + */ + ( x: number, mu: number, sigma: number ): number; + + /** + * Returns a function for evaluating the probability density function (PDF) for Anglit distribution with location `mu` and scale `sigma`. + * + * @param mu - location parameter + * @param sigma - scale parameter + * @returns PDF + * + * @example + * var mypdf = pdf.factory( 0.0, 1.0 ); + * var y = mypdf( 0.0 ); + * // returns 1.0 + * + * y = mypdf( 2.0 ); + * // returns 0.0 + */ + factory( mu: number, sigma: number ): Unary; +} + +/** +* Anglit distribution probability density function (PDF). +* +* @param x - input value +* @param mu - location parameter +* @param sigma - scale parameter +* @returns evaluated PDF +* +* @example +* var y = pdf( 0.5, 0.0, 1.0 ); +* // returns ~0.540 +* +* var mypdf = pdf.factory( 0.0, 1.0 ); +* y = mypdf( 0.0 ); +* // returns 1.0 +* +* y = mypdf( 1.0 ); +* // returns 0.0 +*/ +declare var pdf: PDF; + + +// EXPORTS // + +export = pdf; diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/docs/types/test.ts b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/docs/types/test.ts new file mode 100644 index 000000000000..51bb96e2d2f3 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/docs/types/test.ts @@ -0,0 +1,119 @@ +/* +* @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 pdf = require( './index' ); + + +// TESTS // + +// The function returns a number... +{ + pdf( 2, 0, 4 ); // $ExpectType number + pdf( 1, 2, 8 ); // $ExpectType number +} + +// The compiler throws an error if the function is provided values other than three numbers... +{ + pdf( true, 3, 6 ); // $ExpectError + pdf( false, 2, 4 ); // $ExpectError + pdf( '5', 1, 2 ); // $ExpectError + pdf( [], 1, 2 ); // $ExpectError + pdf( {}, 2, 4 ); // $ExpectError + pdf( ( x: number ): number => x, 2, 4 ); // $ExpectError + + pdf( 9, true, 12 ); // $ExpectError + pdf( 9, false, 12 ); // $ExpectError + pdf( 5, '5', 10 ); // $ExpectError + pdf( 8, [], 16 ); // $ExpectError + pdf( 9, {}, 18 ); // $ExpectError + pdf( 8, ( x: number ): number => x, 16 ); // $ExpectError + + pdf( 9, 5, true ); // $ExpectError + pdf( 9, 5, false ); // $ExpectError + pdf( 5, 2, '5' ); // $ExpectError + pdf( 8, 4, [] ); // $ExpectError + pdf( 9, 4, {} ); // $ExpectError + pdf( 8, 5, ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the function is provided an unsupported number of arguments... +{ + pdf(); // $ExpectError + pdf( 2 ); // $ExpectError + pdf( 2, 0 ); // $ExpectError + pdf( 2, 0, 4, 1 ); // $ExpectError +} + +// Attached to main export is a `factory` method which returns a function... +{ + pdf.factory( 0, 4 ); // $ExpectType Unary +} + +// The `factory` method returns a function which returns a number... +{ + const fcn = pdf.factory( 0, 4 ); + fcn( 2 ); // $ExpectType number +} + +// The compiler throws an error if the function returned by the `factory` method is provided invalid arguments... +{ + const fcn = pdf.factory( 0, 4 ); + fcn( true ); // $ExpectError + fcn( false ); // $ExpectError + fcn( '5' ); // $ExpectError + fcn( [] ); // $ExpectError + fcn( {} ); // $ExpectError + fcn( ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the function returned by the `factory` method is provided an unsupported number of arguments... +{ + const fcn = pdf.factory( 0, 4 ); + fcn(); // $ExpectError + fcn( 2, 0 ); // $ExpectError + fcn( 2, 0, 1 ); // $ExpectError +} + +// The compiler throws an error if the `factory` method is provided values other than two numbers... +{ + pdf.factory( true, 3 ); // $ExpectError + pdf.factory( false, 2 ); // $ExpectError + pdf.factory( '5', 1 ); // $ExpectError + pdf.factory( [], 1 ); // $ExpectError + pdf.factory( {}, 2 ); // $ExpectError + pdf.factory( ( x: number ): number => x, 2 ); // $ExpectError + + pdf.factory( 9, true ); // $ExpectError + pdf.factory( 9, false ); // $ExpectError + pdf.factory( 5, '5' ); // $ExpectError + pdf.factory( 8, [] ); // $ExpectError + pdf.factory( 9, {} ); // $ExpectError + pdf.factory( 8, ( x: number ): number => x ); // $ExpectError + + pdf.factory( [], true ); // $ExpectError + pdf.factory( {}, false ); // $ExpectError + pdf.factory( false, '5' ); // $ExpectError + pdf.factory( {}, [] ); // $ExpectError + pdf.factory( '5', ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the `factory` method is provided an unsupported number of arguments... +{ + pdf.factory( 0 ); // $ExpectError + pdf.factory( 0, 4, 8 ); // $ExpectError +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/examples/c/Makefile b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/examples/c/Makefile new file mode 100644 index 000000000000..c8f8e9a1517b --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/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/stats/base/dists/anglit/pdf/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/examples/c/example.c new file mode 100644 index 000000000000..ffe8f3f46bc3 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/examples/c/example.c @@ -0,0 +1,42 @@ +/** +* @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/stats/base/dists/anglit/pdf.h" +#include +#include + +static double random_uniform( const double min, const double max ) { + double v = (double)rand() / ( (double)RAND_MAX + 1.0 ); + return min + ( v*(max-min) ); +} + +int main( void ) { + double sigma; + double mu; + double x; + double y; + int i; + + for ( i = 0; i < 25; i++ ) { + x = random_uniform( -10.0, 10.0 ); + mu = random_uniform( -5.0, 5.0 ); + sigma = random_uniform( 0.1, 5.0 ); + y = stdlib_base_dists_anglit_pdf( x, mu, sigma ); + printf( "x: %lf, μ: %lf, σ: %lf, f(x;μ,σ): %lf\n", x, mu, sigma, y ); + } +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/examples/index.js b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/examples/index.js new file mode 100644 index 000000000000..66ad0f5f4bdd --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/examples/index.js @@ -0,0 +1,32 @@ +/** +* @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 uniform = require( '@stdlib/random/array/uniform' ); +var logEachMap = require( '@stdlib/console/log-each-map' ); +var pdf = require( './../lib' ); + +var opts = { + 'dtype': 'float64' +}; +var x = uniform( 10, -5.0, 5.0, opts ); +var mu = uniform( x.length, -5.0, 5.0, opts ); +var sigma = uniform( x.length, 0.5, 5.0, opts ); + +logEachMap( 'x: %lf, μ: %lf, σ: %lf, f(x;μ,σ): %lf', x, mu, sigma, pdf ); diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/include.gypi b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/include.gypi new file mode 100644 index 000000000000..bee8d41a2caf --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/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': [ + ' + */ + function pdf( x ) { + var y; + if ( isnan( x ) ) { + return NaN; + } + // Transform x to the standardized form: + y = ( x - mu ) / sigma; + if ( y < -PI_OVER_4 || y > PI_OVER_4 ) { + return 0.0; + } + return cos( 2.0 * y ) / sigma; + } +} + + +// EXPORTS // + +module.exports = factory; diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/lib/index.js b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/lib/index.js new file mode 100644 index 000000000000..bdfa2a48bb7d --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/lib/index.js @@ -0,0 +1,54 @@ +/** +* @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'; + +/** +* Anglit distribution probability density function (PDF). +* +* @module @stdlib/stats/base/dists/anglit/pdf +* +* @example +* var pdf = require( '@stdlib/stats/base/dists/anglit/pdf' ); +* +* var y = pdf( 0.5, 0.0, 1.0 ); +* // returns ~0.540 +* +* var myPDF = pdf.factory( 0.0, 1.0 ); +* y = myPDF( 0.0 ); +* // returns 1.0 +* +* y = myPDF( 1.0 ); +* // returns 0.0 +*/ + +// MODULES // + +var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' ); +var main = require( './main.js' ); +var factory = require( './factory.js' ); + + +// MAIN // + +setReadOnly( main, 'factory', factory ); + + +// EXPORTS // + +module.exports = main; diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/lib/main.js b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/lib/main.js new file mode 100644 index 000000000000..c612146f621d --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/lib/main.js @@ -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. +*/ + +'use strict'; + +// MODULES // + +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var cos = require( '@stdlib/math/base/special/cos' ); +var PI = require( '@stdlib/constants/float64/pi' ); + + +// VARIABLES // + +var PI_OVER_4 = PI / 4.0; + + +// MAIN // + +/** +* Evaluates the probability density function (PDF) for Anglit distribution with location parameter `mu` and scale parameter `sigma` at a value `x`. +* +* @param {number} x - input value +* @param {number} mu - location parameter +* @param {number} sigma - scale parameter +* @returns {number} evaluated PDF +* +* @example +* var y = pdf( 0.0, 0.0, 1.0 ); +* // returns 1.0 +* +* @example +* var y = pdf( 0.5, 0.0, 1.0 ); +* // returns ~0.54 +* +* @example +* var y = pdf( 2.0, 0.0, 1.0 ); +* // returns 0.0 +* +* @example +* var y = pdf( NaN, 0.0, 1.0 ); +* // returns NaN +* +* @example +* var y = pdf( 0.0, NaN, 1.0 ); +* // returns NaN +* +* @example +* var y = pdf( 0.0, 0.0, -1.0 ); +* // returns NaN +*/ +function pdf( x, mu, sigma ) { + var y; + if ( + isnan( x ) || + isnan( mu ) || + isnan( sigma ) || + sigma <= 0.0 + ) { + return NaN; + } + + // Transform x to the standardized form: + y = ( x - mu ) / sigma; + if ( y < -PI_OVER_4 || y > PI_OVER_4 ) { + return 0.0; + } + + return cos( 2.0 * y ) / sigma; +} + + +// EXPORTS // + +module.exports = pdf; diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/lib/native.js b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/lib/native.js new file mode 100644 index 000000000000..5de717d3d83c --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/lib/native.js @@ -0,0 +1,72 @@ +/** +* @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 addon = require( './../src/addon.node' ); + + +// MAIN // + +/** +* Evaluates the probability density function (PDF) for Anglit distribution with location parameter `mu` and scale parameter `sigma` at a value `x`. +* +* @private +* @param {number} x - input value +* @param {number} mu - location parameter +* @param {number} sigma - scale parameter +* @returns {number} evaluated PDF +* +* @example +* var y = pdf( 0.0, 0.0, 1.0 ); +* // returns 1.0 +* +* @example +* var y = pdf( 0.5, 0.0, 1.0 ); +* // returns ~0.540 +* +* @example +* var y = pdf( 2.0, 0.0, 1.0 ); +* // returns 0.0 +* +* @example +* var y = pdf( NaN, 0.0, 1.0 ); +* // returns NaN +* +* @example +* var y = pdf( 0.0, NaN, 1.0 ); +* // returns NaN +* +* @example +* var y = pdf( 0.0, 0.0, NaN ); +* // returns NaN +* +* @example +* var y = pdf( 0.0, 0.0, -1.0 ); +* // returns NaN +*/ +function pdf( x, mu, sigma ) { + return addon( x, mu, sigma ); +} + + +// EXPORTS // + +module.exports = pdf; diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/manifest.json b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/manifest.json new file mode 100644 index 000000000000..4bedefbd2042 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/manifest.json @@ -0,0 +1,82 @@ +{ + "options": { + "task": "build", + "wasm": false + }, + "fields": [ + { + "field": "src", + "resolve": true, + "relative": true + }, + { + "field": "include", + "resolve": true, + "relative": true + }, + { + "field": "libraries", + "resolve": false, + "relative": false + }, + { + "field": "libpath", + "resolve": true, + "relative": false + } + ], + "confs": [ + { + "task": "build", + "wasm": false, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/napi/ternary", + "@stdlib/math/base/assert/is-nan", + "@stdlib/math/base/special/cos", + "@stdlib/constants/float64/pi" + ] + }, + { + "task": "benchmark", + "wasm": false, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/assert/is-nan", + "@stdlib/math/base/special/cos", + "@stdlib/constants/float64/pi" + ] + }, + { + "task": "examples", + "wasm": false, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/assert/is-nan", + "@stdlib/math/base/special/cos", + "@stdlib/constants/float64/pi" + ] + } + ] +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/package.json b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/package.json new file mode 100644 index 000000000000..465e99ea853e --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/package.json @@ -0,0 +1,66 @@ +{ + "name": "@stdlib/stats/base/dists/anglit/pdf", + "version": "0.0.0", + "description": "Anglit distribution probability density function (PDF).", + "license": "Apache-2.0", + "author": { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + }, + "contributors": [ + { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + } + ], + "main": "./lib", + "gypfile": true, + "directories": { + "benchmark": "./benchmark", + "doc": "./docs", + "example": "./examples", + "include": "./include", + "lib": "./lib", + "src": "./src", + "test": "./test" + }, + "types": "./docs/types", + "scripts": {}, + "homepage": "https://github.com/stdlib-js/stdlib", + "repository": { + "type": "git", + "url": "git://github.com/stdlib-js/stdlib.git" + }, + "bugs": { + "url": "https://github.com/stdlib-js/stdlib/issues" + }, + "dependencies": {}, + "devDependencies": {}, + "engines": { + "node": ">=0.10.0", + "npm": ">2.7.0" + }, + "os": [ + "aix", + "darwin", + "freebsd", + "linux", + "macos", + "openbsd", + "sunos", + "win32", + "windows" + ], + "keywords": [ + "stdlib", + "stdmath", + "statistics", + "stats", + "distribution", + "dist", + "probability", + "pdf", + "anglit", + "continuous" + ] +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/src/Makefile b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/src/Makefile new file mode 100644 index 000000000000..2caf905cedbe --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/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/stats/base/dists/anglit/pdf/src/addon.c b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/src/addon.c new file mode 100644 index 000000000000..aab0f6147c84 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/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/stats/base/dists/anglit/pdf.h" +#include "stdlib/math/base/napi/ternary.h" + +STDLIB_MATH_BASE_NAPI_MODULE_DDD_D( stdlib_base_dists_anglit_pdf ) diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/src/main.c b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/src/main.c new file mode 100644 index 000000000000..7699019dd910 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/src/main.c @@ -0,0 +1,56 @@ +/** +* @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/stats/base/dists/anglit/pdf.h" +#include "stdlib/math/base/assert/is_nan.h" +#include "stdlib/math/base/special/cos.h" +#include "stdlib/constants/float64/pi.h" + +static const double PI_OVER_4 = STDLIB_CONSTANT_FLOAT64_PI / 4.0; + +/** +* Evaluates the probability density function (PDF) for an anglit distribution with location `loc` and scale `scale` at a value `x`. +* +* @param x input value +* @param mu location parameter +* @param sigma scale parameter +* @return evaluated PDF +* +* @example +* double y = stdlib_base_dists_anglit_pdf( 0.5, 0.0, 1.0 ); +* // returns ~0.540 +*/ +double stdlib_base_dists_anglit_pdf( const double x, const double mu, const double sigma ) { + double y; + + if ( + stdlib_base_is_nan( x ) || + stdlib_base_is_nan( mu ) || + stdlib_base_is_nan( sigma ) || + sigma <= 0.0 + ) { + return 0.0 / 0.0; // NaN + } + y = ( x - mu ) / sigma; + + if ( y < -PI_OVER_4 || y > PI_OVER_4 ) { + return 0.0; + } + + return stdlib_base_cos( 2.0 * y ) / sigma; +} diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/fixtures/python/data.json b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/fixtures/python/data.json new file mode 100644 index 000000000000..6feb8b26cd02 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/fixtures/python/data.json @@ -0,0 +1 @@ +{"x": [1.819070266563314, 3.23059266160174, 1.3262877768476466, 5.485283470907875, 7.372125864571619, 1.015576463885819, 7.28175908076941, 0.883544734623789, 9.619707005265038, 0.09447078812281351, 0.4556101703665216, 6.61545520785332, 8.413482428572314, 7.599213382349307, 0.9419566136968527, 9.499604218477902, 7.917539520614163, 0.05447314488390087, 5.319954539513018, 0.5792312357620433, 3.042575258987261, 8.05523296308474, 1.6508506968307646, 8.096340754054424, 6.175998666386092, 6.217213795520767, 0.888410627880809, 6.120046480474324, 2.952693991509739, 9.930903642448026, 7.3783391703397285, 5.1748631571132755, 5.667215200915479, 1.4805445761871772, 0.18075574529422767, 6.937493146498358, 4.93849591227482, 3.811129509214344, 4.002441482136179, 3.893953541637524, 0.3794979017226219, 6.28338795521207, 4.344237152616428, 8.219553623846013, 3.4446248703302706, 4.104054478061492, 3.444312807413523, 2.4548997890606685, 1.1574415163131324, 8.292532148914244, 7.08073535953718, 8.470490289835281, 4.930113141739543, 2.508546325389477, 5.0708438408234, 0.55574577189338, 2.202856562348253, 2.594337785199966, 1.5546025892396187, 0.23992171340839197, 4.659316857573011, 2.5053232619096866, 7.470847048302875, 2.4775716398730596, 0.44968956287261497, 2.212322332187523, 4.163581610503815, 6.730541702230287, 7.328133474849642, 3.9038295906635745, 3.4154377794948076, 5.831785323506673, 6.338455501896448, 8.12652879521752, 5.585086169042853, 9.964337296467875, 8.390690107165572, 6.88124173307537, 0.1226683974647036, 7.101162847926041, 5.80178670487125, 4.314303778500788, 1.4456055148321478, 9.622071933740767, 6.680637151843285, 7.070720406363826, 8.79972533592483, 8.013824600637463, 9.30334135498452, 8.194126366967993, 8.583693676350316, 4.715347131281591, 1.8613048541842314, 1.9500606950189237, 0.8567077932381473, 8.60470930685063, 1.0213183961933003, 7.7679732926566185, 2.5444897904049113, 0.8969281944851348, 8.467212623933532, 0.08212589168439588, 6.678596515367054, 8.290391087635298, 6.756292531372655, 4.415185420243884, 4.899972428126241, 9.296947563831665, 6.366607714697246, 7.558256710099735, 2.846636243683255, 2.6341532813230497, 1.8586831904784018, 0.2874880172541072, 2.9788464381086266, 5.658067280329961, 3.66212706530563, 5.103661656875381, 4.674011409596821, 7.762234259528608, 5.321559426364705, 7.98501290057693, 5.607117477308222, 5.6456154068436195, 9.365643672495196, 1.9807734534391497, 2.2832377435782525, 5.128606744458021, 7.965006087773655, 3.0233076205268707, 5.274920182822375, 5.007331627573378, 1.4318039246766012, 9.368706750600866, 8.764336266952839, 4.28315465868435, 1.6143771977192023, 9.966934765694123, 7.377509240858306, 7.890610642642047, 1.4777877578942644, 0.8447732827736831, 3.224841780434856, 2.8583587074868024, 3.430831838145857, 8.712160607792761, 5.850696815762689, 9.311857115002288, 1.5977280071805167, 7.274194538928204, 2.84024257148078, 6.427292412574332, 6.9029647544442705, 8.65618453928066, 1.701915763233297, 6.6272798980862655, 9.966750639644143, 0.06905347772408144, 2.7748111526095567, 1.5664355534212815, 5.371246630385515, 7.395557768568382, 7.175098778304484, 6.150238617841812, 4.340731704496407, 0.5569180015733921, 5.922935731292424, 0.47906323780381066, 9.639653933790111, 6.369980475755883, 5.228364629297611, 2.6845734358494386, 3.054851636842991, 4.731495138663488, 9.112076074703802, 0.48912929566447105, 0.5668278272707106, 8.176090805899936, 3.7515897938198806, 1.2922864408103119, 1.4532234926091636, 9.925251272132098, 7.285678858731006, 8.625278876071933, 8.54209351369785, 6.1583749548518565, 6.846635001928041, 0.11460352731605172, 9.89307653147993, 7.791905924047597, 5.075761015190129, 2.331997389675875, 6.727449689787893, 3.3255540229739458, 3.3014109838479255, 0.21738193567489517, 6.452671459933751, 7.210306071796025, 9.459965409109508, 4.396614397882609, 6.638322309218151, 9.936957762765736, 5.409139505091517, 0.766271678332543, 9.630496390787368, 6.594767796481539, 1.9012642959061277, 6.383435254628762, 8.132591190558498, 0.7719767876307548, 5.919932985769323, 0.19091383344594814, 9.893603638497106, 4.570016291555813, 1.2400066353658956, 5.889734198784305, 9.688032068740283, 1.777792977874867, 2.7732236583822023, 9.793843354986258, 2.580939617277388, 8.312440616248262, 7.571554494190603, 0.16434403076002746, 8.853543004452751, 9.535491786691688, 0.5973079731070363, 8.229398004938437, 7.101981830677294, 8.238224428609813, 2.5177709451284613, 9.706388913660417, 6.734171235632945, 7.188121125564258, 2.910313919450828, 3.8491249214228764, 0.732026438969654, 6.195755564689639, 1.5560309994520838, 6.779493697181895, 9.827052574363423, 6.458199087987817, 3.098418847759582, 8.966983177537877, 8.134049884275925, 8.896281223973293, 2.6147526084827044, 0.6766894297732906, 4.2701655910549725, 1.6306890916536088, 1.6734352649278583, 8.02682355442148, 3.477563930365406, 2.8844575852664, 8.625478615750353, 1.04359034813567, 4.041412897644089, 0.8104108991543613, 8.396514869585118, 3.3577673107636063, 8.715045161027357, 6.4795112605146, 1.626043373257322, 9.995467003977431, 1.2502576790505882, 9.40645114805868, 1.4510491797832237, 6.61024951191484, 1.2620388408538064, 7.787775588988906, 7.665200052440654, 0.6467297835543973, 6.88646974372679, 6.596173161571516, 1.8385364622194356, 7.192626953803179, 6.68051448434746, 4.31676025916008, 4.593997183095037, 0.30972198069444623, 8.849510715271702, 9.30077720081486, 8.40090201777459, 5.156761235899552, 7.547191807870841, 5.633007008281703, 6.9886154530388716, 4.683584663108389, 9.247632805688411, 5.090381201042938, 4.39481215971202, 2.221642096706231, 4.630413377200876, 8.736943757367401, 8.762178256212996, 3.765073711088803, 9.368165605668095, 0.18408151010142149, 9.450026989752423, 1.9859755764895315, 1.6711807589302885, 10.127019167797375, 8.451946265612317, 1.5441422096109214, 17.78351450220481, 18.712503123921376, 0.431255139544231, 10.933586461115894, 11.1504147509047, 17.313320451611798, 13.814315208030356, 2.915930274646763, 18.965766970870177, 1.4331448549694614, 2.812064399021048, 17.74057021401772, 1.6954716904938816, 16.295430216874614, 8.77152160073588, 15.206890379464078, 13.599947287235494, 1.1900827925581203, 14.756167879596418, 9.334456553279736, 16.088415585857827, 11.74176049731069, 6.474290491869259, 1.3763204712635524, 19.83061635407876, 9.508983704437844, 1.333946793916292, 13.430593353610845, 9.515932987438596, 5.7516938763889165, 7.836561391130186, 2.514162836731668, 2.008286837964668, 16.171127641734135, 13.789080236384319, 3.8864446504358385, 8.811385295074043, 8.125606817892763, 10.967584856968537, 8.919227453797278, 2.4481254921500906, 3.6981703993108073, 6.477346611984345, 16.623064363385954, 8.476343465594649, 1.3597321693514353, 3.993373182901714, 12.781675996784934, 11.151842461123149, 12.735852730577065, 5.804536809987595, 19.1098186985244, 14.011730677516079, 19.573299148723194, 14.448466875706828, 12.05812326688245, 2.7846708873292636, 19.53462787704644, 2.7347895011532364, 17.020264117863103, 13.601535466343964, 0.7335396619500378, 6.268197592129203, 3.2639809831979494, 7.536495985723429, 3.134761883926762, 10.83631602104936, 9.592301733604815, 5.7098412523229936, 12.031331538033953, 2.633582276005768, 3.0128128123687636, 15.947821620064698, 19.916104758683865, 11.839379608315665, 10.020190951052601, 15.07716544487648, 19.356629140275633, 4.247538196284184, 11.876306299117918, 13.709543669217766, 13.544104058013698, 17.216885974186734, 4.4617098515193305, 12.9417857998674, 14.921831565463247, 16.432173734942094, 14.435784649135599, 10.65889711435215, 2.7600354691218243, 14.388962880424014, 9.365770098772328, 8.663548452547351, 2.2840913814207275, 17.830724121275402, 13.538740148417581, 14.748179403392745, 5.443968105074692, 5.874260533730413, 3.0909274863881975, 4.679863293176572, 2.9549628683048157, 7.82957588511906, 17.9169945727432, 0.5204656540681984, 9.764571310785893, 5.954394460097998, 7.533045349256751, 13.98871999885657, 3.646859632368966, 6.548517642148348, 1.9405405782913299, 10.783902401069899, 1.180166917660126, 13.623311697952163, 5.783590723506748, 15.889671406875456, 3.799618993521374, 2.413945403961364, 3.5250758740747146, 18.77122404749804, 14.657678378194863, 0.3452889025259309, 1.3435807418508694, 0.3415700545296718, 18.264342735485165, 15.104540116625817, 2.182997016315764, 10.661438176287737, 10.998127929249748, 9.796410130430635, 0.549561532776619, 0.4831349555443154, 9.843483780005112, 18.80562890671465, 10.478691015584547, 15.516498604425376, 10.29583410052286, 14.637373221369606, 4.320923602103219, 16.57360249580399, 13.407081646168233, 0.8930956313400973, 14.959724767727883, 8.985544058191586, 7.287816273096521, 0.03192307922014859, 5.596417381421242, 10.917364654543643, 1.8283508364207868, 9.253088366274456, 12.40789496142173, 18.340736435512582, 6.47262820568423, 13.727916166279373, 12.316025031851051, 1.464557760817382, 12.446568831080473, 0.16123309719949352, 0.19943586905123656, 13.727190060142489, 18.92497265146069, 6.050473791073765, 5.150661766476992, 12.641771681931655, 5.689412094244322, 17.351155109419842, 8.961484492701139, 19.510195036531925, 18.291278967629342, 10.956985001205348, 7.403818344780502, 11.939113311903409, 11.700074162426503, 3.4092336128372303, 10.262219797665978, 4.942577361476772, 19.61268134853374, 19.587512872216625, 9.000771507211706, 17.874992479346474, 12.945077216520355, 1.56099349572425, 11.920057914547037, 17.920865713414194, 2.556372953242838, 1.7960883677305106, 2.1139248005914513, 18.23148121915442, 7.029228988105332, 6.364514808899193, 11.05792148928976, 1.2085918357926495, 18.785842146564274, 6.986580998219576, 15.355865756848818, 1.681575612798334, 13.499860702796704, 18.68565212953175, 0.872330250359763, 13.235799303633538, 4.543253215657616, 16.29952492011022, 6.866919418886404, 8.441866498645956, 17.266610643774534, 1.702659419283421, 18.981135390979688, 10.255505315454432, 7.767039260037151, 19.186173652143932, 18.040251439617933, 7.55763483583187, 0.3459075546541679, 3.9346601889509025, 12.565319233546644, 5.634901391153622, 3.5515657447066307, 18.637811967968563, 15.02755144386278, 4.6256944745011825, 14.152353310096956, 17.98159024467791, 17.488006126059524, 0.6596169942103525, 7.942968875577097, 0.530743155075859, 4.60456703711877, 16.807823105800484, 7.623117421733787, 12.451611580311141, 10.007503084486892, 7.610518430838296, 15.583511150919445, 18.49418944530995, 15.553391827985568, 18.618400109900634, 9.700246814737955, 10.437865260131172, 16.788810246533536, 18.856524782139086, 10.008097714933532, 13.387545575116736, 18.553917195422937, 18.137955299599366, 2.306525235466026, 19.850146188014868, 18.783512554032335, 18.73332485427509, 14.829705943393652, 6.481044567731315, 12.443914056713691, 15.111005455633599, 8.57000079950058, 7.677677140137331, 2.450564297132838, 12.079027201932037, 12.208910118969044, 6.200424376875966, 18.308928582484082, 8.25109818157973, 9.292890431476245, 5.629528574861677, 0.25563098521130145, 7.782817744639368, 9.164307285001815, 0.2718544943908907, 0.7592232421764544, 9.396581239357504, 2.0007336767612283, 5.880937276098079, 19.36949990624163, 10.64017072075201, 15.972458136620977, 11.684539889278799, 3.658073782580755, 15.971779073416146, 19.71886442554393, 3.3337455835060603, 14.742130397017482, 7.150367753418843, 11.93407999347264, 9.154006264443577, 11.726106451503586, 11.802974792664207, 10.681490949013154, 17.573435569293572, 1.424412465515259, 13.226575139701637, 17.204837364984296, 6.332272126978147, 12.651172571383778, 14.110647669918393, 5.641322795545316, 3.665851321149347, 2.8626393328689193, 6.8149390290547505, 20.16781158641504, 31.623574698887147, 0.5984541317196057, 28.37396438060168, 24.57223968341645, 45.2480856209421, 35.03114910600632, 9.760698820814662, 33.17070883905874, 39.08991091278236, 29.72244269778163, 6.82497520574985, 2.1897876572154527, 40.08177486691355, 13.463513216204891, 27.34426308940136, 38.261785165573635, 35.426119856688324, 2.2050820427231055, 24.48688805017825, 41.85949169966076, 40.14382619796985, 3.0183062274173498, 48.14001929952456, 37.200333747061656, 16.7542352308072, 8.385030481868672, 31.57216295480087, 8.70156587840718, 16.02195347298797, 42.98054221431715, 40.162147033494534, 17.073335113345056, 30.113454536816896, 18.27943184012873, 0.7673067505587716, 47.33862234632826, 7.546014876256157, 29.87541854975143, 43.88719385669527, 19.707332172547904, 44.243470608618786, 4.535761309691511, 25.97078498991611, 42.395337681377136, 6.0405706598465, 21.499252215815122, 29.7976907923671, 8.740643468107583, 39.60686960152266, 43.4629411001479, 17.48630776844216, 48.70498378639514, 11.586677402531492, 13.460081037933952, 20.91516523997391, 12.52534387550417, 13.248247260470453, 27.01960747887224, 28.524553736835976, 25.571002328323416, 30.200148490626034, 46.24481086673081, 33.66958530836777, 22.13076883225293, 31.076011479770777, 12.829611550848346, 32.276001805557776, 44.48697509786717, 31.644028492838643, 26.128398619072446, 22.829472654609145, 23.04612684713253, 0.34452186651011685, 43.88619326002249, 1.9142644979500412, 9.199447686646229, 10.570805386629583, 12.943696675523958, 21.66200145044175, 23.582998361221065, 37.33440504443176, 5.485203100128827, 42.48688362395035, 30.671809509407495, 34.14462366010396, 39.46460394277157, 6.44079396981957, 17.319131258276567, 19.26880624331823, 24.557555203426944, 0.09855018630033907, 21.767778331251176, 5.017467938634523, 34.48006738487347, 42.04792885699141, 8.652553590159501, 37.42806552104331, 22.529450348514573, 2.2122976482025383, 7.504310356468452, 40.254957183848646, 0.22787633002895324, 1.4046703726495746, 8.435332217172936, 4.2740681673370196, 9.688845634098724, 41.992213140753805, 8.79079855109864, 21.117120912883923, 11.17776590765695, 43.46694741107373, 41.55709729267098, 16.048595704374346, 14.995797154321833, 46.86363232665006, 46.42305178387316, 9.464485612877988, 40.68014445814501, 29.771156540931887, 3.94079890318722, 38.611839451576216, 32.53900618611452, 22.620097818476463, 25.990820745680104, 0.172464973880343, 14.633121649287755, 1.7537546426685813, 47.62586100538291, 8.96389090501833, 37.857084995487185, 37.598895000560994, 1.3067433619270263, 34.19772691322442, 12.931287010218002, 44.11619970172408, 3.8509614204334577, 33.77686252376789, 17.685912456270525, 2.678079715056553, 5.860637342285696, 13.49343092256403, 48.426491611884146, 10.90227301713309, 46.31236283445859, 49.84774934185128, 16.917519040592353, 47.13233131725122, 31.310428529058075, 26.732241799624195, 11.926209224713535, 23.55183698249142, 24.124609011802477, 2.7673471114359307, 29.193588603098735, 49.39731355480138, 28.480302289796484, 5.624753673998062, 16.732262739534892, 4.544871944014462, 40.98270354640194, 28.829929321370106, 13.673645566451043, 26.155849697974322, 10.264632026885707, 12.566347993670302, 37.52919643458855, 42.19610654795873, 8.113372610245706, 24.037485715371087, 1.6473448368562515, 11.250684443001113, 38.430408476790355, 16.46058153185073, 38.40222410615065, 14.082574546224935, 29.96892306515045, 33.08505971568534, 11.178893535937572, 27.425471847566268, 13.674134745154499, 1.8952505476119097, 1.0613704187910022, 37.792669781520004, 7.1463319581936915, 29.972837156361432, 23.65453619520944, 2.8847250938334126, 15.911321473372148, 36.85192370706788, 30.63337580049921, 44.68131323065451, 3.000833021347271, 0.9333515971484185, 36.849828497588064, 36.18474030727861, 24.89316130522768, 43.61375469987769, 26.499085803800586, 24.926278930749667, 17.15287703624219, 12.038049090858078, 12.391799958097865, 28.24215470727439, 23.31531736289425, 20.5186497081375, 8.066840382780583, 28.565472506245698, 26.915226757582328, 31.741806881366408, 7.549503909344645, 7.895577012123168, 39.604855505046736, 29.080772840353802, 21.50412689971582, 22.322092840014925, 44.29025172234881, 13.103119443383632, 30.606931379526536, 26.09176712586294, 5.924781317355321, 48.660358571880984, 36.56562089518388, 3.7904912060653215, 11.83844180987692, 17.683318023649285, 49.6138502173888, 18.074698171646336, 16.963189506286742, 26.231678551526187, 43.5948199406854, 40.961634053806705, 17.423811665651527, 6.325649871179634, 5.0429575961885575, 36.68938576881774, 10.419334727336139, 14.436207288765985, 13.105460163753785, 29.34882468541964, 23.25730488504158, 38.043911862434776, 13.033972577424802, 18.7526060430445, 12.611952589202652, 4.034891073618035, 43.76042285275739, 22.956845290850946, 11.851659146027337, 26.216095316913286, 6.20268185646839, 9.187830289188437, 49.37570902953062, 17.340694798890098, 0.4607740301306784, 36.67092050025907, 27.043351209322825, 43.50673092198596, 17.562243883476715, 33.40346962275337, 43.781083920423264, 48.317141270120615, 47.30048285658721, 30.392745442467344, 22.64417887042401, 49.445909296873076, 43.02679344484175, 15.114155141719987, 6.796517378109596, 36.03764308841479, 33.8279038893425, 32.31034447732895, 10.30900780539526, 16.834240100845815, 15.395852002541915, 23.81952544791132, 12.664115992963708, 32.17102001116031, 23.55215952844209, 17.773205412491965, 7.74196388969155, 3.804636677834461, 37.60917288734607, 29.32288523100153, 23.728275718028012, 13.544468206865101, 48.033337971539055, 31.27386758165261, 1.1356058061399277, 35.17343469333621, 27.135078847328813, 6.4393800163609125, 46.323769482919005, 39.13899576472772, 23.610395709149934, 47.82209603159441, 23.177847534035205, 38.694839294778696, 45.84311402922171, 17.255600064368615, 20.436287577622693, 47.007023948655394, 33.731395521818456, 10.45408154368561, 2.0253515998049143, 15.736859995784009, 48.18409535859368, 2.606319683869973, 39.039351810591974, 15.30689929538095, 43.42553193362498, 22.529604885140213, 5.289742194778279, 10.19143990312717, 0.2400194802340827, 12.665513526737316, 13.448487288189725, 11.19679881531574, 12.829805467519206, 9.618079866631446, 1.2456265563737468, 37.387848079841646, 36.11412367299071, 1.4413547064957455, 3.6720150863334853, 37.01311113705798, 15.866934961847857, 41.586773844334154, 32.23598309966733, 29.341690528082083, 16.57518376513552, 30.723916164098874, 16.931344678066555, 27.495601972206813, 6.278210108208137, 40.1890920214692, 4.788993125539109, 31.448785255353478, 8.153779383351267, 6.35164101079066, 39.801421069082544, 37.44217154256164, 21.881028227002375, 15.240295008471733, 36.74993850237823, 21.44753694679447, 21.026647950672245, 16.745288616649624, 29.643108732253186, 43.68207671787767, 10.656320038822187, 13.939019935440761, 22.88069902405215, 40.8007803839831, 7.939815634625708, 26.094134369641093, 18.00224773001578, 21.3616865239782, 10.277339271926312, 29.332703695523755, 43.26789066065114, 31.704413612122718, 12.290313852022166, 13.739748259397777, 21.73274387836513, 28.168777272493756, 29.442594932647175, 9.574127689746536, 31.57617481654773, 41.06363251594644, 42.6598578577213, 4.826823291017851, 27.66532788823322, 0.890691476169958, 42.5285076610296, 8.881030093762849, 48.58702649552973, 10.517853239295949, 25.07598192843174, 25.235185100015322, 33.01635492649507, 31.166580550877725, 48.13480374982706, 2.764836777981899, 19.61693342489758, 24.764580525161055, 37.326492797596, 10.356917772294317, 33.189647194630886, 6.634048047029573, 29.848843003812025, 29.244325155364663, 17.899132932916405, 5.375193654078464, 0.8592569908941694, 8.303412961327, 5.454936538248834, 40.37339611628547, 34.38501433414392], "mu": [7.714242172522976, 5.610079066119793, 3.803611725101951, 0.6749597021722631, 9.425761990124997, 8.372592551399194, 1.9112786127606363, 2.0742020577891926, 1.8680937464654834, 2.3502577982698067, 1.4616785910966679, 9.191477755675159, 4.155319302003793, 2.2387147837759835, 1.1122522139825708, 7.6957958355011655, 4.326497400588649, 5.842516952871314, 6.515689713395318, 0.5262525319297084, 8.364049917170526, 1.5264938236848102, 2.4282824130792493, 2.0118646674969587, 8.613932932569334, 1.7826081328256516, 6.388353241688776, 0.5479477827159751, 3.1003300141795433, 0.6167832047854971, 7.7265572495951105, 5.750547259498467, 7.484253383949971, 5.858868075872267, 3.405881736442308, 0.7351958567457506, 1.1820793871928614, 6.836787165881756, 9.461574368271698, 8.748878451259746, 8.810006273768359, 3.505784646088821, 5.588196588303289, 9.914012277226483, 4.204233309021079, 1.9266728824407042, 2.0175074899156864, 7.877712592086846, 5.569888600194757, 7.142265677052722, 6.836102617995299, 5.653778568714344, 0.9564651169475502, 1.267231553798741, 1.8988225522574875, 0.23381927145259107, 0.48333282950672296, 0.9219336660118282, 9.589048665652562, 4.240508074084231, 7.643563324926687, 2.822540644079208, 9.172534236840757, 0.8595106988979873, 2.8799908313958844, 6.716174101529101, 1.813049605775523, 1.8731707232288153, 1.197165305278991, 8.577492114254378, 5.895661559522764, 1.7708156088573823, 3.563347380109285, 3.375017984850052, 8.622752393128406, 8.349251370640888, 5.748209614388494, 9.327010454160593, 0.4769653321187317, 0.8926804083229378, 4.449093452611317, 5.078049492557764, 2.7550938005603767, 4.79894365270872, 7.781592446528592, 0.9231757231224169, 8.544603303997093, 2.0994108086790044, 8.492227844532707, 6.778444238105962, 5.041463693121613, 2.4673674108226775, 3.139103020150263, 2.210641171143558, 4.570368774870023, 3.3062478945665505, 8.063794932138975, 8.475319152608872, 0.8051997886281403, 2.560599393018604, 0.3921762797781514, 3.7868345242457293, 4.061035295234215, 6.964267881339564, 6.534855868719484, 5.449145205655785, 8.698089013721406, 3.048009840233634, 8.349425367668292, 9.31550816171744, 6.6148684240834745, 3.6659437215083535, 5.739154824016831, 9.070728311398792, 8.37069662815398, 9.515981819077235, 7.31532118952519, 1.136674890167001, 5.9914386971055675, 0.3815340018528024, 8.368541429365763, 6.990669031002977, 4.730983012069057, 8.766776907431879, 8.598490423534992, 2.641988520424068, 5.2252754989452255, 2.3735828866476236, 2.8765105191257376, 7.656114739649946, 8.994850127787156, 0.8306709964712922, 8.891361730304563, 5.283343869913814, 2.3857476934537427, 1.9755776625130683, 6.387945095381707, 2.4641409556605307, 3.665821977647271, 8.579486180379964, 2.13091415413302, 4.147019667138539, 4.5403651300284364, 6.326148464589845, 1.9232175416069819, 1.3555199827702824, 9.897020353591172, 6.017132247563344, 8.228950072761364, 2.0778341891447316, 6.221273092448928, 7.241573780020793, 1.0806764587129158, 2.8210542247181913, 6.129098574406515, 2.5171748515201084, 6.3455823364075625, 6.4611341113313125, 7.890796892455211, 5.456234869376929, 7.815034513659124, 2.852813355548421, 9.502844716985614, 6.95281466718577, 7.23464680313465, 4.906241541557973, 7.0235976853726445, 5.242756259181821, 5.2587336274504946, 9.51642043724123, 6.918266019800775, 7.293765842548926, 5.813833521405765, 9.312488522764385, 2.839511669341964, 5.380228641199034, 2.0949549071884253, 6.123843982195596, 7.756452632484492, 0.7776598749523345, 9.93158839836324, 7.498495213490468, 6.276013923957401, 1.1806019225829278, 2.0062480888380154, 5.796455050227189, 2.593851803313777, 1.045520038724922, 1.1022603891581861, 5.823132664289695, 9.037564761308376, 9.698574910603735, 0.7301390052896584, 7.986783836080095, 7.257682733987081, 9.563259969877096, 0.8682767868025898, 3.5519473184476045, 0.4541563647223379, 6.261105145174373, 8.898775342316469, 6.360091917611619, 0.8692156159824116, 3.3369515419039852, 0.2829792222747951, 8.742805065687456, 9.331082670026511, 1.1375411026880766, 8.299152059814388, 4.460892923739772, 4.226876892549746, 9.44687516484731, 7.0031146938222575, 8.666566653223079, 8.503402568676862, 4.923574347961182, 7.100743850716754, 6.647868282252491, 7.6245488309734855, 6.042643844818335, 7.283492853759871, 2.609473722978104, 6.283529128855867, 5.487471984867107, 5.78163965050957, 4.1543252877019565, 0.5773798655261764, 1.5282567090449617, 9.496620400705748, 8.766612810420364, 1.7489419651941585, 2.6328873596052818, 1.7274763170239726, 3.177361258661059, 9.246322727588819, 3.9792264734152605, 6.292891908297441, 3.644776194827503, 1.5228920351075703, 1.104983381400153, 5.64564215180566, 5.309035575545033, 7.433853020671691, 6.114077670014616, 5.947172044794879, 8.564206788804752, 2.45541463164598, 6.33624710826707, 0.3390932561103044, 2.4265263955945904, 0.36528862372611215, 4.389736669188792, 2.8862256271135314, 7.813723937792414, 4.675056499238767, 2.1422326283739457, 8.245821695140709, 8.681817648397484, 2.9816790975596827, 5.634134248692208, 6.966828988585899, 7.638179961764539, 6.673374552626925, 3.7025097743400726, 9.280966323663865, 0.8459268951599586, 3.3864784974372197, 2.2507960246680616, 6.622260182584839, 8.032888504934077, 2.37721843577839, 2.0208504061619346, 7.068167173322509, 6.48580244059599, 3.3358349190377536, 2.760936404081411, 5.084744942508788, 7.817564153852489, 6.824486908070916, 3.2326485349621636, 1.7107261131889273, 5.345577731718584, 8.408619367460089, 4.138867677141182, 9.780303166531684, 2.8141970532186056, 8.807656147611361, 4.486944489414277, 2.591951629207002, 1.6063742981025564, 6.733620480994269, 0.4655682330143729, 7.782124222823245, 6.631591268666441, 0.1890959376222423, 5.951565671355237, 6.678538827840983, 5.602000673865376, 7.804993293301114, 4.6365192243393345, 8.158789225559776, 9.861856887220831, 9.493093525131735, 5.771299240625107, 6.935285191681437, 13.009014394000797, 7.289357615950387, 7.216009772787597, 1.1907156005089514, 9.41001712830002, 17.136194790272587, 5.189123505261282, 7.004105243008116, 3.0989486931143295, 3.2083512766863342, 9.530472406513248, 14.620120329543582, 12.704754037800418, 16.607042299869057, 17.670307471875056, 11.523198455577703, 17.055187474592653, 13.81456358322675, 5.366997186367413, 17.846229313829674, 0.7576626611109738, 6.215730911901862, 6.389384672720679, 13.65847820818615, 8.540800842649293, 5.540082695454895, 0.622358261631013, 15.850368265071392, 4.451844498540867, 8.485257959547898, 10.852617047696192, 5.486735114946586, 17.307780258836, 18.549087714666726, 1.1021435441714922, 18.167091032620736, 0.1629004220483421, 18.19268475079781, 11.62176287663742, 13.219228715800702, 18.045077133681936, 15.772830022172712, 10.986123183252579, 12.023271459703103, 11.047755131301278, 10.62132784247192, 12.993228065456503, 5.589388464672698, 10.906938728929278, 18.286144342225064, 15.102141059556779, 2.2438635502381277, 18.797493686419344, 6.476473052154241, 0.9394640042676072, 13.070171472769578, 0.04608561114322951, 10.222711882198762, 0.8004131540245174, 9.355136986427578, 14.909106814372258, 7.458606240658505, 17.160435013062564, 15.899530466032399, 18.822564634844984, 4.690185629060597, 19.85433807218228, 4.779471560648664, 19.69697780093218, 1.5773872393490618, 17.873794254787665, 14.121278108732369, 1.4439527788960982, 8.451055014999927, 14.154638566591712, 2.6955219466513114, 16.47289633201741, 8.461159513676865, 12.660252051978611, 4.4134867146231045, 0.3492475608649048, 10.146342897243253, 0.04621012465990226, 12.86800711071596, 16.87374364092417, 4.7783332745124145, 9.020381589807611, 4.643383032351023, 10.210622630933093, 16.863661940498616, 7.24689310160297, 19.801346316480004, 16.806750849388962, 3.790372585404025, 17.74209364181303, 18.332511858622645, 17.12603164470071, 10.29585812585568, 16.895779840095237, 0.6693802419524886, 14.783300379122394, 1.5446185975864957, 9.31290258694208, 16.87203175255378, 18.220441292853206, 8.479438808097065, 15.36671484361497, 7.7386653548911895, 16.99499817130067, 17.122072908638103, 8.37306347607809, 1.1206691499129207, 14.125773104299741, 6.469987398542168, 0.28674879092017713, 11.73565811411594, 1.7387626126047007, 2.3093020550948506, 7.058511599068876, 19.805515532877283, 14.472715162735177, 14.112769808085535, 15.934230333198972, 16.161362372167503, 16.999652657757963, 13.32681014152367, 1.7618475586995608, 13.568158908021996, 1.698900793529614, 12.589473503192512, 5.517373953210692, 17.99470717203568, 14.911270083459607, 8.56834279201462, 3.942885934957381, 19.641915869616543, 18.47233924504171, 4.281107294311335, 10.55711586986545, 4.511405597712299, 8.74976822826137, 13.491040302516561, 9.405594886071908, 16.026926423628705, 9.327733704338215, 1.5678432612145476, 10.988428088515377, 18.00810596145836, 8.76390113123325, 0.6169298586348271, 19.254372549539866, 5.698896179468694, 1.7134323192679801, 3.890362499522213, 4.117990781688601, 18.05500157515292, 2.1308712010639774, 4.303457199303626, 7.5642298071890295, 9.978623511243342, 6.363274904112776, 3.832370183189755, 15.232882438340917, 3.2398968696194075, 11.896663846179798, 10.016918110853561, 16.564815110092187, 4.724253986122619, 17.854382193301, 10.76268083560137, 16.16244938341957, 0.6434561981470122, 7.128841025184638, 19.828540017517227, 16.779821002146594, 6.28411394119464, 18.06460805896075, 19.425936714330085, 9.707063775799348, 6.33820401391837, 10.509157902766367, 0.12174524590351066, 0.36655660560241055, 3.814133479462709, 0.6337151564053367, 12.521150464671532, 5.372276157040927, 0.3952763561549366, 6.087025482412858, 16.934733422853057, 14.192038924571248, 15.275352863615918, 14.918677590142996, 9.268889818845008, 12.993605949898726, 16.529372524331855, 15.92024615819316, 16.85570568768209, 2.9044104216450495, 3.1112894840327643, 4.845658316907522, 9.940130189868155, 18.281179870683168, 10.782253958039325, 2.054208650988665, 10.811242924145779, 14.719038105070787, 9.072543701634764, 6.24824032000493, 10.7117317384041, 5.741483998116723, 1.6808641116724776, 6.897320660733389, 9.674793027336303, 10.903521994880492, 10.256197476427268, 10.219575807780695, 2.6735564239784693, 7.850507428685369, 5.311476085835444, 9.845898181628783, 12.419783127952947, 8.329985037641043, 11.6376253264049, 9.944132854291986, 8.27732171272107, 15.901404557338337, 3.5599868259498946, 16.540127467709002, 15.026541798904589, 2.494249453449773, 1.869534745832886, 12.689133262114192, 12.990661836417907, 0.8254166552300912, 5.191837093922739, 17.85321132593949, 14.808431071868984, 2.512099274841928, 12.289500126558739, 4.684093552982862, 2.0872935808976223, 1.4229884291089112, 4.945496219682108, 10.909083898035075, 14.521229836878078, 9.916500807636394, 1.8388785614810899, 9.113735936875166, 16.688137755956134, 2.269023204697911, 7.803930913050854, 12.210786021590911, 16.433036706977205, 14.650468010121996, 10.631219307426127, 5.987029467379044, 9.419993019048174, 14.551696801609218, 15.586444385124631, 9.613547555353087, 17.861250808847835, 16.070405873884628, 5.793530840320722, 19.778230600491014, 7.761636412308183, 15.536076340118854, 13.201207696462019, 8.516943240167159, 2.3232412929910784, 4.61098937749491, 5.552741858621475, 19.22162850643542, 13.027470303299758, 14.264608468302933, 17.538813132224938, 16.70490623322592, 14.041433434328972, 13.732139791292553, 5.42016771521558, 15.002138777010021, 14.895329644320107, 6.421581176671156, 14.986693180661455, 6.71857509136961, 16.470501287356083, 14.468297749831907, 1.751727096271174, 11.315286814799547, 2.854259054416637, 7.1616116123109546, 19.605709996661854, 3.2418191747476843, 11.127144071347761, 14.21071562996885, 33.60886337035129, 5.3042025425003745, 42.063642074702486, 24.005439125053936, 11.712299089977563, 5.281156357194483, 33.97206816403172, 4.163987303423222, 37.425830767836715, 30.419509034543307, 36.391485215431125, 6.3071996850940835, 28.401590030057044, 19.653587284969383, 25.293219273800876, 5.480628513565261, 22.22595709640841, 16.65172532736154, 8.185917669119725, 8.328896064643471, 48.96045707293581, 14.103714345409907, 33.36078266022108, 33.03012891423298, 19.54266902783029, 27.13092160984048, 13.60319245474324, 27.942656571045593, 18.54183217919007, 33.36440486584383, 14.146968701201468, 35.77320867024654, 16.597339386571182, 34.3131984739115, 22.479049113594172, 22.306556398422718, 17.106330352552373, 4.710968362192086, 1.154424042044977, 14.346894214500699, 8.445054114647121, 36.86570994575356, 31.080407864731924, 16.69763635316157, 35.817511463513355, 2.4374933591008343, 39.642805586312875, 5.817386394146606, 19.451936723833906, 7.7927381393227355, 7.369109207486435, 14.565859993792513, 24.951876358457366, 38.178172960132386, 8.766243445973505, 15.1425326711366, 27.426835927571393, 13.407518446133343, 24.330185922456693, 16.889675460960234, 38.09244908232812, 3.2447405492745274, 24.971810812814788, 39.89468536253919, 26.353580653684915, 23.55732897888329, 18.069521399336224, 31.586115306932566, 26.285034711113465, 22.34874660663041, 13.966271751084758, 1.015421294808888, 28.030983839979722, 23.98110486753029, 41.486070529489325, 41.80749471159724, 28.6908898855896, 9.616795931954607, 37.588026099448946, 47.45160642997152, 10.783148234509548, 44.982784077903034, 8.125812143831396, 15.401769168836921, 33.154637901992096, 35.8755566549257, 6.7514261213355855, 41.308246102432165, 22.174936093268727, 26.353169005946732, 43.72669656086236, 21.20224010056572, 7.434828096582541, 2.3431425281755516, 22.329602941818795, 41.05976227530735, 2.2802135093149323, 46.63281006563491, 3.550808183940596, 42.763565671879284, 10.161919457638657, 12.18846578495037, 38.519960047242805, 35.664207946112455, 48.75879551488768, 19.874612771899486, 25.096445379152225, 37.79315867928633, 28.906703248039488, 43.58308922562541, 8.396082841559943, 2.448227465035208, 42.129401300070846, 16.450856815668452, 19.800796770026434, 49.50387942917906, 3.542087892180912, 49.6616188475327, 26.539743485247126, 45.09476171610566, 33.05269200457345, 14.414273322960746, 1.2889417064162245, 4.974685833238429, 36.44704051320844, 17.36323310395576, 7.797198570635322, 31.537789819054307, 9.53289025070302, 25.63990155146727, 7.562697603726792, 39.96655922258747, 24.9360665106306, 27.10806853841486, 29.584329990533924, 5.250125014883849, 16.226311102702468, 48.08896878090458, 46.858648704217146, 16.33854039503777, 5.8057152124365965, 21.69726649519019, 39.38645820937437, 28.94935836129481, 41.38527651346089, 36.20216375986329, 31.340553281639217, 45.672339926911185, 45.45838261438101, 20.80576813466537, 37.62432437469677, 39.610569109181604, 41.001410202902, 9.967469258386451, 20.19721018552323, 15.155688119124761, 1.4375471957852048, 31.662016550220233, 8.937850198091962, 12.23524413549823, 13.344821501803516, 42.996111973394136, 24.84684277987631, 23.837450503862133, 33.55997831630065, 34.26232509552873, 40.509369834170634, 5.447467611764739, 26.628534414667786, 12.067542486624294, 25.600112190798992, 49.01440099378417, 35.935151913564525, 6.219150364307902, 10.124760701886675, 9.632496960555088, 5.994905537082673, 49.9590971774271, 48.905422660968355, 41.10378865556374, 28.810009923917402, 9.01307595686285, 6.122637267699798, 43.73293216795203, 0.026030372510055244, 40.687945309698605, 18.661318300041085, 7.561399960324911, 11.134823172186238, 34.797335700392594, 7.423894499364575, 36.87128437649277, 2.61813649322975, 19.821794371131947, 27.832284935138663, 19.052673656572306, 13.555526692924612, 38.58091527431929, 30.256953703298468, 38.78457444173716, 35.64284380168546, 37.501214070600746, 35.236295458577146, 29.35407440352381, 2.4330321934347143, 47.59282254421272, 2.468528283494437, 28.986006603411425, 46.67674739451437, 7.666789594727269, 26.966737474883523, 18.107284574095168, 32.79007560509537, 41.17978810128072, 41.36541341341498, 35.55506414269756, 15.719369203310846, 35.65178795391829, 2.329125809592403, 11.341260013992244, 7.549442849009552, 43.71705966012305, 48.75090308777985, 28.569484699390557, 12.46508977244697, 49.56633873918069, 32.50617122986013, 46.10596463483796, 17.542297924812907, 16.310668640185373, 37.2509953352262, 14.164508197243798, 1.8991358290769589, 45.38988859338718, 16.61507322550786, 2.7178264555360707, 18.668796693803053, 35.068049117006325, 4.931099524795118, 35.02188136640422, 33.07360406858518, 8.237500587594305, 35.05682884946038, 34.187112862817784, 17.699563723911705, 28.67852510536073, 36.2631841439063, 33.27521067048186, 17.338464215093396, 21.14296412764344, 5.5760472815090445, 36.362145537552244, 28.38001890374458, 14.198278258404013, 31.04317062579921, 46.36882402804135, 23.751493117446326, 27.249359915956884, 30.38344714059366, 13.642787634949705, 24.256838380596417, 36.97945493846646, 5.162924687113496, 37.628648478626644, 22.553994046593523, 11.700081630365528, 5.384671033202798, 22.86081134311859, 38.32611529009447, 22.4795275674576, 21.206076143493934, 37.128788132076664, 40.99072690333718, 13.021483513215305, 19.25891070914874, 18.68380882969685, 39.58209147621476, 25.96196944625424, 29.055710664739674, 39.933715150533075, 20.736708046472195, 49.910960893128994, 40.24231468889851, 39.91972155374859, 49.23838831059902, 46.862841779006764, 35.74573479924854, 35.30034580627876, 31.69251618488289, 20.17185013471926, 17.948323321274007, 3.853005845571922, 4.699509954502129, 42.37974631481951, 3.337442395149648, 47.24447713998664, 4.564648074579713, 47.90821813314158, 27.73184437632873, 42.36874516531015, 38.28506796707374, 12.281734228795122, 15.376053204422785, 44.06516093150219, 23.498768198268987, 31.478007415258453, 19.818927171039753, 9.184997884492569, 18.64773748106638, 26.93124036884644, 33.75577864146824, 24.60769455114314, 0.8660280577721557, 43.11492081809877, 40.347737245932876, 38.67304416424732, 40.694335857488575, 0.4054036617534118, 48.22959467805114, 18.911299044125858, 15.089704273864129, 39.44269410405572, 4.492812364361626, 38.0765125447154, 30.337367271736344, 2.195886966720728, 32.50272879391693, 30.3650276673863, 1.3753516299838786, 48.80577222695564, 25.466937709527077, 19.330995770177594, 48.15913159439837, 11.371230638190216, 12.267275668043098, 39.63813660025995, 35.1980809259667, 44.26370312764299, 44.596535395504375, 27.72984553763733, 10.401759301620356, 29.984638657979527, 21.906560596820473, 39.937880674085754, 1.529759848814738, 11.563718131494326, 43.17130493899219, 13.075863111808495, 14.756144673446258, 31.709177562433045, 30.857758126131053, 17.16484775407941, 44.32371933464933, 41.30766795210364, 22.477212932278274, 38.16696231397521, 6.691907560572075, 26.707132419507257, 41.51872261600504, 8.738578867853263, 17.803617548368884, 20.456615330249278, 32.31683326150725, 29.345522404665044, 0.9688471932676246, 35.081539423393195, 19.300611821952035, 24.781431302372543, 27.49467695895848, 39.019968949212064, 37.40980369003532, 0.3589014295772852, 13.22310738312103, 19.79411961856542, 30.86994321143913, 43.363851752215176, 36.685404365252985, 6.668550176595856, 27.550108967905985, 24.914535681674828, 30.946610138951318, 4.492711868732563, 12.199484788689318, 1.9646683329557046, 4.535151538376808, 23.385928631615528, 15.197670233711547, 28.431179113815453, 19.171250436197656, 0.4456804187745911, 32.48716303460636, 5.973515667964002, 46.99785763500716, 42.475236120751426, 49.61405028568313, 12.187723014677971, 7.240874461258523, 36.04280785705335, 7.531139371004258, 49.07429616832433], "sigma": [3.856632912544532, 1.7441945432215844, 2.2922703014255563, 4.157138241206322, 0.894373013711918, 1.76818966231643, 4.15651361140432, 1.98793500013077, 2.802617848289078, 2.7195869906021275, 0.8663649746406468, 4.720930170685998, 3.471597145649977, 3.168466009275196, 4.970193102366512, 1.194768129952813, 2.8170275513789615, 0.16440750719279407, 1.9288661693934372, 0.1732355776089948, 4.82844358045215, 0.3414376919235621, 1.5057358394491864, 1.4242016966702793, 1.9759426520893468, 3.3887289964179144, 0.28587310443587327, 3.6789008012000495, 3.766473025135317, 2.672691446340385, 1.3853701521742074, 4.499759843063209, 2.202064292905391, 2.452508886445828, 1.6678089366382927, 3.9880661133864557, 4.445191757439675, 3.239859152211954, 2.8424845946917916, 1.0532607279640753, 3.0370056127040903, 2.3518958278859707, 2.042313685577001, 4.00233559677752, 1.332475382597026, 0.9849228891079735, 2.706452656475256, 4.729068159478658, 0.9308525606206399, 3.579918418447298, 3.707124581329154, 1.9572252157262104, 4.371719508577407, 3.260082437447734, 4.700705903746054, 2.1319573154181874, 4.161321762785032, 0.41443117118906375, 0.6034168846152214, 0.956130022855073, 3.2379835648011035, 1.2303789323153094, 1.129385215373073, 4.162942966156895, 1.4154329406184578, 3.3788955469566955, 4.0362711090894, 2.2022199797871482, 1.942450112897543, 3.984156393155894, 4.54434527540465, 2.0296996414999353, 0.13737607943505006, 1.9185752702238905, 1.9361782360040427, 1.09217181982563, 0.8296851802497077, 3.1482133543529707, 0.34509134778804684, 3.3480802319099694, 1.1078935734002082, 2.7552541518661524, 0.8147713841376285, 4.509722103959296, 3.9520016629033248, 2.8621760645182426, 2.634180128376391, 4.045006398294805, 2.2339423917923984, 4.895704495255507, 3.947209414602201, 4.404214497895229, 1.9917633413984104, 4.616295983990325, 0.7287478002805133, 2.6846404509589927, 2.768405403916144, 2.3081721226484744, 1.2056415740325195, 4.276061772173445, 0.38608167863926557, 3.127834801738098, 0.34733171194649737, 4.032776605508736, 3.018735782252432, 0.958056303151363, 1.9440346576414418, 3.5447284545575206, 3.416361188534777, 2.8309846909525125, 4.688893925766011, 4.336980852997743, 2.890433830474783, 0.483773475277106, 1.5377445714383127, 1.607669334587628, 4.86506332649064, 0.9116093803931797, 1.0202108208615073, 0.12793447547392275, 3.61287728323144, 1.514123068911644, 3.457223668125272, 2.941604169085093, 0.4945370195895873, 4.9781305162838905, 1.5131255148075629, 1.8449295627025688, 4.860386196085528, 4.183459076470115, 3.988294500977622, 2.9717959466400514, 0.7654798140672716, 1.4709425263733702, 0.41357114135848805, 1.8721119520165577, 3.4000008591548285, 2.909011099817953, 1.3844228946685857, 1.6537472567364553, 2.2191008895502917, 3.2416117014327765, 3.57615220926095, 2.0309693539631217, 3.2068577753433316, 1.27743854282053, 2.7017490767018417, 3.7135331556132463, 3.3833884830467476, 3.2089617343037453, 3.8300690341058186, 1.2116861030437085, 4.093706751465051, 2.0518018677523457, 1.3621688091055522, 3.764439142459866, 3.571272619978203, 4.086043637762303, 2.9030996488041945, 4.132851315560688, 1.0527077358558408, 0.20765125865383044, 1.3340274063381075, 1.6920763505996734, 4.866398311337427, 0.5838723305254645, 0.624657552564573, 2.755672305925562, 1.670720038158474, 0.9597058226192126, 4.101416996673128, 0.8686689621800692, 3.4893803493351108, 3.9927061441603793, 4.5974371110299055, 3.890544934598083, 2.980876024040234, 1.2871695579243008, 4.2004899842665475, 3.099814874817672, 1.6851987532354062, 2.1828693241626635, 2.4754176348471066, 4.354830075034204, 0.5078217449731308, 1.875543143739347, 3.943437669771931, 2.398776733489535, 1.6439987988857732, 3.6542014080834866, 2.6780474654324937, 4.81464132768398, 2.230561092328039, 4.197879414677828, 2.8146463689779524, 1.0797798276548485, 3.368245341240017, 2.9924985472893675, 1.5764445340559499, 0.8077668460985579, 0.5091687245454327, 4.432911403422256, 0.1924343781227077, 1.6929681757307122, 2.287119578132105, 4.728529013077249, 2.8222398939455196, 2.806476185244594, 0.5290226056103666, 1.8881759018392854, 0.4900939067863139, 3.3065747873402787, 2.0803824810734586, 1.7787575203317199, 4.11980676524544, 0.22444119690383055, 3.550620066917863, 2.3313656818867283, 3.03589272795927, 4.529423414883472, 3.5624160954858337, 0.20401285715525494, 2.3050941940035106, 2.601467214558397, 3.1017487752446167, 3.023548344004421, 2.770986729596427, 4.769946986807909, 1.236998894240379, 2.7488099511593287, 1.1176201528987224, 2.091717088304415, 2.4827895488461635, 0.18922848158810124, 4.770645696033098, 2.9010762398481185, 0.4362518313823823, 0.27311276202547824, 0.720240684461033, 1.1792149393671665, 3.9182685722132216, 4.623781710361255, 1.8426781788336506, 4.574008275153989, 2.272289081006674, 2.7106393859993427, 4.836886522673629, 2.4784081774321822, 2.039340199490963, 4.704185570116163, 3.1451116457835453, 3.9456946808665765, 3.806693689467229, 3.053719884538551, 3.2003348794867787, 0.804547751350647, 4.56490441757884, 3.1870076453847425, 4.616407621996272, 3.6688487430897885, 0.9144441932500731, 4.520318176982046, 4.172689147499621, 3.936287184009885, 3.51833511111574, 3.260525129981629, 0.583622596621965, 0.8088530819227582, 4.622563766436512, 4.195513995670356, 2.2818878323308995, 3.763685215548184, 1.875302704624645, 3.994333654032276, 0.8435935908669122, 4.224163548662894, 1.2882339569363042, 1.651320321566642, 4.31570045960744, 2.5014730937003393, 4.378563060614963, 2.5104978190512606, 4.2489529642281, 2.3299703226096997, 4.634273164183229, 4.400466121861999, 4.422461279244911, 1.430932956442115, 3.8014195643045237, 4.5916472952279515, 1.5676836273899855, 1.6816016350824254, 2.9805200040365056, 0.24115458482847896, 4.611421768000086, 4.598066059787913, 2.2510691255985353, 1.4015063044059188, 0.23731498155691153, 4.819048308733043, 7.232733616880444, 2.5331118005695914, 6.603811270172831, 5.147099727396481, 9.293573078690427, 3.184806175238288, 8.115912044141679, 4.74089753643215, 3.0616118889759982, 4.677962865474118, 9.009207397803484, 8.28494696070872, 1.0313615205349562, 8.570262009976709, 6.991385188217314, 8.739600196893877, 1.9386013520433374, 1.1273431869410493, 4.070537847867977, 8.768612894710687, 5.155920655205335, 2.83418863138213, 1.1372129804623297, 9.825416183262181, 2.9267405943651634, 4.034675253931404, 9.028141790979074, 6.08774640051291, 4.794663738154241, 8.431060905223806, 5.290579512494004, 6.764046536651739, 6.904776404422832, 8.033312145826923, 0.985139243243188, 8.787284859085732, 6.443915573475753, 1.7428536425034622, 9.804481615278672, 5.999957159407063, 5.554971899352367, 2.8471916807229194, 9.519295489184119, 0.5704425057916753, 6.45406475643271, 2.4567987579004344, 4.20568121246208, 0.8335415161269649, 5.723197521844991, 2.8795271887607217, 4.732833040396941, 9.60060920389196, 2.208652099904017, 5.658204166693852, 6.863061621350478, 8.224031538483809, 3.2382765804207607, 0.14284294690095176, 3.039611980013721, 7.434901821160155, 0.12747358090665342, 8.194768394677691, 8.612338824349528, 1.3671064142527765, 6.124028005979855, 1.7161828526275404, 5.791224469026452, 6.430997270239341, 8.146533530406629, 9.702983880402888, 5.459659462015683, 1.606653689447649, 4.347607039062062, 0.946180947218237, 9.464191900834974, 9.60657454865114, 2.6374972914223442, 7.774890230463125, 1.9667559772552692, 6.115118159291711, 9.501884248905267, 3.504169656184825, 2.7865315678653366, 4.804728469816138, 2.751517178752255, 0.837971406923037, 7.7663832355664635, 0.9386162170187087, 6.422167718597173, 5.007532807918477, 5.460178486156502, 0.6130181813568533, 5.100513990700611, 2.371732619102601, 6.871059593905312, 7.278206322094659, 1.7212019714769142, 7.263377062619344, 5.231325907537357, 5.614613891670651, 9.705919991086185, 7.966271403255718, 3.1138288593538497, 7.806377164803534, 5.136408739142448, 0.8555404640320902, 6.653387966882222, 9.813007749114062, 7.130694505751762, 6.2546763135369305, 5.215105707520267, 6.226109726421377, 3.932054232299519, 5.346218247814574, 7.9681830356159855, 3.4934481546166722, 3.8733570420015053, 7.060098438448211, 6.581658631945784, 1.4002759420187443, 1.5590463808429715, 2.7707317146117414, 9.528699568712069, 0.8150290676077246, 5.842990438028331, 1.6654147342277226, 4.912318266075021, 6.308885050813022, 3.2188638754975516, 9.528707123594913, 5.149937309427563, 7.913482924322635, 8.692902847789046, 0.3102091148840749, 6.003210334635903, 6.15177505346721, 2.6383981135014323, 1.2944672879217924, 8.98930412392451, 2.812998366274963, 6.890243875651178, 6.5665151524343175, 4.844754727240413, 6.01051598685416, 9.809394867810946, 0.228036872442257, 5.09226192893818, 6.420570417002566, 6.928004194243381, 5.832508050004096, 7.259486366933047, 3.7562519424007337, 7.65273943918663, 4.084823721700166, 1.0004632924119004, 5.709985736790666, 1.1172895986295794, 8.73360099254275, 3.9718258257010968, 4.89833478041313, 4.356786553960766, 2.3261756075407325, 4.96850872214664, 0.5973107179000369, 5.6430841313888385, 2.9002022911982888, 5.774864747490652, 4.5498502194492625, 8.94572027966873, 5.797162694202455, 1.900013769290238, 0.5087193025478022, 0.5298451351445583, 9.618664007864357, 6.230120935646064, 5.78828085882897, 3.732446860225619, 5.365161110287825, 5.938149754998237, 6.508765201269617, 9.695369774355834, 3.495003610354574, 7.000555788298966, 1.4653663886942263, 8.17727987711245, 1.2133087019114572, 3.6530015459003096, 5.457932039196004, 4.066024137037369, 1.3862087663601872, 3.524805829273854, 2.9376823328948998, 7.491693383515092, 9.328405571120616, 2.608505710528546, 3.991993074504725, 1.6854081338796225, 9.463366145374465, 3.493233574338961, 5.639555449328146, 6.814532330837989, 5.745244556935825, 3.7892133796800636, 3.2876391586924325, 7.184706442389636, 2.167547989836429, 3.4128787631732638, 4.841870300386654, 6.588786873455123, 6.349341876811372, 3.6049256714213054, 7.624767658342103, 7.317804309421253, 3.936052515699401, 7.31897586812677, 9.135590911520866, 4.234176281891258, 8.206366044691146, 6.185098671768452, 8.480453035150758, 3.1839361619327913, 6.366998057041344, 0.6370602734667474, 5.137942899753513, 1.8118353712978943, 0.37776852433478775, 6.170684720688213, 7.0884082258942, 4.673928865042752, 7.576123727581204, 6.333105751419155, 7.078398782413203, 8.714759839851032, 4.116186559768982, 5.570659774589502, 8.457038355575456, 7.615160629293224, 5.660965227654608, 1.0231946302825237, 1.3290210518149135, 0.36507860653529645, 2.3229429195737663, 9.450806656601932, 9.842189137947067, 7.104382623952974, 9.305623408856244, 4.922466511252683, 8.808157989438273, 1.5086911623911632, 8.3396862226982, 1.7092352902634969, 5.770382285928762, 7.62679449039418, 6.220151239423588, 7.952076190910719, 9.855454494856694, 3.6682477029758656, 5.198654144229501, 2.5600720913995874, 8.872719005206637, 7.63555898591385, 9.24314107840022, 7.432145981648243, 0.18485028073403253, 1.645928048526362, 7.891931834261916, 1.770344010465757, 9.30860445888164, 6.174282691685972, 0.9185872023219771, 3.1399983262474023, 2.6859662720743143, 7.821525057641827, 5.298542457829873, 1.8634265301235826, 5.413931471530989, 5.048688803107006, 3.4083537240355324, 0.16392958277103947, 1.1709465103704846, 0.4965553208620864, 1.862818025606515, 1.490988154549092, 4.113256249284564, 2.3920698577476043, 4.40859110491529, 8.314777635251332, 7.315321012326411, 2.312745146252764, 1.8338419004821624, 3.5320980108446745, 1.1380786660719255, 1.0606225038121837, 5.264834860146313, 4.61922901002393, 1.9869839111407939, 4.798355733070753, 7.12889923733184, 4.461394560610919, 5.306204411858702, 15.105568588259036, 18.172730649785112, 17.69768942511893, 10.87231278625956, 17.695909497055194, 2.0858595540482328, 16.952807551034844, 17.827820269004555, 12.685703649034117, 19.35201659169696, 5.803088244094236, 18.45098635009139, 15.680537675592356, 4.798434519567318, 13.854708588528274, 16.19721574184127, 5.562025313397972, 2.2955148131949823, 19.184863182072863, 11.19256893738697, 6.670904934584414, 9.15369978061187, 2.17712187631296, 17.456334247887256, 13.140852007638877, 16.070413723564993, 17.63296670003273, 14.404483350733067, 10.380287032430362, 18.870088392293113, 18.699766413870112, 18.200916705548813, 14.800311313845075, 9.50198483661651, 0.1961024420208615, 0.26432353445974055, 10.8711689207053, 13.49853440353579, 5.932578378230339, 9.151546373904537, 2.236578164473952, 18.48048031072377, 6.730074568848324, 11.464650459789624, 15.43428817621833, 12.773278618717134, 14.668226955209382, 4.7525623788037725, 10.16399243414975, 13.889551679602093, 12.799284214198833, 17.356176243621633, 12.730815022490226, 12.20469419910514, 13.890335797747209, 13.064842074143359, 17.728922416446522, 14.735194770704856, 5.55745556190767, 16.43954301502304, 17.817899604909094, 2.445114038149453, 13.847852197392939, 7.41215111874882, 11.13281856851167, 7.874387217083625, 7.738056261559883, 15.901013816465122, 13.190505570747064, 16.38239974889527, 19.295327482823335, 14.444186871301104, 19.203700764650115, 3.9567726770535154, 1.7520686483029762, 17.649925492560175, 2.4272225051361795, 4.660145435263661, 6.8214760230114795, 8.007939183989047, 2.217783782020216, 1.090312807809837, 4.449533021909361, 12.767000075722871, 12.696972066297144, 18.56233508112174, 7.491436369916796, 15.869789206731, 8.814864168681968, 19.99830435301551, 6.329829792758884, 6.675215655409496, 8.765518298388688, 19.634402589986514, 11.484019535718438, 18.462459274728694, 17.787357857499128, 6.7179010540824144, 15.788087536575894, 9.641448771081281, 18.18294539934914, 16.896813127653967, 7.3582623061700945, 2.1406948803061554, 3.8061518944266726, 17.33518657891634, 11.820425177417292, 1.4980073780921068, 9.227321163546216, 0.9097719244366365, 5.221500937482935, 4.920988783692519, 5.193496784307989, 7.349999614932051, 15.373819275912576, 0.14430214695835317, 3.7834825369710647, 7.274382959731691, 7.68222656461907, 10.304560125508257, 16.540592584426427, 13.057440983104053, 7.860302744867465, 2.8591546052951475, 4.719548539652956, 16.991654006955116, 8.376704567902468, 15.952274725658256, 5.460657253338596, 4.967584663418242, 9.637664753188119, 4.540072403974982, 18.474062711722603, 11.61346745095405, 9.329595095287777, 19.29995095582403, 8.967770655554506, 14.30385322128642, 15.084082741538358, 19.139095424774112, 19.868913501248528, 2.6730809571240672, 19.618879884190374, 1.2244666915503633, 16.136242259367503, 4.403033130424545, 3.587001269590104, 5.657743864342574, 1.3836807068579218, 4.466253753600737, 4.439713104002937, 4.261982610767225, 13.158166649258712, 6.4354180905887715, 4.3882083355289625, 17.17928212889017, 13.856975042936826, 8.041777579955303, 12.753861920620187, 12.139884567850773, 2.3239973882877756, 14.316658210289939, 5.784817250335109, 9.464319521987587, 4.480768085582559, 19.791709944872537, 6.0602739637808885, 8.493703424419424, 7.354541201425219, 15.185860563368808, 13.874989261363918, 18.762767842359644, 10.234846883402618, 0.1863665494070389, 2.546966723741308, 9.626524227599097, 16.927549088581966, 15.810228462712676, 2.413934446638952, 13.979354507553607, 5.306796075040392, 17.459426429605543, 13.829105264384246, 16.370077783510705, 19.042672541548566, 14.251034195845076, 1.825587326935453, 13.280163508806869, 2.579792087958739, 19.25669156333175, 19.732158644831788, 10.678709702689064, 12.910281727386613, 12.133095093472338, 8.005939452743183, 2.8554155655392703, 15.572709363533647, 3.013115274326549, 7.694088103417795, 1.8414308140993643, 7.240614990570591, 1.4312368833715676, 1.5269277890660922, 1.16679608092624, 16.108115809090727, 5.010743541832444, 16.46615746263868, 7.188803327050216, 12.768807908096056, 2.4942341839418356, 5.320535215810324, 12.01977481764481, 7.014850217148018, 14.218833136678153, 13.896578040869839, 8.705588740768718, 0.6974443868487711, 7.418017479545384, 16.632959425996614, 11.863033816330244, 4.0176594238926535, 15.018834046876412, 13.99557009823033, 12.823250768211029, 19.771139068268518, 8.859349682937493, 10.25495994518089, 15.186914529539832, 9.755126085534645, 4.787187531361194, 2.5478788763733764, 9.967864006888048, 5.921388130085179, 9.196194151751307, 7.145258425523293, 12.745708338437689, 8.500141521390834, 12.74405011932035, 9.090139011315456, 3.373476968839762, 7.3267470424629995, 0.6489255959856487, 13.39708031886753, 13.90667897017658, 3.1847598132621555, 1.7999923324699305, 0.7013543404674435, 18.669377484853904, 3.7080857212791973, 3.201165217462127, 4.207811944371463, 14.115168796764374, 13.091005032499341, 12.705904142396582, 11.373243299302196, 4.131280874751754, 11.39115250838506, 17.81390450110983, 1.260006165605522, 6.406614509059932, 17.555907328784453, 4.324657793400211, 10.28059080769103, 4.351972682348412, 0.33940779847077857, 19.535233791411194, 19.887792552400956, 12.641744456873733, 15.85089110931164, 19.976785147422035, 9.67491692897622, 1.4361375965553422, 2.5227697627852623, 16.09237263791612, 5.41838704552988, 0.8835480350254572, 17.72939271439133, 13.889161485515743, 13.574169248501619, 0.3712903240334786, 9.022568886094696, 6.673569504905864, 9.86279002407333, 4.788868413330839, 5.065822267905727, 6.680260169274107, 11.929625070548957, 14.42127204024335, 13.737797369799601, 19.024989472284933, 14.32395258229531, 2.407261324685206, 17.475952361358374, 19.667023304943076, 17.95787989335606, 10.640284868152406, 6.854165282988872, 12.979075057976049, 16.325253750045295, 12.496576834569309, 4.965915491499896, 19.812896130426232, 15.826630242092172, 14.697046137924833, 10.190520443897253, 17.5101620303848, 1.786143251293094, 4.1284146103775345, 15.718249870319305, 5.75091979908366, 11.745646066118077, 13.600109525407268, 17.735167639697636, 2.5997145164175484, 7.476458917094796, 19.305054077432725, 15.8474171211322, 9.135957712758023, 5.2194232790407185, 6.88434775248527, 14.428177561207406, 4.822403766595188, 8.839283217801258, 2.43363574252294, 12.111709224126285, 0.6949458871240134, 13.520071627450191, 1.6724517230328868, 5.5969988904706796, 4.184208179351567, 9.462014938712345, 11.4088596264527, 11.094724427638113, 2.4066574796082145, 5.91252503916684, 6.673125769546815, 6.3255426041160066, 10.743176189447711, 19.030833487462594, 4.0463164936588, 16.78855882431905, 10.472779340290462, 2.10782820418945, 10.247038530226142, 15.82264473108352, 4.614667257175094, 1.5123997210557332, 17.2596932678515, 11.532048518811454, 12.340108940279313, 12.882175923891653, 1.4806565374492298, 0.6517827895692471, 4.797869197030258, 8.692019697104323, 2.869695543048978, 6.494946381446895, 7.463879740315233, 16.778325043751817, 15.097767739279162, 6.139394324067882, 0.8058231966920184, 19.87773634736269, 8.410746492101906, 15.096619579743768, 16.011690462926627, 14.634106217330707, 0.9696911743961956, 11.289805417281972, 10.125073328873482, 10.201892116441607, 15.627630847211067, 10.95247480224291, 15.61092563852357, 18.784885196394317, 4.568752245772708, 6.296600313880429, 17.489359640972495, 0.5205467321157897, 15.28371503244666, 18.820119134350218, 11.964444632978417, 7.899272180022858, 11.15001475879651, 10.43099911483914, 11.623015963416064, 7.0737472629967435, 11.426116559528602, 19.632398643728763, 19.447185210094343, 2.74483122614334, 18.65604797510561, 10.879276863367618, 8.585409008195695, 14.958880361863306, 7.775743722573201, 7.647446123733779, 19.42656404736823, 7.973916340976721, 13.310117824391675], "expected": [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.18327035355242627, 0.0, 0.0, 0.0, 0.09771690719682191, 0.0, 0.0, 0.20072720317144407, 0.0, 0.0, 0.0, 0.16846943892380914, 4.725985007881525, 0.0, 0.0, 0.3404082986476738, 0.0, 0.0, 0.0, 0.0, 0.0, 0.26468494302665246, 0.0, 0.632524985357654, 0.21499873005266315, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.16909658902415556, 0.1655115117259208, 0.3132977703094307, 0.0, 0.18244289269726388, 0.0, 0.0, 0.2236160982531169, 0.26740490694396635, 0.0, 0.0, 0.22201449785019262, 0.046674523543389676, 0.44782476066827653, 0.16280998927594237, 0.0, 0.0, 0.0, 0.0, 0.7070800942684946, 0.0, 0.1712165973751805, 0.0, 0.0, 0.09786801077741991, 0.0, 0.0, 0.0, 0.10146618199972698, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.005414246081967026, 0.0, 0.0, 0.0, 0.30858153521733495, 0.0, 0.0, 0.2147667304055074, 0.0, 0.3725252043086198, 0.0, 0.3347094069447212, 0.17104242490820395, 0.0, 0.118672131834129, 0.142467196550137, 0.21524486851978744, 0.0, 0.0, 0.0, 0.3543843214998888, 0.0, 0.16656113714267626, 0.0, 0.0, 0.0, 0.1962462035485136, 0.32770593723622604, 0.0, 0.0, 0.0, 0.11668162228561758, 0.11424726500630945, 0.0, 0.20496326880339547, 0.0, 0.0, 0.0, 0.0, 0.014169226723688802, 0.0, 0.0, 0.0, 0.0, 0.16810980225635966, 0.25288541084243077, 0.0, 0.0, 0.1938323217419504, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.40669973036063756, 0.374789742031131, 0.0, 0.20730318393041752, 0.0, 0.18385318346402937, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.18553610845263352, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.34442773337883087, 0.07649747207740212, 0.0, 0.0, 0.0, 0.0, 0.0, 0.16561327300057996, 0.0, 0.0, 0.0, 0.0, 0.0, 0.17406069582731862, 0.0, 0.0, 0.30497960500155236, 0.0, 0.0, 0.27685028283229435, 0.0, 0.0, 0.4939616922402261, 0.0, 0.2974888487164803, 0.0, 0.1295760415681616, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.13003937772573235, 0.0, 0.0, 1.5277389523977571, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03187631238481975, 0.0, 0.0, 0.0, 0.0, 0.0, 0.18996587007791757, 0.0, 0.0, 0.0, 0.3608449473906788, 0.0, 0.0, 0.3372387709975117, 0.17334261123581218, 0.0, 0.0, 0.0, 0.0, 0.3433140989767122, 0.0, 0.0, 1.3825502944031387, 0.0, 0.0, 0.1901009239471506, 0.0, 0.06946831042270014, 0.0, 0.3578984175549273, 0.20629601195382222, 0.0, 0.0, 0.20052412360226973, 0.21414243765948404, 0.0, 0.2501185917827432, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.08829189605330387, 0.0, 0.19278469868208253, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.2367245966259613, 0.0, 0.1979552677715985, 0.5232667054387021, 0.2499724459159899, 0.0, 0.0, 0.0, 0.0, 0.11856133279549036, 0.0, 0.0, 0.0, 0.23535054649348586, 0.2755244677645317, 0.1230946600226852, 0.06490708236762745, 0.15382624721544963, 0.6726157993760599, 0.0, 0.011592786773475526, 0.0, 0.0, 0.0, 0.0, 0.0, 0.12632047473438376, 0.0, 0.0, 0.0, 0.09412051553077384, 0.0, 0.3861514372843082, 0.1439620325435996, 0.0, 0.0, 0.0, 0.0, 0.0005263832260570433, 0.0, 0.0, 0.08215946089490182, 0.10297903479193037, 0.0, 0.10797653373557779, 0.14211506237923552, 0.0, 0.0, 0.0, 0.0, 0.09651007432502796, 0.13434929040192403, 0.0, 0.0, 0.07035235434883026, 0.12342417894440899, 0.0, 0.11058311373281805, 0.0, 0.0, 0.11549453603058987, 0.0, 0.0, 0.0, 0.11802007744218693, 0.2549031616212206, 0.0, 0.07319694141519975, 0.1510723733231469, 0.057574269938584395, 0.09989244330082371, 0.0, 0.0, 0.005544667034326236, 0.0, 0.0, 0.0, 0.0, 0.0, 0.056834829898374704, 0.0, 0.0, 0.1040589226387612, 0.0, 0.14107275823119123, 0.0, 0.06827278572652441, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.02987587735100627, 0.0, 0.15828754405023257, 0.0, 0.0, 0.0, 0.0, 0.0, 0.2092648076013674, 0.0, 0.10304044339213606, 0.0, 0.07001792972708114, 0.0, 0.0, 0.03546201080307845, 0.0, 0.0, 0.0, 0.19714487188225366, 0.0, 0.0, 0.0, 0.0, 0.1556695653580716, 0.14202634681990814, 0.0, 0.0, 0.0, 0.0, 0.10938504837618852, 0.11475713319998827, 0.0, 0.0, 0.0, 0.013080122707201642, 0.07618244216595708, 0.0, 0.0, 0.0, 0.001792281284059727, 0.0, 0.0, 0.0, 0.0, 0.0, 0.10865615448690671, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.029995816594105596, 0.10492080243107219, 0.0, 0.1605314725391094, 0.5856991614129341, 0.0, 0.0, 0.0, 0.10127029287432855, 0.14971024151692636, 0.12565147139112762, 0.09025671041018862, 0.0, 0.09951582996556607, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05293254414296308, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.13118965164451674, 0.1256012560127703, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.052969052172869714, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.04250895595993205, 0.0, 0.12453414250971861, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.10353585083630103, 0.16019837740293189, 0.0, 0.0, 0.0, 0.0, 0.0, 0.09921766091596636, 0.18100405076349219, 0.1426936378570003, 0.3891550346873983, 0.0, 0.0, 0.19306806769544183, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.10501963984740459, 0.0, 0.00014891980567864117, 0.0, 0.07229974721134678, 0.07472573744527829, 0.11865100983170478, 0.0, 0.11584308687160816, 0.15008400387252918, 0.09988730699914151, 0.0, 0.44350720087340856, 0.0, 0.0, 0.11505776813902108, 0.0, 0.0, 0.05738098047703322, 0.1295728601900119, 0.1820860634691632, 0.0, 0.0, 0.0, 0.0, 0.0, 0.10897770472233188, 0.0, 0.0, 0.0, 0.0, 0.401194315855871, 0.0, 0.0, 0.0, 0.0, 0.11399951431437541, 0.0, 0.0, 0.0, 0.0, 0.10811368751509441, 0.0, 0.008269974847110197, 0.09206396521038011, 0.0, 0.0, 0.0, 0.0, 0.0636884942600362, 0.0, 0.0, 0.04445959882703627, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.032568477049477934, 0.0, 0.0, 0.12166113926196884, 0.02408757339805167, 0.0871325099148501, 0.0, 0.0, 0.1262935464076655, 0.0, 0.0, 0.1294844138404101, 0.0, 0.0, 0.0, 0.07197287198090023, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.21438467400785208, 0.05178113432236628, 0.040373068785752066, 0.022690457186928296, 0.016438614077750987, 0.0, 0.0, 0.22967749506637616, 0.0, 0.0, 0.0, 0.014042621308893131, 0.0, 0.0, 0.10967263582676567, 0.0, 0.0, 0.0, 0.0, 0.021533251744168418, 0.0, 0.0, 0.0, 0.014050830420909153, 0.0, 0.0, 0.04959944188614635, 0.10302046302632933, 0.0, 0.0, 0.0, 0.0, 0.03409664657765042, 0.0, 0.0, 0.0, 0.04098258210451027, 0.0, 0.0, 0.0, 0.0, 0.0, 0.019434505675202378, 0.05312823039353783, 0.0, 0.0, 0.0, 0.05240493390109607, 0.0385937236446699, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.01882745003006501, 0.0, 0.0, 0.0, 0.0, 0.07132695585667137, 0.0, 0.0, 0.030441129589722875, 0.0, 0.06613723106051016, 0.0740166533960396, 0.04185878871116312, 0.06737882275450348, 0.0, 0.020194473788912773, 0.02675078774875195, 0.0, 0.052585409961922375, 0.0, 0.022786717582038044, 0.1187487671790857, 0.0, 0.0, 0.04662651930435136, 0.0, 0.0, 0.0, 0.044617612079596185, 0.22502921914608112, 0.0, 0.05509296328974116, 0.0, 0.21115262839391458, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0013571696263927713, 0.007925530434127899, 0.045624364001630856, 0.0, 0.06131958245512901, 0.0, 0.030946620648283938, 0.0, 0.08027162384999616, 0.0, 0.0, 0.0, 0.04097811808123134, 0.0, 0.0, 0.0, 0.0, 0.03527778593717684, 0.05031683115855165, 0.12106833221609493, 0.0, 0.0, 0.014328807678232613, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.16756549515528968, 0.0, 0.05744909500562959, 0.0, 0.0, 0.08648137382688491, 0.0, 0.06354838374192863, 0.05576027152382385, 0.0, 0.0, 0.0, 0.0, 0.0310262668778252, 0.0, 0.0, 0.0, 0.0, 0.0994039696281902, 0.13180843146643553, 0.010749710833170343, 0.0, 0.05890302214247436, 0.03624810924477717, 0.0, 0.0, 0.0, 0.05173203434888445, 0.047856162546939876, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.08740902309550734, 0.0680613957109029, 0.06078178231049099, 0.0, 0.06710394642976913, 0.03327748284547668, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05712131683723938, 0.07025156591623558, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03409525823853203, 0.0, 0.0, 0.0, 0.16336115439967766, 0.04930365235818126, 0.05921275166767775, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03668296018292485, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.06861061770045186, 0.0, 0.0, 0.0, 0.0, 0.016343020238995233, 0.0, 0.0, 0.0, 0.07144519908552686, 0.02480211839230755, 0.03232363683642515, 0.0, 0.0, 0.05894416069688833, 0.10190014190163223, 0.20117680190249684, 0.0, 0.0, 0.0, 0.0, 0.0, 0.07329349860826415, 0.0, 0.0, 0.0, 0.0, 0.0718084347416087, 0.0, 0.0673437595862607, 0.0, 0.2478032834785879, 0.0, 0.0, 0.02317980865073521, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05612091032818168, 0.0, 0.053116229807774, 0.028736429364018067, 0.0, 0.0, 0.0, 0.0, 0.022281793236055704, 0.0, 0.0, 0.0, 3.367901754468115e-05, 0.0, 0.0, 0.0, 0.05864352410841891, 0.11544679208192134, 0.0, 0.0, 0.0, 0.06604408217568601, 0.0, 0.0877669880925386, 0.0, 0.0, 0.1189100242344106, 0.0, 0.0, 0.0, 0.0, 0.07265671434444111, 0.0, 0.0, 0.0, 0.05609089659227832, 0.046811352066582594, 0.0, 0.0, 0.0, 0.0, 0.013828959333832784, 0.06796629263103786, 0.0, 0.034316549438058674, 0.0, 0.011146292600039985, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05031082738850309, 0.0, 0.0, 0.0, 0.0, 0.0, 0.05810246130512669, 0.0, 0.0, 0.0, 0.07217734354176637, 0.07762028956640249, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.009875582305056656, 0.0, 0.052764768596691186, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.03503975804970918, 0.08160591202585235, 0.07653564723802564, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.06486366660856965, 0.0, 0.0, 0.050212064094143874, 0.0, 0.0, 0.0, 0.023996042455748174, 0.0, 0.0, 0.009902067311228113, 0.07393114608927329, 0.0, 0.0034299222638891824, 0.03415527827224945, 0.0, 0.20440632317034976, 0.0, 0.0, 0.0, 0.06541451735870582, 0.04364436197569675, 0.0, 0.0, 0.0, 0.0903981926150711, 0.06856796403888264, 0.0, 0.001206144490544889, 0.03174468555550203, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.12574642162151725, 0.0, 0.0, 0.0]} \ No newline at end of file diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/fixtures/python/runner.py b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/fixtures/python/runner.py new file mode 100644 index 000000000000..a706004989c8 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/fixtures/python/runner.py @@ -0,0 +1,98 @@ +#!/usr/bin/env python +# +# @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. + +"""Generate fixtures.""" + +import os +import json +import numpy as np +from scipy.stats import anglit + + +# Get the file path: +FILE = os.path.realpath(__file__) + +# Extract the directory in which this file resides: +DIR = os.path.dirname(FILE) + + +def gen(x, loc, sigma, name): + """Generate fixture data and write to file. + + # Arguments + + * `x`: input values + * `loc`: location parameter + * `sigma`: scale parameter + * `name::str`: output filename + + # Examples + + ``` python + python> x = np.random.uniform(0.0, 10.0, 300) + python> loc = np.random.uniform(0.0, 10.0, 300) + python> sigma = np.random.uniform(0.1, 5.0, 300) + python> gen(x, loc, sigma, './data.json') + ``` + """ + z = anglit.pdf(x, loc, sigma) + + # Store data to be written to file as a dictionary: + data = { + "x": x.tolist(), + "mu": loc.tolist(), + "sigma": sigma.tolist(), + "expected": z.tolist() + } + + # Based on the script directory, create an output filepath: + filepath = os.path.join(DIR, name) + + # Write the data to the output filepath as JSON: + with open(filepath, "w", encoding="utf-8") as outfile: + json.dump(data, outfile) + + +def main(): + """Generate fixture data.""" + # Set seed for reproducibility: + np.random.seed(457) + + # Generate data covering small, medium, and large ranges: + x_small = np.random.uniform(0.0, 10.0, 300) + loc_small = np.random.uniform(0.0, 10.0, 300) + sigma_small = np.random.uniform(0.1, 5.0, 300) + + x_med = np.random.uniform(0.0, 20.0, 300) + loc_med = np.random.uniform(0.0, 20.0, 300) + sigma_med = np.random.uniform(0.1, 10.0, 300) + + x_large = np.random.uniform(0.0, 50.0, 400) + loc_large = np.random.uniform(0.0, 50.0, 400) + sigma_large = np.random.uniform(0.1, 20.0, 400) + + # Concatenate arrays: + x = np.concatenate([x_small, x_med, x_large]) + loc = np.concatenate([loc_small, loc_med, loc_large]) + sigma = np.concatenate([sigma_small, sigma_med, sigma_large]) + + gen(x, loc, sigma, "data.json") + + +if __name__ == "__main__": + main() diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.factory.js b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.factory.js new file mode 100644 index 000000000000..7db8c2cfabd1 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.factory.js @@ -0,0 +1,151 @@ +/** +* @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 isnan = require( '@stdlib/math/base/assert/is-nan' ); +var abs = require( '@stdlib/math/base/special/abs' ); +var PINF = require( '@stdlib/constants/float64/pinf' ); +var NINF = require( '@stdlib/constants/float64/ninf' ); +var EPS = require( '@stdlib/constants/float64/eps' ); +var factory = require( './../lib/factory.js' ); + + +// FIXTURES // + +var data = require( './fixtures/python/data.json' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof factory, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function returns a function', function test( t ) { + var pdf = factory( 0.0, 1.0 ); + t.strictEqual( typeof pdf, 'function', 'returns expected value' ); + t.end(); +}); + +tape( 'if provided `NaN` for any parameter, the created function returns `NaN`', function test( t ) { + var pdf; + var y; + + pdf = factory( 0.0, 1.0 ); + y = pdf( NaN ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + pdf = factory( NaN, 1.0 ); + y = pdf( 0.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + pdf = factory( 1.0, NaN ); + y = pdf( 0.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + pdf = factory( NaN, NaN ); + y = pdf( 0.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + pdf = factory( NaN, NaN ); + y = pdf( NaN ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided a finite `mu` and `sigma`, the function returns a function which returns `0` when provided `+infinity` for `x`', function test( t ) { + var pdf; + var y; + + pdf = factory( 0.5, 1.0 ); + y = pdf( PINF ); + t.strictEqual( y, 0.0, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided a finite `mu` and `sigma`, the function returns a function which returns `0` when provided `-infinity` for `x`', function test( t ) { + var pdf; + var y; + + pdf = factory( 0.5, 1.0 ); + y = pdf( NINF ); + t.strictEqual( y, 0.0, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided `sigma <= 0`, the created function always returns `NaN`', function test( t ) { + var pdf; + var y; + + pdf = factory( 0.0, -1.0 ); + + y = pdf( 2.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = pdf( 0.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + pdf = factory( 0.0, 0.0 ); + y = pdf( 2.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + pdf = factory( 0.0, NINF ); + y = pdf( 2.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the created function evaluates the pdf for `x` given `mu` and `sigma`', function test( t ) { + var expected; + var delta; + var sigma; + var pdf; + var tol; + var mu; + var x; + var y; + var i; + + expected = data.expected; + x = data.x; + mu = data.mu; + sigma = data.sigma; + + for ( i = 0; i < x.length; i++ ) { + pdf = factory( mu[i], sigma[i] ); + y = pdf( x[i] ); + if ( y === expected[i] ) { + t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', σ: '+sigma[i]+', y: '+y+', expected: '+expected[i] ); + } else { + delta = abs( y - expected[ i ] ); + tol = 2.0 * EPS * abs( expected[ i ] ); + t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. σ: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' ); + } + } + t.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.js b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.js new file mode 100644 index 000000000000..e2e51b8221dd --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.js @@ -0,0 +1,38 @@ +/** +* @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 pdf = require( './../lib' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof pdf, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'attached to the main export is a factory method for generating `pdf` functions', function test( t ) { + t.strictEqual( typeof pdf.factory, 'function', 'exports a factory method' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.native.js new file mode 100644 index 000000000000..8739b02e2189 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.native.js @@ -0,0 +1,117 @@ +/** +* @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 tryRequire = require( '@stdlib/utils/try-require' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var abs = require( '@stdlib/math/base/special/abs' ); +var PINF = require( '@stdlib/constants/float64/pinf' ); +var NINF = require( '@stdlib/constants/float64/ninf' ); +var EPS = require( '@stdlib/constants/float64/eps' ); + + +// FIXTURES // + +var data = require( './fixtures/python/data.json' ); + + +// VARIABLES // + +var pdf = tryRequire( resolve( __dirname, './../lib/native.js' ) ); +var opts = { + 'skip': ( pdf instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof pdf, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'if provided `NaN` for any parameter, the function returns `NaN`', opts, function test( t ) { + var y = pdf( NaN, 0.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + y = pdf( 0.0, NaN, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + y = pdf( 0.0, 0.0, NaN ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'if provided `+infinity` for `x` and a finite `mu` and `sigma`, the function returns `0`', opts, function test( t ) { + var y = pdf( PINF, 0.0, 1.0 ); + t.strictEqual( y, 0.0, 'returns expected value' ); + t.end(); +}); + +tape( 'if provided `-infinity` for `x` and a finite `mu` and `sigma`, the function returns `0`', opts, function test( t ) { + var y = pdf( NINF, 0.0, 1.0 ); + t.strictEqual( y, 0.0, 'returns expected value' ); + t.end(); +}); + +tape( 'if provided `sigma <= 0`, the function returns `NaN`', opts, function test( t ) { + var y; + + y = pdf( 0.0, 0.0, -1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = pdf( 0.0, 0.0, 0.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = pdf( 0.0, 0.0, NINF ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function evaluates the pdf for `x` given `mu` and `sigma`', opts, function test( t ) { + var expected; + var delta; + var sigma; + var tol; + var mu; + var x; + var y; + var i; + + expected = data.expected; + x = data.x; + mu = data.mu; + sigma = data.sigma; + + for ( i = 0; i < x.length; i++ ) { + y = pdf( x[i], mu[i], sigma[i] ); + if ( y === expected[i] ) { + t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', σ: '+sigma[i]+', y: '+y+', expected: '+expected[i] ); + } else { + delta = abs( y - expected[ i ] ); + tol = 2.0 * EPS * abs( expected[ i ] ); + t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. σ: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' ); + } + } + t.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.pdf.js b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.pdf.js new file mode 100644 index 000000000000..e864748dd2af --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dists/anglit/pdf/test/test.pdf.js @@ -0,0 +1,108 @@ +/** +* @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 isnan = require( '@stdlib/math/base/assert/is-nan' ); +var abs = require( '@stdlib/math/base/special/abs' ); +var PINF = require( '@stdlib/constants/float64/pinf' ); +var NINF = require( '@stdlib/constants/float64/ninf' ); +var EPS = require( '@stdlib/constants/float64/eps' ); +var pdf = require( './../lib' ); + + +// FIXTURES // + +var data = require( './fixtures/python/data.json' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof pdf, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'if provided `NaN` for any parameter, the function returns `NaN`', function test( t ) { + var y = pdf( NaN, 0.0, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + y = pdf( 0.0, NaN, 1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + y = pdf( 0.0, 0.0, NaN ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + t.end(); +}); + +tape( 'if provided `+infinity` for `x` and a finite `mu` and `sigma`, the function returns `0`', function test( t ) { + var y = pdf( PINF, 0.0, 1.0 ); + t.strictEqual( y, 0.0, 'returns expected value' ); + t.end(); +}); + +tape( 'if provided `-infinity` for `x` and a finite `mu` and `sigma`, the function returns `0`', function test( t ) { + var y = pdf( NINF, 0.0, 1.0 ); + t.strictEqual( y, 0.0, 'returns expected value' ); + t.end(); +}); + +tape( 'if provided `sigma <= 0`, the function returns `NaN`', function test( t ) { + var y; + + y = pdf( 0.0, 0.0, -1.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = pdf( 0.0, 0.0, 0.0 ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + y = pdf( 0.0, 0.0, NINF ); + t.strictEqual( isnan( y ), true, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function evaluates the pdf for `x` given `mu` and `sigma`', function test( t ) { + var expected; + var delta; + var sigma; + var tol; + var mu; + var x; + var y; + var i; + + expected = data.expected; + x = data.x; + mu = data.mu; + sigma = data.sigma; + + for ( i = 0; i < x.length; i++ ) { + y = pdf( x[i], mu[i], sigma[i] ); + if ( y === expected[i] ) { + t.strictEqual( y, expected[i], 'x: '+x[i]+', mu: '+mu[i]+', σ: '+sigma[i]+', y: '+y+', expected: '+expected[i] ); + } else { + delta = abs( y - expected[ i ] ); + tol = 2.0 * EPS * abs( expected[ i ] ); + t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. mu: '+mu[i]+'. σ: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' ); + } + } + t.end(); +});