public class oadd.org.apache.commons.math.special.Beta extends java.lang.Object
{
private static final double DEFAULT_EPSILON;
private static final double HALF_LOG_TWO_PI;
private static final double[] DELTA;
private void <init>()
{
oadd.org.apache.commons.math.special.Beta v;
v := @this: oadd.org.apache.commons.math.special.Beta;
specialinvoke v.<java.lang.Object: void <init>()>();
return;
}
public static double regularizedBeta(double, double, double)
{
double v, v, v, v;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double regularizedBeta(double,double,double,double,int)>(v, v, v, 1.0E-14, 2147483647);
return v;
}
public static double regularizedBeta(double, double, double, double)
{
double v, v, v, v, v;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double regularizedBeta(double,double,double,double,int)>(v, v, v, v, 2147483647);
return v;
}
public static double regularizedBeta(double, double, double, int)
{
int v;
double v, v, v, v;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double;
v := @parameter: int;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double regularizedBeta(double,double,double,double,int)>(v, v, v, 1.0E-14, v);
return v;
}
public static double regularizedBeta(double, double, double, double, int)
{
oadd.org.apache.commons.math.special.Beta$1 v;
byte v, v, v, v, 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, v, v, v;
v := @parameter: double;
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 = 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 cmpl 1.0;
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 + 1.0;
v = v + v;
v = v + 2.0;
v = v / v;
v = v cmpl v;
if v <= 0 goto label;
v = 1.0 - v;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double regularizedBeta(double,double,double,double,int)>(v, v, v, v, v);
v = 1.0 - v;
goto label;
label:
v = new oadd.org.apache.commons.math.special.Beta$1;
specialinvoke v.<oadd.org.apache.commons.math.special.Beta$1: void <init>(double,double)>(v, v);
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = v * v;
v = 1.0 - 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.util.FastMath: double log(double)>(v);
v = v - v;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double logBeta(double,double)>(v, v);
v = v - v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double exp(double)>(v);
v = v * 1.0;
v = virtualinvoke v.<oadd.org.apache.commons.math.util.ContinuedFraction: double evaluate(double,double,int)>(v, v, v);
v = v / v;
label:
return v;
}
public static double logBeta(double, double, double, int)
{
int v;
double v, v, v, v;
v := @parameter: double;
v := @parameter: double;
v := @parameter: double;
v := @parameter: int;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double logBeta(double,double)>(v, v);
return v;
}
private static double logGammaSum(double, double) throws oadd.org.apache.commons.math.exception.OutOfRangeException
{
java.lang.Double v, v, v, v, v, v;
byte v, v, v, v, v, v;
oadd.org.apache.commons.math.exception.OutOfRangeException v, v;
double v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v;
v := @parameter: double;
v := @parameter: double;
v = v cmpg 1.0;
if v < 0 goto label;
v = v cmpl 2.0;
if v <= 0 goto label;
label:
v = new oadd.org.apache.commons.math.exception.OutOfRangeException;
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(v);
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(1.0);
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(2.0);
specialinvoke v.<oadd.org.apache.commons.math.exception.OutOfRangeException: void <init>(java.lang.Number,java.lang.Number,java.lang.Number)>(v, v, v);
throw v;
label:
v = v cmpg 1.0;
if v < 0 goto label;
v = v cmpl 2.0;
if v <= 0 goto label;
label:
v = new oadd.org.apache.commons.math.exception.OutOfRangeException;
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(v);
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(1.0);
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(2.0);
specialinvoke v.<oadd.org.apache.commons.math.exception.OutOfRangeException: void <init>(java.lang.Number,java.lang.Number,java.lang.Number)>(v, v, v);
throw v;
label:
v = v - 1.0;
v = v - 1.0;
v = v + v;
v = v cmpg 0.5;
if v > 0 goto label;
v = 1.0 + v;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma1p(double)>(v);
return v;
label:
v = v cmpg 1.5;
if v > 0 goto label;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma1p(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log1p(double)>(v);
v = v + v;
return v;
label:
v = v - 1.0;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma1p(double)>(v);
v = 1.0 + v;
v = v * v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = v + v;
return v;
}
private static double logGammaMinusLogGammaSum(double, double) throws oadd.org.apache.commons.math.exception.NumberIsTooSmallException
{
java.lang.Double v, v, v, v;
byte v, v, v, v;
double v, v, v, v, v, v, v, v, v, v, v, v, v, v, v;
oadd.org.apache.commons.math.exception.NumberIsTooSmallException v, v;
v := @parameter: double;
v := @parameter: double;
v = v cmpg 0.0;
if v >= 0 goto label;
v = new oadd.org.apache.commons.math.exception.NumberIsTooSmallException;
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(v);
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(0.0);
specialinvoke v.<oadd.org.apache.commons.math.exception.NumberIsTooSmallException: void <init>(java.lang.Number,java.lang.Number,boolean)>(v, v, 1);
throw v;
label:
v = v cmpg 10.0;
if v >= 0 goto label;
v = new oadd.org.apache.commons.math.exception.NumberIsTooSmallException;
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(v);
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(10.0);
specialinvoke v.<oadd.org.apache.commons.math.exception.NumberIsTooSmallException: void <init>(java.lang.Number,java.lang.Number,boolean)>(v, v, 1);
throw v;
label:
v = v cmpg v;
if v > 0 goto label;
v = v - 0.5;
v = v + v;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double deltaMinusDeltaSum(double,double)>(v, v);
goto label;
label:
v = v - 0.5;
v = v + v;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double deltaMinusDeltaSum(double,double)>(v, v);
label:
v = v / v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log1p(double)>(v);
v = v * v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = v - 1.0;
v = v * v;
v = v cmpg v;
if v > 0 goto label;
v = v - v;
v = v - v;
goto label;
label:
v = v - v;
v = v - v;
label:
return v;
}
private static double deltaMinusDeltaSum(double, double) throws oadd.org.apache.commons.math.exception.OutOfRangeException, oadd.org.apache.commons.math.exception.NumberIsTooSmallException
{
java.lang.Integer v, v;
byte v, v, v;
oadd.org.apache.commons.math.exception.NumberIsTooSmallException v;
java.lang.Double v, v, v;
double[] v, v, v, v, v, v;
int v, v, v, v, v, v, v, v, v, v;
oadd.org.apache.commons.math.exception.OutOfRangeException v;
double v, v, v, 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 = v cmpg 0.0;
if v < 0 goto label;
v = v cmpl v;
if v <= 0 goto label;
label:
v = new oadd.org.apache.commons.math.exception.OutOfRangeException;
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(v);
v = staticinvoke <java.lang.Integer: java.lang.Integer valueOf(int)>(0);
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(v);
specialinvoke v.<oadd.org.apache.commons.math.exception.OutOfRangeException: void <init>(java.lang.Number,java.lang.Number,java.lang.Number)>(v, v, v);
throw v;
label:
v = v cmpg 10.0;
if v >= 0 goto label;
v = new oadd.org.apache.commons.math.exception.NumberIsTooSmallException;
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(v);
v = staticinvoke <java.lang.Integer: java.lang.Integer valueOf(int)>(10);
specialinvoke v.<oadd.org.apache.commons.math.exception.NumberIsTooSmallException: void <init>(java.lang.Number,java.lang.Number,boolean)>(v, v, 1);
throw v;
label:
v = v / v;
v = 1.0 + v;
v = v / v;
v = 1.0 + v;
v = 1.0 / v;
v = v * v;
v = <oadd.org.apache.commons.math.special.Beta: double[] DELTA>;
v = lengthof v;
v = newarray (double)[v];
v[0] = 1.0;
v = 1;
label:
v = lengthof v;
if v >= v goto label;
v = v - 1;
v = v[v];
v = v * v;
v = v + v;
v = 1.0 + v;
v[v] = v;
v = v + 1;
goto label;
label:
v = 10.0 / v;
v = v * v;
v = <oadd.org.apache.commons.math.special.Beta: double[] DELTA>;
v = <oadd.org.apache.commons.math.special.Beta: double[] DELTA>;
v = lengthof v;
v = v - 1;
v = v[v];
v = lengthof v;
v = v - 1;
v = v[v];
v = v * v;
v = <oadd.org.apache.commons.math.special.Beta: double[] DELTA>;
v = lengthof v;
v = v - 2;
label:
if v < 0 goto label;
v = v * v;
v = <oadd.org.apache.commons.math.special.Beta: double[] DELTA>;
v = v[v];
v = v[v];
v = v * v;
v = v + v;
v = v - 1;
goto label;
label:
v = v * v;
v = v / v;
return v;
}
private static double sumDeltaMinusDeltaSum(double, double)
{
java.lang.Double v, v, v, v;
double[] v, v, v, v;
byte v, v;
int v, v, v, v;
double v, v, v, v, v, v, v, v, v, v, v, v;
oadd.org.apache.commons.math.exception.NumberIsTooSmallException v, v;
v := @parameter: double;
v := @parameter: double;
v = v cmpg 10.0;
if v >= 0 goto label;
v = new oadd.org.apache.commons.math.exception.NumberIsTooSmallException;
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(v);
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(10.0);
specialinvoke v.<oadd.org.apache.commons.math.exception.NumberIsTooSmallException: void <init>(java.lang.Number,java.lang.Number,boolean)>(v, v, 1);
throw v;
label:
v = v cmpg 10.0;
if v >= 0 goto label;
v = new oadd.org.apache.commons.math.exception.NumberIsTooSmallException;
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(v);
v = staticinvoke <java.lang.Double: java.lang.Double valueOf(double)>(10.0);
specialinvoke v.<oadd.org.apache.commons.math.exception.NumberIsTooSmallException: void <init>(java.lang.Number,java.lang.Number,boolean)>(v, v, 1);
throw v;
label:
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double min(double,double)>(v, v);
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double max(double,double)>(v, v);
v = 10.0 / v;
v = v * v;
v = <oadd.org.apache.commons.math.special.Beta: double[] DELTA>;
v = <oadd.org.apache.commons.math.special.Beta: double[] DELTA>;
v = lengthof v;
v = v - 1;
v = v[v];
v = <oadd.org.apache.commons.math.special.Beta: double[] DELTA>;
v = lengthof v;
v = v - 2;
label:
if v < 0 goto label;
v = v * v;
v = <oadd.org.apache.commons.math.special.Beta: double[] DELTA>;
v = v[v];
v = v + v;
v = v - 1;
goto label;
label:
v = v / v;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double deltaMinusDeltaSum(double,double)>(v, v);
v = v + v;
return v;
}
public static double logBeta(double, double)
{
byte v, v, v, v, v, v, v, v, v, v, v, v, v, v;
boolean v, v;
int v, v;
double v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, v, 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 = 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:
return #NaN;
label:
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double min(double,double)>(v, v);
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double max(double,double)>(v, v);
v = v cmpl 10.0;
if v < 0 goto label;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double sumDeltaMinusDeltaSum(double,double)>(v, v);
v = v / v;
v = 1.0 + v;
v = v / v;
v = v - 0.5;
v = neg v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = v * v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log1p(double)>(v);
v = v * v;
v = v cmpg v;
if v > 0 goto label;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = -0.5 * v;
v = v + 0.9189385332046727;
v = v + v;
v = v - v;
v = v - v;
return v;
label:
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = -0.5 * v;
v = v + 0.9189385332046727;
v = v + v;
v = v - v;
v = v - v;
return v;
label:
v = v cmpl 2.0;
if v <= 0 goto label;
v = v cmpl 1000.0;
if v <= 0 goto label;
v = v - 1.0;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double floor(double)>(v);
v = 1.0;
v = v;
v = 0;
label:
if v >= v goto label;
v = v - 1.0;
v = v / v;
v = 1.0 + v;
v = v / v;
v = v * v;
v = v + 1;
goto label;
label:
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(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 = staticinvoke <oadd.org.apache.commons.math.special.Beta: double logGammaMinusLogGammaSum(double,double)>(v, v);
v = v + v;
v = v + v;
return v;
label:
v = 1.0;
v = v;
label:
v = v cmpl 2.0;
if v <= 0 goto label;
v = v - 1.0;
v = v / v;
v = 1.0 + v;
v = v / v;
v = v * v;
goto label;
label:
v = v cmpg 10.0;
if v >= 0 goto label;
v = 1.0;
v = v;
label:
v = v cmpl 2.0;
if v <= 0 goto label;
v = v - 1.0;
v = v + v;
v = v / v;
v = v * v;
goto label;
label:
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = v + v;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double logGammaSum(double,double)>(v, v);
v = v - v;
v = v + v;
v = v + v;
return v;
label:
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma(double)>(v);
v = v + v;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double logGammaMinusLogGammaSum(double,double)>(v, v);
v = v + v;
return v;
label:
v = v cmpl 1.0;
if v < 0 goto label;
v = v cmpl 2.0;
if v <= 0 goto label;
v = v cmpg 10.0;
if v >= 0 goto label;
v = 1.0;
v = v;
label:
v = v cmpl 2.0;
if v <= 0 goto label;
v = v - 1.0;
v = v + v;
v = v / v;
v = v * v;
goto label;
label:
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double logGammaSum(double,double)>(v, v);
v = v - v;
v = v + v;
v = v + v;
return v;
label:
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double logGammaMinusLogGammaSum(double,double)>(v, v);
v = v + v;
return v;
label:
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma(double)>(v);
v = v + v;
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double logGammaSum(double,double)>(v, v);
v = v - v;
return v;
label:
v = v cmpl 10.0;
if v < 0 goto label;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double logGamma(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.special.Beta: double logGammaMinusLogGammaSum(double,double)>(v, v);
v = v + v;
return v;
label:
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double gamma(double)>(v);
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double gamma(double)>(v);
v = v * v;
v = v + v;
v = staticinvoke <oadd.org.apache.commons.math.special.Gamma: double gamma(double)>(v);
v = v / v;
v = staticinvoke <oadd.org.apache.commons.math.util.FastMath: double log(double)>(v);
return v;
}
static void <clinit>()
{
double[] v;
v = newarray (double)[15];
v[0] = 0.08333333333333333;
v[1] = -2.777777777777778E-5;
v[2] = 7.936507936507937E-8;
v[3] = -5.952380952380953E-10;
v[4] = 8.417508417508329E-12;
v[5] = -1.917526917518546E-13;
v[6] = 6.410256405103255E-15;
v[7] = -2.955065141253382E-16;
v[8] = 1.7964371635940225E-17;
v[9] = -1.3922896466162779E-18;
v[10] = 1.338028550140209E-19;
v[11] = -1.542460098679661E-20;
v[12] = 1.9770199298095743E-21;
v[13] = -2.3406566479399704E-22;
v[14] = 1.713480149663986E-23;
<oadd.org.apache.commons.math.special.Beta: double[] DELTA> = v;
return;
}
}