From richard.roger at agric.nsw.gov.au  Fri Jul 18 02:31:15 2003
From: richard.roger at agric.nsw.gov.au (richard.roger@agric.nsw.gov.au)
Date: Fri Jul 18 01:31:09 2003
Subject: [Rd] (PR#2867)
Message-ID: <200307172331.h6HNVFJH023196@pubhealth.ku.dk>

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). &nbsp;But the error messages on lines on lines 809 and 815 refer to &quot;dqeqp3&quot;, i.e., they currently have a &quot;q&quot; where they should have a &quot;g&quot;.</font>
<br>
<br><font size=2 face="Courier New">Secondly, the non-trivial one</font>
<br>
<br><font size=2 face="Courier New">(2) &nbsp;I have built R-1.7.1 under Solaris 7 on a SunUltrasparc 2 using gcc-2.95.3. &nbsp;On running make check-all, bug #2867 still rears its head. &nbsp;Here is the tail of the output from reg-tests-1.Rout.fail</font>
<br>
<br><font size=2 face="Courier New">&gt; ## as.hclust.twins() lost labels and more till (incl) 1.6.2</font>
<br><font size=2 face="Courier New">&gt; </font>
<br><font size=2 face="Courier New">&gt; </font>
<br><font size=2 face="Courier New">&gt; ## PR#2867 qr(LAPACK=TRUE) didn't always pivot in 1.7.0</font>
<br><font size=2 face="Courier New">&gt; set.seed(1)</font>
<br><font size=2 face="Courier New">&gt; X &lt;- matrix(rnorm(40),10,4)</font>
<br><font size=2 face="Courier New">&gt; X[,1] &lt;- X[,2]</font>
<br><font size=2 face="Courier New">&gt; (qrx &lt;- 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. &nbsp;The call has an extra argument &quot;rwork&quot;. &nbsp;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">&nbsp; &nbsp; F77_CALL(dgeqp3)(&amp;m, &amp;n, REAL(A), &amp;m, INTEGER(jpvt), REAL(tau),</font>
<br><font size=2 face="Courier New">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&amp;tmp, &amp;lwork, rwork, &amp;info);</font>
<br>
<br><font size=2 face="Courier New">The subroutine:</font>
<br><font size=2 face="Courier New">&nbsp; &nbsp; &nbsp; 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 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;NSW Agriculture<br>
Spatial Information Officer Systems &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 161 Kite St<br>
Resource Information Unit &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;Locked Bag 21<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;ORANGE NSW 2800<br>
ph: (02) 6391 3697 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;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_=--


[Rd] (PR#2867)
richard.roger@agric.nsw.gov.au