(PR#2867)
richard.roger@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@biostat.ku.dk) FAX: (+45) 35327907