search for: lapack

Prof. Ripley, > From: Prof Brian Ripley [mailto:ripley at stats.ox.ac.uk] [snipped...] > R cannot currently make use of anyone else's LAPACK routines (because > we've found too many problems with them). > > If you are getting no useful diagnostics, try > dyn.load("RHOME/modules/lapack.so"). That ought to tell you > why you can't > load it. My guess would be to do with libraries, so do an ldd on &g...

Prof. Ripley, > From: Prof Brian Ripley [mailto:ripley at stats.ox.ac.uk] [snipped...] > R cannot currently make use of anyone else's LAPACK routines (because > we've found too many problems with them). > > If you are getting no useful diagnostics, try > dyn.load("RHOME/modules/lapack.so"). That ought to tell you > why you can't > load it. My guess would be to do with libraries, so do an ldd on &g...

...gt;>> - figure out Fortran2003 specification for C/Fortran interoperability > >>> -- this _sounds_ like the right solution, but I don't think many > >>> understand how to use it and what is implied (in particular, will > >>> it require making changes to LAPACK itself?) > >> That would actually be fairly easy. If you declare the subroutines > >> BIND(C), as in > >> > >> subroutine foo(a,b) BIND(C,name="foo_") > >> real a > >> character*1 b > >> end &g...

Full_Name: Ravi Varadhan Version: 2.8.1 OS: Windows Submission from: (NULL) (162.129.251.19) Inverting a matrix with solve(), but using LAPACK=TRUE, gives erroneous results: Here is an example: hilbert <- function(n) { i <- 1:n; 1 / outer(i - 1, i, "+") } h5 <- hilbert(5) hinv1 <- solve(qr(h5)) hinv2 <- solve(qr(h5, LAPACK=TRUE)) all.equal(hinv1, hinv2) # They are not equal Here is a function that I...

...U 3.20GHz R-version: R-2.1.0 I've started using a new Linux server, upgraded at the same time to R-2.1.0 (see above) and have problems with some elementary analysis that ran without a problem on my previous configuration. anova.glm gives the following error: Error in La.chol2inv(x, size) : lapack routines cannot be loaded This was with the after the following configure ./configure --with-readline --prefix=/opt/R-2.1.0 --with-libpng --with-jpeglib --with-pcre --without-x The obvious thing to try next was: ./configure --with-readline --prefix=/opt/R-2.1.0 --with-libpng --with-jpeglib --...

Hi, Why doesn't R include a full Lapack but only a subset? My cda package (now archived) relying on RcppArmadillo has broken multiple times on CRAN over the past few years following updates in the underlying Armadillo library, Every time it follows the same pattern: Armadillo adds a function to solve a specialised linear system more eff...

Professor Ripley commented on LAPACK error codes: https://stat.ethz.ch/pipermail/r-help/2007-March/127702.html and says "Internal LAPACK errors are usually problems with arithmetic accuracy, and as such are compiler- and CPU-specific." Is there a listing for the error codes from Lapack routine 'dsyevr'? Especially I...

...t; - figure out Fortran2003 specification for C/Fortran interoperability >>>>> -- this _sounds_ like the right solution, but I don't think many >>>>> understand how to use it and what is implied (in particular, will >>>>> it require making changes to LAPACK itself?) >>>> That would actually be fairly easy. If you declare the subroutines >>>> BIND(C), as in >>>> >>>> subroutine foo(a,b) BIND(C,name="foo_") >>>> real a >>>> character*1 b >&gt...

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 the parameter correlation...

...Thomas K?nig wrote: > > > - figure out Fortran2003 specification for C/Fortran interoperability > > -- this _sounds_ like the right solution, but I don't think many > > understand how to use it and what is implied (in particular, will > > it require making changes to LAPACK itself?) > > That would actually be fairly easy. If you declare the subroutines > BIND(C), as in > > subroutine foo(a,b) BIND(C,name="foo_") > real a > character*1 b > end > > you will get the calling signature that you already...

R developers, I am trying to understand how symbols are resolved, so that I can configure a package that I contributed to, and so that I can provide guidance to (linux / OSX) users of the package. To be concrete, my package uses the LAPACK Fortran symbol zsysv. This is not in libRlapack, but is defined on my system in the library /usr/lib64/liblapack.so. * I suspect that the reason the symbol is not in libRlapack is just one of economy, i.e., no use for the symbol in R routines, rather than for other nefarious reasons (?? some f...

...he new GotoBLAS2 library. I think the symbolic links and everything are set up correctly, because the BLAS itself works: I'm using Dirk Eddelbuettel's gbcd package to benchmark the difference, and I can still run the matrix multiplication benchmarks (I guess crossprod doesn't depend on lapack?), and I'm getting a 15x speedup over the default BLAS with 8 cores used. The problem, though, is that when I try to run anything that depends on lapack--like qr()--I get the following: Error in qr.default(a, LAPACK = TRUE) : lapack routines cannot be loaded In addition: Warning message: In q...

The following snippet suggests that there is either a bug in qr(,LAPACK=T), or some bug in my understanding. Note that the detected rank is correct (= 2) using the default LINPACK qr, but incorrect (=3) using LAPACK. This is running on Linux Redhat 9.0, using the lapack library that comes with the Redhat distribution. I'm running R 1.7.1 compiled from the sou...

? Wed, 8 May 2024 16:59:25 +0000 Simon Andrews <simon.andrews at babraham.ac.uk> ?????: > The lapack libraries are: > > $ rpm -qa | grep lapack > lapack-3.9.0-10.el9.x86_64 > lapack64_-3.9.0-10.el9.x86_64 > lapack64-3.9.0-10.el9.x86_64 > lapack-devel-3.9.0-10.el9.x86_64 Thanks for this information! I figured out I needed to enable the "PowerTools/CRB" repo, installed...

I am interested in tapping into LAPACK functionality directly from R. Using R-2.4.1 for Windows, I was able to do so ala: dyn.load("bin/Rlapack.dll") is.loaded("dstebz") # returns TRUE N <- 100 NW <- 4 n.tapers <- 5 tpW <- (2 * pi * NW)/N otNmo <- 1:N D <- as.double(cos(tpW) * ((N - 1 - 2 * (0:(N...

I am having trouble using R routines from the RExcel add-in, that use lapack.dll. As an example if I start the R kernel from within Excel and execute "x<-rbind(c(1,2),c(1,-1)) z<-solve(x)", I get the following error: "Error in solve.default(x): lapack routines could not be loaded. In addition: Warning message: unable to load shared library "C...

R-devel has the --with-lapack flag for configure. If you build R with --with-blas="-framework vecLib" --with-lapack="--framework vecLib" then the build goes through, using the native optimized BLAS and Lapack in /System/Library/Frameworks/vecLib.framework. Also, it works, in the sense that it does eigenval...

Le 22/01/2018 ? 17:40, Keith O'Hara a ?crit?: > This behavior is noted in the qr documentation, no? > > rank - the rank of x as computed by the decomposition(*): always full rank in the LAPACK case. For a me a "full rank matrix" is a matrix the rank of which is indeed min(nrow(A), ncol(A)) but here the meaning of "always is full rank" is somewhat confusing. Does it mean that only full rank matrices must be submitted to qr() when LAPACK=TRUE? May be there is a jargon w...

hello, i want to compute the top k eigenvalues+eigenvectors of a (large) real symmetric matrix. since it doesn't look like any top-level R function does this, i'll call LAPACK from a C shlib and then use .Call. the only LAPACK function i see to do this in R_ext/Lapack.h is dsyevx. however, i know that in LAPACK dsyevr can also return a partial eigendecomposition. why is dsyevr not exported in R_ext/Lapack.h? my superficial understanding is that dsyevr is "better&quo...

hello, i want to compute the top k eigenvalues+eigenvectors of a (large) real symmetric matrix. since it doesn't look like any top-level R function does this, i'll call LAPACK from a C shlib and then use .Call. the only LAPACK function i see to do this in R_ext/Lapack.h is dsyevx. however, i know that in LAPACK dsyevr can also return a partial eigendecomposition. why is dsyevr not exported in R_ext/Lapack.h? my superficial understanding is that dsyevr is "better&quo...