From e27e42552984783e2a7370a3482a967b3f22968e Mon Sep 17 00:00:00 2001 From: Ujjwal Verma Date: Fri, 24 Jul 2026 15:36:35 +0530 Subject: [PATCH 01/11] feat: add C implementation for invgamma/quantile --- .../base/dists/invgamma/quantile/README.md | 91 ++++++++++++ .../quantile/benchmark/benchmark.native.js | 64 ++++++++ .../invgamma/quantile/benchmark/c/Makefile | 60 ++++++++ .../invgamma/quantile/benchmark/c/benchmark.c | 138 ++++++++++++++++++ .../base/dists/invgamma/quantile/binding.gyp | 90 ++++++++++++ .../invgamma/quantile/examples/c/Makefile | 60 ++++++++ .../invgamma/quantile/examples/c/example.c | 42 ++++++ .../base/dists/invgamma/quantile/include.gypi | 33 +++++ .../stats/base/dists/invgamma/quantile.h | 43 ++++++ .../base/dists/invgamma/quantile/lib/index.js | 15 +- .../dists/invgamma/quantile/lib/native.js | 80 ++++++++++ .../dists/invgamma/quantile/manifest.json | 74 ++++++++++ .../base/dists/invgamma/quantile/package.json | 8 +- .../base/dists/invgamma/quantile/src/Makefile | 70 +++++++++ .../base/dists/invgamma/quantile/src/addon.c | 22 +++ .../base/dists/invgamma/quantile/src/main.c | 49 +++++++ .../invgamma/quantile/test/test.native.js | 128 ++++++++++++++++ 17 files changed, 1062 insertions(+), 5 deletions(-) create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/benchmark/benchmark.native.js create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/benchmark/c/Makefile create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/benchmark/c/benchmark.c create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/binding.gyp create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/examples/c/Makefile create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/examples/c/example.c create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/include.gypi create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/include/stdlib/stats/base/dists/invgamma/quantile.h create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/lib/native.js create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/manifest.json create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/src/Makefile create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/src/addon.c create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/src/main.c create mode 100644 lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/test/test.native.js diff --git a/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/README.md b/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/README.md index feafe97f0993..d98bf66895ab 100644 --- a/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/README.md +++ b/lib/node_modules/@stdlib/stats/base/dists/invgamma/quantile/README.md @@ -144,6 +144,97 @@ logEachMap( 'p: %0.4f, α: %0.4f, β: %0.4f, Q(p;α,β): %0.4f', p, alpha, beta, + + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/stats/base/dists/invgamma/quantile.h" +``` + +#### stdlib_base_dists_invgamma_quantile( p, alpha, beta ) + +Evaluates the quantile function for an inverse gamma distribution. + +```c +double y = stdlib_base_dists_invgamma_quantile( 0.8, 2.0, 1.0 ); +// returns ~1.213 +``` + +The function accepts the following arguments: + +- **p**: `[in] double` input value. +- **alpha**: `[in] double` shape parameter. +- **beta**: `[in] double` scale parameter. + +```c +double stdlib_base_dists_invgamma_quantile( const double p, const double alpha, const double beta ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/stats/base/dists/invgamma/quantile.h" +#include +#include + +int main( void ) { + double alpha; + double beta; + double p; + double y; + int i; + + for ( i = 0; i < 10; i++ ) { + p = (double)rand() / ( (double)RAND_MAX + 1.0 ); + alpha = ( (double)rand() / ( (double)RAND_MAX + 1.0 ) ) * 10.0; + beta = ( (double)rand() / ( (double)RAND_MAX + 1.0 ) ) * 10.0; + y = stdlib_base_dists_invgamma_quantile( p, alpha, beta ); + printf( "p: %lf, alpha: %lf, beta: %lf, y: %lf\n", p, alpha, beta, y ); + } +} +``` + +
+ + + +
+ + +