search for: dgeqp3

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...

...eturned which is only occasionally a true rank. Would not it be more consistent to replace the rank in the latter case by something based on the following pseudo code ? d=abs(diag(qr)) rank=sum(d >= d[1]*tol) Here, we rely on the fact column pivoting is activated in the called lapack routine (dgeqp3) and diagonal term in qr matrix are put in decreasing order (according to their absolute values). Serguei. How to reproduce: a=diag(2) a[2,2]=0 qaf=qr(a, LAPACK=FALSE) qaf$rank # shows 1. OK it's the true rank value qat=qr(a, LAPACK=TRUE) qat$rank #shows 2. Bad, it's not the expected val...

...ld not it be more consistent to replace the rank in the latter case by something >> based on the following pseudo code ? >> >> d=abs(diag(qr)) >> rank=sum(d >= d[1]*tol) >> >> Here, we rely on the fact column pivoting is activated in the called lapack routine (dgeqp3) >> and diagonal term in qr matrix are put in decreasing order (according to their absolute values). >> >> Serguei. >> >> How to reproduce: >> >> a=diag(2) >> a[2,2]=0 >> qaf=qr(a, LAPACK=FALSE) >> qaf$rank # shows 1. OK it's the true...

...; >>> based on the following pseudo code ? > >>> > >>> d=abs(diag(qr)) > >>> rank=sum(d >= d[1]*tol) > >>> > >>> Here, we rely on the fact column pivoting is activated in the called lapack routine (dgeqp3) > >>> and diagonal term in qr matrix are put in decreasing order (according to their absolute values). > >>> > >>> Serguei. > >>> > >>> How to reproduce: > >>> > >>> a=diag(2) &gt...

...e rank. > > Would not it be more consistent to replace the rank in the latter case by something > based on the following pseudo code ? > > d=abs(diag(qr)) > rank=sum(d >= d[1]*tol) > > Here, we rely on the fact column pivoting is activated in the called lapack routine (dgeqp3) > and diagonal term in qr matrix are put in decreasing order (according to their absolute values). > > Serguei. > > How to reproduce: > > a=diag(2) > a[2,2]=0 > qaf=qr(a, LAPACK=FALSE) > qaf$rank # shows 1. OK it's the true rank value > qat=qr(a, LAPACK=TRUE...

...t.seed(0) X<-matrix(rnorm(40),10,4);X[,1]<-X[,2] qrx<-qr(X,LAPACK=TRUE) qrx$pivot # note, no pivoting qr.R(qrx) # note zero on leading diagonal is not in trailing block 3. ?qr includes the statement: Using LAPACK (including in the complex case) uses column pivoting 4. If I call F77_NAME(dgeqp3) directly I get correct pivoting. platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 7.0 year 2003 month 04 day 16...

...n the latter case by something >>> based on the following pseudo code ? >>> >>> d=abs(diag(qr)) >>> rank=sum(d >= d[1]*tol) >>> >>> Here, we rely on the fact column pivoting is activated in the called lapack routine (dgeqp3) >>> and diagonal term in qr matrix are put in decreasing order (according to their absolute values). >>> >>> Serguei. >>> >>> How to reproduce: >>> >>> a=diag(2) >>> a[2,2]=0 >>>...

man for `qr` says that the function uses LINPACK's DQRDC, while it in fact uses DQRDC2. ``` The QR decomposition of the matrix as computed by LINPACK or LAPACK. The components in the returned value correspond directly to the values returned by DQRDC/DGEQP3/ZGEQP3 ```

Hello. I am trying to invert a matrix, and I am finding that I can get different answers depending on whether I set LAPACK true or false using "qr". I had understood that LAPACK is, in general more robust and faster than LINPACK, so I am confused as to why I am getting what seems to be invalid answers. The matrix is ostensibly the Hessian for a function I am optimizing. I want to get

...(related to PR#3324) o Removed warning about names in persistent strings when a namespace is saved. o Fixed some malformed error messages in the methods package. o pipes were not opening properly when profiling on a Mac. (PR#1140) o Lapack error messages (PR#3494) and call to DGEQP3 (PR#2867) are corrected. o Rd conversion was limiting a file to 1000 pairs of braces, without any warning. Now the limit is 10000, with a warning. (PR#3400) o In the tcltk package, the tkimage.*() commands were defined nonsensically as widget commands. They have been redefined to be...

...(related to PR#3324) o Removed warning about names in persistent strings when a namespace is saved. o Fixed some malformed error messages in the methods package. o pipes were not opening properly when profiling on a Mac. (PR#1140) o Lapack error messages (PR#3494) and call to DGEQP3 (PR#2867) are corrected. o Rd conversion was limiting a file to 1000 pairs of braces, without any warning. Now the limit is 10000, with a warning. (PR#3400) o In the tcltk package, the tkimage.*() commands were defined nonsensically as widget commands. They have been redefined to be...