Displaying 6 results from an estimated 6 matches for "la_r".
Did you mean:
lar
2011 Feb 21
2
Segfaults of eigen
...Diag")= logi FALSE
..- attr(*, "Upper")= logi FALSE
..- attr(*, "method")= chr "euclidean"
..- attr(*, "call")= language dist(x = iris[, -5])
> eigen(D)
*** caught segfault ***
address (nil), cause 'unknown'
Traceback:
1: .Call("La_rs", x, only.values, PACKAGE = "base")
2: eigen(D)
All systems are Gentoo systems, i.e. R-2.12.1 and blas-atlas-3.9.23-r4
is installed by compiling the sources. Recompiling R and blas-atlas did
not solve the issue. This issue seems the reason that example(svm)
creates segfaults, too....
2001 Jun 19
5
core dump on 64-bit Solaris (PR#990)
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 erro...
2007 Mar 18
1
eigen returns NAs from a real matrix
...sets" "methods"
[7] "base"
Using EISPACK does NOT return NAs:
> sum(is.na(eigen(eigenBug, EISPACK=TRUE)$vectors))
[1] 0
I traced the problem to the following line in 'eigen':
z <- if (!complex.x)
.Call("La_rs", x, only.values, PACKAGE = "base")
Comments?
Best Wishes,
Spencer Graves
2011 Feb 18
1
segfault during example(svm)
...m> # visualize (classes by color, SV by crosses):
svm> plot(cmdscale(dist(iris[,-5])),
svm+ col = as.integer(iris[,5]),
svm+ pch = c("o","+")[1:150 %in% model$index + 1])
*** caught segfault ***
address (nil), cause 'unknown'
Traceback:
1: .Call("La_rs", x, only.values, PACKAGE = "base")
2: eigen(-x/2, symmetric = TRUE)
3: cmdscale(dist(iris[, -5]))
4: plot(cmdscale(dist(iris[, -5])), col = as.integer(iris[, 5]),
pch = c("o", "+")[1:150 %in% model$index + 1])
5: eval.with.vis(expr, envir, enclos)
6: eval.w...
2004 Dec 10
0
strange gee behavior
...rts execution and never finishes). I was able to track down the
offending 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. Runni...
2005 May 02
14
eigenvalues of a circulant matrix
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