From 5a772d46bc93d48241ac5c8f3b62537bac1fbbc0 Mon Sep 17 00:00:00 2001
From: Vinit Pandit <vinitsharma2480@gmail.com>
Date: Mon, 16 Dec 2024 21:37:59 +0530
Subject: [PATCH 1/9] feat: add c implementation of
 stats/base/dists/rayleigh/logpdf

---
 .../base/dists/rayleigh/logpdf/README.md      |  95 +++++++++
 .../rayleigh/logpdf/benchmark/benchmark.js    |  22 +-
 .../logpdf/benchmark/benchmark.native.js      |  70 +++++++
 .../rayleigh/logpdf/benchmark/c/Makefile      | 146 +++++++++++++
 .../rayleigh/logpdf/benchmark/c/benchmark.c   | 141 +++++++++++++
 .../base/dists/rayleigh/logpdf/binding.gyp    | 170 ++++++++++++++++
 .../dists/rayleigh/logpdf/examples/c/Makefile | 146 +++++++++++++
 .../rayleigh/logpdf/examples/c/example.c      |  40 ++++
 .../base/dists/rayleigh/logpdf/include.gypi   |  53 +++++
 .../stdlib/stats/base/dists/rayleigh/logpdf.h |  38 ++++
 .../base/dists/rayleigh/logpdf/lib/main.js    |   2 +-
 .../base/dists/rayleigh/logpdf/lib/native.js  |  67 ++++++
 .../base/dists/rayleigh/logpdf/manifest.json  |  89 ++++++++
 .../base/dists/rayleigh/logpdf/package.json   |   3 +
 .../base/dists/rayleigh/logpdf/src/Makefile   |  70 +++++++
 .../base/dists/rayleigh/logpdf/src/addon.c    |  23 +++
 .../base/dists/rayleigh/logpdf/src/main.c     |  57 ++++++
 .../dists/rayleigh/logpdf/test/test.native.js | 191 ++++++++++++++++++
 18 files changed, 1417 insertions(+), 6 deletions(-)
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.native.js
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/Makefile
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/binding.gyp
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/examples/c/Makefile
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/examples/c/example.c
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/include.gypi
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/include/stdlib/stats/base/dists/rayleigh/logpdf.h
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/native.js
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/manifest.json
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/Makefile
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/addon.c
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/main.c
 create mode 100644 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/test/test.native.js

diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/README.md b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/README.md
index ead9028b0838..6adb817687f5 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/README.md
@@ -153,6 +153,101 @@ for ( i = 0; i < 10; i++ ) {
 
 <!-- /.examples -->
 
+<!-- C interface documentation. -->
+
+* * *
+
+<section class="c">
+
+## C APIs
+
+<!-- Section to include introductory text. Make sure to keep an empty line after the intro `section` element and another before the `/section` close. -->
+
+<section class="intro">
+
+</section>
+
+<!-- /.intro -->
+
+<!-- C usage documentation. -->
+
+<section class="usage">
+
+### Usage
+
+```c
+#include "stdlib/stats/base/dists/rayleigh/logpdf.h"
+```
+
+#### stdlib_base_dists_rayleigh_logpdf( x, sigma )
+
+Evaluates the logarithm of the probability density function (PDF) for a Rayleigh distribution.
+
+```c
+double out = stdlib_base_dists_rayleigh_logpdf( 0.3, 1.0 );
+// returns ~-1.249
+```
+
+The function accepts the following arguments:
+
+-   **x**: `[in] double` input value.
+-   **sigma**: `[in] double` scale parameter.
+
+```c
+double stdlib_base_dists_rayleigh_logpdf( const double x, const double sigma );
+```
+
+</section>
+
+<!-- /.usage -->
+
+<!-- C API usage notes. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
+
+<section class="notes">
+
+</section>
+
+<!-- /.notes -->
+
+<!-- C API usage examples. -->
+
+<section class="examples">
+
+### Examples
+
+```c
+#include "stdlib/stats/base/dists/rayleigh/logpdf.h"
+#include <stdlib.h>
+#include <stdio.h>
+
+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 x;
+    double y;
+    int i;
+
+    for ( i = 0; i < 25; i++ ) {
+        x = random_uniform( 0.0, 10.0 );
+        sigma = random_uniform( 0.0, 10.0 );
+        y = stdlib_base_dists_rayleigh_logpdf( x, sigma );
+        printf( "x: %lf, σ: %lf, f(x;σ): %lf\n", x, sigma, y );
+    }
+}
+```
+
+</section>
+
+<!-- /.examples -->
+
+</section>
+
+<!-- /.c -->
+
 <!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. -->
 
 <section class="related">
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js
index 4164fff25eae..5b6b93e25ead 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js
@@ -21,6 +21,7 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
+var Float64Array = require( '@stdlib/array/float64' );
 var randu = require( '@stdlib/random/base/randu' );
 var isnan = require( '@stdlib/math/base/assert/is-nan' );
 var EPS = require( '@stdlib/constants/float64/eps' );
@@ -32,15 +33,21 @@ var logpdf = require( './../lib' );
 
 bench( pkg, function benchmark( b ) {
 	var sigma;
+	var len;
 	var x;
 	var y;
 	var i;
 
+	len = 100;
 	b.tic();
+	x = new Float64Array( len );
+	sigma = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = ( randu() * 100.0 ) - 100.0;
+		sigma[ i ] = ( randu() * 20.0 ) + EPS;
+	}
 	for ( i = 0; i < b.iterations; i++ ) {
-		x = ( randu()*100.0 ) - 100;
-		sigma = ( randu()*20.0 ) + EPS;
-		y = logpdf( x, sigma );
+		y = logpdf( x[ i % 100 ], sigma[ i % 100 ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
@@ -56,17 +63,22 @@ bench( pkg, function benchmark( b ) {
 bench( pkg+':factory', function benchmark( b ) {
 	var mylogpdf;
 	var sigma;
+	var len;
 	var x;
 	var y;
 	var i;
 
+	len = 100;
 	sigma = 4.0;
+	x = new Float64Array( len );
 	mylogpdf = logpdf.factory( sigma );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = ( randu() * 50.0 ) - 25.0;
+	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		x = ( randu()*50.0 ) - 25.0;
-		y = mylogpdf( x );
+		y = mylogpdf( x[ i % len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.native.js
new file mode 100644
index 000000000000..c85b037d4064
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.native.js
@@ -0,0 +1,70 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+// MODULES //
+
+var resolve = require( 'path' ).resolve;
+var bench = require( '@stdlib/bench' );
+var Float64Array = require( '@stdlib/array/float64' );
+var randu = require( '@stdlib/random/base/randu' );
+var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var EPS = require( '@stdlib/constants/float64/eps' );
+var tryRequire = require( '@stdlib/utils/try-require' );
+var pkg = require( './../package.json' ).name;
+
+
+// VARIABLES //
+
+var logpdf = tryRequire( resolve( __dirname, './../lib/native.js' ) );
+var opts = {
+	'skip': ( logpdf instanceof Error )
+};
+
+
+// MAIN //
+
+bench( pkg+'::native', opts, function benchmark( b ) {
+	var sigma;
+	var len;
+	var x;
+	var y;
+	var i;
+
+	len = 100;
+	b.tic();
+	x = new Float64Array( len );
+	sigma = new Float64Array( len );
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = ( randu() * 100.0 ) - 100.0;
+		sigma[ i ] = ( randu() * 20.0 ) + EPS;
+	}
+	for ( i = 0; i < b.iterations; i++ ) {
+		y = logpdf( x[ i % 100 ], sigma[ i % 100 ] );
+		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/rayleigh/logpdf/benchmark/c/Makefile b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/Makefile
new file mode 100644
index 000000000000..f69e9da2b4d3
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/Makefile
@@ -0,0 +1,146 @@
+#/
+# @license Apache-2.0
+#
+# Copyright (c) 2024 The Stdlib Authors.
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#    http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+#/
+
+# 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/rayleigh/logpdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
new file mode 100644
index 000000000000..8026c7364a3c
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
@@ -0,0 +1,141 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+#include "stdlib/stats/base/dists/rayleigh/logpdf.h"
+#include "stdlib/constants/float64/eps.h"
+#include <stdlib.h>
+#include <stdio.h>
+#include <math.h>
+#include <time.h>
+#include <sys/time.h>
+
+#define NAME "rayleigh-cdf"
+#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 x[ 100 ];
+	double sigma;
+	double y;
+	double t;
+	int i;
+
+	sigma = 4.0;
+	for ( i = 0; i < 100; i++ ) {
+		x[ i ] = random_uniform( 0.0, 25.0 );
+	}
+
+	t = tic();
+	for ( i = 0; i < ITERATIONS; i++ ) {
+		y = stdlib_base_dists_rayleigh_logpdf( x[ i%100 ], sigma );
+		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/rayleigh/logpdf/binding.gyp b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/binding.gyp
new file mode 100644
index 000000000000..ec3992233442
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/binding.gyp
@@ -0,0 +1,170 @@
+# @license Apache-2.0
+#
+# Copyright (c) 2024 The Stdlib Authors.
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#    http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+# 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/rayleigh/logpdf/examples/c/Makefile b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/examples/c/Makefile
new file mode 100644
index 000000000000..6aed70daf167
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/examples/c/Makefile
@@ -0,0 +1,146 @@
+#/
+# @license Apache-2.0
+#
+# Copyright (c) 2024 The Stdlib Authors.
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#    http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+#/
+
+# 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/rayleigh/logpdf/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/examples/c/example.c
new file mode 100644
index 000000000000..932429c532b4
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/examples/c/example.c
@@ -0,0 +1,40 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+#include "stdlib/stats/base/dists/rayleigh/logpdf.h"
+#include <stdlib.h>
+#include <stdio.h>
+
+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 x;
+	double y;
+	int i;
+
+	for ( i = 0; i < 25; i++ ) {
+		x = random_uniform( 0.0, 10.0 );
+		sigma = random_uniform( 0.0, 10.0 );
+		y = stdlib_base_dists_rayleigh_logpdf( x, sigma );
+		printf( "x: %lf, σ: %lf, f(x;σ): %lf\n", x, sigma, y );
+	}
+}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/include.gypi b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/include.gypi
new file mode 100644
index 000000000000..575cb043c0bf
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/include.gypi
@@ -0,0 +1,53 @@
+# @license Apache-2.0
+#
+# Copyright (c) 2024 The Stdlib Authors.
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#    http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+# 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': [
+      '<!@(node -e "var arr = require(\'@stdlib/utils/library-manifest\')(\'./manifest.json\',{},{\'basedir\':process.cwd(),\'paths\':\'posix\'}).include; for ( var i = 0; i < arr.length; i++ ) { console.log( arr[ i ] ); }")',
+    ],
+
+    # Add-on destination directory:
+    'addon_output_dir': './src',
+
+    # Source files:
+    'src_files': [
+      '<(src_dir)/addon.c',
+      '<!@(node -e "var arr = require(\'@stdlib/utils/library-manifest\')(\'./manifest.json\',{},{\'basedir\':process.cwd(),\'paths\':\'posix\'}).src; for ( var i = 0; i < arr.length; i++ ) { console.log( arr[ i ] ); }")',
+    ],
+
+    # Library dependencies:
+    'libraries': [
+      '<!@(node -e "var arr = require(\'@stdlib/utils/library-manifest\')(\'./manifest.json\',{},{\'basedir\':process.cwd(),\'paths\':\'posix\'}).libraries; for ( var i = 0; i < arr.length; i++ ) { console.log( arr[ i ] ); }")',
+    ],
+
+    # Library directories:
+    'library_dirs': [
+      '<!@(node -e "var arr = require(\'@stdlib/utils/library-manifest\')(\'./manifest.json\',{},{\'basedir\':process.cwd(),\'paths\':\'posix\'}).libpath; for ( var i = 0; i < arr.length; i++ ) { console.log( arr[ i ] ); }")',
+    ],
+  }, # end variables
+}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/include/stdlib/stats/base/dists/rayleigh/logpdf.h b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/include/stdlib/stats/base/dists/rayleigh/logpdf.h
new file mode 100644
index 000000000000..0e2f7ba91f6c
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/include/stdlib/stats/base/dists/rayleigh/logpdf.h
@@ -0,0 +1,38 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+#ifndef STDLIB_STATS_BASE_DISTS_RAYLEIGH_LOGPDF_H
+#define STDLIB_STATS_BASE_DISTS_RAYLEIGH_LOGPDF_H
+
+/*
+* If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
+*/
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+/**
+* Evaluates the logarithm of the probability density function (PDF) for a Rayleigh distribution with scale parameter `sigma` at a value `x`.
+*/
+double stdlib_base_dists_rayleigh_logpdf( const double x, const double sigma );
+
+#ifdef __cplusplus
+}
+#endif
+
+#endif // !STDLIB_STATS_BASE_DISTS_RAYLEIGH_LOGPDF_H
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/main.js b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/main.js
index 7e035a993c63..29ad4c033c3c 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/main.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/main.js
@@ -79,7 +79,7 @@ function logpdf( x, sigma ) {
 	}
 	s2 = pow( sigma, 2.0 );
 	s2i = 1.0 / s2;
-	return ln( s2i * x ) - (pow( x, 2.0 ) / ( 2.0 * s2 ));
+	return ln( s2i * x ) - ( pow( x, 2.0 ) / ( 2.0 * s2 ) );
 }
 
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/native.js b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/native.js
new file mode 100644
index 000000000000..2d3167779a3c
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/native.js
@@ -0,0 +1,67 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+// MODULES //
+
+var addon = require( './../src/addon.node' );
+
+
+// MAIN //
+
+/**
+* Evaluates the logarithm of the probability density function (PDF) for a Rayleigh distribution with scale parameter `sigma` at a value `x`.
+*
+* @param {number} x - input value
+* @param {NonNegativeNumber} sigma - scale parameter
+* @returns {number} evaluated logPDF
+*
+* @example
+* var y = logpdf( 0.3, 1.0 );
+* // returns ~-1.249
+*
+* @example
+* var y = logpdf( 2.0, 0.8 );
+* // returns ~-1.986
+*
+* @example
+* var y = logpdf( -1.0, 0.5 );
+* // returns -Infinity
+*
+* @example
+* var y = logpdf( 0.0, NaN );
+* // returns NaN
+*
+* @example
+* var y = logpdf( NaN, 2.0 );
+* // returns NaN
+*
+* @example
+* // Negative scale parameter:
+* var y = logpdf( 2.0, -1.0 );
+* // returns NaN
+*/
+function logpdf( x, sigma ) {
+	return addon( x, sigma );
+}
+
+
+// EXPORTS //
+
+module.exports = logpdf;
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/manifest.json b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/manifest.json
new file mode 100644
index 000000000000..4b9f8e03256e
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/manifest.json
@@ -0,0 +1,89 @@
+{
+  "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/binary",
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/math/base/special/pow",
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/ninf"
+      ]
+    },
+    {
+      "task": "benchmark",
+      "wasm": false,
+      "src": [
+        "./src/main.c"
+      ],
+      "include": [
+        "./include"
+      ],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/constants/float64/eps",
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/math/base/special/pow",
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/ninf"
+      ]
+    },
+    {
+      "task": "examples",
+      "wasm": false,
+      "src": [
+        "./src/main.c"
+      ],
+      "include": [
+        "./include"
+      ],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/math/base/special/pow",
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/ninf"
+      ]
+    }
+  ]
+}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/package.json b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/package.json
index b17482301bea..7ee4a451ed3f 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/package.json
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/package.json
@@ -14,11 +14,14 @@
     }
   ],
   "main": "./lib",
+  "gypfile": true,
   "directories": {
     "benchmark": "./benchmark",
     "doc": "./docs",
     "example": "./examples",
+    "include": "./include",
     "lib": "./lib",
+    "src": "./src",
     "test": "./test"
   },
   "types": "./docs/types",
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/Makefile b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/Makefile
new file mode 100644
index 000000000000..bcf18aa46655
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/Makefile
@@ -0,0 +1,70 @@
+#/
+# @license Apache-2.0
+#
+# Copyright (c) 2024 The Stdlib Authors.
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#    http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+#/
+
+# 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/rayleigh/logpdf/src/addon.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/addon.c
new file mode 100644
index 000000000000..958f7301d557
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/addon.c
@@ -0,0 +1,23 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+#include "stdlib/stats/base/dists/rayleigh/logpdf.h"
+#include "stdlib/math/base/napi/binary.h"
+
+// cppcheck-suppress shadowFunction
+STDLIB_MATH_BASE_NAPI_MODULE_DD_D( stdlib_base_dists_rayleigh_logpdf )
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/main.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/main.c
new file mode 100644
index 000000000000..a3d02a71984f
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/main.c
@@ -0,0 +1,57 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+#include "stdlib/stats/base/dists/rayleigh/logpdf.h"
+#include "stdlib/math/base/assert/is_nan.h"
+#include "stdlib/math/base/special/pow.h"
+#include "stdlib/math/base/special/ln.h"
+#include "stdlib/constants/float64/pinf.h"
+#include "stdlib/constants/float64/ninf.h"
+
+
+/**
+* Evaluates the logarithm of the probability density function (PDF) for a Rayleigh distribution with scale parameter `sigma` at a value `x`.
+*
+* @param x        input value
+* @param sigma    scale parameter
+* @return         evaluated LOGPDF
+*
+* @example
+* double y = stdlib_base_arcsine_cdf( 0.3, 1.0 );
+* // returns ~-1.249
+*/
+double stdlib_base_dists_rayleigh_logpdf( const double x, const double sigma ) {
+	double s2i;
+	double s2;
+	if (
+		stdlib_base_is_nan( x ) ||
+		stdlib_base_is_nan( sigma ) ||
+		sigma < 0.0
+	) {
+		return 0.0/0.0; // NaN
+	}
+	if ( sigma == 0.0 ) {
+		return ( x == 0.0 ) ? STDLIB_CONSTANT_FLOAT64_PINF : STDLIB_CONSTANT_FLOAT64_NINF;
+	}
+	if ( x < 0.0 || x == STDLIB_CONSTANT_FLOAT64_PINF ) {
+		return STDLIB_CONSTANT_FLOAT64_NINF;
+	}
+	s2 = stdlib_base_pow( sigma, 2.0 );
+	s2i = 1.0 / s2;
+	return stdlib_base_ln( s2i * x ) - stdlib_base_pow( x, 2.0 ) / ( 2.0 * s2 );
+}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/test/test.native.js b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/test/test.native.js
new file mode 100644
index 000000000000..38f8758d2c0c
--- /dev/null
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/test/test.native.js
@@ -0,0 +1,191 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+// MODULES //
+
+var 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 smallScale = require( './fixtures/julia/small_scale.json' );
+var mediumScale = require( './fixtures/julia/medium_scale.json' );
+var largeScale = require( './fixtures/julia/large_scale.json' );
+
+
+// VARIABLES //
+
+var logpdf = tryRequire( resolve( __dirname, './../lib/native.js' ) );
+var opts = {
+	'skip': ( logpdf instanceof Error )
+};
+
+
+// TESTS //
+
+tape( 'main export is a function', opts, function test( t ) {
+	t.ok( true, __filename );
+	t.strictEqual( typeof logpdf, '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 = logpdf( NaN, 1.0 );
+	t.equal( isnan( y ), true, 'returns NaN' );
+	y = logpdf( 0.0, NaN );
+	t.equal( isnan( y ), true, 'returns NaN' );
+	t.end();
+});
+
+tape( 'if provided `+infinity` for `x` and a valid `sigma`, the function returns `-Infinity`', opts, function test( t ) {
+	var y = logpdf( PINF, 1.0 );
+	t.equal( y, NINF, 'returns -Infinity' );
+	t.end();
+});
+
+tape( 'if provided `-infinity` for `x` and a valid `sigma`, the function returns `-Infinity`', opts, function test( t ) {
+	var y = logpdf( NINF, 1.0 );
+	t.equal( y, NINF, 'returns -Infinity' );
+	t.end();
+});
+
+tape( 'if provided a negative `sigma`, the function returns `NaN`', opts, function test( t ) {
+	var y;
+
+	y = logpdf( 2.0, -1.0 );
+	t.equal( isnan( y ), true, 'returns NaN' );
+
+	y = logpdf( 0.0, -1.0 );
+	t.equal( isnan( y ), true, 'returns NaN' );
+
+	y = logpdf( 2.0, NINF );
+	t.equal( isnan( y ), true, 'returns NaN' );
+
+	t.end();
+});
+
+tape( 'if provided `sigma` equal to `0`, the function evaluates a degenerate distribution centered at `0`', opts, function test( t ) {
+	var y;
+
+	y = logpdf( 0.0, 0.0 );
+	t.equal( y, PINF, 'returns +infinity for x equal to 0' );
+
+	y = logpdf( 1.0, 0.0 );
+	t.equal( y, NINF, 'returns -infinity ' );
+
+	y = logpdf( PINF, 0.0 );
+	t.equal( y, NINF, 'returns -infinity ' );
+
+	y = logpdf( NINF, 0.0 );
+	t.equal( y, NINF, 'returns -infinity ' );
+
+	y = logpdf( NaN, 0.0 );
+	t.equal( isnan( y ), true, 'returns NaN' );
+
+	t.end();
+});
+
+tape( 'the function evaluates the logpdf for `x` given small scale parameter `sigma`', opts, function test( t ) {
+	var expected;
+	var delta;
+	var sigma;
+	var tol;
+	var i;
+	var x;
+	var y;
+
+	expected = smallScale.expected;
+	x = smallScale.x;
+	sigma = smallScale.sigma;
+	for ( i = 0; i < x.length; i++ ) {
+		y = logpdf( x[i], sigma[i] );
+		if ( expected[i] !== null) {
+			if ( y === expected[i] ) {
+				t.equal( y, expected[i], 'x: '+x[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
+			} else {
+				delta = abs( y - expected[ i ] );
+				tol = 15.0 * EPS * abs( expected[ i ] );
+				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
+			}
+		}
+	}
+	t.end();
+});
+
+tape( 'the function evaluates the logpdf for `x` given medium scale parameter `sigma`', opts, function test( t ) {
+	var expected;
+	var delta;
+	var sigma;
+	var tol;
+	var i;
+	var x;
+	var y;
+
+	expected = mediumScale.expected;
+	x = mediumScale.x;
+	sigma = mediumScale.sigma;
+	for ( i = 0; i < x.length; i++ ) {
+		y = logpdf( x[i], sigma[i] );
+		if ( expected[i] !== null ) {
+			if ( y === expected[i] ) {
+				t.equal( y, expected[i], 'x: '+x[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
+			} else {
+				delta = abs( y - expected[ i ] );
+				tol = 3.0 * EPS * abs( expected[ i ] );
+				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
+			}
+		}
+	}
+	t.end();
+});
+
+tape( 'the function evaluates the logpdf for `x` given large scale parameter `sigma`', opts, function test( t ) {
+	var expected;
+	var delta;
+	var sigma;
+	var tol;
+	var i;
+	var x;
+	var y;
+
+	expected = largeScale.expected;
+	x = largeScale.x;
+	sigma = largeScale.sigma;
+	for ( i = 0; i < x.length; i++ ) {
+		y = logpdf( x[i], sigma[i] );
+		if ( expected[i] !== null ) {
+			if ( y === expected[i] ) {
+				t.equal( y, expected[i], 'x: '+x[i]+', sigma: '+sigma[i]+', y: '+y+', expected: '+expected[i] );
+			} else {
+				delta = abs( y - expected[ i ] );
+				tol = 3.0 * EPS * abs( expected[ i ] );
+				t.ok( delta <= tol, 'within tolerance. x: '+x[ i ]+'. sigma: '+sigma[i]+'. y: '+y+'. E: '+expected[ i ]+'. Δ: '+delta+'. tol: '+tol+'.' );
+			}
+		}
+	}
+	t.end();
+});

From 21bfea299bff8b9cc30598b928c68a35bf32933d Mon Sep 17 00:00:00 2001
From: Vinit Pandit <vinitsharma2480@gmail.com>
Date: Tue, 17 Dec 2024 10:51:45 +0530
Subject: [PATCH 2/9] fix: fixing c benchmark

---
 .../base/dists/rayleigh/logpdf/benchmark/c/benchmark.c    | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
index 8026c7364a3c..ed4c22d9a422 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
@@ -95,19 +95,19 @@ static double random_uniform( const double min, const double max ) {
 static double benchmark( void ) {
 	double elapsed;
 	double x[ 100 ];
-	double sigma;
+	double sigma[ 100 ];
 	double y;
 	double t;
 	int i;
 
-	sigma = 4.0;
 	for ( i = 0; i < 100; i++ ) {
-		x[ i ] = random_uniform( 0.0, 25.0 );
+		x[ i ] = random_uniform( 20, 20 + STDLIB_CONSTANT_FLOAT64_EPS );
+		sigma[ i ] = random_uniform( -100, 0 );
 	}
 
 	t = tic();
 	for ( i = 0; i < ITERATIONS; i++ ) {
-		y = stdlib_base_dists_rayleigh_logpdf( x[ i%100 ], sigma );
+		y = stdlib_base_dists_rayleigh_logpdf( x[ i % 100 ], sigma[ i % 100 ] );
 		if ( y != y ) {
 			printf( "should not return NaN\n" );
 			break;

From cdb3ff7ff9fd44ce0ffc090010d93851b99427f4 Mon Sep 17 00:00:00 2001
From: Vinit Pandit <106718914+MeastroZI@users.noreply.github.com>
Date: Tue, 17 Dec 2024 10:58:46 +0530
Subject: [PATCH 3/9] fix: clean ups

Signed-off-by: Vinit Pandit <106718914+MeastroZI@users.noreply.github.com>
---
 .../stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js     | 2 +-
 .../base/dists/rayleigh/logpdf/benchmark/benchmark.native.js    | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js
index 5b6b93e25ead..759acb0120b2 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js
@@ -47,7 +47,7 @@ bench( pkg, function benchmark( b ) {
 		sigma[ i ] = ( randu() * 20.0 ) + EPS;
 	}
 	for ( i = 0; i < b.iterations; i++ ) {
-		y = logpdf( x[ i % 100 ], sigma[ i % 100 ] );
+		y = logpdf( x[ i % len ], sigma[ i % len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.native.js
index c85b037d4064..8d0e52019bf2 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.native.js
@@ -56,7 +56,7 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 		sigma[ i ] = ( randu() * 20.0 ) + EPS;
 	}
 	for ( i = 0; i < b.iterations; i++ ) {
-		y = logpdf( x[ i % 100 ], sigma[ i % 100 ] );
+		y = logpdf( x[ i % len ], sigma[ i % len ] );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}

From 03d5538cef8f6ec903c5473f2ea82ce2b43d6a9c Mon Sep 17 00:00:00 2001
From: Vinit Pandit <106718914+MeastroZI@users.noreply.github.com>
Date: Tue, 17 Dec 2024 11:41:32 +0530
Subject: [PATCH 4/9] fix: add private annotation

Signed-off-by: Vinit Pandit <106718914+MeastroZI@users.noreply.github.com>
---
 .../@stdlib/stats/base/dists/rayleigh/logpdf/lib/native.js       | 1 +
 1 file changed, 1 insertion(+)

diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/native.js b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/native.js
index 2d3167779a3c..cde8f69b8195 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/native.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/lib/native.js
@@ -28,6 +28,7 @@ var addon = require( './../src/addon.node' );
 /**
 * Evaluates the logarithm of the probability density function (PDF) for a Rayleigh distribution with scale parameter `sigma` at a value `x`.
 *
+* @private
 * @param {number} x - input value
 * @param {NonNegativeNumber} sigma - scale parameter
 * @returns {number} evaluated logPDF

From 20ebafa193e645e150bbdaba0fd76159e6768187 Mon Sep 17 00:00:00 2001
From: Vinit Pandit <106718914+MeastroZI@users.noreply.github.com>
Date: Tue, 17 Dec 2024 20:36:21 +0530
Subject: [PATCH 5/9] fix: logical fix and clean up

Signed-off-by: Vinit Pandit <106718914+MeastroZI@users.noreply.github.com>
---
 .../stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c    | 2 +-
 .../@stdlib/stats/base/dists/rayleigh/logpdf/src/main.c         | 1 -
 2 files changed, 1 insertion(+), 2 deletions(-)

diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
index ed4c22d9a422..4dfb190c2be0 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
@@ -101,7 +101,7 @@ static double benchmark( void ) {
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		x[ i ] = random_uniform( 20, 20 + STDLIB_CONSTANT_FLOAT64_EPS );
+		x[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, STDLIB_CONSTANT_FLOAT64_EPS  + 20.0 );
 		sigma[ i ] = random_uniform( -100, 0 );
 	}
 
diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/main.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/main.c
index a3d02a71984f..444f27e7ee30 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/main.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/src/main.c
@@ -23,7 +23,6 @@
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/ninf.h"
 
-
 /**
 * Evaluates the logarithm of the probability density function (PDF) for a Rayleigh distribution with scale parameter `sigma` at a value `x`.
 *

From 50f321c7bb32c1aaf3db1bf097bcbb7b864c4535 Mon Sep 17 00:00:00 2001
From: Vinit Pandit <106718914+MeastroZI@users.noreply.github.com>
Date: Tue, 17 Dec 2024 20:37:26 +0530
Subject: [PATCH 6/9] fix: undoing by mistake changes on benchmark:factory

Signed-off-by: Vinit Pandit <106718914+MeastroZI@users.noreply.github.com>
---
 .../base/dists/rayleigh/logpdf/benchmark/benchmark.js    | 9 ++-------
 1 file changed, 2 insertions(+), 7 deletions(-)

diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js
index 759acb0120b2..99689ad9967a 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/benchmark.js
@@ -63,22 +63,17 @@ bench( pkg, function benchmark( b ) {
 bench( pkg+':factory', function benchmark( b ) {
 	var mylogpdf;
 	var sigma;
-	var len;
 	var x;
 	var y;
 	var i;
 
-	len = 100;
 	sigma = 4.0;
-	x = new Float64Array( len );
 	mylogpdf = logpdf.factory( sigma );
-	for ( i = 0; i < len; i++ ) {
-		x[ i ] = ( randu() * 50.0 ) - 25.0;
-	}
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		y = mylogpdf( x[ i % len ] );
+		x = ( randu()*50.0 ) - 25.0;
+		y = mylogpdf( x );
 		if ( isnan( y ) ) {
 			b.fail( 'should not return NaN' );
 		}

From ca4392bf307bf00c4382b67c98b578caf00a1837 Mon Sep 17 00:00:00 2001
From: Vinit Pandit <106718914+MeastroZI@users.noreply.github.com>
Date: Tue, 17 Dec 2024 20:44:10 +0530
Subject: [PATCH 7/9] fix: clean ups

Signed-off-by: Vinit Pandit <106718914+MeastroZI@users.noreply.github.com>
---
 .../stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c    | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
index 4dfb190c2be0..03173d15ec40 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/benchmark/c/benchmark.c
@@ -24,7 +24,7 @@
 #include <time.h>
 #include <sys/time.h>
 
-#define NAME "rayleigh-cdf"
+#define NAME "rayleigh-logpdf"
 #define ITERATIONS 1000000
 #define REPEATS 3
 

From 1b4e6f5733d965e265547fba45cecdc25774a684 Mon Sep 17 00:00:00 2001
From: Philipp Burckhardt <pburckhardt@outlook.com>
Date: Sat, 28 Dec 2024 20:19:37 -0500
Subject: [PATCH 8/9] Update
 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/README.md

Signed-off-by: Philipp Burckhardt <pburckhardt@outlook.com>
---
 .../@stdlib/stats/base/dists/rayleigh/logpdf/README.md          | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/README.md b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/README.md
index 6adb817687f5..09a23396b958 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/README.md
@@ -235,7 +235,7 @@ int main( void ) {
         x = random_uniform( 0.0, 10.0 );
         sigma = random_uniform( 0.0, 10.0 );
         y = stdlib_base_dists_rayleigh_logpdf( x, sigma );
-        printf( "x: %lf, σ: %lf, f(x;σ): %lf\n", x, sigma, y );
+        printf( "x: %lf, σ: %lf, ln(f(x;σ)): %lf\n", x, sigma, y );
     }
 }
 ```

From 6c21df07ef32ab4a1a5d358501e4e009a2ff8170 Mon Sep 17 00:00:00 2001
From: Philipp Burckhardt <pburckhardt@outlook.com>
Date: Sat, 28 Dec 2024 20:19:43 -0500
Subject: [PATCH 9/9] Update
 lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/examples/c/example.c

Signed-off-by: Philipp Burckhardt <pburckhardt@outlook.com>
---
 .../stats/base/dists/rayleigh/logpdf/examples/c/example.c       | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/examples/c/example.c
index 932429c532b4..1dda8e50ec8a 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/examples/c/example.c
+++ b/lib/node_modules/@stdlib/stats/base/dists/rayleigh/logpdf/examples/c/example.c
@@ -35,6 +35,6 @@ int main( void ) {
 		x = random_uniform( 0.0, 10.0 );
 		sigma = random_uniform( 0.0, 10.0 );
 		y = stdlib_base_dists_rayleigh_logpdf( x, sigma );
-		printf( "x: %lf, σ: %lf, f(x;σ): %lf\n", x, sigma, y );
+		printf( "x: %lf, σ: %lf, ln(f(x;σ)): %lf\n", x, sigma, y );
 	}
 }