Lines Matching refs:fB
837 double fB; in GetBeta() local
840 fA = fAlpha; fB = fBeta; in GetBeta()
844 fA = fBeta; fB = fAlpha; in GetBeta()
846 if (fA+fB < fMaxGammaArgument) // simple case in GetBeta()
847 return GetGamma(fA)/GetGamma(fA+fB)*GetGamma(fB); in GetBeta()
854 fLanczos /= lcl_getLanczosSum(fA+fB); in GetBeta()
855 fLanczos *= lcl_getLanczosSum(fB); in GetBeta()
856 double fABgm = fA+fB+fgm; in GetBeta()
857 fLanczos *= sqrt((fABgm/(fA+fgm))/(fB+fgm)); in GetBeta()
858 double fTempA = fB/(fA+fgm); // (fA+fgm)/fABgm = 1 / ( 1 + fB/(fA+fgm)) in GetBeta()
859 double fTempB = fA/(fB+fgm); in GetBeta()
861 -fB * ::boost::math::log1p(fTempB)-fgm); in GetBeta()
870 double fB; in GetLogBeta() local
873 fA = fAlpha; fB = fBeta; in GetLogBeta()
877 fA = fBeta; fB = fAlpha; in GetLogBeta()
882 fLanczos /= lcl_getLanczosSum(fA+fB); in GetLogBeta()
883 fLanczos *= lcl_getLanczosSum(fB); in GetLogBeta()
885 double fABgm = fA+fB+fgm; in GetLogBeta()
886 fLogLanczos += 0.5*(log(fABgm)-log(fA+fgm)-log(fB+fgm)); in GetLogBeta()
887 double fTempA = fB/(fA+fgm); // (fA+fgm)/fABgm = 1 / ( 1 + fB/(fA+fgm)) in GetLogBeta()
888 double fTempB = fA/(fB+fgm); in GetLogBeta()
890 -fB * ::boost::math::log1p(fTempB)-fgm; in GetLogBeta()
896 double ScInterpreter::GetBetaDistPDF(double fX, double fA, double fB) in GetBetaDistPDF() argument
901 if (fB == 1.0) in GetBetaDistPDF()
903 if (fB == 2.0) in GetBetaDistPDF()
905 if (fX == 1.0 && fB < 1.0) in GetBetaDistPDF()
911 return fB + fB * ::boost::math::expm1((fB-1.0) * ::boost::math::log1p(-fX)); in GetBetaDistPDF()
913 return fB * pow(0.5-fX+0.5,fB-1.0); in GetBetaDistPDF()
915 if (fB == 1.0) // result a*x^(a-1) in GetBetaDistPDF()
938 if (fB < 1.0 && fX == 1.0) in GetBetaDistPDF()
953 double fBm1LogY = (fB-1.0) * fLogY; in GetBetaDistPDF()
954 double fLogBeta = GetLogBeta(fA,fB); in GetBetaDistPDF()
960 return pow(fX,fA-1.0) * pow(0.5-fX+0.5,fB-1.0) / GetBeta(fA,fB); in GetBetaDistPDF()
972 double lcl_GetBetaHelperContFrac(double fX, double fA, double fB) in lcl_GetBetaHelperContFrac() argument
976 b2 = 1.0 - (fA+fB)/(fA+1.0)*fX; in lcl_GetBetaHelperContFrac()
1000 d2m = rm*(fB-rm)*fX/((apl2m-1.0)*apl2m); in lcl_GetBetaHelperContFrac()
1001 d2m1 = -(fA+rm)*(fA+fB+rm)*fX/(apl2m*(apl2m+1.0)); in lcl_GetBetaHelperContFrac()
1041 double fB = fBeta; in GetBetaDist() local
1046 fB = fAlpha; in GetBetaDist()
1052 fResult = lcl_GetBetaHelperContFrac(fX,fA,fB); in GetBetaDist()
1054 double fP = fA/(fA+fB); in GetBetaDist()
1055 double fQ = fB/(fA+fB); in GetBetaDist()
1057 if (fA > 1.0 && fB > 1.0 && fP < 0.97 && fQ < 0.97) //found experimental in GetBetaDist()
1058 fTemp = GetBetaDistPDF(fX,fA,fB)*fX*fY; in GetBetaDist()
1060 fTemp = exp(fA*flnX + fB*flnY - GetLogBeta(fA,fB)); in GetBetaDist()
2075 double fP, fA, fB, fAlpha, fBeta; in ScBetaInv() local
2077 fB = GetDouble(); in ScBetaInv()
2079 fB = 1.0; in ScBetaInv()
2087 if (fP < 0.0 || fP >= 1.0 || fA == fB || fAlpha <= 0.0 || fBeta <= 0.0) in ScBetaInv()
2103 PushDouble(fA + fVal*(fB-fA)); // scale to (A,B) in ScBetaInv()