@stdlib/stats-base-dists-rayleigh-quantile 0.2.2 → 0.3.0
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/NOTICE +1 -1
- package/README.md +108 -16
- package/include/stdlib/stats/base/dists/rayleigh/quantile.h +38 -0
- package/lib/native.js +72 -0
- package/manifest.json +83 -0
- package/package.json +10 -5
- package/src/addon.c +23 -0
- package/src/main.c +51 -0
package/NOTICE
CHANGED
|
@@ -1 +1 @@
|
|
|
1
|
-
Copyright (c) 2016-
|
|
1
|
+
Copyright (c) 2016-2026 The Stdlib Authors.
|
package/README.md
CHANGED
|
@@ -128,7 +128,7 @@ Returns a function for evaluating the [quantile function][quantile-function] of
|
|
|
128
128
|
```javascript
|
|
129
129
|
var myQuantile = quantile.factory( 0.4 );
|
|
130
130
|
|
|
131
|
-
y = myQuantile( 0.4 );
|
|
131
|
+
var y = myQuantile( 0.4 );
|
|
132
132
|
// returns ~0.404
|
|
133
133
|
|
|
134
134
|
y = myQuantile( 1.0 );
|
|
@@ -146,19 +146,107 @@ y = myQuantile( 1.0 );
|
|
|
146
146
|
<!-- eslint no-undef: "error" -->
|
|
147
147
|
|
|
148
148
|
```javascript
|
|
149
|
-
var
|
|
149
|
+
var uniform = require( '@stdlib/random-array-uniform' );
|
|
150
|
+
var logEachMap = require( '@stdlib/console-log-each-map' );
|
|
150
151
|
var quantile = require( '@stdlib/stats-base-dists-rayleigh-quantile' );
|
|
151
152
|
|
|
152
|
-
var
|
|
153
|
-
|
|
154
|
-
|
|
155
|
-
var
|
|
153
|
+
var opts = {
|
|
154
|
+
'dtype': 'float64'
|
|
155
|
+
};
|
|
156
|
+
var p = uniform( 10, 0.0, 1.0, opts );
|
|
157
|
+
var sigma = uniform( 10, 0.0, 10.0, opts );
|
|
156
158
|
|
|
157
|
-
|
|
158
|
-
|
|
159
|
-
|
|
160
|
-
|
|
161
|
-
|
|
159
|
+
logEachMap( 'p: %0.4f, σ: %0.4f, Q(p;σ): %0.4f', p, sigma, quantile );
|
|
160
|
+
```
|
|
161
|
+
|
|
162
|
+
</section>
|
|
163
|
+
|
|
164
|
+
<!-- /.examples -->
|
|
165
|
+
|
|
166
|
+
<!-- C interface documentation. -->
|
|
167
|
+
|
|
168
|
+
* * *
|
|
169
|
+
|
|
170
|
+
<section class="c">
|
|
171
|
+
|
|
172
|
+
## C APIs
|
|
173
|
+
|
|
174
|
+
<!-- Section to include introductory text. Make sure to keep an empty line after the intro `section` element and another before the `/section` close. -->
|
|
175
|
+
|
|
176
|
+
<section class="intro">
|
|
177
|
+
|
|
178
|
+
</section>
|
|
179
|
+
|
|
180
|
+
<!-- /.intro -->
|
|
181
|
+
|
|
182
|
+
<!-- C usage documentation. -->
|
|
183
|
+
|
|
184
|
+
<section class="usage">
|
|
185
|
+
|
|
186
|
+
### Usage
|
|
187
|
+
|
|
188
|
+
```c
|
|
189
|
+
#include "stdlib/stats/base/dists/rayleigh/quantile.h"
|
|
190
|
+
```
|
|
191
|
+
|
|
192
|
+
#### stdlib_base_dists_rayleigh_quantile( p, sigma )
|
|
193
|
+
|
|
194
|
+
Evaluates the [quantile function][quantile-function] for a [Rayleigh][rayleigh-distribution] distribution with parameter `sigma` (scale parameter).
|
|
195
|
+
|
|
196
|
+
```c
|
|
197
|
+
double out = stdlib_base_dists_rayleigh_quantile( 0.8, 1.0 );
|
|
198
|
+
// returns ~1.794
|
|
199
|
+
```
|
|
200
|
+
|
|
201
|
+
The function accepts the following arguments:
|
|
202
|
+
|
|
203
|
+
- **p**: `[in] double` input probability.
|
|
204
|
+
- **sigma**: `[in] double` scale parameter.
|
|
205
|
+
|
|
206
|
+
```c
|
|
207
|
+
double stdlib_base_dists_rayleigh_quantile( const double p, const double sigma );
|
|
208
|
+
```
|
|
209
|
+
|
|
210
|
+
</section>
|
|
211
|
+
|
|
212
|
+
<!-- /.usage -->
|
|
213
|
+
|
|
214
|
+
<!-- C API usage notes. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
|
|
215
|
+
|
|
216
|
+
<section class="notes">
|
|
217
|
+
|
|
218
|
+
</section>
|
|
219
|
+
|
|
220
|
+
<!-- /.notes -->
|
|
221
|
+
|
|
222
|
+
<!-- C API usage examples. -->
|
|
223
|
+
|
|
224
|
+
<section class="examples">
|
|
225
|
+
|
|
226
|
+
### Examples
|
|
227
|
+
|
|
228
|
+
```c
|
|
229
|
+
#include "stdlib/stats/base/dists/rayleigh/quantile.h"
|
|
230
|
+
#include <stdlib.h>
|
|
231
|
+
#include <stdio.h>
|
|
232
|
+
|
|
233
|
+
static double random_uniform( const double min, const double max ) {
|
|
234
|
+
double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
|
|
235
|
+
return min + ( v*(max-min) );
|
|
236
|
+
}
|
|
237
|
+
|
|
238
|
+
int main( void ) {
|
|
239
|
+
double sigma;
|
|
240
|
+
double p;
|
|
241
|
+
double y;
|
|
242
|
+
int i;
|
|
243
|
+
|
|
244
|
+
for ( i = 0; i < 25; i++ ) {
|
|
245
|
+
p = random_uniform( 0.0, 1.0 );
|
|
246
|
+
sigma = random_uniform( 0.1, 10.0 );
|
|
247
|
+
y = stdlib_base_dists_rayleigh_quantile( p, sigma );
|
|
248
|
+
printf( "p: %lf, σ: %lf, Q(p;σ): %lf\n", p, sigma, y );
|
|
249
|
+
}
|
|
162
250
|
}
|
|
163
251
|
```
|
|
164
252
|
|
|
@@ -166,6 +254,10 @@ for ( i = 0; i < 10; i++ ) {
|
|
|
166
254
|
|
|
167
255
|
<!-- /.examples -->
|
|
168
256
|
|
|
257
|
+
</section>
|
|
258
|
+
|
|
259
|
+
<!-- /.c -->
|
|
260
|
+
|
|
169
261
|
<!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. -->
|
|
170
262
|
|
|
171
263
|
<section class="related">
|
|
@@ -200,7 +292,7 @@ See [LICENSE][stdlib-license].
|
|
|
200
292
|
|
|
201
293
|
## Copyright
|
|
202
294
|
|
|
203
|
-
Copyright © 2016-
|
|
295
|
+
Copyright © 2016-2026. The Stdlib [Authors][stdlib-authors].
|
|
204
296
|
|
|
205
297
|
</section>
|
|
206
298
|
|
|
@@ -213,8 +305,8 @@ Copyright © 2016-2024. The Stdlib [Authors][stdlib-authors].
|
|
|
213
305
|
[npm-image]: http://img.shields.io/npm/v/@stdlib/stats-base-dists-rayleigh-quantile.svg
|
|
214
306
|
[npm-url]: https://npmjs.org/package/@stdlib/stats-base-dists-rayleigh-quantile
|
|
215
307
|
|
|
216
|
-
[test-image]: https://github.com/stdlib-js/stats-base-dists-rayleigh-quantile/actions/workflows/test.yml/badge.svg?branch=v0.
|
|
217
|
-
[test-url]: https://github.com/stdlib-js/stats-base-dists-rayleigh-quantile/actions/workflows/test.yml?query=branch:v0.
|
|
308
|
+
[test-image]: https://github.com/stdlib-js/stats-base-dists-rayleigh-quantile/actions/workflows/test.yml/badge.svg?branch=v0.3.0
|
|
309
|
+
[test-url]: https://github.com/stdlib-js/stats-base-dists-rayleigh-quantile/actions/workflows/test.yml?query=branch:v0.3.0
|
|
218
310
|
|
|
219
311
|
[coverage-image]: https://img.shields.io/codecov/c/github/stdlib-js/stats-base-dists-rayleigh-quantile/main.svg
|
|
220
312
|
[coverage-url]: https://codecov.io/github/stdlib-js/stats-base-dists-rayleigh-quantile?branch=main
|
|
@@ -226,8 +318,8 @@ Copyright © 2016-2024. The Stdlib [Authors][stdlib-authors].
|
|
|
226
318
|
|
|
227
319
|
-->
|
|
228
320
|
|
|
229
|
-
[chat-image]: https://img.shields.io/
|
|
230
|
-
[chat-url]: https://
|
|
321
|
+
[chat-image]: https://img.shields.io/badge/zulip-join_chat-brightgreen.svg
|
|
322
|
+
[chat-url]: https://stdlib.zulipchat.com
|
|
231
323
|
|
|
232
324
|
[stdlib]: https://github.com/stdlib-js/stdlib
|
|
233
325
|
|
|
@@ -0,0 +1,38 @@
|
|
|
1
|
+
/**
|
|
2
|
+
* @license Apache-2.0
|
|
3
|
+
*
|
|
4
|
+
* Copyright (c) 2025 The Stdlib Authors.
|
|
5
|
+
*
|
|
6
|
+
* Licensed under the Apache License, Version 2.0 (the "License");
|
|
7
|
+
* you may not use this file except in compliance with the License.
|
|
8
|
+
* You may obtain a copy of the License at
|
|
9
|
+
*
|
|
10
|
+
* http://www.apache.org/licenses/LICENSE-2.0
|
|
11
|
+
*
|
|
12
|
+
* Unless required by applicable law or agreed to in writing, software
|
|
13
|
+
* distributed under the License is distributed on an "AS IS" BASIS,
|
|
14
|
+
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
15
|
+
* See the License for the specific language governing permissions and
|
|
16
|
+
* limitations under the License.
|
|
17
|
+
*/
|
|
18
|
+
|
|
19
|
+
#ifndef STDLIB_STATS_BASE_DISTS_RAYLEIGH_QUANTILE_H
|
|
20
|
+
#define STDLIB_STATS_BASE_DISTS_RAYLEIGH_QUANTILE_H
|
|
21
|
+
|
|
22
|
+
/*
|
|
23
|
+
* If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
|
|
24
|
+
*/
|
|
25
|
+
#ifdef __cplusplus
|
|
26
|
+
extern "C" {
|
|
27
|
+
#endif
|
|
28
|
+
|
|
29
|
+
/**
|
|
30
|
+
* Evaluates the quantile function for a Rayleigh distribution.
|
|
31
|
+
*/
|
|
32
|
+
double stdlib_base_dists_rayleigh_quantile( const double p, const double sigma );
|
|
33
|
+
|
|
34
|
+
#ifdef __cplusplus
|
|
35
|
+
}
|
|
36
|
+
#endif
|
|
37
|
+
|
|
38
|
+
#endif // !STDLIB_STATS_BASE_DISTS_RAYLEIGH_QUANTILE_H
|
package/lib/native.js
ADDED
|
@@ -0,0 +1,72 @@
|
|
|
1
|
+
/**
|
|
2
|
+
* @license Apache-2.0
|
|
3
|
+
*
|
|
4
|
+
* Copyright (c) 2025 The Stdlib Authors.
|
|
5
|
+
*
|
|
6
|
+
* Licensed under the Apache License, Version 2.0 (the "License");
|
|
7
|
+
* you may not use this file except in compliance with the License.
|
|
8
|
+
* You may obtain a copy of the License at
|
|
9
|
+
*
|
|
10
|
+
* http://www.apache.org/licenses/LICENSE-2.0
|
|
11
|
+
*
|
|
12
|
+
* Unless required by applicable law or agreed to in writing, software
|
|
13
|
+
* distributed under the License is distributed on an "AS IS" BASIS,
|
|
14
|
+
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
15
|
+
* See the License for the specific language governing permissions and
|
|
16
|
+
* limitations under the License.
|
|
17
|
+
*/
|
|
18
|
+
|
|
19
|
+
'use strict';
|
|
20
|
+
|
|
21
|
+
// MODULES //
|
|
22
|
+
|
|
23
|
+
var addon = require( './../src/addon.node' );
|
|
24
|
+
|
|
25
|
+
|
|
26
|
+
// MAIN //
|
|
27
|
+
|
|
28
|
+
/**
|
|
29
|
+
* Evaluates the quantile function for a Rayleigh distribution with scale parameter `sigma` at a probability `p`.
|
|
30
|
+
*
|
|
31
|
+
* @private
|
|
32
|
+
* @param {Probability} p - input value
|
|
33
|
+
* @param {NonNegativeNumber} sigma - scale parameter
|
|
34
|
+
* @returns {number} evaluated quantile function
|
|
35
|
+
*
|
|
36
|
+
* @example
|
|
37
|
+
* var y = quantile( 0.8, 1.0 );
|
|
38
|
+
* // returns ~1.794
|
|
39
|
+
*
|
|
40
|
+
* @example
|
|
41
|
+
* var y = quantile( 0.5, 4.0 );
|
|
42
|
+
* // returns ~4.71
|
|
43
|
+
*
|
|
44
|
+
* @example
|
|
45
|
+
* var y = quantile( 1.1, 1.0 );
|
|
46
|
+
* // returns NaN
|
|
47
|
+
*
|
|
48
|
+
* @example
|
|
49
|
+
* var y = quantile( -0.2, 1.0 );
|
|
50
|
+
* // returns NaN
|
|
51
|
+
*
|
|
52
|
+
* @example
|
|
53
|
+
* var y = quantile( NaN, 1.0 );
|
|
54
|
+
* // returns NaN
|
|
55
|
+
*
|
|
56
|
+
* @example
|
|
57
|
+
* var y = quantile( 0.0, NaN );
|
|
58
|
+
* // returns NaN
|
|
59
|
+
*
|
|
60
|
+
* @example
|
|
61
|
+
* // Negative scale parameter:
|
|
62
|
+
* var y = quantile( 0.5, -1.0 );
|
|
63
|
+
* // returns NaN
|
|
64
|
+
*/
|
|
65
|
+
function quantile( p, sigma ) {
|
|
66
|
+
return addon( p, sigma );
|
|
67
|
+
}
|
|
68
|
+
|
|
69
|
+
|
|
70
|
+
// EXPORTS //
|
|
71
|
+
|
|
72
|
+
module.exports = quantile;
|
package/manifest.json
ADDED
|
@@ -0,0 +1,83 @@
|
|
|
1
|
+
{
|
|
2
|
+
"options": {
|
|
3
|
+
"task": "build",
|
|
4
|
+
"wasm": false
|
|
5
|
+
},
|
|
6
|
+
"fields": [
|
|
7
|
+
{
|
|
8
|
+
"field": "src",
|
|
9
|
+
"resolve": true,
|
|
10
|
+
"relative": true
|
|
11
|
+
},
|
|
12
|
+
{
|
|
13
|
+
"field": "include",
|
|
14
|
+
"resolve": true,
|
|
15
|
+
"relative": true
|
|
16
|
+
},
|
|
17
|
+
{
|
|
18
|
+
"field": "libraries",
|
|
19
|
+
"resolve": false,
|
|
20
|
+
"relative": false
|
|
21
|
+
},
|
|
22
|
+
{
|
|
23
|
+
"field": "libpath",
|
|
24
|
+
"resolve": true,
|
|
25
|
+
"relative": false
|
|
26
|
+
}
|
|
27
|
+
],
|
|
28
|
+
"confs": [
|
|
29
|
+
{
|
|
30
|
+
"task": "build",
|
|
31
|
+
"wasm": false,
|
|
32
|
+
"src": [
|
|
33
|
+
"./src/main.c"
|
|
34
|
+
],
|
|
35
|
+
"include": [
|
|
36
|
+
"./include"
|
|
37
|
+
],
|
|
38
|
+
"libraries": [],
|
|
39
|
+
"libpath": [],
|
|
40
|
+
"dependencies": [
|
|
41
|
+
"@stdlib/math-base-napi-binary",
|
|
42
|
+
"@stdlib/math-base-assert-is-nan",
|
|
43
|
+
"@stdlib/math-base-special-log1p",
|
|
44
|
+
"@stdlib/math-base-special-sqrt"
|
|
45
|
+
]
|
|
46
|
+
},
|
|
47
|
+
{
|
|
48
|
+
"task": "benchmark",
|
|
49
|
+
"wasm": false,
|
|
50
|
+
"src": [
|
|
51
|
+
"./src/main.c"
|
|
52
|
+
],
|
|
53
|
+
"include": [
|
|
54
|
+
"./include"
|
|
55
|
+
],
|
|
56
|
+
"libraries": [],
|
|
57
|
+
"libpath": [],
|
|
58
|
+
"dependencies": [
|
|
59
|
+
"@stdlib/math-base-assert-is-nan",
|
|
60
|
+
"@stdlib/math-base-special-log1p",
|
|
61
|
+
"@stdlib/constants-float64-eps",
|
|
62
|
+
"@stdlib/math-base-special-sqrt"
|
|
63
|
+
]
|
|
64
|
+
},
|
|
65
|
+
{
|
|
66
|
+
"task": "examples",
|
|
67
|
+
"wasm": false,
|
|
68
|
+
"src": [
|
|
69
|
+
"./src/main.c"
|
|
70
|
+
],
|
|
71
|
+
"include": [
|
|
72
|
+
"./include"
|
|
73
|
+
],
|
|
74
|
+
"libraries": [],
|
|
75
|
+
"libpath": [],
|
|
76
|
+
"dependencies": [
|
|
77
|
+
"@stdlib/math-base-assert-is-nan",
|
|
78
|
+
"@stdlib/math-base-special-log1p",
|
|
79
|
+
"@stdlib/math-base-special-sqrt"
|
|
80
|
+
]
|
|
81
|
+
}
|
|
82
|
+
]
|
|
83
|
+
}
|
package/package.json
CHANGED
|
@@ -1,6 +1,6 @@
|
|
|
1
1
|
{
|
|
2
2
|
"name": "@stdlib/stats-base-dists-rayleigh-quantile",
|
|
3
|
-
"version": "0.
|
|
3
|
+
"version": "0.3.0",
|
|
4
4
|
"description": "Rayleigh distribution quantile function.",
|
|
5
5
|
"license": "Apache-2.0",
|
|
6
6
|
"author": {
|
|
@@ -14,9 +14,12 @@
|
|
|
14
14
|
}
|
|
15
15
|
],
|
|
16
16
|
"main": "./lib",
|
|
17
|
+
"gypfile": false,
|
|
17
18
|
"directories": {
|
|
18
19
|
"doc": "./docs",
|
|
20
|
+
"include": "./include",
|
|
19
21
|
"lib": "./lib",
|
|
22
|
+
"src": "./src",
|
|
20
23
|
"dist": "./dist"
|
|
21
24
|
},
|
|
22
25
|
"types": "./docs/types",
|
|
@@ -31,11 +34,13 @@
|
|
|
31
34
|
},
|
|
32
35
|
"dependencies": {
|
|
33
36
|
"@stdlib/math-base-assert-is-nan": "^0.2.2",
|
|
34
|
-
"@stdlib/math-base-
|
|
35
|
-
"@stdlib/math-base-special-
|
|
36
|
-
"@stdlib/
|
|
37
|
+
"@stdlib/math-base-napi-binary": "^0.3.1",
|
|
38
|
+
"@stdlib/math-base-special-log1p": "^0.2.3",
|
|
39
|
+
"@stdlib/math-base-special-sqrt": "^0.2.2",
|
|
40
|
+
"@stdlib/stats-base-dists-degenerate-quantile": "^0.3.0",
|
|
37
41
|
"@stdlib/utils-constant-function": "^0.2.2",
|
|
38
|
-
"@stdlib/utils-define-nonenumerable-read-only-property": "^0.2.2"
|
|
42
|
+
"@stdlib/utils-define-nonenumerable-read-only-property": "^0.2.2",
|
|
43
|
+
"@stdlib/utils-library-manifest": "^0.2.3"
|
|
39
44
|
},
|
|
40
45
|
"devDependencies": {},
|
|
41
46
|
"engines": {
|
package/src/addon.c
ADDED
|
@@ -0,0 +1,23 @@
|
|
|
1
|
+
/**
|
|
2
|
+
* @license Apache-2.0
|
|
3
|
+
*
|
|
4
|
+
* Copyright (c) 2025 The Stdlib Authors.
|
|
5
|
+
*
|
|
6
|
+
* Licensed under the Apache License, Version 2.0 (the "License");
|
|
7
|
+
* you may not use this file except in compliance with the License.
|
|
8
|
+
* You may obtain a copy of the License at
|
|
9
|
+
*
|
|
10
|
+
* http://www.apache.org/licenses/LICENSE-2.0
|
|
11
|
+
*
|
|
12
|
+
* Unless required by applicable law or agreed to in writing, software
|
|
13
|
+
* distributed under the License is distributed on an "AS IS" BASIS,
|
|
14
|
+
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
15
|
+
* See the License for the specific language governing permissions and
|
|
16
|
+
* limitations under the License.
|
|
17
|
+
*/
|
|
18
|
+
|
|
19
|
+
#include "stdlib/stats/base/dists/rayleigh/quantile.h"
|
|
20
|
+
#include "stdlib/math/base/napi/binary.h"
|
|
21
|
+
|
|
22
|
+
// cppcheck-suppress shadowFunction
|
|
23
|
+
STDLIB_MATH_BASE_NAPI_MODULE_DD_D( stdlib_base_dists_rayleigh_quantile )
|
package/src/main.c
ADDED
|
@@ -0,0 +1,51 @@
|
|
|
1
|
+
/**
|
|
2
|
+
* @license Apache-2.0
|
|
3
|
+
*
|
|
4
|
+
* Copyright (c) 2025 The Stdlib Authors.
|
|
5
|
+
*
|
|
6
|
+
* Licensed under the Apache License, Version 2.0 (the "License");
|
|
7
|
+
* you may not use this file except in compliance with the License.
|
|
8
|
+
* You may obtain a copy of the License at
|
|
9
|
+
*
|
|
10
|
+
* http://www.apache.org/licenses/LICENSE-2.0
|
|
11
|
+
*
|
|
12
|
+
* Unless required by applicable law or agreed to in writing, software
|
|
13
|
+
* distributed under the License is distributed on an "AS IS" BASIS,
|
|
14
|
+
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
15
|
+
* See the License for the specific language governing permissions and
|
|
16
|
+
* limitations under the License.
|
|
17
|
+
*/
|
|
18
|
+
|
|
19
|
+
#include "stdlib/stats/base/dists/rayleigh/quantile.h"
|
|
20
|
+
#include "stdlib/math/base/assert/is_nan.h"
|
|
21
|
+
#include "stdlib/math/base/special/log1p.h"
|
|
22
|
+
#include "stdlib/math/base/special/sqrt.h"
|
|
23
|
+
|
|
24
|
+
/**
|
|
25
|
+
* Evaluates the quantile function for a Rayleigh distribution with scale parameter `sigma` at a probability `p`.
|
|
26
|
+
*
|
|
27
|
+
* @param p input probability
|
|
28
|
+
* @param sigma scale parameter
|
|
29
|
+
* @return evaluated quantile function
|
|
30
|
+
*
|
|
31
|
+
* @example
|
|
32
|
+
* double y = stdlib_base_dists_rayleigh_quantile( 0.8, 1.0 );
|
|
33
|
+
* // returns ~1.794
|
|
34
|
+
*/
|
|
35
|
+
double stdlib_base_dists_rayleigh_quantile( const double p, const double sigma ) {
|
|
36
|
+
double s2;
|
|
37
|
+
if (
|
|
38
|
+
stdlib_base_is_nan( p ) ||
|
|
39
|
+
stdlib_base_is_nan( sigma ) ||
|
|
40
|
+
sigma < 0.0 ||
|
|
41
|
+
p < 0.0 ||
|
|
42
|
+
p > 1.0
|
|
43
|
+
) {
|
|
44
|
+
return 0.0 / 0.0; // NaN
|
|
45
|
+
}
|
|
46
|
+
if ( sigma == 0.0 ) {
|
|
47
|
+
return 0.0;
|
|
48
|
+
}
|
|
49
|
+
s2 = sigma * sigma;
|
|
50
|
+
return stdlib_base_sqrt( -2.0 * s2 * stdlib_base_log1p( -p ) );
|
|
51
|
+
}
|