qsys_gig1_approx_gelenbe.h File Reference

Gelenbe diffusion approximation with instantaneous-return boundary. More...

#include "line/api/qsys/qsys_types.h"
#include "line/num/number.h"

Include dependency graph for qsys_gig1_approx_gelenbe.h:

Go to the source code of this file.

Namespaces
namespace	line
namespace	line::qsys

Functions
template<class T>
QsysResult< T >	line::qsys::qsys_gig1_approx_gelenbe (const T &lambda, const T &mu, const T &ca, const T &cs)
	Gelenbe diffusion approximation with instantaneous-return boundary.

Detailed Description

Gelenbe diffusion approximation with instantaneous-return boundary.

Templated port of matlab/src/api/qsys/qsys_gig1_approx_gelenbe.m.

rhat = exp(-2(1-rho)/(rho ca^2 + cs^2)) W = 1/(mu (1-rhat))

DIVERGENCE: jar/src/main/java/jline/api/qsys/Qsys_gig1_approx_gelenbe.java computes rhat = exp(-2(1-rho)/(rho*ca + cs)), treating its arguments as already-squared coefficients of variation. MATLAB is ground truth and squares, so this port squares. The two disagree unless ca = cs = 1.

Reference: Gelenbe, E. (1975). On approximate computer system models. Journal of the ACM 22(2), 261-269.

The geometric fit carries an exp, so this requires transcendental arithmetic and cannot be instantiated at T = Rational.

Definition in file qsys_gig1_approx_gelenbe.h.