| # | Line 374 | Line 374 | init_ptable( | |
|---|---|---|
| 374 | "code error 2 in init_ptable()"); | |
| 375 | frac1 = 1. + prob/((1.-frac)*rowp[p]/rowsum[t] | |
| 376 | + frac*rowp1[p]/rowsum[t+1]); | |
| 377 | – | if (p <= 0 || frac1 > 0.5) |
| 378 | – | frac1 -= 0.5; |
| 379 | – | else if (p >= sntp[1]-1 || frac1 < 0.5) { |
| 380 | – | frac1 += 0.5; |
| 381 | – | --p; |
| 382 | – | } |
| 377 | pvals[i*(nphi+1) + j] = (1.-frac1)*phdiv[p] + | |
| 378 | frac1*phdiv[p+1]; | |
| 379 | } | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines (old) |
| > | Changed lines (new) |