public class oadd.org.apache.commons.math.special.Gamma extends java.lang.Object
{
public static final double GAMMA;
private static final double DEFAULT_EPSILON;
private static final double[] LANCZOS;
private static final double HALF_LOG_2_PI;
private static final double C_LIMIT;
private static final double S_LIMIT;
private void <init>()
{
oadd.org.apache.commons.math.special.Gamma v;
v := @this: oadd.org.apache.commons.math.special.Gamma;
specialinvoke v.<java.lang.Object: void <init>()>();
return;
}
public static double logGamma(double)
{
double[] v, v, v;
byte v;
int v, v;
boolean v;
double v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v;
v := @parameter: double;
v = staticinvoke <java.lang.Double: boolean isNaN(double)>(v);
if v != 0 goto label;
v = v cmpg 0.0;
if v > 0 goto label;
label:
v = #NaN;
goto label;
label:
v = 0.0;
v = <oadd.org.apache.commons.math.special.Gamma: double[] LANCZOS>;
v = lengthof v;
v = v - 1;
label:
if v <= 0 goto label;
v = <oadd.org.apache.commons.math.special.Gamma: double[] LANCZOS>;
v = v[v];
v = v + v;
v = v / v;
v = v + v;
v = v - 1;
goto label;
label:
v = <oadd.org.apache.commons.math.special.Gamma: double[] LANCZOS>;
v = v[0];
v = v + v;
v = v + 4.7421875;
v = v + 0.5;
v = v + 0.5;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = v * v;
v = v - v;
v = <oadd.org.apache.commons.math.special.Gamma: double HALF_LOG_2_PI>;
v = v + v;
v = v / v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = v + v;
label:
return v;
}
public static double regularizedGammaP(double, double) throws oadd.org.apache.commons.math.MathException
{
double v, v, v;
v := @parameter: double;
v := @parameter: double;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double regularizedGammaP(double,double,double,int)>(v, v, 1.0E-14, 2147483647);
return v;
}
public static double regularizedGammaP(double, double, double, int) throws oadd.org.apache.commons.math.MathException
{
byte v, v, v, v, v, v, v, v;
oadd.org.apache.commons.math.MaxIterationsExceededException v;
boolean v, v, v;
int v;
double v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double;
v := @parameter: int;
v = staticinvoke <java.lang.Double: boolean isNaN(double)>(v);
if v != 0 goto label;
v = staticinvoke <java.lang.Double: boolean isNaN(double)>(v);
if v != 0 goto label;
v = v cmpg 0.0;
if v <= 0 goto label;
v = v cmpg 0.0;
if v >= 0 goto label;
label:
v = #NaN;
goto label;
label:
v = v cmpl 0.0;
if v != 0 goto label;
v = 0.0;
goto label;
label:
v = v + 1.0;
v = v cmpl v;
if v < 0 goto label;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double regularizedGammaQ(double,double,double,int)>(v, v, v, v);
v = 1.0 - v;
goto label;
label:
v = 0.0;
v = 1.0 / v;
v = v;
label:
v = v / v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double abs(double)>(v);
v = v cmpl v;
if v <= 0 goto label;
v = v cmpg v;
if v >= 0 goto label;
v = v cmpg #Infinity;
if v >= 0 goto label;
v = v + 1.0;
v = v + v;
v = v / v;
v = v * v;
v = v + v;
goto label;
label:
v = v cmpl v;
if v < 0 goto label;
v = new oadd.org.apache.commons.math.MaxIterationsExceededException;
specialinvoke v.<oadd.org.apache.commons.math.MaxIterationsExceededException: void <init>(int)>(v);
throw v;
label:
v = staticinvoke <java.lang.Double: boolean isInfinite(double)>(v);
if v == 0 goto label;
v = 1.0;
goto label;
label:
v = neg v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = v * v;
v = v + v;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma(double)>(v);
v = v - v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double exp(double)>(v);
v = v * v;
label:
return v;
}
public static double regularizedGammaQ(double, double) throws oadd.org.apache.commons.math.MathException
{
double v, v, v;
v := @parameter: double;
v := @parameter: double;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double regularizedGammaQ(double,double,double,int)>(v, v, 1.0E-14, 2147483647);
return v;
}
public static double regularizedGammaQ(double, double, double, int) throws oadd.org.apache.commons.math.MathException
{
byte v, v, v, v;
int v;
boolean v, v;
double v, v, v, v, v, v, v, v, v, v, v, v, v, v, v;
oadd.org.apache.commons.math.special.Gamma$1 v;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double;
v := @parameter: int;
v = staticinvoke <java.lang.Double: boolean isNaN(double)>(v);
if v != 0 goto label;
v = staticinvoke <java.lang.Double: boolean isNaN(double)>(v);
if v != 0 goto label;
v = v cmpg 0.0;
if v <= 0 goto label;
v = v cmpg 0.0;
if v >= 0 goto label;
label:
v = #NaN;
goto label;
label:
v = v cmpl 0.0;
if v != 0 goto label;
v = 1.0;
goto label;
label:
v = v + 1.0;
v = v cmpg v;
if v >= 0 goto label;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double regularizedGammaP(double,double,double,int)>(v, v, v, v);
v = 1.0 - v;
goto label;
label:
v = new oadd.org.apache.commons.math.special.Gamma$1;
specialinvoke v.<oadd.org.apache.commons.math.special.Gamma$1: void <init>(double)>(v);
v = virtualinvoke v.<oadd.org.apache.commons.math.util.ContinuedFraction: double evaluate(double,double,int)>(v, v, v);
v = 1.0 / v;
v = neg v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = v * v;
v = v + v;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma(double)>(v);
v = v - v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double exp(double)>(v);
v = v * v;
label:
return v;
}
public static double digamma(double)
{
byte v, v, v;
double v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v;
v := @parameter: double;
v = v cmpl 0.0;
if v <= 0 goto label;
v = v cmpg 1.0E-5;
if v > 0 goto label;
v = 1.0 / v;
v = -0.5772156649015329 - v;
return v;
label:
v = v cmpl 49.0;
if v < 0 goto label;
v = v * v;
v = 1.0 / v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = 0.5 / v;
v = v - v;
v = v / 252.0;
v = 0.008333333333333333 - v;
v = v * v;
v = 0.08333333333333333 + v;
v = v * v;
v = v - v;
return v;
label:
v = v + 1.0;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double digamma(double)>(v);
v = 1.0 / v;
v = v - v;
return v;
}
public static double trigamma(double)
{
byte v, v, v;
double v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v;
v := @parameter: double;
v = v cmpl 0.0;
if v <= 0 goto label;
v = v cmpg 1.0E-5;
if v > 0 goto label;
v = v * v;
v = 1.0 / v;
return v;
label:
v = v cmpl 49.0;
if v < 0 goto label;
v = v * v;
v = 1.0 / v;
v = 1.0 / v;
v = v / 2.0;
v = v + v;
v = v / v;
v = v / 42.0;
v = 0.03333333333333333 + v;
v = v * v;
v = 0.16666666666666666 - v;
v = v * v;
v = v + v;
return v;
label:
v = v + 1.0;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double trigamma(double)>(v);
v = v * v;
v = 1.0 / v;
v = v + v;
return v;
}
static void <clinit>()
{
double[] v;
double v, v;
v = newarray (double)[15];
v[0] = 0.9999999999999971;
v[1] = 57.15623566586292;
v[2] = -59.59796035547549;
v[3] = 14.136097974741746;
v[4] = -0.4919138160976202;
v[5] = 3.399464998481189E-5;
v[6] = 4.652362892704858E-5;
v[7] = -9.837447530487956E-5;
v[8] = 1.580887032249125E-4;
v[9] = -2.1026444172410488E-4;
v[10] = 2.1743961811521265E-4;
v[11] = -1.643181065367639E-4;
v[12] = 8.441822398385275E-5;
v[13] = -2.6190838401581408E-5;
v[14] = 3.6899182659531625E-6;
<oadd.org.apache.commons.math.special.Gamma: double[] LANCZOS> = v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(6.283185307179586);
v = 0.5 * v;
<oadd.org.apache.commons.math.special.Gamma: double HALF_LOG_2_PI> = v;
return;
}
}