19 using namespace shogun;
30 return sum/values.
vlen;
44 return sum_squared_diff/(values.
vlen-1);
62 int32_t deg=values.
vlen-1;
73 conf_int_low=mean-interval;
74 conf_int_up=mean+interval;
168 big=4.503599627370496e15;
169 biginv=2.22044604925031308085e-16;
170 maxgam=171.624376956302725;
186 if (b*x<=1.0&&x<=0.95)
205 if ((flag==1&&b*x<=1.0)&&x<=0.95)
218 y=x*(a+b-2.0)-(a-1.0);
608 int32_t breaknewtcycle;
609 int32_t breakihalvecycle;
704 x=2.0/(1.0/(2.0*aaa-1.0)+1.0/(2.0*bbb-1.0));
706 -(1.0/(2.0*bbb-1.0)-1.0/(2.0*aaa-1.0))
707 *(lgm+5.0/6.0-2.0/(3.0*x));
729 if (mainlooppos==ihalve)
734 mainlooppos=ihalvecycle;
741 if (mainlooppos==ihalvecycle)
788 di=1.0-(1.0-di)*(1.0-di);
864 mainlooppos=ihalvecycle;
869 mainlooppos=breakihalvecycle;
877 if (mainlooppos==breakihalvecycle)
896 if (mainlooppos==newt)
905 mainlooppos=newtcycle;
912 if (mainlooppos==newtcycle)
948 mainlooppos=breaknewtcycle;
958 mainlooppos=breaknewtcycle;
967 xt=x0+0.5*yyy*(x-x0);
970 mainlooppos=breaknewtcycle;
977 xt=x1-0.5*yyy*(x1-x);
980 mainlooppos=breaknewtcycle;
990 mainlooppos=newtcycle;
995 mainlooppos=breaknewtcycle;
1003 if (mainlooppos==breaknewtcycle)
1048 expm2=0.13533528323661269189;
1049 s2pi=2.50662827463100050242;
1071 p0=-59.9633501014107895267;
1072 p0=98.0010754185999661536+y2*p0;
1073 p0=-56.6762857469070293439+y2*p0;
1074 p0=13.9312609387279679503+y2*p0;
1075 p0=-1.23916583867381258016+y2*p0;
1077 q0=1.95448858338141759834+y2*q0;
1078 q0=4.67627912898881538453+y2*q0;
1079 q0=86.3602421390890590575+y2*q0;
1080 q0=-225.462687854119370527+y2*q0;
1081 q0=200.260212380060660359+y2*q0;
1082 q0=-82.0372256168333339912+y2*q0;
1083 q0=15.9056225126211695515+y2*q0;
1084 q0=-1.18331621121330003142+y2*q0;
1095 p1=4.05544892305962419923;
1096 p1=31.5251094599893866154+z*p1;
1097 p1=57.1628192246421288162+z*p1;
1098 p1=44.0805073893200834700+z*p1;
1099 p1=14.6849561928858024014+z*p1;
1100 p1=2.18663306850790267539+z*p1;
1101 p1=-1.40256079171354495875*0.1+z*p1;
1102 p1=-3.50424626827848203418*0.01+z*p1;
1103 p1=-8.57456785154685413611*0.0001+z*p1;
1105 q1=15.7799883256466749731+z*q1;
1106 q1=45.3907635128879210584+z*q1;
1107 q1=41.3172038254672030440+z*q1;
1108 q1=15.0425385692907503408+z*q1;
1109 q1=2.50464946208309415979+z*q1;
1110 q1=-1.42182922854787788574*0.1+z*q1;
1111 q1=-3.80806407691578277194*0.01+z*q1;
1112 q1=-9.33259480895457427372*0.0001+z*q1;
1117 p2=3.23774891776946035970;
1118 p2=6.91522889068984211695+z*p2;
1119 p2=3.93881025292474443415+z*p2;
1120 p2=1.33303460815807542389+z*p2;
1121 p2=2.01485389549179081538*0.1+z*p2;
1122 p2=1.23716634817820021358*0.01+z*p2;
1123 p2=3.01581553508235416007*0.0001+z*p2;
1124 p2=2.65806974686737550832*0.000001+z*p2;
1125 p2=6.23974539184983293730*0.000000001+z*p2;
1127 q2=6.02427039364742014255+z*q2;
1128 q2=3.67983563856160859403+z*q2;
1129 q2=1.37702099489081330271+z*q2;
1130 q2=2.16236993594496635890*0.1+z*q2;
1131 q2=1.34204006088543189037*0.01+z*q2;
1132 q2=3.28014464682127739104*0.0001+z*q2;
1133 q2=2.89247864745380683936*0.000001+z*q2;
1134 q2=6.79019408009981274425*0.000000001+z*q2;