From 2825e0acc1bc3f816bdea8c8c3b2c7ddd0b7ef3a Mon Sep 17 00:00:00 2001 From: yuvi-mittal Date: Tue, 4 Feb 2025 21:11:17 +0530 Subject: [PATCH 1/3] feat : add C implementation for stats/base/dists/truncated-normal/pdf --- type: pre_push_report description: Results of running various checks prior to pushing changes. report: - task: run_javascript_examples status: na - task: run_c_examples status: na - task: run_cpp_examples status: na - task: run_javascript_readme_examples status: na - task: run_c_benchmarks status: na - task: run_cpp_benchmarks status: na - task: run_fortran_benchmarks status: na - task: run_javascript_benchmarks status: na - task: run_julia_benchmarks status: na - task: run_python_benchmarks status: na - task: run_r_benchmarks status: na - task: run_javascript_tests status: na --- --- .../base/dists/truncated-normal/pdf/README.md | 109 +++++++++- .../pdf/benchmark/benchmark.js | 111 ++++++++++ .../pdf/benchmark/benchmark.native.js | 80 +++++++ .../truncated-normal/pdf/benchmark/c/Makefile | 146 +++++++++++++ .../pdf/benchmark/c/benchmark.c | 0 .../dists/truncated-normal/pdf/binding.gyp | 170 +++++++++++++++ .../truncated-normal/pdf/examples/c/Makefile | 146 +++++++++++++ .../truncated-normal/pdf/examples/c/example.c | 49 +++++ .../dists/truncated-normal/pdf/include.gypi | 53 +++++ .../stats/base/dists/truncated-normal/pdf.h | 38 ++++ .../dists/truncated-normal/pdf/lib/native.js | 82 +++++++ .../dists/truncated-normal/pdf/manifest.json | 91 ++++++++ .../dists/truncated-normal/pdf/package.json | 3 + .../dists/truncated-normal/pdf/src/Makefile | 70 ++++++ .../dists/truncated-normal/pdf/src/addon.c | 23 ++ .../dists/truncated-normal/pdf/src/main.c | 95 +++++++++ .../pdf/test/fixtures/python/data.json | 1 + .../pdf/test/fixtures/python/runner.py | 200 ++++++++++++++++++ .../truncated-normal/pdf/test/test.native.js | 138 ++++++++++++ .../truncated-normal/pdf/test/test.pdf.js | 130 ++++++++++++ 20 files changed, 1733 insertions(+), 2 deletions(-) create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/benchmark/benchmark.js create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/benchmark/benchmark.native.js create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/benchmark/c/Makefile create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/benchmark/c/benchmark.c create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/binding.gyp create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/examples/c/Makefile create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/examples/c/example.c create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/include.gypi create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/include/stdlib/stats/base/dists/truncated-normal/pdf.h create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/lib/native.js create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/manifest.json create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/src/Makefile create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/src/addon.c create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/src/main.c create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/test/fixtures/python/data.json create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/test/fixtures/python/runner.py create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/test/test.native.js create mode 100644 lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/test/test.pdf.js diff --git a/lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/README.md b/lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/README.md index 26232640a149..45c76a6c2a28 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/README.md +++ b/lib/node_modules/@stdlib/stats/base/dists/truncated-normal/pdf/README.md @@ -40,7 +40,7 @@ f(x;\mu,\sigma,a,b) = \begin{cases} \frac{\frac{1}{\sigma}\phi(\frac{x - \mu}{\ -where `Phi` and `phi` denote the [cumulative distribution function][cdf] and [density function][pdf] of the [normal][normal-distribution] distribution, respectively, `mu` is the location and `sigma > 0` is the scale parameter of the distribution. `a` and `b` are the minimum and maximum support. +where `Phi` and `phi` denote the [cumulative distribution function][pdf] and [density function][pdf] of the [normal][normal-distribution] distribution, respectively, `mu` is the location and `sigma > 0` is the scale parameter of the distribution. `a` and `b` are the minimum and maximum support. @@ -145,6 +145,111 @@ for ( i = 0; i < 25; i++ ) { + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/stats/base/dists/truncated-normal/pdf.h" +``` + +#### stdlib_base_dists_truncated_normal_pdf( x, a, b , mu , sigma ) + +Evaluates the r distribution function (pdf) for an truncated-normal distribution. + +```c +double out = stdlib_base_dists_truncated_normal_pdf( 0.9, 0.0, 1.0, 0.0, 1.0 ); +// returns ~0.7795 +``` + +The function accepts the following arguments: + +- **x**: `[in] double` input value. +- **a**: `[in] double` lower limit. +- **b**: `[in] double` upper limit. +- **mu**: `[in] double` location parameter. +- **sigma**: `[in] double` scale parameter. + + +```c +double stdlib_base_dists_truncated-normal_pdf( const double x, const double a, const double b , const double mu , const double sigma ); +``` +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/stats/base/dists/truncated-normal/pdf.h" +#include +#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 x; + double a; + double b; + double mu; + double sigma; + double y; + int i; + + for ( i = 0; i < 25; i++ ) { + x = random_uniform( -10.0, 10.0 ); + a = random_uniform( -20.0, 0.0 ); + b = random_uniform( a, a + 40.0 ); + mu = random_uniform( a, b ); + sigma = random_uniform( 0.1, 5.0 ); + y = stdlib_base_dists_truncated_normal_pdf( x, a, b, mu, sigma ); + printf( "x: %lf, a: %lf, b: %lf, mu: %lf, sigma: %lf, f(x;a,b,mu,sigma): %lf\n", x, a, b, mu, sigma, y ); + } +} +``` + +
+ + + +
+ + +