- case 1 : PMAT_JC69(l1, coef_gamma[i], P1);
- PMAT_JC69(l2, coef_gamma[i], P2);
- PMAT_JC69(l3, coef_gamma[i], P3); break;
- case 2 : PMAT_K80(l1, coef_gamma[i], u[0], P1);
- PMAT_K80(l2, coef_gamma[i], u[0], P2);
- PMAT_K80(l3, coef_gamma[i], u[0], P3); break;
- case 3 : PMAT_F81(l1, coef_gamma[i], BF, P1);
- PMAT_F81(l2, coef_gamma[i], BF, P3);
- PMAT_F81(l3, coef_gamma[i], BF, P3); break;
- case 4 : PMAT_F84(l1, coef_gamma[i], u[0], BF, P1);
- PMAT_F84(l2, coef_gamma[i], u[0], BF, P2);
- PMAT_F84(l3, coef_gamma[i], u[0], BF, P3); break;
- case 5 : PMAT_HKY85(l1, coef_gamma[i], u[0], BF, P1);
- PMAT_HKY85(l2, coef_gamma[i], u[0], BF, P2);
- PMAT_HKY85(l3, coef_gamma[i], u[0], BF, P3); break;
- case 6 : PMAT_T92(l1, coef_gamma[i], u[0], BF, P1);
- PMAT_T92(l2, coef_gamma[i], u[0], BF, P2);
- PMAT_T92(l3, coef_gamma[i], u[0], BF, P3); break;
- case 7 : PMAT_TN93(l1, coef_gamma[i], u[0], u[1], BF, P1);
- PMAT_TN93(l2, coef_gamma[i], u[0], u[1], BF, P2);
- PMAT_TN93(l3, coef_gamma[i], u[0], u[1], BF, P3);break;
- case 8 : PMAT_GTR(l1, coef_gamma[i], u[0], u[1], u[2], u[3], u[4], BF, P1);
- PMAT_GTR(l2, coef_gamma[i], u[0], u[1], u[2], u[3], u[4], BF, P2);
- PMAT_GTR(l3, coef_gamma[i], u[0], u[1], u[2], u[3], u[4], BF, P3); break;
+ case 1 : PMAT_JC69(el, coef_gamma[i], P); break;
+ case 2 : PMAT_K80(el, coef_gamma[i], u[0], P); break;
+ case 3 : PMAT_F81(el, coef_gamma[i], BF, P); break;
+ case 4 : PMAT_F84(el, coef_gamma[i], u[0], BF, P); break;
+ case 5 : PMAT_HKY85(el, coef_gamma[i], u[0], BF, P); break;
+ case 6 : PMAT_T92(el, coef_gamma[i], u[0], BF, P); break;
+ case 7 : PMAT_TN93(el, coef_gamma[i], u[0], u[1], BF, P); break;
+ case 8 : PMAT_GTR(el, coef_gamma[i], u[0], u[1], u[2], u[3], u[4], BF, P); break;