Peter Dalgaard BSA
p.dalgaard at biostat.ku.dk
Fri Jul 18 01:09:15 MEST 2003
richard.roger at agric.nsw.gov.au writes:
> As an R novice, I am puzzled that the call to dgeqp3 has more arguments
> than dgeqp3 itself. The call has an extra argument "rwork". But this
> could explain the strange error code!
> The call:
> F77_CALL(dgeqp3)(&m, &n, REAL(A), &m, INTEGER(jpvt), REAL(tau),
> &tmp, &lwork, rwork, &info);
> The subroutine:
> SUBROUTINE DGEQP3( M, N, A, LDA, JPVT, TAU, WORK, LWORK, INFO )
I don't think you need to be a novice to be puzzled by that one...
The extra argument would likely cause info not to be modified, and
being uninitialized it can contain random garbage, quite possibly < 0.
But what was the allocation of rwork supposed to be good for then??
O__ ---- Peter Dalgaard Blegdamsvej 3
c/ /'_ --- Dept. of Biostatistics 2200 Cph. N
(*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk) FAX: (+45) 35327907
More information about the R-devel