@stdlib/stats-base-dists-gamma-logpdf 0.0.7 → 0.2.0

package/NOTICE

@@ -1 +1 @@
-Copyright (c) 2016-2022 The Stdlib Authors.
+Copyright (c) 2016-2024 The Stdlib Authors.

package/README.md

@@ -18,6 +18,17 @@ limitations under the License.
 
 -->
 
+
+<details>
+  <summary>
+    About stdlib...
+  </summary>
+  <p>We believe in a future in which the web is a preferred environment for numerical computation. To help realize this future, we've built stdlib. stdlib is a standard library, with an emphasis on numerical and scientific computation, written in JavaScript (and C) for execution in browsers and in Node.js.</p>
+  <p>The library is fully decomposable, being architected in such a way that you can swap out and mix and match APIs and functionality to cater to your exact preferences and use cases.</p>
+  <p>When you use stdlib, you can be absolutely certain that you are using the most thorough, rigorous, well-written, studied, documented, tested, measured, and high-quality code out there.</p>
+  <p>To join us in bringing numerical computing to the web, get started by checking us out on <a href="https://github.com/stdlib-js/stdlib">GitHub</a>, and please consider <a href="https://opencollective.com/stdlib">financially supporting stdlib</a>. We greatly appreciate your continued support!</p>
+</details>
+
 # Logarithm of Probability Density Function
 
 [![NPM version][npm-image]][npm-url] [![Build Status][test-image]][test-url] [![Coverage Status][coverage-image]][coverage-url] <!-- [![dependencies][dependencies-image]][dependencies-url] -->
@@ -197,7 +208,7 @@ See [LICENSE][stdlib-license].
 
 ## Copyright
 
-Copyright &copy; 2016-2022. The Stdlib [Authors][stdlib-authors].
+Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
 
 </section>
 
@@ -210,8 +221,8 @@ Copyright &copy; 2016-2022. The Stdlib [Authors][stdlib-authors].
 [npm-image]: http://img.shields.io/npm/v/@stdlib/stats-base-dists-gamma-logpdf.svg
 [npm-url]: https://npmjs.org/package/@stdlib/stats-base-dists-gamma-logpdf
 
-[test-image]: https://github.com/stdlib-js/stats-base-dists-gamma-logpdf/actions/workflows/test.yml/badge.svg?branch=v0.0.7
-[test-url]: https://github.com/stdlib-js/stats-base-dists-gamma-logpdf/actions/workflows/test.yml?query=branch:v0.0.7
+[test-image]: https://github.com/stdlib-js/stats-base-dists-gamma-logpdf/actions/workflows/test.yml/badge.svg?branch=v0.2.0
+[test-url]: https://github.com/stdlib-js/stats-base-dists-gamma-logpdf/actions/workflows/test.yml?query=branch:v0.2.0
 
 [coverage-image]: https://img.shields.io/codecov/c/github/stdlib-js/stats-base-dists-gamma-logpdf/main.svg
 [coverage-url]: https://codecov.io/github/stdlib-js/stats-base-dists-gamma-logpdf?branch=main
@@ -224,7 +235,7 @@ Copyright &copy; 2016-2022. The Stdlib [Authors][stdlib-authors].
 -->
 
 [chat-image]: https://img.shields.io/gitter/room/stdlib-js/stdlib.svg
-[chat-url]: https://gitter.im/stdlib-js/stdlib/
+[chat-url]: https://app.gitter.im/#/room/#stdlib-js_stdlib:gitter.im
 
 [stdlib]: https://github.com/stdlib-js/stdlib
 
@@ -234,8 +245,11 @@ Copyright &copy; 2016-2022. The Stdlib [Authors][stdlib-authors].
 [es-module]: https://developer.mozilla.org/en-US/docs/Web/JavaScript/Guide/Modules
 
 [deno-url]: https://github.com/stdlib-js/stats-base-dists-gamma-logpdf/tree/deno
+[deno-readme]: https://github.com/stdlib-js/stats-base-dists-gamma-logpdf/blob/deno/README.md
 [umd-url]: https://github.com/stdlib-js/stats-base-dists-gamma-logpdf/tree/umd
+[umd-readme]: https://github.com/stdlib-js/stats-base-dists-gamma-logpdf/blob/umd/README.md
 [esm-url]: https://github.com/stdlib-js/stats-base-dists-gamma-logpdf/tree/esm
+[esm-readme]: https://github.com/stdlib-js/stats-base-dists-gamma-logpdf/blob/esm/README.md
 [branches-url]: https://github.com/stdlib-js/stats-base-dists-gamma-logpdf/blob/main/branches.md
 
 [stdlib-license]: https://raw.githubusercontent.com/stdlib-js/stats-base-dists-gamma-logpdf/main/LICENSE

package/SECURITY.md

@@ -0,0 +1,5 @@
+# Security
+
+> Policy for reporting security vulnerabilities.
+
+See the security policy [in the main project repository](https://github.com/stdlib-js/stdlib/security).

package/dist/index.d.ts

@@ -0,0 +1,3 @@
+/// <reference path="../docs/types/index.d.ts" />
+import logPDF from '../docs/types/index';
+export = logPDF;

package/dist/index.js

@@ -0,0 +1,13 @@
+"use strict";var o=function(r,e){return function(){return e||r((e={exports:{}}).exports,e),e.exports}};var y=o(function(vr,P){
+var T=require('@stdlib/math-base-special-gamma-lanczos-sum-expg-scaled/dist'),j=require('@stdlib/math-base-special-gammaln/dist'),k=require('@stdlib/math-base-special-gamma/dist'),B=require('@stdlib/math-base-special-log1p/dist'),C=require('@stdlib/math-base-special-sqrt/dist'),H=require('@stdlib/math-base-special-abs/dist'),t=require('@stdlib/math-base-special-exp/dist'),s=require('@stdlib/math-base-special-pow/dist'),c=require('@stdlib/math-base-special-max/dist'),l=require('@stdlib/math-base-special-min/dist'),I=require('@stdlib/math-base-special-ln/dist'),N=require('@stdlib/constants-float64-max-ln/dist'),m=require('@stdlib/constants-float64-min-ln/dist'),g=require('@stdlib/constants-float64-gamma-lanczos-g/dist'),J=require('@stdlib/constants-float64-e/dist');function K(r,e){var i,n,v,a,u,q,f;return v=r+g-.5,f=(e-r-g+.5)/v,r<1?e<=m?t(r*I(e)-e-j(r)):s(e,r)*t(-e)/k(r):(H(f*f*r)<=100&&r>150?(i=r*(B(f)-f)+e*(.5-g)/v,i=t(i)):(a=r*I(e/v),u=r-e,l(a,u)<=m||c(a,u)>=N?(n=u/r,l(a,u)/2>m&&c(a,u)/2<N?(q=s(e/v,r/2)*t(u/2),i=q*q):l(a,u)/4>m&&c(a,u)/4<N&&e>r?(q=s(e/v,r/4)*t(u/4),i=q*q,i*=i):n>m&&n<N?i=s(e*t(n)/v,r):i=t(a+u)):i=s(e/v,r)*t(u)),i*=C(v/J)/T(r),i)}P.exports=K
+});var d=o(function(ar,D){
+var Q=require('@stdlib/math-base-special-gammaln/dist'),U=require('@stdlib/math-base-special-exp/dist'),G=require('@stdlib/math-base-special-ln/dist'),A=require('@stdlib/constants-float64-pinf/dist'),V=require('@stdlib/constants-float64-max/dist'),W=y();function Y(r,e){var i;return r<=0?NaN:e<0?NaN:e===0?r>1?0:r===1?1:A:(i=W(r,e),e<1&&V*e<i?A:(i===0?(i=r*G(e)-e-Q(r)-G(e),i=U(i)):i/=e,i))}D.exports=Y
+});var O=o(function(tr,E){
+var p=require('@stdlib/math-base-assert-is-nan/dist'),L=require('@stdlib/math-base-special-ln/dist'),M=require('@stdlib/constants-float64-pinf/dist'),_=require('@stdlib/constants-float64-ninf/dist'),Z=d();function $(r,e,i){return p(r)||p(e)||p(i)||e<0||i<=0?NaN:r<0||r===M?_:e===0?r===0?M:_:L(Z(e,r*i))+L(i)}E.exports=$
+});var w=o(function(qr,X){
+var x=require('@stdlib/utils-constant-function/dist'),b=require('@stdlib/stats-base-dists-degenerate-logpdf/dist').factory,F=require('@stdlib/math-base-assert-is-nan/dist'),S=require('@stdlib/math-base-special-ln/dist'),h=require('@stdlib/constants-float64-ninf/dist'),z=require('@stdlib/constants-float64-pinf/dist'),rr=d();function er(r,e){if(F(r)||F(e)||r<0||e<=0)return x(NaN);if(r===0)return b(0);return i;function i(n){return F(n)?NaN:n<0||n===z?h:S(rr(r,n*e))+S(e)}}X.exports=er
+});var ir=require('@stdlib/utils-define-nonenumerable-read-only-property/dist'),R=O(),ur=w();ir(R,"factory",ur);module.exports=R;
+/** @license Apache-2.0 */
+/** @license Apache-2.0 */
+/** @license Apache-2.0 */
+//# sourceMappingURL=index.js.map

package/dist/index.js.map

@@ -0,0 +1,7 @@
+{
+  "version": 3,
+  "sources": ["../lib/regularised_gamma_prefix.js", "../lib/gamma_p_derivative.js", "../lib/main.js", "../lib/factory.js", "../lib/index.js"],
+  "sourcesContent": ["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*\n*\n* ## Notice\n*\n* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_64_0/boost/math/special_functions/gamma.hpp}. The implementation has been modified for JavaScript.\n*\n* ```text\n* Copyright John Maddock 2006-7, 2013-14.\n* Copyright Paul A. Bristow 2007, 2013-14.\n* Copyright Nikhar Agrawal 2013-14.\n* Copyright Christopher Kormanyos 2013-14.\n*\n* Use, modification and distribution are subject to the\n* Boost Software License, Version 1.0. (See accompanying file\n* LICENSE or copy at http://www.boost.org/LICENSE_1_0.txt)\n* ```\n*/\n\n'use strict';\n\n// MODULES //\n\nvar lanczosSumExpGScaled = require( '@stdlib/math-base-special-gamma-lanczos-sum-expg-scaled' );\nvar gammaln = require( '@stdlib/math-base-special-gammaln' );\nvar gamma = require( '@stdlib/math-base-special-gamma' );\nvar log1p = require( '@stdlib/math-base-special-log1p' );\nvar sqrt = require( '@stdlib/math-base-special-sqrt' );\nvar abs = require( '@stdlib/math-base-special-abs' );\nvar exp = require( '@stdlib/math-base-special-exp' );\nvar pow = require( '@stdlib/math-base-special-pow' );\nvar max = require( '@stdlib/math-base-special-max' );\nvar min = require( '@stdlib/math-base-special-min' );\nvar ln = require( '@stdlib/math-base-special-ln' );\nvar MAX_LN = require( '@stdlib/constants-float64-max-ln' );\nvar MIN_LN = require( '@stdlib/constants-float64-min-ln' );\nvar G = require( '@stdlib/constants-float64-gamma-lanczos-g' );\nvar E = require( '@stdlib/constants-float64-e' );\n\n\n// MAIN //\n\n/**\n* Computes `(z^a)*(e^-z) / gamma(a)`.\n*\n* @private\n* @param {number} a - input value\n* @param {number} z - input value\n* @returns {number} function value\n*/\nfunction regularisedGammaPrefix( a, z ) {\n\tvar prefix;\n\tvar amza;\n\tvar agh;\n\tvar alz;\n\tvar amz;\n\tvar sq;\n\tvar d;\n\n\tagh = a + G - 0.5;\n\td = ( (z - a) - G + 0.5 ) / agh;\n\tif ( a < 1.0 ) {\n\t\t// Treat a < 1 as a special case because our Lanczos approximations are optimized against the factorials with a > 1, and for high precision types very small values of `a` can give rather erroneous results for gamma:\n\t\tif ( z <= MIN_LN ) {\n\t\t\t// Use logs, so should be free of cancellation errors:\n\t\t\treturn exp( ( a * ln(z) ) - z - gammaln( a ) );\n\t\t}\n\t\t// No danger of overflow as gamma(a) < 1/a for small a, so direct calculation:\n\t\treturn pow( z, a ) * exp( -z ) / gamma( a );\n\t}\n\tif ( abs(d*d*a) <= 100.0 && a > 150.0 ) {\n\t\t// Special case for large a and a ~ z:\n\t\tprefix = ( a * ( log1p( d ) - d ) ) + ( z * ( 0.5-G ) / agh );\n\t\tprefix = exp( prefix );\n\t}\n\telse {\n\t\t// General case. Direct computation is most accurate, but use various fallbacks for different parts of the problem domain:\n\t\talz = a * ln(z / agh);\n\t\tamz = a - z;\n\t\tif (\n\t\t\tmin(alz, amz) <= MIN_LN ||\n\t\t\tmax(alz, amz) >= MAX_LN\n\t\t) {\n\t\t\tamza = amz / a;\n\t\t\tif (\n\t\t\t\tmin(alz, amz)/2.0 > MIN_LN &&\n\t\t\t\tmax(alz, amz)/2.0 < MAX_LN\n\t\t\t) {\n\t\t\t\t// Compute square root of the result and then square it:\n\t\t\t\tsq = pow( z / agh, a / 2.0 ) * exp( amz / 2.0 );\n\t\t\t\tprefix = sq * sq;\n\t\t\t}\n\t\t\telse if (\n\t\t\t\tmin(alz, amz)/4.0 > MIN_LN &&\n\t\t\t\tmax(alz, amz)/4.0 < MAX_LN &&\n\t\t\t\tz > a\n\t\t\t) {\n\t\t\t\t// Compute the 4th root of the result then square it twice:\n\t\t\t\tsq = pow( z / agh, a / 4.0 ) * exp( amz / 4.0 );\n\t\t\t\tprefix = sq * sq;\n\t\t\t\tprefix *= prefix;\n\t\t\t}\n\t\t\telse if (\n\t\t\t\tamza > MIN_LN &&\n\t\t\t\tamza < MAX_LN\n\t\t\t) {\n\t\t\t\tprefix = pow( (z * exp(amza)) / agh, a );\n\t\t\t}\n\t\t\telse {\n\t\t\t\tprefix = exp( alz + amz );\n\t\t\t}\n\t\t}\n\t\telse\n\t\t{\n\t\t\tprefix = pow( z / agh, a ) * exp( amz );\n\t\t}\n\t}\n\tprefix *= sqrt( agh / E ) / lanczosSumExpGScaled( a );\n\treturn prefix;\n}\n\n\n// EXPORTS //\n\nmodule.exports = regularisedGammaPrefix;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*\n*\n* ## Notice\n*\n* The original C++ code and copyright notice are from the [Boost library]{@link http://www.boost.org/doc/libs/1_58_0/boost/math/special_functions/gamma.hpp}. The implementation has been modified for JavaScript.\n*\n* ```text\n* Copyright John Maddock 2006-7, 2013-14.\n* Copyright Paul A. Bristow 2007, 2013-14.\n* Copyright Nikhar Agrawal 2013-14.\n* Copyright Christopher Kormanyos 2013-14.\n*\n* Use, modification and distribution are subject to the\n* Boost Software License, Version 1.0. (See accompanying file\n* LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)\n* ```\n*/\n\n'use strict';\n\n// MODULES //\n\nvar gammaln = require( '@stdlib/math-base-special-gammaln' );\nvar exp = require( '@stdlib/math-base-special-exp' );\nvar ln = require( '@stdlib/math-base-special-ln' );\nvar PINF = require( '@stdlib/constants-float64-pinf' );\nvar FLOAT64_MAX = require( '@stdlib/constants-float64-max' );\nvar regularisedGammaPrefix = require( './regularised_gamma_prefix.js' );\n\n\n// MAIN //\n\n/**\n* Calculates the partial derivative with respect to x of the incomplete gamma function.\n*\n* @private\n* @param {PositiveNumber} a - function parameter\n* @param {NonNegativeNumber} x - function parameter\n* @returns {number} function value\n*/\nfunction gammaPDerivative( a, x ) {\n\tvar f1;\n\tif ( a <= 0.0 ) {\n\t\treturn NaN;\n\t}\n\tif ( x < 0.0 ) {\n\t\treturn NaN;\n\t}\n\tif ( x === 0.0 ) {\n\t\tif ( a > 1.0 ) {\n\t\t\treturn 0.0;\n\t\t}\n\t\treturn ( a === 1.0 ) ? 1.0 : PINF;\n\t}\n\tf1 = regularisedGammaPrefix( a, x );\n\tif ( x < 1.0 && ( FLOAT64_MAX * x < f1 ) ) {\n\t\treturn PINF;\n\t}\n\tif ( f1 === 0.0 ) {\n\t\t// Underflow in calculation, use logs instead:\n\t\tf1 = (a * ln( x )) - x - gammaln( a ) - ln( x );\n\t\tf1 = exp( f1 );\n\t} else {\n\t\tf1 /= x;\n\t}\n\treturn f1;\n}\n\n\n// EXPORTS //\n\nmodule.exports = gammaPDerivative;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar isnan = require( '@stdlib/math-base-assert-is-nan' );\nvar ln = require( '@stdlib/math-base-special-ln' );\nvar PINF = require( '@stdlib/constants-float64-pinf' );\nvar NINF = require( '@stdlib/constants-float64-ninf' );\nvar gammaDeriv = require( './gamma_p_derivative.js' );\n\n\n// MAIN //\n\n/**\n* Evaluates the logarithm of the probability density function (PDF) for a gamma distribution with shape parameter `alpha` and rate parameter `beta` at a value `x`.\n*\n* @param {number} x - input value\n* @param {NonNegativeNumber} alpha - shape parameter\n* @param {PositiveNumber} beta - rate parameter\n* @returns {number} evaluated logPDF\n*\n* @example\n* var y = logpdf( 2.0, 0.5, 1.0 );\n* // returns ~-2.919\n*\n* @example\n* var y = logpdf( 0.1, 1.0, 1.0 );\n* // returns ~-0.1\n*\n* @example\n* var y = logpdf( -1.0, 4.0, 2.0 );\n* // returns -Infinity\n*\n* @example\n* var y = logpdf( NaN, 0.6, 1.0 );\n* // returns NaN\n*\n* @example\n* var y = logpdf( 0.0, NaN, 1.0 );\n* // returns NaN\n*\n* @example\n* var y = logpdf( 0.0, 1.0, NaN );\n* // returns NaN\n*\n* @example\n* // Negative shape parameter:\n* var y = logpdf( 2.0, -1.0, 1.0 );\n* // returns NaN\n*\n* @example\n* // Negative rate parameter:\n* var y = logpdf( 2.0, 1.0, -1.0 );\n* // returns NaN\n*/\nfunction logpdf( x, alpha, beta ) {\n\tif (\n\t\tisnan( x ) ||\n\t\tisnan( alpha ) ||\n\t\tisnan( beta ) ||\n\t\talpha < 0.0 ||\n\t\tbeta <= 0.0\n\t) {\n\t\treturn NaN;\n\t}\n\tif ( x < 0.0 || x === PINF ) {\n\t\treturn NINF;\n\t}\n\tif ( alpha === 0.0 ) {\n\t\t// Point mass at 0...\n\t\treturn ( x === 0.0 ) ? PINF : NINF;\n\t}\n\treturn ln( gammaDeriv( alpha, x * beta ) ) + ln( beta );\n}\n\n\n// EXPORTS //\n\nmodule.exports = logpdf;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar constantFunction = require( '@stdlib/utils-constant-function' );\nvar degenerate = require( '@stdlib/stats-base-dists-degenerate-logpdf' ).factory;\nvar isnan = require( '@stdlib/math-base-assert-is-nan' );\nvar ln = require( '@stdlib/math-base-special-ln' );\nvar NINF = require( '@stdlib/constants-float64-ninf' );\nvar PINF = require( '@stdlib/constants-float64-pinf' );\nvar gammaDeriv = require( './gamma_p_derivative.js' );\n\n\n// MAIN //\n\n/**\n* Returns a function for evaluating the logarithm of the probability density function (PDF) for a gamma distribution with shape parameter `alpha` and rate parameter `beta`.\n*\n* @param {NonNegativeNumber} alpha - shape parameter\n* @param {PositiveNumber} beta - rate parameter\n* @returns {Function} logPDF\n*\n* @example\n* var logpdf = factory( 3.0, 1.5 );\n*\n* var y = logpdf( 1.0 );\n* // returns ~-0.977\n*\n* y = logpdf( 4.0 );\n* // returns ~-2.704\n*/\nfunction factory( alpha, beta ) {\n\tif (\n\t\tisnan( alpha ) ||\n\t\tisnan( beta ) ||\n\t\talpha < 0.0 ||\n\t\tbeta <= 0.0\n\t) {\n\t\treturn constantFunction( NaN );\n\t}\n\tif ( alpha === 0.0 ) {\n\t\treturn degenerate( 0.0 );\n\t}\n\treturn logpdf;\n\n\t/**\n\t* Evaluates the logarithm of the probability density function (PDF) for a gamma distribution.\n\t*\n\t* @private\n\t* @param {number} x - input value\n\t* @returns {number} evaluated logPDF\n\t*\n\t* @example\n\t* var y = logpdf( -1.2 );\n\t* // returns <number>\n\t*/\n\tfunction logpdf( x ) {\n\t\tif ( isnan( x ) ) {\n\t\t\treturn NaN;\n\t\t}\n\t\tif ( x < 0.0 || x === PINF ) {\n\t\t\treturn NINF;\n\t\t}\n\t\treturn ln( gammaDeriv( alpha, x * beta ) ) + ln( beta );\n\t}\n}\n\n\n// EXPORTS //\n\nmodule.exports = factory;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n/**\n* Gamma distribution logarithm of probability density function (PDF).\n*\n* @module @stdlib/stats-base-dists-gamma-logpdf\n*\n* @example\n* var logpdf = require( '@stdlib/stats-base-dists-gamma-logpdf' );\n*\n* var y = logpdf( 2.0, 0.5, 1.0 );\n* // returns ~-2.919\n*\n* @example\n* var factory = require( '@stdlib/stats-base-dists-gamma-logpdf' ).factory;\n*\n* var logpdf = factory( 6.0, 7.0 );\n* var y = logpdf( 2.0 );\n* // returns ~-3.646\n*/\n\n// MODULES //\n\nvar setReadOnly = require( '@stdlib/utils-define-nonenumerable-read-only-property' );\nvar main = require( './main.js' );\nvar factory = require( './factory.js' );\n\n\n// MAIN //\n\nsetReadOnly( main, 'factory', factory );\n\n\n// EXPORTS //\n\nmodule.exports = main;\n"],
+  "mappings": "uGAAA,IAAAA,EAAAC,EAAA,SAAAC,GAAAC,EAAA,cAsCA,IAAIC,EAAuB,QAAS,yDAA0D,EAC1FC,EAAU,QAAS,mCAAoC,EACvDC,EAAQ,QAAS,iCAAkC,EACnDC,EAAQ,QAAS,iCAAkC,EACnDC,EAAO,QAAS,gCAAiC,EACjDC,EAAM,QAAS,+BAAgC,EAC/CC,EAAM,QAAS,+BAAgC,EAC/CC,EAAM,QAAS,+BAAgC,EAC/CC,EAAM,QAAS,+BAAgC,EAC/CC,EAAM,QAAS,+BAAgC,EAC/CC,EAAK,QAAS,8BAA+B,EAC7CC,EAAS,QAAS,kCAAmC,EACrDC,EAAS,QAAS,kCAAmC,EACrDC,EAAI,QAAS,2CAA4C,EACzDC,EAAI,QAAS,6BAA8B,EAa/C,SAASC,EAAwBC,EAAGC,EAAI,CACvC,IAAIC,EACAC,EACAC,EACAC,EACAC,EACAC,EACAC,EAIJ,OAFAJ,EAAMJ,EAAIH,EAAI,GACdW,GAAOP,EAAID,EAAKH,EAAI,IAAQO,EACvBJ,EAAI,EAEHC,GAAKL,EAEFN,EAAOU,EAAIN,EAAGO,CAAC,EAAMA,EAAIhB,EAASe,CAAE,CAAE,EAGvCT,EAAKU,EAAGD,CAAE,EAAIV,EAAK,CAACW,CAAE,EAAIf,EAAOc,CAAE,GAEtCX,EAAImB,EAAEA,EAAER,CAAC,GAAK,KAASA,EAAI,KAE/BE,EAAWF,GAAMb,EAAOqB,CAAE,EAAIA,GAAUP,GAAM,GAAIJ,GAAMO,EACxDF,EAASZ,EAAKY,CAAO,IAIrBG,EAAML,EAAIN,EAAGO,EAAIG,CAAG,EACpBE,EAAMN,EAAIC,EAETR,EAAIY,EAAKC,CAAG,GAAKV,GACjBJ,EAAIa,EAAKC,CAAG,GAAKX,GAEjBQ,EAAOG,EAAMN,EAEZP,EAAIY,EAAKC,CAAG,EAAE,EAAMV,GACpBJ,EAAIa,EAAKC,CAAG,EAAE,EAAMX,GAGpBY,EAAKhB,EAAKU,EAAIG,EAAKJ,EAAI,CAAI,EAAIV,EAAKgB,EAAM,CAAI,EAC9CJ,EAASK,EAAKA,GAGdd,EAAIY,EAAKC,CAAG,EAAE,EAAMV,GACpBJ,EAAIa,EAAKC,CAAG,EAAE,EAAMX,GACpBM,EAAID,GAGJO,EAAKhB,EAAKU,EAAIG,EAAKJ,EAAI,CAAI,EAAIV,EAAKgB,EAAM,CAAI,EAC9CJ,EAASK,EAAKA,EACdL,GAAUA,GAGVC,EAAOP,GACPO,EAAOR,EAEPO,EAASX,EAAMU,EAAIX,EAAIa,CAAI,EAAKC,EAAKJ,CAAE,EAGvCE,EAASZ,EAAKe,EAAMC,CAAI,GAKzBJ,EAASX,EAAKU,EAAIG,EAAKJ,CAAE,EAAIV,EAAKgB,CAAI,GAGxCJ,GAAUd,EAAMgB,EAAMN,CAAE,EAAId,EAAsBgB,CAAE,EAC7CE,EACR,CAKAnB,EAAO,QAAUgB,IC3IjB,IAAAU,EAAAC,EAAA,SAAAC,GAAAC,EAAA,cAsCA,IAAIC,EAAU,QAAS,mCAAoC,EACvDC,EAAM,QAAS,+BAAgC,EAC/CC,EAAK,QAAS,8BAA+B,EAC7CC,EAAO,QAAS,gCAAiC,EACjDC,EAAc,QAAS,+BAAgC,EACvDC,EAAyB,IAa7B,SAASC,EAAkBC,EAAGC,EAAI,CACjC,IAAIC,EACJ,OAAKF,GAAK,EACF,IAEHC,EAAI,EACD,IAEHA,IAAM,EACLD,EAAI,EACD,EAECA,IAAM,EAAQ,EAAMJ,GAE9BM,EAAKJ,EAAwBE,EAAGC,CAAE,EAC7BA,EAAI,GAASJ,EAAcI,EAAIC,EAC5BN,GAEHM,IAAO,GAEXA,EAAMF,EAAIL,EAAIM,CAAE,EAAKA,EAAIR,EAASO,CAAE,EAAIL,EAAIM,CAAE,EAC9CC,EAAKR,EAAKQ,CAAG,GAEbA,GAAMD,EAEAC,GACR,CAKAV,EAAO,QAAUO,ICvFjB,IAAAI,EAAAC,EAAA,SAAAC,GAAAC,EAAA,cAsBA,IAAIC,EAAQ,QAAS,iCAAkC,EACnDC,EAAK,QAAS,8BAA+B,EAC7CC,EAAO,QAAS,gCAAiC,EACjDC,EAAO,QAAS,gCAAiC,EACjDC,EAAa,IA+CjB,SAASC,EAAQC,EAAGC,EAAOC,EAAO,CACjC,OACCR,EAAOM,CAAE,GACTN,EAAOO,CAAM,GACbP,EAAOQ,CAAK,GACZD,EAAQ,GACRC,GAAQ,EAED,IAEHF,EAAI,GAAOA,IAAMJ,EACdC,EAEHI,IAAU,EAELD,IAAM,EAAQJ,EAAOC,EAExBF,EAAIG,EAAYG,EAAOD,EAAIE,CAAK,CAAE,EAAIP,EAAIO,CAAK,CACvD,CAKAT,EAAO,QAAUM,IChGjB,IAAAI,EAAAC,EAAA,SAAAC,GAAAC,EAAA,cAsBA,IAAIC,EAAmB,QAAS,iCAAkC,EAC9DC,EAAa,QAAS,4CAA6C,EAAE,QACrEC,EAAQ,QAAS,iCAAkC,EACnDC,EAAK,QAAS,8BAA+B,EAC7CC,EAAO,QAAS,gCAAiC,EACjDC,EAAO,QAAS,gCAAiC,EACjDC,GAAa,IAqBjB,SAASC,GAASC,EAAOC,EAAO,CAC/B,GACCP,EAAOM,CAAM,GACbN,EAAOO,CAAK,GACZD,EAAQ,GACRC,GAAQ,EAER,OAAOT,EAAkB,GAAI,EAE9B,GAAKQ,IAAU,EACd,OAAOP,EAAY,CAAI,EAExB,OAAOS,EAaP,SAASA,EAAQC,EAAI,CACpB,OAAKT,EAAOS,CAAE,EACN,IAEHA,EAAI,GAAOA,IAAMN,EACdD,EAEDD,EAAIG,GAAYE,EAAOG,EAAIF,CAAK,CAAE,EAAIN,EAAIM,CAAK,CACvD,CACD,CAKAV,EAAO,QAAUQ,KC/CjB,IAAIK,GAAc,QAAS,uDAAwD,EAC/EC,EAAO,IACPC,GAAU,IAKdF,GAAaC,EAAM,UAAWC,EAAQ,EAKtC,OAAO,QAAUD",
+  "names": ["require_regularised_gamma_prefix", "__commonJSMin", "exports", "module", "lanczosSumExpGScaled", "gammaln", "gamma", "log1p", "sqrt", "abs", "exp", "pow", "max", "min", "ln", "MAX_LN", "MIN_LN", "G", "E", "regularisedGammaPrefix", "a", "z", "prefix", "amza", "agh", "alz", "amz", "sq", "d", "require_gamma_p_derivative", "__commonJSMin", "exports", "module", "gammaln", "exp", "ln", "PINF", "FLOAT64_MAX", "regularisedGammaPrefix", "gammaPDerivative", "a", "x", "f1", "require_main", "__commonJSMin", "exports", "module", "isnan", "ln", "PINF", "NINF", "gammaDeriv", "logpdf", "x", "alpha", "beta", "require_factory", "__commonJSMin", "exports", "module", "constantFunction", "degenerate", "isnan", "ln", "NINF", "PINF", "gammaDeriv", "factory", "alpha", "beta", "logpdf", "x", "setReadOnly", "main", "factory"]
+}

package/docs/types/index.d.ts

@@ -16,7 +16,7 @@
 * limitations under the License.
 */
 
-// TypeScript Version: 2.0
+// TypeScript Version: 4.1
 
 /**
 * Evaluates the logarithm of the probability density function (PDF) for a gamma distribution.

package/lib/index.js

@@ -40,15 +40,15 @@
 // MODULES //
 
 var setReadOnly = require( '@stdlib/utils-define-nonenumerable-read-only-property' );
-var logpdf = require( './logpdf.js' );
+var main = require( './main.js' );
 var factory = require( './factory.js' );
 
 
 // MAIN //
 
-setReadOnly( logpdf, 'factory', factory );
+setReadOnly( main, 'factory', factory );
 
 
 // EXPORTS //
 
-module.exports = logpdf;
+module.exports = main;

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@stdlib/stats-base-dists-gamma-logpdf",
-  "version": "0.0.7",
+  "version": "0.2.0",
   "description": "Gamma distribution logarithm of probability density function (PDF).",
   "license": "Apache-2.0",
   "author": {
@@ -37,36 +37,37 @@
     "url": "https://github.com/stdlib-js/stdlib/issues"
   },
   "dependencies": {
-    "@stdlib/constants-float64-e": "^0.0.x",
-    "@stdlib/constants-float64-gamma-lanczos-g": "^0.0.x",
-    "@stdlib/constants-float64-max": "^0.0.x",
-    "@stdlib/constants-float64-max-ln": "^0.0.x",
-    "@stdlib/constants-float64-min-ln": "^0.0.x",
-    "@stdlib/constants-float64-ninf": "^0.0.x",
-    "@stdlib/constants-float64-pinf": "^0.0.x",
-    "@stdlib/math-base-assert-is-nan": "^0.0.x",
-    "@stdlib/math-base-special-abs": "^0.0.x",
-    "@stdlib/math-base-special-exp": "^0.0.x",
-    "@stdlib/math-base-special-gamma": "^0.0.x",
-    "@stdlib/math-base-special-gamma-lanczos-sum-expg-scaled": "^0.0.x",
-    "@stdlib/math-base-special-gammaln": "^0.0.x",
-    "@stdlib/math-base-special-ln": "^0.0.x",
-    "@stdlib/math-base-special-log1p": "^0.0.x",
-    "@stdlib/math-base-special-max": "^0.0.x",
-    "@stdlib/math-base-special-min": "^0.0.x",
-    "@stdlib/math-base-special-pow": "^0.0.x",
-    "@stdlib/math-base-special-sqrt": "^0.0.x",
-    "@stdlib/stats-base-dists-degenerate-logpdf": "^0.0.x",
-    "@stdlib/utils-constant-function": "^0.0.x",
-    "@stdlib/utils-define-nonenumerable-read-only-property": "^0.0.x"
+    "@stdlib/constants-float64-e": "^0.2.0",
+    "@stdlib/constants-float64-gamma-lanczos-g": "^0.2.0",
+    "@stdlib/constants-float64-max": "^0.2.0",
+    "@stdlib/constants-float64-max-ln": "^0.2.0",
+    "@stdlib/constants-float64-min-ln": "^0.2.0",
+    "@stdlib/constants-float64-ninf": "^0.2.0",
+    "@stdlib/constants-float64-pinf": "^0.2.0",
+    "@stdlib/math-base-assert-is-nan": "^0.2.0",
+    "@stdlib/math-base-special-abs": "^0.2.0",
+    "@stdlib/math-base-special-exp": "^0.2.0",
+    "@stdlib/math-base-special-gamma": "^0.2.0",
+    "@stdlib/math-base-special-gamma-lanczos-sum-expg-scaled": "^0.2.0",
+    "@stdlib/math-base-special-gammaln": "^0.2.0",
+    "@stdlib/math-base-special-ln": "^0.2.0",
+    "@stdlib/math-base-special-log1p": "^0.2.0",
+    "@stdlib/math-base-special-max": "^0.2.0",
+    "@stdlib/math-base-special-min": "^0.2.0",
+    "@stdlib/math-base-special-pow": "^0.2.0",
+    "@stdlib/math-base-special-sqrt": "^0.2.0",
+    "@stdlib/stats-base-dists-degenerate-logpdf": "^0.2.0",
+    "@stdlib/utils-constant-function": "^0.2.0",
+    "@stdlib/utils-define-nonenumerable-read-only-property": "^0.2.0"
   },
   "devDependencies": {
-    "@stdlib/bench": "^0.0.x",
-    "@stdlib/constants-float64-eps": "^0.0.x",
-    "@stdlib/random-base-randu": "^0.0.x",
+    "@stdlib/constants-float64-eps": "^0.2.0",
+    "@stdlib/random-base-randu": "^0.1.0",
     "tape": "git+https://github.com/kgryte/tape.git#fix/globby",
     "istanbul": "^0.4.1",
-    "tap-spec": "5.x.x"
+    "tap-min": "git+https://github.com/Planeshifter/tap-min.git",
+    "@stdlib/bench-harness": "^0.2.0",
+    "@stdlib/bench": "^0.3.1"
   },
   "engines": {
     "node": ">=0.10.0",
@@ -99,7 +100,7 @@
     "continuous"
   ],
   "funding": {
-    "type": "patreon",
-    "url": "https://www.patreon.com/athan"
+    "type": "opencollective",
+    "url": "https://opencollective.com/stdlib"
   }
 }

package/docs/repl.txt

@@ -1,83 +0,0 @@
-
-{{alias}}( x, α, β )
-    Evaluates the logarithm of the probability density function (PDF) for a
-    gamma distribution with shape parameter `α` and rate parameter `β` at a
-    value `x`.
-
-    If `α < 0` or `β <= 0`, the function returns `NaN`.
-
-    If provided `NaN` as any argument, the function returns `NaN`.
-
-    Parameters
-    ----------
-    x: number
-        Input value.
-
-    α: number
-        Shape parameter.
-
-    β: number
-        Rate parameter.
-
-    Returns
-    -------
-    out: number
-        Evaluated logPDF.
-
-    Examples
-    --------
-    > var y = {{alias}}( 2.0, 0.5, 1.0 )
-    ~-2.919
-    > y = {{alias}}( 0.1, 1.0, 1.0 )
-    ~-0.1
-    > y = {{alias}}( -1.0, 4.0, 2.0 )
-    -Infinity
-
-    > y = {{alias}}( NaN, 0.6, 1.0 )
-    NaN
-    > y = {{alias}}( 0.0, NaN, 1.0 )
-    NaN
-    > y = {{alias}}( 0.0, 1.0, NaN )
-    NaN
-
-    // Negative shape parameter:
-    > y = {{alias}}( 2.0, -1.0, 1.0 )
-    NaN
-    // Non-positive rate parameter:
-    > y = {{alias}}( 2.0, 1.0, -1.0 )
-    NaN
-
-    // Degenerate distribution centered at `0.0` when `α = 0.0`:
-    > y = {{alias}}( 2.0, 0.0, 2.0 )
-    -Infinity
-    > y = {{alias}}( 0.0, 0.0, 2.0 )
-    Infinity
-
-
-{{alias}}.factory( α, β )
-    Returns a function for evaluating the logarithm of the probability density
-    function (PDF) of a gamma distribution with shape parameter `α` and rate
-    parameter `β`.
-
-    Parameters
-    ----------
-    α: number
-        Shape parameter.
-
-    β: number
-        Rate parameter.
-
-    Returns
-    -------
-    logpdf: Function
-        Logarithm of probability density function (PDF).
-
-    Examples
-    --------
-    > var mylogPDF = {{alias}}.factory( 6.0, 7.0 );
-    > var y = mylogPDF( 2.0 )
-    ~-3.646
-
-    See Also
-    --------
-

package/docs/types/test.ts

@@ -1,119 +0,0 @@
-/*
-* @license Apache-2.0
-*
-* Copyright (c) 2019 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 logpdf = require( './index' );
-
-
-// TESTS //
-
-// The function returns a number...
-{
-	logpdf( 2, 2, 4 ); // $ExpectType number
-	logpdf( 1, 2, 8 ); // $ExpectType number
-}
-
-// The compiler throws an error if the function is provided values other than three numbers...
-{
-	logpdf( true, 3, 6 ); // $ExpectError
-	logpdf( false, 2, 4 ); // $ExpectError
-	logpdf( '5', 1, 2 ); // $ExpectError
-	logpdf( [], 1, 2 ); // $ExpectError
-	logpdf( {}, 2, 4 ); // $ExpectError
-	logpdf( ( x: number ): number => x, 2, 4 ); // $ExpectError
-
-	logpdf( 9, true, 12 ); // $ExpectError
-	logpdf( 9, false, 12 ); // $ExpectError
-	logpdf( 5, '5', 10 ); // $ExpectError
-	logpdf( 8, [], 16 ); // $ExpectError
-	logpdf( 9, {}, 18 ); // $ExpectError
-	logpdf( 8, ( x: number ): number => x, 16 ); // $ExpectError
-
-	logpdf( 9, 5, true ); // $ExpectError
-	logpdf( 9, 5, false ); // $ExpectError
-	logpdf( 5, 2, '5' ); // $ExpectError
-	logpdf( 8, 4, [] ); // $ExpectError
-	logpdf( 9, 4, {} ); // $ExpectError
-	logpdf( 8, 5, ( x: number ): number => x ); // $ExpectError
-}
-
-// The compiler throws an error if the function is provided an unsupported number of arguments...
-{
-	logpdf(); // $ExpectError
-	logpdf( 2 ); // $ExpectError
-	logpdf( 2, 0 ); // $ExpectError
-	logpdf( 2, 0, 4, 1 ); // $ExpectError
-}
-
-// Attached to main export is a `factory` method which returns a function...
-{
-	logpdf.factory( 3, 4 ); // $ExpectType Unary
-}
-
-// The `factory` method returns a function which returns a number...
-{
-	const fcn = logpdf.factory( 3, 4 );
-	fcn( 2 ); // $ExpectType number
-}
-
-// The compiler throws an error if the function returned by the `factory` method is provided invalid arguments...
-{
-	const fcn = logpdf.factory( 3, 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 = logpdf.factory( 3, 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...
-{
-	logpdf.factory( true, 3 ); // $ExpectError
-	logpdf.factory( false, 2 ); // $ExpectError
-	logpdf.factory( '5', 1 ); // $ExpectError
-	logpdf.factory( [], 1 ); // $ExpectError
-	logpdf.factory( {}, 2 ); // $ExpectError
-	logpdf.factory( ( x: number ): number => x, 2 ); // $ExpectError
-
-	logpdf.factory( 9, true ); // $ExpectError
-	logpdf.factory( 9, false ); // $ExpectError
-	logpdf.factory( 5, '5' ); // $ExpectError
-	logpdf.factory( 8, [] ); // $ExpectError
-	logpdf.factory( 9, {} ); // $ExpectError
-	logpdf.factory( 8, ( x: number ): number => x ); // $ExpectError
-
-	logpdf.factory( [], true ); // $ExpectError
-	logpdf.factory( {}, false ); // $ExpectError
-	logpdf.factory( false, '5' ); // $ExpectError
-	logpdf.factory( {}, [] ); // $ExpectError
-	logpdf.factory( '5', ( x: number ): number => x ); // $ExpectError
-}
-
-// The compiler throws an error if the `factory` method is provided an unsupported number of arguments...
-{
-	logpdf.factory( 0 ); // $ExpectError
-	logpdf.factory( 0, 4, 8 ); // $ExpectError
-}