| line |
true |
false |
branch |
|
41
|
5 |
1 |
qsame = olda == aa && oldb == bb; |
|
|
5 |
0 |
qsame = olda == aa && oldb == bb; |
|
42
|
5 |
1 |
if(qsame) goto S20; |
|
43
|
1 |
0 |
if(!(aa < minlog || bb < minlog)) goto S10; |
|
|
1 |
0 |
if(!(aa < minlog || bb < minlog)) goto S10; |
|
51
|
6 |
0 |
if(!(min(aa,bb) > 1.0)) goto S100; |
|
|
0 |
6 |
if(!(min(aa,bb) > 1.0)) goto S100; |
|
56
|
5 |
1 |
if(qsame) goto S30; |
|
57
|
1 |
0 |
a = min(aa,bb); |
|
58
|
0 |
1 |
b = max(aa,bb); |
|
70
|
0 |
6 |
if(v > expmax) goto S55; |
|
76
|
0 |
6 |
if(w > infnty/a) goto S55; |
|
88
|
4 |
2 |
if(s+2.60943791243410 >= 5.0*z) goto S70; |
|
93
|
2 |
0 |
if(s > t) goto S70; |
|
103
|
0 |
0 |
if(alpha/(b+w) < minlog) goto S30; |
|
104
|
0 |
0 |
if(r+alpha*log(alpha/(b+w)) < t) goto S30; |
|
109
|
6 |
0 |
if(aa == a) { |
|
120
|
0 |
0 |
if(qsame) goto S110; |
|
121
|
0 |
0 |
a = max(aa,bb); |
|
122
|
0 |
0 |
b = min(aa,bb); |
|
136
|
0 |
0 |
if(u1 >= 0.5) goto S130; |
|
142
|
0 |
0 |
if(0.25*u2+z-y >= k1) goto S120; |
|
149
|
0 |
0 |
if(!(z <= 0.25)) goto S160; |
|
155
|
0 |
0 |
if(a > 1.0) goto S135; |
|
157
|
0 |
0 |
if(v > expmax) goto S132; |
|
162
|
0 |
0 |
if(w > expmax) goto S140; |
|
167
|
0 |
0 |
if(v > expmax) goto S140; |
|
169
|
0 |
0 |
if(w > infnty/a) goto S140; |
|
186
|
0 |
0 |
if(z >= k2) goto S120; |
|
194
|
0 |
0 |
if(a > 1.0) goto S175; |
|
196
|
0 |
0 |
if(v > expmax) goto S172; |
|
201
|
0 |
0 |
if(w > expmax) goto S180; |
|
206
|
0 |
0 |
if(v > expmax) goto S180; |
|
208
|
0 |
0 |
if(w > infnty/a) goto S180; |
|
226
|
0 |
0 |
if(alpha/(b+w) < minlog) goto S120; |
|
227
|
0 |
0 |
if(alpha*(log(alpha/(b+w))+v)-1.38629436111989 < log(z)) goto S120; |
|
232
|
0 |
0 |
if(a == aa) { |
|
263
|
6 |
0 |
if(!(df <= 0.0)) goto S10; |
|
302
|
6 |
0 |
if(av >= 0.0) goto S10; |
|
332
|
6 |
0 |
if(!(dfn <= 0.0 || dfd <= 0.0)) goto S10; |
|
|
6 |
0 |
if(!(dfn <= 0.0 || dfd <= 0.0)) goto S10; |
|
350
|
6 |
0 |
if(!(xden <= 1.0E-37*xnum)) goto S20; |
|
404
|
6 |
0 |
if(a > 0.0 && r > 0.0) goto S10; |
|
|
6 |
0 |
if(a > 0.0 && r > 0.0) goto S10; |
|
443
|
8 |
4 |
for(i=1; i<=p; i++) *(work+i-1) = snorm(); |
|
444
|
8 |
4 |
for(i=1,D3=1,D4=(p-i+D3)/D3; D4>0; D4--,i+=D3) { |
|
458
|
12 |
8 |
for(j=1,D1=1,D2=(i-j+D1)/D1; D2>0; D2--,j+=D1) { |
|
496
|
0 |
2 |
if(n < 0) ftnstop("N < 0 in GENMUL"); |
|
497
|
0 |
2 |
if(ncat <= 1) ftnstop("NCAT <= 1 in GENMUL"); |
|
499
|
4 |
2 |
for(i=0; i
|
|
500
|
0 |
4 |
if(*(p+i) < 0.0F) ftnstop("Some P(i) < 0 in GENMUL"); |
|
501
|
0 |
4 |
if(*(p+i) > 1.0F) ftnstop("Some P(i) > 1 in GENMUL"); |
|
504
|
0 |
2 |
if(ptot > 0.99999F) ftnstop("Sum of P(i) > 1 in GENMUL"); |
|
510
|
6 |
2 |
for(i=0; i
|
|
514
|
4 |
2 |
for(icat=0; icat
|
|
518
|
0 |
4 |
if(ntot <= 0) return; |
|
551
|
6 |
0 |
if(!(df < 1.0 || xnonc < 0.0)) goto S10; |
|
|
6 |
0 |
if(!(df < 1.0 || xnonc < 0.0)) goto S10; |
|
557
|
6 |
0 |
if(df >= 1.000000001) goto S20; |
|
614
|
6 |
0 |
qcond = dfn < 1.0 || dfd <= 0.0 || xnonc < 0.0; |
|
|
6 |
0 |
qcond = dfn < 1.0 || dfd <= 0.0 || xnonc < 0.0; |
|
|
0 |
6 |
qcond = dfn < 1.0 || dfd <= 0.0 || xnonc < 0.0; |
|
615
|
6 |
0 |
if(!qcond) goto S10; |
|
630
|
6 |
0 |
if(dfn >= 1.000001) goto S20; |
|
656
|
6 |
0 |
if(!(xden <= 1.0E-37*xnum)) goto S40; |
|
700
|
6 |
0 |
if(sd >= 0.0) goto S10; |
|
723
|
34 |
4 |
for(i=1,D1=1,D2=(larray-i+D1)/D1; D2>0; D2--,i+=D1) { |
|
746
|
14 |
0 |
if(!(low > high)) goto S10; |
|
770
|
1022 |
513 |
if(getset == 0) *g = curntg; |
|
772
|
513 |
0 |
if(*g < 0 || *g > numg) { |
|
|
0 |
513 |
if(*g < 0 || *g > numg) { |
|
795
|
1019 |
3 |
if(getset == 0) *qvalue = qinit; |
|
812
|
16 |
469 |
if(getset != 0) qstate = 1; |
|
931
|
6 |
4 |
if(pp != psave) goto S10; |
|
932
|
0 |
4 |
if(n != nsave) goto S20; |
|
933
|
4 |
0 |
if(xnp < 30.0) goto S150; |
|
940
|
0 |
6 |
if(pp < 0.0F) ftnstop("PP < 0.0 in IGNBIN"); |
|
941
|
0 |
6 |
if(pp > 1.0F) ftnstop("PP > 1.0 in IGNBIN"); |
|
943
|
4 |
2 |
p = min(psave,1.0-psave); |
|
949
|
0 |
6 |
if(n < 0L) ftnstop("N < 0 in IGNBIN"); |
|
952
|
6 |
0 |
if(xnp < 30.0) goto S140; |
|
978
|
0 |
0 |
if(u > p1) goto S40; |
|
985
|
0 |
0 |
if(u > p2) goto S50; |
|
987
|
0 |
0 |
v = v*c+1.0-ABS(xm-x)/p1; |
|
988
|
0 |
0 |
if(v > 1.0 || v <= 0.0) goto S30; |
|
|
0 |
0 |
if(v > 1.0 || v <= 0.0) goto S30; |
|
995
|
0 |
0 |
if(u > p3) goto S60; |
|
997
|
0 |
0 |
if(ix < 0) goto S30; |
|
1005
|
0 |
0 |
if(ix > n) goto S30; |
|
1012
|
0 |
0 |
if(k > 20 && k < xnpq/2-1) goto S130; |
|
|
0 |
0 |
if(k > 20 && k < xnpq/2-1) goto S130; |
|
1020
|
0 |
0 |
if(T1 < 0) goto S80; |
|
1021
|
0 |
0 |
else if(T1 == 0) goto S120; |
|
1025
|
0 |
0 |
for(i=mp; i<=ix; i++) f *= (g/i-r); |
|
1029
|
0 |
0 |
for(i=ix1; i<=m; i++) f /= (g/i-r); |
|
1031
|
0 |
0 |
if(v <= f) goto S170; |
|
1040
|
0 |
0 |
if(alv < ynorm-amaxp) goto S170; |
|
1041
|
0 |
0 |
if(alv > ynorm+amaxp) goto S30; |
|
1058
|
0 |
0 |
-140.0/w2)/w2)/w2)/w2)/w/166320.0) goto S170; |
|
1075
|
10 |
76 |
if(u < f) goto S170; |
|
1076
|
0 |
76 |
if(ix > 110) goto S150; |
|
1082
|
2 |
8 |
if(psave > 0.5) ix = n-ix; |
|
1121
|
0 |
6 |
if(n <= 0L) ftnstop("N <= 0 in IGNNBN"); |
|
1122
|
0 |
6 |
if(p <= 0.0F) ftnstop("P <= 0.0 in IGNNBN"); |
|
1123
|
0 |
6 |
if(p >= 1.0F) ftnstop("P >= 1.0 in IGNNBN"); |
|
1217
|
4 |
8 |
if(mu == muprev) goto S10; |
|
1218
|
6 |
2 |
if(mu < 10.0) goto S120; |
|
1237
|
0 |
6 |
if(g < 0.0) goto S20; |
|
1242
|
2 |
4 |
if(ignpoi >= ll) return ignpoi; |
|
1249
|
4 |
0 |
if(d*u >= difmuk*difmuk*difmuk) return ignpoi; |
|
1259
|
0 |
0 |
if(mu == muold) goto S30; |
|
1270
|
0 |
0 |
if(g < 0.0) goto S50; |
|
1280
|
0 |
0 |
if(fy-u*fy <= py*exp(px-fx)) return ignpoi; |
|
1291
|
0 |
0 |
if(t <= -0.6744) goto S50; |
|
1304
|
0 |
0 |
if(c*fabs(u) > py*exp(px+e)-fy*exp(fx+e)) goto S50; |
|
1311
|
0 |
0 |
if(ignpoi >= 10) goto S80; |
|
1324
|
0 |
0 |
if(fabs(v) <= 0.25) goto S90; |
|
1336
|
0 |
0 |
if(kflag <= 0) goto S40; |
|
1344
|
0 |
6 |
if(mu == muold) goto S130; |
|
1346
|
6 |
0 |
if(mu >= 0.0) goto S125; |
|
1362
|
2 |
4 |
if(u <= p0) return ignpoi; |
|
1368
|
4 |
0 |
if(l == 0) goto S150; |
|
1370
|
0 |
0 |
if(u > 0.458) j = min(l,m); |
|
1371
|
0 |
0 |
for(k=j; k<=l; k++) { |
|
1372
|
0 |
0 |
if(u <= *(pp+k-1)) goto S180; |
|
1374
|
0 |
0 |
if(l == 35) goto S130; |
|
1381
|
14 |
0 |
for(k=l; k<=35; k++) { |
|
1385
|
4 |
10 |
if(u <= q) goto S170; |
|
1417
|
40 |
0 |
if(!(low > high)) goto S10; |
|
1423
|
40 |
0 |
if(!(range > maxnum)) goto S20; |
|
1428
|
36 |
4 |
if(!(low == high)) goto S30; |
|
1442
|
0 |
36 |
if(!(ign <= maxnow)) goto S40; |
|
1458
|
205 |
16 |
for (i=0, i_nb= -1L; *(str+i); i++) |
|
1459
|
180 |
25 |
if ( *(str+i) != ' ' ) i_nb = i; |
|
1486
|
1062 |
0 |
if (a <= 0 || a >= m || s <= 0 || s >= m) { |
|
|
1062 |
0 |
if (a <= 0 || a >= m || s <= 0 || s >= m) { |
|
|
1062 |
0 |
if (a <= 0 || a >= m || s <= 0 || s >= m) { |
|
|
0 |
1062 |
if (a <= 0 || a >= m || s <= 0 || s >= m) { |
|
1493
|
0 |
1062 |
if (a < h) { |
|
1501
|
529 |
533 |
if (a1 >= h) { /* A2=1 */ |
|
1505
|
127 |
529 |
while (p < 0) { p += m; } |
|
1512
|
1062 |
0 |
if (a1 != 0) { |
|
1516
|
388 |
674 |
if (p > 0) { p -= m; } |
|
1518
|
276 |
1062 |
while (p < 0) { p += m; } |
|
1525
|
279 |
1062 |
while (p < 0) { p += m; } |
|
1530
|
1062 |
0 |
if (a0 != 0) { |
|
1534
|
1060 |
2 |
if (p > 0) { p -= m; } |
|
1536
|
579 |
1062 |
while (p < 0) { p += m; } |
|
1616
|
0 |
16 |
if(lphr < 1) return; |
|
1617
|
189 |
16 |
for(i=0; i<(lphr-1); i++) { |
|
1683
|
3 |
0 |
if(!(p <= 0)) goto S10; |
|
1692
|
6 |
3 |
for(i=2,D2=1,D3=(p+1-i+D2)/D2; D3>0; D3--,i+=D2) *(parm+i-1) = *(meanv+i-2); |
|
1697
|
3 |
0 |
if(!(info != 0)) goto S30; |
|
1710
|
6 |
3 |
for(i=1,D4=1,D5=(p-i+D4)/D4; D5>0; D5--,i+=D4) { |
|
1711
|
9 |
6 |
for(j=i-1; j
|
|
1766
|
6 |
24 |
if(u < 1.0) goto S20; |
|
1768
|
6 |
18 |
if(u > *q1) goto S60; |
|
1777
|
2 |
4 |
if(ustar < umin) umin = ustar; |
|
1779
|
0 |
6 |
if(u > *(q+i-1)) goto S70; |
|
1868
|
18 |
30 |
if(a == aa) goto S10; |
|
1869
|
0 |
30 |
if(a < 1.0) goto S120; |
|
1886
|
20 |
28 |
if(t >= 0.0) return sgamma; |
|
1891
|
22 |
6 |
if(d*u <= t*t*t) return sgamma; |
|
1895
|
2 |
4 |
if(a == aaa) goto S40; |
|
1904
|
4 |
0 |
if(a <= 3.686) goto S30; |
|
1905
|
0 |
0 |
if(a <= 13.022) goto S20; |
|
1932
|
2 |
4 |
if(x <= 0.0) goto S70; |
|
1937
|
0 |
4 |
if(fabs(v) <= 0.25) goto S50; |
|
1946
|
2 |
2 |
if(log(1.0-u) <= q) return sgamma; |
|
1960
|
0 |
18 |
if(t < -0.71874483771719) goto S70; |
|
1965
|
2 |
16 |
if(fabs(v) <= 0.25) goto S80; |
|
1974
|
0 |
18 |
if(q <= 0.0) goto S70; |
|
1975
|
10 |
8 |
if(q <= 0.5) goto S100; |
|
1979
|
8 |
0 |
if(q < 15.0) goto S95; |
|
1986
|
0 |
0 |
if((q+e-0.5*t*t) > 87.4982335337737) goto S115; |
|
1987
|
0 |
0 |
if(c*fabs(u) > exp(q+e-0.5*t*t)) goto S70; |
|
1998
|
14 |
4 |
if(c*fabs(u) > w*exp(e-0.5*t*t)) goto S70; |
|
2020
|
0 |
0 |
if(p >= 1.0) goto S140; |
|
2022
|
0 |
0 |
if(sexpo() < sgamma) goto S130; |
|
2026
|
0 |
0 |
if(sexpo() < (1.0-a)*log(sgamma)) goto S130; |
|
2103
|
44 |
36 |
if(u > 0.5) s = 1.0; |
|
2107
|
0 |
80 |
if(i == 32) i = 31; |
|
2108
|
2 |
78 |
if(i == 0) goto S100; |
|
2115
|
6 |
76 |
if(ustar <= *(t+i-1)) goto S60; |
|
2123
|
44 |
36 |
if(s == 1.0) snorm = -y; |
|
2137
|
2 |
4 |
if(ustar > tt) goto S50; |
|
2139
|
0 |
4 |
if(ustar >= u) goto S70; |
|
2154
|
10 |
2 |
if(u < 1.0) goto S110; |
|
2164
|
2 |
0 |
if(ustar > tt) goto S50; |
|
2166
|
0 |
0 |
if(ustar >= u) goto S150; |
|
2174
|
10 |
8 |
if ( ( sign>0.0f && num<0.0f ) || ( sign<0.0f && num>0.0f ) ) |
|
|
10 |
0 |
if ( ( sign>0.0f && num<0.0f ) || ( sign<0.0f && num>0.0f ) ) |
|
|
8 |
10 |
if ( ( sign>0.0f && num<0.0f ) || ( sign<0.0f && num>0.0f ) ) |
|
|
8 |
0 |
if ( ( sign>0.0f && num<0.0f ) || ( sign<0.0f && num>0.0f ) ) |
|
2186
|
0 |
0 |
if (msg != NULL) fprintf(stderr,"%s\n",msg); |