The R Project SVN R

Rev

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

Rev 77052 Rev 77395
Line 1144... Line 1144...
1144 
	Memcpy(REAL(A), REAL(Ain), (size_t)m * n);
 1144 
	Memcpy(REAL(A), REAL(Ain), (size_t)m * n);
1145 
    }
 1145 
    }
1146 
 
 1146 
 
1147 
    SEXP jpvt = PROTECT(allocVector(INTSXP, n));
 1147 
    SEXP jpvt = PROTECT(allocVector(INTSXP, n));
1148 
    for (int i = 0; i < n; i++) INTEGER(jpvt)[i] = 0;
 1148 
    for (int i = 0; i < n; i++) INTEGER(jpvt)[i] = 0;
1149 
    SEXP tau = PROTECT(allocVector(REALSXP, m < n ? m : n));
 1149 
    SEXP tau = PROTECT(allocVector(REALSXP, m < n ? m : n)); // qraux
1150 
    int info, lwork = -1;
 1150 
    int info, lwork = -1;
1151 
    double tmp;
 1151 
    double tmp;
1152 
    F77_CALL(dgeqp3)(&m, &n, REAL(A), &m, INTEGER(jpvt), REAL(tau),
 1152 
    F77_CALL(dgeqp3)(&m, &n, REAL(A), &m, INTEGER(jpvt), REAL(tau),
1153 
		     &tmp, &lwork, &info);
 1153 
		     &tmp, &lwork, &info);
1154 
    if (info < 0)
 1154 
    if (info < 0)