The R Project SVN R-packages

Rev

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

Rev 1535 Rev 1576
Line 59... Line 59...
59 
    int piv = asLogical(pivot);
 59 
    int piv = asLogical(pivot);
60 
    SEXP ans = PROTECT(NEW_OBJECT(MAKE_CLASS("lCholCMatrix")));
 60 
    SEXP ans = PROTECT(NEW_OBJECT(MAKE_CLASS("lCholCMatrix")));
61 
    int j, n = INTEGER(GET_SLOT(x, Matrix_DimSym))[0];
 61 
    int j, n = INTEGER(GET_SLOT(x, Matrix_DimSym))[0];
62 
    int *Xi = INTEGER(GET_SLOT(x, Matrix_iSym)),
 62 
    int *Xi = INTEGER(GET_SLOT(x, Matrix_iSym)),
63 
	*Xp = INTEGER(GET_SLOT(x, Matrix_pSym)),
 63 
	*Xp = INTEGER(GET_SLOT(x, Matrix_pSym)),
64 
	*P, *Pinv, *Parent, *Lp;
 64 
	*P, *Pinv = (int *) NULL, *Parent, *Lp;
65 
    double *D = Calloc(n, double), *Tx, *Xx = Calloc(Xp[n], double);
 65 
    double *D = Calloc(n, double), *Tx, *Xx = Calloc(Xp[n], double);
66 
 
 66 
 
67 
    if (CHAR(asChar(GET_SLOT(x, Matrix_uploSym)))[0] != 'U')
 67 
    if (CHAR(asChar(GET_SLOT(x, Matrix_uploSym)))[0] != 'U')
68 
	error(_("Must have uplo == 'U' in x argument to lsCMatrix_chol"));
 68 
	error(_("Must have uplo == 'U' in x argument to lsCMatrix_chol"));
69 
    SET_SLOT(ans, Matrix_uploSym, mkString("L"));
 69 
    SET_SLOT(ans, Matrix_uploSym, mkString("L"));