log Gamma - Java java.lang

Description

log Gamma

Demo Code

//package com.java2s;

import static java.lang.Math.*;

public class Main {
    private static final double sqrt2Pi = sqrt(2 * PI);

    private static double logGamma(double x) {
        double tmp = (x - 0.5) * log(x + 4.5) - (x + 4.5);
        double ser = 1.0 + 76.18009173 / (x + 0) - 86.50532033 / (x + 1)
                + 24.01409822 / (x + 2) - 1.231739516 / (x + 3)
                + 0.00120858003 / (x + 4) - 0.00000536382 / (x + 5);
        return tmp + log(ser * sqrt2Pi);
    }/*from w  w w .ja v a2  s  .  c o m*/
}

• Previous
• Next

Related Tutorials

• Returns natural logarithm of Gamma function of x Algorithm is from Numerical recipes 6.1
• Calculates the natural logarithm of a number
• Computes log2 (ceil (a)).
• Returns the natural log of the Binomial Coefficient, "n choose k", the number of k-element subsets that can be selected from an n-element set.
• Given log(a) and log(b) calculate log(a + b) boils down to log( exp(log_a) + exp(log_b) ) but this might overflow, so we turn this into log([exp(log_a - log_c) + exp(log_b - log_c)]exp(log_c)) and we set log_c == max(log_a,log_b) and so it becomes: LARGE + log(1 + exp(SMALL - LARGE)) == LARGE + log(1 + SMALL) ~= large the whole idea being to avoid an overflow (exp(LARGE) == VERY LARGE == overflow)