[Rd] (PR#2867)
richard.roger at agric.nsw.gov.au
richard.roger at agric.nsw.gov.au
Fri Jul 18 02:31:15 MEST 2003
This is a multipart message in MIME format.
--=_alternative 00812CFCCA256D66_=
Content-Type: text/plain; charset="us-ascii"
Two points in respect to PR#2867.
First, the trivial one:
(1) In Lapack.c, on lines 806 and 812, there are calls of
F77_CALL(dgeqp3). But the error messages on lines on lines 809 and 815
refer to "dqeqp3", i.e., they currently have a "q" where they should have
a "g".
Secondly, the non-trivial one
(2) I have built R-1.7.1 under Solaris 7 on a SunUltrasparc 2 using
gcc-2.95.3. On running make check-all, bug #2867 still rears its head.
Here is the tail of the output from reg-tests-1.Rout.fail
> ## as.hclust.twins() lost labels and more till (incl) 1.6.2
>
>
> ## PR#2867 qr(LAPACK=TRUE) didn't always pivot in 1.7.0
> set.seed(1)
> X <- matrix(rnorm(40),10,4)
> X[,1] <- X[,2]
> (qrx <- qr(X, LAPACK=TRUE))
Error in qr(X, LAPACK = TRUE) : error code -14983444 from Lapack routine
dqeqp3
Execution halted
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 )
Thanks for R!
Regards
Richard E. Roger
Dr. R. E. Roger NSW Agriculture
Spatial Information Officer Systems 161 Kite St
Resource Information Unit Locked Bag 21
ORANGE NSW 2800
ph: (02) 6391 3697 fax: (02) 6391 3740
This message is intended for the addressee named and may contain
confidential information. If you are not the intended recipient or
received it in error, please delete the message and notify sender. Views
expressed are those of the individual sender and are not necessarily the
views of their organisation.
--=_alternative 00812CFCCA256D66_=
Content-Type: text/html; charset="us-ascii"
<br><font size=2 face="Courier New">Two points in respect to PR#2867.</font>
<br>
<br><font size=2 face="Courier New">First, the trivial one:</font>
<br>
<br><font size=2 face="Courier New">(1) In Lapack.c, on lines 806 and 812, there are calls of F77_CALL(dgeqp3). But the error messages on lines on lines 809 and 815 refer to "dqeqp3", i.e., they currently have a "q" where they should have a "g".</font>
<br>
<br><font size=2 face="Courier New">Secondly, the non-trivial one</font>
<br>
<br><font size=2 face="Courier New">(2) I have built R-1.7.1 under Solaris 7 on a SunUltrasparc 2 using gcc-2.95.3. On running make check-all, bug #2867 still rears its head. Here is the tail of the output from reg-tests-1.Rout.fail</font>
<br>
<br><font size=2 face="Courier New">> ## as.hclust.twins() lost labels and more till (incl) 1.6.2</font>
<br><font size=2 face="Courier New">> </font>
<br><font size=2 face="Courier New">> </font>
<br><font size=2 face="Courier New">> ## PR#2867 qr(LAPACK=TRUE) didn't always pivot in 1.7.0</font>
<br><font size=2 face="Courier New">> set.seed(1)</font>
<br><font size=2 face="Courier New">> X <- matrix(rnorm(40),10,4)</font>
<br><font size=2 face="Courier New">> X[,1] <- X[,2]</font>
<br><font size=2 face="Courier New">> (qrx <- qr(X, LAPACK=TRUE))</font>
<br><font size=2 face="Courier New">Error in qr(X, LAPACK = TRUE) : error code -14983444 from Lapack routine dqeqp3</font>
<br><font size=2 face="Courier New">Execution halted</font>
<br>
<br>
<br><font size=2 face="Courier New">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!</font>
<br>
<br><font size=2 face="Courier New">The call:</font>
<br>
<br><font size=2 face="Courier New"> F77_CALL(dgeqp3)(&m, &n, REAL(A), &m, INTEGER(jpvt), REAL(tau),</font>
<br><font size=2 face="Courier New"> &tmp, &lwork, rwork, &info);</font>
<br>
<br><font size=2 face="Courier New">The subroutine:</font>
<br><font size=2 face="Courier New"> SUBROUTINE DGEQP3( M, N, A, LDA, JPVT, TAU, WORK, LWORK, INFO )</font>
<br>
<br>
<br><font size=2 face="sans-serif">Thanks for R!</font>
<br>
<br><font size=2 face="sans-serif">Regards</font>
<br>
<br><font size=2 face="sans-serif">Richard E. Roger</font>
<br>
<br><font size=2 face="sans-serif">Dr. R. E. Roger NSW Agriculture<br>
Spatial Information Officer Systems 161 Kite St<br>
Resource Information Unit Locked Bag 21<br>
ORANGE NSW 2800<br>
ph: (02) 6391 3697 fax: (02) 6391 3740<br>
<br>
This message is intended for the addressee named and may contain confidential information. If you are not the intended recipient or received it in error, please delete the message and notify sender. Views expressed are those of the individual sender and are not necessarily the views of their organisation.</font>
--=_alternative 00812CFCCA256D66_=--
More information about the R-devel
mailing list