The R Project SVN R-packages

Rev

Rev 4529 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 4529 Rev 4560
Line 147... Line 147...
147 
    AZERO(vx, n * n);
 147 
    AZERO(vx, n * n);
148 
    F77_CALL(dlacpy)(uplo, &n, &n, REAL(GET_SLOT(x, Matrix_xSym)), &n, vx, &n);
 148 
    F77_CALL(dlacpy)(uplo, &n, &n, REAL(GET_SLOT(x, Matrix_xSym)), &n, vx, &n);
149 
    perm = INTEGER(ALLOC_SLOT(val, Matrix_permSym, INTSXP, n));
 149 
    perm = INTEGER(ALLOC_SLOT(val, Matrix_permSym, INTSXP, n));
150 
    F77_CALL(dsytrf)(uplo, &n, vx, &n, perm, &tmp, &lwork, &info);
 150 
    F77_CALL(dsytrf)(uplo, &n, vx, &n, perm, &tmp, &lwork, &info);
151 
    lwork = (int) tmp;
 151 
    lwork = (int) tmp;
152 
    work = Calloc(lwork, double);
 152 
    work = Alloca(lwork, double);
- 
 
 153 
    R_CheckStack();
153 
    F77_CALL(dsytrf)(uplo, &n, vx, &n, perm, work, &lwork, &info);
 154 
    F77_CALL(dsytrf)(uplo, &n, vx, &n, perm, work, &lwork, &info);
154 
    if (info) error(_("Lapack routine dsytrf returned error code %d"), info);
 155 
    if (info) error(_("Lapack routine dsytrf returned error code %d"), info);
155 
    UNPROTECT(1);
 156 
    UNPROTECT(1);
156 
    Free(work);
 - 
 
157 
    return set_factors(x, val, "BunchKaufman");
 157 
    return set_factors(x, val, "BunchKaufman");
158 
}
 158 
}
159 
 
 159 
 
160 
SEXP dsyMatrix_as_dspMatrix(SEXP from)
 160 
SEXP dsyMatrix_as_dspMatrix(SEXP from)
161 
{
 161 
{