fix test.constructionModelesLassoMLE; TODO: det(rho[,,r]) ? numerical errors ? ....

[valse.git] / src / sources / constructionModelesLassoRank.c

diff --git a/src/sources/constructionModelesLassoRank.c b/src/sources/constructionModelesLassoRank.c
index 59be2f7..1115311 100644 (file)
--- a/src/sources/constructionModelesLassoRank.c
+++ b/src/sources/constructionModelesLassoRank.c
@@ -55,6 +55,8 @@ void constructionModelesLassoRank_core(
        //Initialize phi to zero, because unactive variables won't be assigned
        for (int i=0; i<p*m*k*L*Size; i++)
                phi[i] = 0.0;
+       for (int i=0; i<L*Size*2; i++)
+               llh[i] = INFINITY;
 
        //initiate parallel section
        int lambdaIndex;