# | 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 |
> | Changed lines |