similar to: Listing of LAPACK error codes

Using the 64-bit Solaris compilers make check dumps core in La_rs at 73 F77_CALL(dsyev)(jobv, uplo, &n, rx, &n, rvalues, work, &lwork, &info); I can't reproduce it easily, but example(eigen) occasionally stops with Error: abs(sm - V %*% diag(lam) %*% t(V)) < 60 * Meps is not TRUE which might be connected. The tiny rounding errors in example(eigen) aren't

dear list, since matrix manipulations is often of interest in statistical computations, i'd like to get a working example of using Lapack for regression. However, i run into an error. My matrix-lapack-example.c file: #include <R_ext/Lapack.h> void reg(const char* trans, const int* m, const int* n, const int* nrhs, double* a, const int* lda, double* b, const int*

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

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

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.

Dear All, Under Debian GNU/Linux La.svd (with method = "dgesdd") sometimes gives the error "Error in La.svd(data, nu = 0, nv = min(nrow, ncol), method = "dgesdd") : LAPACK routine DGESDD gave error code -12" It seems not to depend on the data per se, but on the relationship between numbers of rows and columns. For example, if the number of columns is 100,

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"

On 2 February 2015 at 10:07, William Dunlap <wdunlap at tibco.com> wrote: <snip> > > If all goes well then > eigen(lastEigenX) > will cause the same error and you or someone on this list can see what > is odd about that matrix (e.g., by looking at its singular values). Preferably *not* this list as this doesn't really seem to be about developing R or with/for

Hi! everybody: We used opus-codec for a VOIP gateway. The GW is running at a UBUNTU server. The opus stream is transcoded to G711 pcmu stream.So there are many opus codecs running simultaneously. We noticed that if there more than 5 streams in. the voice then has notisable glitchs.More streams in, worse voice got. Then we write test code for opus-codec which encode a .pcm file simultaneously.

DeaR list, I happily use eigen() to compute the eigenvalues and eigenvectors of a fairly large matrix (200x200, say), but it seems over-killed as its rank is limited to typically 2 or 3. I sort of remember being taught that numerical techniques can find iteratively decreasing eigenvalues and corresponding orthogonal eigenvectors, which would provide a nice alternative (once I have the

Thank you for 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

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

I want to create a list of functions in a for loop, with the index of the loop appearing explicitly in the function code. After quite a bit of thrashing around I figured out how to do it. Here is a toy example: junk <- vector("list",4) for(i in 1:4) { itmp <- i junk[[i]] <- eval(bquote(function(x){42 + .(itmp)*x})) } So I'm *basically* happy, but there's

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

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

On 02 Feb 2015, at 04:08 , eigen <liguowei1991 at gmail.com> wrote: > Thank you for 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. > How could we? All we know is that you are having trouble running some unspecified code

You can start diagnosing the problem by capturing the matrix that caused eigen() to stop. You can do this in a variety of ways; here is one trace(eigen, quote(lastEigenX <<- x)) After setting the trace, make your offending function call and after the error the global variable 'lastEigenX' will contain the bad matrix. Take the trace off of eigen untrace(eigen) If all goes

(cross-posted on SO: https://stackoverflow.com/questions/78022766) Hi all, I'm trying to compile R as a static library with the -fPIC flag so I can use it within java+JNI (is it only possible ?), but I cannot find the right flags in '.configure' to compile R this way. I tested various flags but I cannot find the correct syntax. for now, my latest attempt was ``` rm -rvf?

Hello List, I have create a small package with C code. In the C code, I call Fortran routines of Lapack. Here is the call in the C function: F77_CALL(dgesvd)(&jobu, &jobvt,&size, &size,A, &size, D,U,&size,V,&size,&work1, &lwork,&error); In the src directory, I have a Makevars File with: PKG_LIBS = $(LAPACK_LIBS) $(BLAS_LIBS) $(FLIBS) And in the

Like kile for LaTeX, Linux/KDE's kate editor is an excellent editor for R, with easy code submission to a running R process. Syntax highlighting is good. I have not been able to figure out two things: - how to automatically reformat a line or region of text using good indentation rules (Emacs/ESS make this so easy by just hitting Tab while the cursor is in a line, or highlighting a