NAG Library Manual, Mark 27.3

Interfaces: FL CL CPP AD

NAG CL Interface Introduction

G01 (Stat) Chapter Contents

G01 (Stat) Chapter Introduction

g01mu: FL CL CPP AD

NAG CL Interface
g01muc (pdf_vavilov)

Keyword Search:

NAG Library Manual, Mark 27.3

Interfaces: FL CL CPP AD

NAG CL Interface Introduction

G01 (Stat) Chapter Contents

G01 (Stat) Chapter Introduction

g01mu: FL CL CPP AD

▸▿ Contents

2 Specification

3 Description

4 References

5 Arguments

6 Error Indicators and Warnings

8 Parallelism and Performance

9 Further Comments

▸▿ 10 Example

10.1 Program Text

10.2 Program Data

10.3 Program Results

© The Numerical Algorithms Group Ltd. 2021

CL Name Style:
Short (a00aac)
Long (impl_details)
Full (nag_info_impl_details)

1 Purpose

g01muc returns the value of the Vavilov density function

ϕ_{V} (λ; κ, β^{2})

.

It is intended to be used after a call to g01zuc.

2 Specification

#include <nag.h>

double

g01muc (double x, const double comm_arr[])

The function may be called by the names: g01muc, nag_stat_pdf_vavilov or nag_prob_density_vavilov.

3 Description

g01muc evaluates an approximation to the Vavilov density function

ϕ_{V} (λ; κ, β^{2})

given by

ϕ_{V} (λ; κ, β^{2}) = \frac{1}{2 π i} \int_{c - i \infty}^{c + i \infty} e^{λ s} f (s; κ, β^{2}) d s,

where

κ > 0

and

0 \leq β^{2} \leq 1

,

c

is an arbitrary real constant and

f (s; κ, β^{2}) = C (κ, β^{2}) \exp {s \ln κ + (s + κ β^{2}) [\ln (\frac{s}{κ}) + E_{1} (\frac{s}{κ})] - κ \exp (- \frac{s}{κ})} .

E_{1} (x) = \int_{0}^{x} t^{- 1} (1 - e^{- t}) d t

is the exponential integral,

C (κ, β^{2}) = \exp {κ (1 + γ β^{2})}

and

γ

is Euler's constant.

The method used is based on Fourier expansions. Further details can be found in Schorr (1974).

For values of

κ \leq 0.01

, the Vavilov distribution can be replaced by the Landau distribution since

λ_{V} = (λ_{L} - \ln κ) / κ

. For values of

κ \geq 10

, the Vavilov distribution can be replaced by a Gaussian distribution with mean

μ = γ - 1 - β^{2} - \ln κ

and variance

σ^{2} = (2 - β^{2}) / 2 κ

.

4 References

Schorr B (1974) Programs for the Landau and the Vavilov distributions and the corresponding random numbers Comp. Phys. Comm. 7 215–224

5 Arguments

1: $x$ – double Input: On entry: the argument $λ$ of the function.
2: $comm_arr [322]$ – const double Communication Array: On entry: this must be the same argument comm_arr as returned by a previous call to g01zuc.

6 Error Indicators and Warnings

None.

7 Accuracy

At least five significant digits are usually correct.

8 Parallelism and Performance

g01muc is not threaded in any implementation.

9 Further Comments

g01muc can be called repeatedly with different values of

λ

provided that the values of

κ

and

β^{2}

remain unchanged between calls. Otherwise, g01zuc must be called again.

10 Example

This example evaluates

ϕ_{V} (λ; κ, β^{2})

at

λ = 2.5

,

κ = 0.4

and

β^{2} = 0.1

, and prints the results.

10.1 Program Text

Program Text (g01muce.c)

10.2 Program Data

Program Data (g01muce.d)

10.3 Program Results

Program Results (g01muce.r)

NAG Library Manual, Mark 27.3

Interfaces: FL CL CPP AD

NAG CL Interface Introduction

G01 (Stat) Chapter Contents

G01 (Stat) Chapter Introduction

g01mu: FL CL CPP AD

© The Numerical Algorithms Group Ltd. 2021