search for: dsyevr

...ipley 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 am interested about the meaning and handling of error codes 1 and 2. In Lapack.c I only see the reference to the variable info in certain Fortran code: F77_CALL(dsyevr)(jobv, range, uplo, &n, rx, &n, &vl, &vu, &il, &iu, &abstol, &m, rvalues,...

...ute 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" (faster? stabler?) for both complete and partial eigenproblems than dsyevd/dsyevx, but only the complete eigenproblem interface to dsyevr app...

...ute 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" (faster? stabler?) for both complete and partial eigenproblems than dsyevd/dsyevx, but only the complete eigenproblem interface to dsyevr app...

Hi, I got an error message in my program saying "Error in eigen(gene_intersection.kernel) : error code 1 from Lapack routine 'dsyevr' Execution halted". As you see, I was trying to compute the eigenvalues of a matrix but got this error. Is there anyone who knows what this error means and how I can fix it? Theoretically the eigenvalues should be nonnegative, if it helps. Thank you! -- View this message in contex...

...d Null function to calculate eigenvectors and nearPD to get > > approximate positive definite matrix first but they also had errors. > > > > > > > > -- > > View this message in context: > > > http://r.789695.n4.nabble.com/error-code-1-from-Lapack-routine-dsyevr-tp4702571p4702639.html > > Sent from the R devel mailing list archive at Nabble.com. > > > > ______________________________________________ > > R-devel at r-project.org mailing list > > https://stat.ethz.ch/mailman/listinfo/r-devel > > > > [[alterna...

...your reply. Do you have any idea of how to get rid of the errors? I tried Null function to calculate eigenvectors and nearPD to get approximate positive definite matrix first but they also had errors. -- View this message in context: http://r.789695.n4.nabble.com/error-code-1-from-Lapack-routine-dsyevr-tp4702571p4702639.html Sent from the R devel mailing list archive at Nabble.com.

On 31/01/2015 14:15, eigen wrote: > Hi, > > I got an error message in my program saying > > "Error in eigen(gene_intersection.kernel) : > error code 1 from Lapack routine 'dsyevr' > Execution halted". > > As you see, I was trying to compute the eigenvalues of a matrix but got this > error. Is there anyone who knows what this error means and how I can fix it? The help page for eigen says Source: By default ?eigen? uses the LAPACK routines ?DSYEV...

While using the rmvnorm function, I get the error: Error in eigen(sigma, sym = TRUE) : error code 5 from Lapack routine 'dsyevr' The same thing happens when I try the eigen() function on my covariance matrix. The matrix is a symmetric 111x111 matrix. Well, it is almost symmetric; there are slight deviations from symmetry (the largest is 3e-18). I have this in an MCMC loop, and it happens about once in every 40 ite...

...pecified platform. It is a bit like someone calling from Finland and telling you that they are lost looking for the University and are standing at a roadsign pointing to "Keskus". > -- > View this message in context: http://r.789695.n4.nabble.com/error-code-1-from-Lapack-routine-dsyevr-tp4702571p4702639.html > Sent from the R devel mailing list archive at Nabble.com. > > ______________________________________________ > R-devel at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-devel -- Peter Dalgaard, Professor, Center for Statistics, Copenh...

...ow to get rid of the > errors? I tried Null function to calculate eigenvectors and nearPD to get > approximate positive definite matrix first but they also had errors. > > > > -- > View this message in context: > http://r.789695.n4.nabble.com/error-code-1-from-Lapack-routine-dsyevr-tp4702571p4702639.html > Sent from the R devel mailing list archive at Nabble.com. > > ______________________________________________ > R-devel at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-devel > [[alternative HTML version deleted]]

...R wrapper. I'm not experienced enough to try and make my own interface with Fortran code, so here are two questions: - is this option (choosing a desired number of eigenvectors) already implemented in some function / package that I missed? - is the "range of indices" option in DSYEVR.f < http:// www.netlib.org/lapack/double/dsyevr.f > what I think, the indices of the desired eigenvalues ordered from the highest to lowest? Many thanks in advance for any piece of advice, Sincerely, Baptiste dummy example if needed: test <- matrix(c(1, 2, 0, 4, 5, 6, 1.00001, 2, 0...

Dear all R users, Can anyone tell me to calculate Eigen value of any real symmetric matrix which algorithm R uses? Is it Jacobi method ? If not is it possible to get explicit algorithm for calculating it? Thanks and regards, Arun [[alternative HTML version deleted]]

I noticed that there is some redundant text in eigen.Rd (both the patched and the devel version). "whereas" appears to be a leftover of previous versions. \source{ By default \code{eigen} uses the LAPACK routines \code{DSYEVR}, \code{DGEEV}, \code{ZHEEV} and \code{ZGEEV} whereas LAPACK is from \url{http://www.netlib.org/lapack} and its guide is listed in the references. } Berend

Hello everyone. I am doing a logistic gam (package mgcv) on a pretty large dataframe (130.000 cases with 100 variables). Because of that, the gam is fitted on a random subset of 10000. Now when I want to predict the values for the rest of the data, I get the following error: > gam.basis_alleakti.1.pr=predict(gam.basis_alleakti.1, +

...data and determine that this particular dataset results in a correlation matrix that is not positive definite and a line within gee with eigen(z$wcor) is where things lock up . Specifically, the line within eigen .Call("La_rs", x, only.values, if (capabilities("IEEE754")) "dsyevr" else "dsyev", PACKAGE = "base") is where it happens. The issue is that I get different behavior when I run the function "interactively" and when it's run by sourcing a script. Running at a command line (i.e. not embedding the call within functions and sourc...

Dear all I am so glad the R can provide the efficient calculate about eigenvector and eigenvalue. However, i have some puzzle about the procedure of eigen. Fristly, what kind of procedue does the R utilize such that the eigen are obtained? For example, A=matrix(c(1,2,4,3),2,2) we can define the eigenvalue lamda, such as det | 1-lamda 4 | =0 | 2 3-lamda | then

Hi, It is my understanding that the eigenvectors of a circulant matrix are given as follows: 1,omega,omega^2,....,omega^{p-1} where the matrix has dimension given by p x p and omega is one of p complex roots of unity. (See Bellman for an excellent discussion on this). The matrix created by the attached row and obtained using the following commands indicates no imaginary parts for the

--Apple-Mail-1-507788993 Content-Transfer-Encoding: 7bit Content-Type: text/plain; charset=US-ASCII; format=flowed Here is the same thing as before, now done a bit more systematically. The rather bad performance of R-patched on eigen is a mystery (to me), the gain in La.eigen when using ATLAS is probably not Altivec, but using threaded level 3 BLAS on a dual processor machine.