The R Project SVN R-packages

Rev

Rev 1862 | Rev 3454 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 1862 Rev 1877
Line 26... Line 26...
26 
    SET_SLOT(val, Matrix_uploSym, duplicate(uploP));
 26 
    SET_SLOT(val, Matrix_uploSym, duplicate(uploP));
27 
    SET_SLOT(val, Matrix_diagSym, mkString("N"));
 27 
    SET_SLOT(val, Matrix_diagSym, mkString("N"));
28 
    SET_SLOT(val, Matrix_DimSym, duplicate(dimP));
 28 
    SET_SLOT(val, Matrix_DimSym, duplicate(dimP));
29 
    SET_SLOT(val, Matrix_xSym, duplicate(GET_SLOT(x, Matrix_xSym)));
 29 
    SET_SLOT(val, Matrix_xSym, duplicate(GET_SLOT(x, Matrix_xSym)));
30 
    F77_CALL(dpptrf)(uplo, dims, REAL(GET_SLOT(val, Matrix_xSym)), &info);
 30 
    F77_CALL(dpptrf)(uplo, dims, REAL(GET_SLOT(val, Matrix_xSym)), &info);
- 
 
 31 
    if (info) {
- 
 
 32 
	if(info > 0) /* e.g. x singular */
- 
 
 33 
	    error(_("the leading minor of order %d is not positive definite"),
- 
 
 34 
		    info);
- 
 
 35 
	else /* should never happen! */
- 
 
 36 
	    error(_("Lapack routine %s returned error code %d"), "dpptrf", info);
- 
 
 37 
    }
31 
    UNPROTECT(1);
 38 
    UNPROTECT(1);
32 
    return set_factors(x, val, "pCholesky");
 39 
    return set_factors(x, val, "pCholesky");
33 
}
 40 
}
34 
 
 41 
 
35 
static
 42 
static