final fixes for ape 2.7-3

[ape.git] / src / rTrait.c

diff --git a/src/rTrait.c b/src/rTrait.c
index 0e543a6d8642f205891e3d0d32a192695f960040..43dc73ad85f3163b65e49333889a2088d30a06a5 100644 (file)
--- a/src/rTrait.c
+++ b/src/rTrait.c
@@ -12,6 +12,7 @@ void rTraitCont(int *model, int *Nedge, int *edge1, int *edge2, double *el,
 {
 /* The tree must be in pruningwise order */
        int i;
+       double alphaT, M, S;
 
        switch(*model) {
        case 1 : for (i = *Nedge - 1; i >= 0; i--) {
@@ -25,7 +26,7 @@ void rTraitCont(int *model, int *Nedge, int *edge1, int *edge2, double *el,
                                alphaT = alpha[i] * el[i];
                                M = exp(-alphaT);
                                S = sigma[i] * sqrt((1 - exp(-2 * alphaT))/(2 * alpha[i]));
-                       } else {
+                       } else { /* same than if (alpha[i] ==0) */
                                M = 1;
                                S = sqrt(el[i]) * sigma[i];
                        }