bose() [1/2]

template<typename T >

T cntr::bose	(	T	beta,
		T	omega
	)

Evaluates the Bose distribution function.

Purpose

Evaluates the Bose distribution function \( n_\mathrm{B}(\omega)\) at given inverse temperature \(\beta\). The calculation is numerically stable for any \(\omega\ne 0 \).

Parameters

beta	[T] inverse temperature \(\beta\)
omega	[T] energy

Definition at line 131 of file cntr_equilibrium_impl.hpp.

Referenced by bose().

                        {
    T arg = omega * beta;
    if (arg < 0)
        return (-1.0 - bose<T>(beta, -omega));
    if (fabs(arg) > EXPMAX) {
        return 0.0;
    } else if (arg < 1e-10) {
        return 1.0 / arg;
    } else {
        return 1.0 / (exp(arg) - 1.0);
    }
}

Here is the caller graph for this function: