search for: dqrdc2

...of symbol `bakslv' changed from 120 to 338 in uncmin.o ../lib/libappl.a(dqrdc.o): In function `dqrdc': dqrdc.o(.text+0x3d6): undefined reference to `dnrm2' dqrdc.o(.text+0x6cf): undefined reference to `dnrm2' dqrdc.o(.text+0xc48): undefined reference to `dnrm2' ../lib/libappl.a(dqrdc2.o): In function `dqrdc2': dqrdc2.o(.text+0x98): undefined reference to `dnrm2' dqrdc2.o(.text+0x516): undefined reference to `dnrm2' ../lib/libappl.a(dqrdc2.o)(.text+0xa58): more undefined references to `dnrm2' follow make[3]: Leaving directory `/home/fredrigl/R-0.62.3/src/main'...

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

...slv.o binning.o chisqsim.o chull.o c poly.o cumsum.o distance.o fft.o fmin.o fortran.o lbfgsb.o loglin.o lowess.o mac har.o massdist.o pretty.o pythag.o rowsum.o sock.o splines.o stem.o strsignif.o tabulate.o uncmin.o zeroin.o ch2inv.o chol.o dpbfa.o dpbsl.o dpoco.o dpodi.o dpo fa.o dposl.o dqrdc.o dqrdc2.o dqrls.o dqrsl.o dqrutl.o dsvdc.o dtrco.o dtrsl.o e igen.o lminfl.o blas.o make[3]: cr: Command not found make[3]: *** [libappl.a] Error 127 make[3]: Leaving directory `/home/ppp/paradis/R/R-1.0.1/src/appl' make[2]: *** [R] Error 2 make[2]: Leaving directory `/home/ppp/paradis/R/R-1.0.1/src/ap...

I found the problem but not a solution. It turns out if I add the following lines to dqrdc2.f I get the error: write(*,300) ldx,n,p 300 format(3i4) I don't get a compile error but I get the seemingly unrelated error in linking R.DLL I guess the question now is, "How do I add a simple print statement?". Or, what is wrong with the above print statement? Thank yo...

...>> >> qr() has a tolerance argument which by default is tol=1e-07. > > > But the tolerance is not about the size of the values! It is about the > reduction of the norm during the householder transformation. From the > sources <https://svn.r-project.org/R/trunk/src/appl/dqrdc2.f> > > > c cycle the columns from l to p left-to-right until one > c with non-negligible norm is located. a column is considered > c to have become negligible if its norm has fallen below > c tol times its original norm.... > > Best, > Uwe Ligges Oh,...

2017-06-22 19:49 GMT+02:00 Uwe Ligges <ligges at statistik.tu-dortmund.de>: > On 22.06.2017 17:11, Bernd Funovits wrote: >> >> Hello, >> >> I experienced some unexpected behaviour while determining the rank of matrices (sometimes 1x1 matrices): >> base::qr(matrix(1e-20))$rank returns 1 (incorrect) >> base::qr(diag(c(1, 1e-20)))$rank returns 2

...o(.text+0x24e8): undefined reference to `__d_sign' ../appl/libappl.a(dpoco.o): In function `dpoco_': dpoco.o(.text+0x534): undefined reference to `__d_sign' ../appl/libappl.a(dqrdc.o): In function `dqrdc_': dqrdc.o(.text+0xce0): undefined reference to `__d_sign' ./appl/libappl.a(dqrdc2.o): In function `dqrdc2_': dqrdc2.o(.text+0xb30): undefined reference to `__d_sign' ../appl/libappl.a(dsvdc.o): In function `dsvdc_': dsvdc.o(.text+0x52c): undefined reference to `__d_sign' ../appl/libappl.a(dsvdc.o)(.text+0xd58): more undefined references to `__d_sign' follow ....

...ing.o chisqsim.o chull.o cpoly.o cumsum.o distance.o fft.o fmin.o fortran.o lbfgsb.o loglin.o lowess.o machar.o maxcol.o massdist.o pretty.o rowsum.o sock.o splines.o stem.o strsignif.o tabulate.o uncmin.o zeroin.o ch2inv.o chol.o dpbfa.o dpbsl.o dpoco.o dpodi.o dpofa.o dposl.o dqrdc.o dqrdc2.o dqrls.o dqrsl.o dqrutl.o dsvdc.o dtrco.o dtrsl.o eigen.o lminfl.o ranlib libappl.a #make: No such file or directory *** Exit 1 Stop. *** Exit 1 Stop. *** Exit 1 Stop. fenchel:/net/prod/appl/R/R-1.2.1/> The same problem appears when I try to compile v.1.2.0. What is wrong? Best regards Kl...

...0 -0.4082483 -0.4082483 [3,] 0.5773503 0.7071068 0.000000 -0.7071068 0.7071068 $rank [1] 2 $qraux [1] 1.5773503 1.7071068 0.0000000 0.7071068 0.7071068 $pivot [1] 1 3 2 4 5 It seems that there is a problem here (and S gets this right). I think the changes to the Linpack pivoting strategy in dqrdc2 fail in this example. -- Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272860 (secr) Oxford OX1 3TG, UK...

...ths.o binning.o chull.o cpoly.o cumsum.o fft.o fmin.o integrate.o interv.o lbfgsb.o loglin.o lowess.o machar.o maxcol.o massdist.o pretty.o rcont.o rowsum.o splines.o stem.o strsignif.o tabulate.o uncmin.o zeroin.o ch2inv.o chol.o dchdc.o dpbfa.o dpbsl.o dpoco.o dpodi.o dpofa.o dposl.o dqrdc.o dqrdc2.o dqrls.o dqrsl.o dqrutl.o dsvdc.o dtrco.o dtrsl.o eigen.o lminfl.o blas.o zgemm.o make[3]: *** [libappl.a] Error 1 make[3]: Leaving directory `/home/gruen/Applications/R-2.0.1/src/appl' make[2]: *** [R] Error 2 make[2]: Leaving directory `/home/gruen/Applications/R-2.0.1/src/appl' make[1...

...cc -g -O2 -I../include -c distance.c f77 -c dpoco.f -o dpoco.o dpoco.f: dpoco: f77 -c dpodi.f -o dpodi.o dpodi.f: dpodi: f77 -c dpofa.f -o dpofa.o dpofa.f: dpofa: f77 -c dposl.f -o dposl.o dposl.f: dposl: f77 -c dqrdc.f -o dqrdc.o dqrdc.f: dqrdc: f77 -c dqrdc2.f -o dqrdc2.o dqrdc2.f: dqrdc2: f77 -c dqrls.f -o dqrls.o dqrls.f: dqrls: f77 -c dqrsl.f -o dqrsl.o dqrsl.f: dqrsl: f77 -c dqrutl.f -o dqrutl.o dqrutl.f: dqrqty: dqrqy: dqrcf: dqrrsd: dqrxb: f77 -c dsvdc.f -o dsvdc.o dsvdc.f:...

...t; Best, >> Uwe Ligges > > qr() has a tolerance argument which by default is tol=1e-07. But the tolerance is not about the size of the values! It is about the reduction of the norm during the householder transformation. From the sources <https://svn.r-project.org/R/trunk/src/appl/dqrdc2.f> c cycle the columns from l to p left-to-right until one c with non-negligible norm is located. a column is considered c to have become negligible if its norm has fallen below c tol times its original norm.... Best, Uwe Ligges > > Regards, > I?aki >

Yesterday I ported FLAME to OS X (http://www.cs.utexas.edu/users/flame/index.html). FLAME is a system to automatically (currently semi-automatically) generate efficient and correct linear algebra routines in C or MATLAB or MPI from descriptions of the algorithms. It uses code that relies heavily on partitioning of matrices, and avoids indices as much as possible. There are a number of

Dear R Developers, I should have had the Design library running in R by now but have kept putting off changing some calls to LINPACK routines to use those builtin to R. Specifically I call dqrsl1 and dqr. Would it be an easy task to put those in the next release of R? If not I'll finally bite the bullet and get back into reading LINPACK documentation (which I have but haven't examined

...act -g -I../../src/include -c cpoly.c cc -ieee_with_inexact -g -I../../src/include -c cumsum.c cc -ieee_with_inexact -g -I../../src/include -c distance.c f77 -c dpoco.f -o dpoco.o f77 -c dpodi.f -o dpodi.o f77 -c dpofa.f -o dpofa.o f77 -c dposl.f -o dposl.o f77 -c dqrdc.f -o dqrdc.o f77 -c dqrdc2.f -o dqrdc2.o f77 -c dqrls.f -o dqrls.o f77 -c dqrsl.f -o dql.o f77 -c dqrutl.f -o dqrutl.o f77 -c dsvdc.f -o dsvdc.o f77 -c dtrsl.f -o dtrsl.o cc -ieee_with_inexact -g -I../../src/include -c eigen.c cc -ieee_with_inexact -g -I../../src/include -c fortran.c cc -ieee_with_inexact -g -I../../s...

...um.o gcc -g -I../include -c distance.c -o distance.o gcc -g -I../include -c dpoco.c -o dpoco.o gcc -g -I../include -c dpodi.c -o dpodi.o gcc -g -I../include -c dpofa.c -o dpofa.o gcc -g -I../include -c dposl.c -o dposl.o gcc -g -I../include -c dqrdc.c -o dqrdc.o gcc -g -I../include -c dqrdc2.c -o dqrdc2.o gcc -g -I../include -c dqrls.c -o dqrls.o gcc -g -I../include -c dqrsl.c -o dqrsl.o gcc -g -I../include -c dqrutl.c -o dqrutl.o gcc -g -I../include -c dsvdc.c -o dsvdc.o gcc -g -I../include -c dtrsl.c -o dtrsl.o gcc -g -I../include -c eigen.c -o eigen.o gcc -g -I../include...

...y.c -o cpoly.o gcc -g -O2 -I../include -I../../src/include -c cumsum.c -o cumsum.o gcc -g -O2 -I../include -I../../src/include -c distance.c -o distance.o f77 -c dpoco.f -o dpoco.o f77 -c dpodi.f -o dpodi.o f77 -c dpofa.f -o dpofa.o f77 -c dposl.f -o dposl.o f77 -c dqrdc.f -o dqrdc.o f77 -c dqrdc2.f -o dqrdc2.o f77 -c dqrls.f -o dqrls.o f77 -c dqrsl.f -o dqrsl.o f77 -c dqrutl.f -o dqrutl.o f77 -c dsvdc.f -o dsvdc.o f77 -c dtrco.f -o dtrco.o f77 -c dtrsl.f -o dtrsl.o gcc -g -O2 -I../include -I../../src/include -c eigen.c -o eigen.o gcc -g -O2 -I../include -I../../src/include -c fortra...

...y.c -o cpoly.o gcc -g -O2 -I../include -I../../src/include -c cumsum.c -o cumsum.o gcc -g -O2 -I../include -I../../src/include -c distance.c -o distance.o f77 -c dpoco.f -o dpoco.o f77 -c dpodi.f -o dpodi.o f77 -c dpofa.f -o dpofa.o f77 -c dposl.f -o dposl.o f77 -c dqrdc.f -o dqrdc.o f77 -c dqrdc2.f -o dqrdc2.o f77 -c dqrls.f -o dqrls.o f77 -c dqrsl.f -o dqrsl.o f77 -c dqrutl.f -o dqrutl.o f77 -c dsvdc.f -o dsvdc.o f77 -c dtrco.f -o dtrco.o f77 -c dtrsl.f -o dtrsl.o gcc -g -O2 -I../include -I../../src/include -c eigen.c -o eigen.o gcc -g -O2 -I../include -I../../src/include -c fortra...

Recent sensible changes to the dynload mechanism have made an old problem resurface: how should we deal with packages which contain Fortran code and may need to be linked against additional libraries such as -lf2c? The current consensus is that extra Fortran libraries maybe needed are handled via the make variable FLIBS, and that `-lf2c' or `-lg2c' are added by default if g77 is used.

...7 -g -O2 -c chol.f -o chol.o g77 -g -O2 -c dpbfa.f -o dpbfa.o g77 -g -O2 -c dpbsl.f -o dpbsl.o g77 -g -O2 -c dpoco.f -o dpoco.o g77 -g -O2 -c dpodi.f -o dpodi.o g77 -g -O2 -c dpofa.f -o dpofa.o g77 -g -O2 -c dposl.f -o dposl.o g77 -g -O2 -c dqrdc.f -o dqrdc.o g77 -g -O2 -c dqrdc2.f -o dqrdc2.o g77 -g -O2 -c dqrls.f -o dqrls.o g77 -g -O2 -c dqrsl.f -o dqrsl.o g77 -g -O2 -c dqrutl.f -o dqrutl.o g77 -g -O2 -c dsvdc.f -o dsvdc.o g77 -g -O2 -c dtrco.f -o dtrco.o g77 -g -O2 -c dtrsl.f -o dtrsl.o g77 -g -O2 -c eigen.f -o eigen.o g77 -g -O2 -c lminfl.f -o lm...