Lines Matching refs:fX
720 double ScInterpreter::GetChiDist(double fX, double fDF) in GetChiDist() argument
723 if (fX <= 0.0) in GetChiDist()
726 return GetUpRegIGamma( fDF/2.0, fX/2.0); in GetChiDist()
733 double ScInterpreter::GetChiSqDistCDF(double fX, double fDF) in GetChiSqDistCDF() argument
736 if (fX <= 0.0) in GetChiSqDistCDF()
739 return GetLowRegIGamma( fDF/2.0, fX/2.0); in GetChiSqDistCDF()
742 double ScInterpreter::GetChiSqDistPDF(double fX, double fDF) in GetChiSqDistPDF() argument
747 if (fX <= 0.0) in GetChiSqDistPDF()
749 if (fDF*fX > 1391000.0) in GetChiSqDistPDF()
752 fValue = exp((0.5*fDF - 1) * log(fX*0.5) - 0.5 * fX - log(2.0) - GetLogGamma(0.5*fDF)); in GetChiSqDistPDF()
764 fValue = 1/sqrt(fX*2*F_PI); in GetChiSqDistPDF()
769 fValue *= (fX / fCount); in GetChiSqDistPDF()
772 if (fX>=1425.0) // underflow in e^(-x/2) in GetChiSqDistPDF()
773 fValue = exp(log(fValue)-fX/2); in GetChiSqDistPDF()
775 fValue *= exp(-fX/2); in GetChiSqDistPDF()
785 double fX; in ScChiSqDist() local
796 fX = GetDouble(); in ScChiSqDist()
798 PushDouble(GetChiSqDistCDF(fX,fDF)); in ScChiSqDist()
800 PushDouble(GetChiSqDistPDF(fX,fDF)); in ScChiSqDist()
896 double ScInterpreter::GetBetaDistPDF(double fX, double fA, double fB) in GetBetaDistPDF() argument
904 return -2.0*fX + 2.0; in GetBetaDistPDF()
905 if (fX == 1.0 && fB < 1.0) in GetBetaDistPDF()
910 if (fX <= 0.01) 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()
918 return fA * fX; in GetBetaDistPDF()
919 if (fX == 0.0 && fA < 1.0) in GetBetaDistPDF()
924 return fA * pow(fX,fA-1); in GetBetaDistPDF()
926 if (fX <= 0.0) in GetBetaDistPDF()
928 if (fA < 1.0 && fX == 0.0) in GetBetaDistPDF()
936 if (fX >= 1.0) in GetBetaDistPDF()
938 if (fB < 1.0 && fX == 1.0) in GetBetaDistPDF()
950 double fLogY = (fX < 0.1) ? ::boost::math::log1p(-fX) : log(0.5-fX+0.5); in GetBetaDistPDF()
951 double fLogX = log(fX); 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()
1038 double fX = fXin; in GetBetaDist() local
1047 fX = fY; in GetBetaDist()
1052 fResult = lcl_GetBetaHelperContFrac(fX,fA,fB); in GetBetaDist()
1058 fTemp = GetBetaDistPDF(fX,fA,fB)*fX*fY; in GetBetaDist()
1962 double fX = GetDouble(); // x in ScGammaDist() local
1968 PushDouble( GetGammaDist( fX, fAlpha, fBeta)); in ScGammaDist()
1970 PushDouble( GetGammaDistPDF( fX, fAlpha, fBeta)); in ScGammaDist()