Skip to content

Commit e211cdb

Browse files
committed
added examples, package.json and fixed if statement flow #209
1 parent cc278e6 commit e211cdb

File tree

7 files changed

+303
-21
lines changed

7 files changed

+303
-21
lines changed
Lines changed: 146 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,146 @@
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+
# VARIABLES #
20+
21+
ifndef VERBOSE
22+
QUIET := @
23+
else
24+
QUIET :=
25+
endif
26+
27+
# Determine the OS ([1][1], [2][2]).
28+
#
29+
# [1]: https://en.wikipedia.org/wiki/Uname#Examples
30+
# [2]: http://stackoverflow.com/a/27776822/2225624
31+
OS ?= $(shell uname)
32+
ifneq (, $(findstring MINGW,$(OS)))
33+
OS := WINNT
34+
else
35+
ifneq (, $(findstring MSYS,$(OS)))
36+
OS := WINNT
37+
else
38+
ifneq (, $(findstring CYGWIN,$(OS)))
39+
OS := WINNT
40+
else
41+
ifneq (, $(findstring Windows_NT,$(OS)))
42+
OS := WINNT
43+
endif
44+
endif
45+
endif
46+
endif
47+
48+
# Define the program used for compiling C source files:
49+
ifdef C_COMPILER
50+
CC := $(C_COMPILER)
51+
else
52+
CC := gcc
53+
endif
54+
55+
# Define the command-line options when compiling C files:
56+
CFLAGS ?= \
57+
-std=c99 \
58+
-O3 \
59+
-Wall \
60+
-pedantic
61+
62+
# Determine whether to generate position independent code ([1][1], [2][2]).
63+
#
64+
# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options
65+
# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option
66+
ifeq ($(OS), WINNT)
67+
fPIC ?=
68+
else
69+
fPIC ?= -fPIC
70+
endif
71+
72+
# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`):
73+
INCLUDE ?=
74+
75+
# List of source files:
76+
SOURCE_FILES ?=
77+
78+
# List of libraries (e.g., `-lopenblas -lpthread`):
79+
LIBRARIES ?=
80+
81+
# List of library paths (e.g., `-L /foo/bar -L /beep/boop`):
82+
LIBPATH ?=
83+
84+
# List of C targets:
85+
c_targets := example.out
86+
87+
88+
# RULES #
89+
90+
#/
91+
# Compiles source files.
92+
#
93+
# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`)
94+
# @param {string} [CFLAGS] - C compiler options
95+
# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`)
96+
# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`)
97+
# @param {string} [SOURCE_FILES] - list of source files
98+
# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`)
99+
# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`)
100+
#
101+
# @example
102+
# make
103+
#
104+
# @example
105+
# make all
106+
#/
107+
all: $(c_targets)
108+
109+
.PHONY: all
110+
111+
#/
112+
# Compiles C source files.
113+
#
114+
# @private
115+
# @param {string} CC - C compiler (e.g., `gcc`)
116+
# @param {string} CFLAGS - C compiler options
117+
# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`)
118+
# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`)
119+
# @param {string} SOURCE_FILES - list of source files
120+
# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`)
121+
# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`)
122+
#/
123+
$(c_targets): %.out: %.c
124+
$(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES)
125+
126+
#/
127+
# Runs compiled examples.
128+
#
129+
# @example
130+
# make run
131+
#/
132+
run: $(c_targets)
133+
$(QUIET) ./$<
134+
135+
.PHONY: run
136+
137+
#/
138+
# Removes generated files.
139+
#
140+
# @example
141+
# make clean
142+
#/
143+
clean:
144+
$(QUIET) -rm -f *.o *.out
145+
146+
.PHONY: clean
Lines changed: 43 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,43 @@
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/wald/pdf.h"
20+
#include "stdlib/constants/float64/eps.h"
21+
#include <stdlib.h>
22+
#include <stdio.h>
23+
24+
static double random_uniform( const double min, const double max ) {
25+
double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
26+
return min + ( v*(max-min) );
27+
}
28+
29+
int main( void ) {
30+
double lambda;
31+
double mu;
32+
double x;
33+
double y;
34+
int i;
35+
36+
for ( i = 0; i < 10; i++ ) {
37+
x = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
38+
mu = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
39+
lambda = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
40+
y = stdlib_base_dists_wald_pdf( x, mu, lambda );
41+
printf( "x: %lf, µ: %lf, λ: %lf, f(x;µ,λ): %lf\n", x, mu, lambda, y );
42+
}
43+
}
Lines changed: 36 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,36 @@
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+
var randu = require( '@stdlib/random/base/randu' );
22+
var pdf = require( './../lib' );
23+
24+
var lambda;
25+
var mu;
26+
var x;
27+
var y;
28+
var i;
29+
30+
for ( i = 0; i < 10; i++ ) {
31+
x = randu() * 10.0;
32+
mu = randu() * 10.0;
33+
lambda = randu() * 10.0;
34+
y = pdf( x, mu, lambda );
35+
console.log( 'x: %d, µ: %d, λ: %d, f(x;µ,λ): %d', x.toFixed( 4 ), mu.toFixed( 4 ), lambda.toFixed( 4 ), y.toFixed( 4 ) );
36+
}

lib/node_modules/@stdlib/stats/base/dists/wald/pdf/lib/factory.js

Lines changed: 3 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -51,16 +51,12 @@ function factory( mu, lambda ) {
5151
var B;
5252
if (
5353
isnan( mu ) ||
54-
isnan( lambda )
54+
isnan( lambda ) ||
55+
mu <= 0.0 ||
56+
lambda < 0.0
5557
) {
5658
return constantFunction( NaN );
5759
}
58-
if ( mu <= 0.0 ) {
59-
return constantFunction( NaN );
60-
}
61-
if ( lambda < 0.0 ) {
62-
return constantFunction( NaN );
63-
}
6460
if ( lambda === 0.0 ) {
6561
return degenerate( mu );
6662
}

lib/node_modules/@stdlib/stats/base/dists/wald/pdf/lib/main.js

Lines changed: 3 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -97,16 +97,12 @@ function pdf( x, mu, lambda ) {
9797
if (
9898
isnan( x ) ||
9999
isnan( mu ) ||
100-
isnan( lambda )
100+
isnan( lambda ) ||
101+
mu <= 0.0 ||
102+
lambda < 0.0
101103
) {
102104
return NaN;
103105
}
104-
if ( mu <= 0.0 ) {
105-
return NaN;
106-
}
107-
if ( lambda < 0.0 ) {
108-
return NaN;
109-
}
110106
if ( lambda === 0.0 ) {
111107
return ( x === mu ) ? PINF : 0.0;
112108
}
Lines changed: 69 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,69 @@
1+
{
2+
"name": "@stdlib/stats/base/dists/wald/pdf",
3+
"version": "0.0.0",
4+
"description": "Wald distribution probability density function (PDF).",
5+
"license": "Apache-2.0",
6+
"author": {
7+
"name": "The Stdlib Authors",
8+
"url": "https://github.com/stdlib-js/stdlib/graphs/contributors"
9+
},
10+
"contributors": [
11+
{
12+
"name": "The Stdlib Authors",
13+
"url": "https://github.com/stdlib-js/stdlib/graphs/contributors"
14+
}
15+
],
16+
"main": "./lib",
17+
"gypfile": true,
18+
"directories": {
19+
"benchmark": "./benchmark",
20+
"doc": "./docs",
21+
"example": "./examples",
22+
"include": "./include",
23+
"lib": "./lib",
24+
"src": "./src",
25+
"test": "./test"
26+
},
27+
"types": "./docs/types",
28+
"scripts": {},
29+
"homepage": "https://github.com/stdlib-js/stdlib",
30+
"repository": {
31+
"type": "git",
32+
"url": "git://github.com/stdlib-js/stdlib.git"
33+
},
34+
"bugs": {
35+
"url": "https://github.com/stdlib-js/stdlib/issues"
36+
},
37+
"dependencies": {},
38+
"devDependencies": {},
39+
"engines": {
40+
"node": ">=0.10.0",
41+
"npm": ">2.7.0"
42+
},
43+
"os": [
44+
"aix",
45+
"darwin",
46+
"freebsd",
47+
"linux",
48+
"macos",
49+
"openbsd",
50+
"sunos",
51+
"win32",
52+
"windows"
53+
],
54+
"keywords": [
55+
"stdlib",
56+
"stdmath",
57+
"statistics",
58+
"stats",
59+
"distribution",
60+
"dist",
61+
"continuous",
62+
"probability",
63+
"pdf",
64+
"wald",
65+
"inverse-gaussian",
66+
"inverse gaussian",
67+
"univariate"
68+
]
69+
}

lib/node_modules/@stdlib/stats/base/dists/wald/pdf/src/main.c

Lines changed: 3 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -42,16 +42,12 @@ double stdlib_base_dists_wald_pdf( const double x, const double mu, const double
4242
if (
4343
stdlib_base_is_nan( x ) ||
4444
stdlib_base_is_nan( mu ) ||
45-
stdlib_base_is_nan( lambda )
45+
stdlib_base_is_nan( lambda ) ||
46+
mu <= 0.0 ||
47+
lambda < 0.0
4648
) {
4749
return 0.0/0.0; // NaN
4850
}
49-
if ( mu <= 0.0 ) {
50-
return 0.0/0.0; // NaN
51-
}
52-
if ( lambda < 0.0 ) {
53-
return 0.0/0.0; // NaN
54-
}
5551
if ( lambda == 0.0 ) {
5652
return ( x == mu ) ? STDLIB_CONSTANT_FLOAT64_PINF : 0.0;
5753
}

0 commit comments

Comments
 (0)