42     double p[] = {25, 62.5, 12.5};  in test1()  local
44     D d(b, b+Np+1, p);  in test1()  local
54     double s = std::accumulate(prob.begin(), prob.end(), 0.0);  in test1()  local
73             for (I j = lb; j != ub; ++j)  in test1()  local
105     double p[] = {0, 62.5, 12.5};  in test2()  local
107     D d(b, b+Np+1, p);  in test2()  local
117     double s = std::accumulate(prob.begin(), prob.end(), 0.0);  in test2()  local
136             for (I j = lb; j != ub; ++j)  in test2()  local
168     double p[] = {25, 0, 12.5};  in test3()  local
170     D d(b, b+Np+1, p);  in test3()  local
180     double s = std::accumulate(prob.begin(), prob.end(), 0.0);  in test3()  local
199             for (I j = lb; j != ub; ++j)  in test3()  local
231     double p[] = {25, 62.5, 0};  in test4()  local
233     D d(b, b+Np+1, p);  in test4()  local
243     double s = std::accumulate(prob.begin(), prob.end(), 0.0);  in test4()  local
262             for (I j = lb; j != ub; ++j)  in test4()  local
294     double p[] = {25, 0, 0};  in test5()  local
296     D d(b, b+Np+1, p);  in test5()  local
306     double s = std::accumulate(prob.begin(), prob.end(), 0.0);  in test5()  local
325             for (I j = lb; j != ub; ++j)  in test5()  local
357     double p[] = {0, 25, 0};  in test6()  local
359     D d(b, b+Np+1, p);  in test6()  local
369     double s = std::accumulate(prob.begin(), prob.end(), 0.0);  in test6()  local
388             for (I j = lb; j != ub; ++j)  in test6()  local
420     double p[] = {0, 0, 1};  in test7()  local
422     D d(b, b+Np+1, p);  in test7()  local
432     double s = std::accumulate(prob.begin(), prob.end(), 0.0);  in test7()  local
451             for (I j = lb; j != ub; ++j)  in test7()  local
483     double p[] = {75, 25};  in test8()  local
485     D d(b, b+Np+1, p);  in test8()  local
495     double s = std::accumulate(prob.begin(), prob.end(), 0.0);  in test8()  local
514             for (I j = lb; j != ub; ++j)  in test8()  local
546     double p[] = {0, 25};  in test9()  local
548     D d(b, b+Np+1, p);  in test9()  local
558     double s = std::accumulate(prob.begin(), prob.end(), 0.0);  in test9()  local
577             for (I j = lb; j != ub; ++j)  in test9()  local
609     double p[] = {1, 0};  in test10()  local
611     D d(b, b+Np+1, p);  in test10()  local
621     double s = std::accumulate(prob.begin(), prob.end(), 0.0);  in test10()  local
640             for (I j = lb; j != ub; ++j)  in test10()  local
672     double p[] = {1};  in test11()  local
674     D d(b, b+Np+1, p);  in test11()  local
684     double s = std::accumulate(prob.begin(), prob.end(), 0.0);  in test11()  local
703             for (I j = lb; j != ub; ++j)  in test11()  local