Paul Gilbert wrote:
> Peter Dalgaard wrote:
>
>
>> Well, there's make check-all...
>>
> Ok, this is what I should be running. It reports
>
> running tests of LAPACK-based functions
> make[3]: Entering directory `/home/mfa/gilp/toolchain/R/src/R-beta/tests'
> running code in 'lapa
Peter Dalgaard wrote:
> Well, there's make check-all...
Ok, this is what I should be running. It reports
running tests of LAPACK-based functions
make[3]: Entering directory `/home/mfa/gilp/toolchain/R/src/R-beta/tests'
running code in 'lapack.R' ...make[3]: *** [lapack.Rout] Error 1
make
Paul Gilbert wrote:
> Prof Brian Ripley wrote:
>
>> On Wed, 11 Apr 2007, Paul Gilbert wrote:
>>
>>> Hmmm. It is a bit disconcerting that make check passes and I can get
>>> a fairly seriously wrong answer. Perhaps this test could be added to
>>> make check.
>>>
>> Well, you have l
> "PaulG" == Paul Gilbert <[EMAIL PROTECTED]>
> on Wed, 11 Apr 2007 10:51:22 -0400 writes:
PaulG> Hmmm. It is a bit disconcerting that make check passes and I can
get a
PaulG> fairly seriously wrong answer. Perhaps this test could be added to
make
PaulG> check.
Yes.
Prof Brian Ripley wrote:
> On Wed, 11 Apr 2007, Paul Gilbert wrote:
>> Hmmm. It is a bit disconcerting that make check passes and I can get
>> a fairly seriously wrong answer. Perhaps this test could be added to
>> make check.
> Well, you have learnt something new about software engineering! 'm
On Wed, 11 Apr 2007, Paul Gilbert wrote:
Hmmm. It is a bit disconcerting that make check passes and I can get a
fairly seriously wrong answer. Perhaps this test could be added to make
check.
Well, you have learnt something new about software engineering! 'make
check' is supposed to test th
Hmmm. It is a bit disconcerting that make check passes and I can get a
fairly seriously wrong answer. Perhaps this test could be added to make
check. Whether or not the one answer is correct may be questionable, but
there is no question that
prod(eigen(z, symmetric = FALSE, only.values = TRU
I've also tried
on Redhat EL4(x86_64) [but with /usr/local/.. gcc/gfortran 4.1.1]
and FC6 (ix86 & x86_64)
and Ubuntu dapper (x86_64)
and never got the wrong results.
> "BDR" == Prof Brian Ripley <[EMAIL PROTECTED]>
> on Wed, 11 Apr 2007 07:43:38 +0100 (BST) write
All the systems I tried this on give the 'correct' answer, including
x86_64 Linux, FC5 (gcc 4.1.1)
i686 Linux, FC5
ix86 Windows (both gcc 3.4.5 and gcc pre-4.3.0)
Sparc Solaris, with gcc3, gcc4 and SunPro compilers.
Mainly with R 2.5.0 beta, some with R-devel (where the code is unchanged).
We h
Peter Dalgaard wrote:
> Paul Gilbert wrote:
>
>> Here is the example. Pehaps others could check on other platforms. It
>> is only the first eigenvalue that is different. I am relatively sure
>> the old values are correct, since I compare with an alternate
>> calculation using the expansion o
Paul Gilbert wrote:
> Here is the example. Pehaps others could check on other platforms. It is
> only the first eigenvalue that is different. I am relatively sure the
> old values are correct, since I compare with an alternate calculation
> using the expansion of a polynomial determinant.
>
>
>
Here is another check that R-2.4.1 is right. This should give 1.
prod(eigen(z, symmetric = FALSE, only.values = TRUE)$values ) *
prod(eigen(solve(z), symmetric = FALSE, only.values = TRUE)$values )
R-2.4.1 gives
[1] 1+0i
R-2.5.0 gives
[1] 1.01677-0i
Paul
Paul Gilbert wrote:
> Here is the exa
Here is the example. Pehaps others could check on other platforms. It is
only the first eigenvalue that is different. I am relatively sure the
old values are correct, since I compare with an alternate calculation
using the expansion of a polynomial determinant.
z <- t(matrix(c(
0, 0, 0, 0, 0
We are only aware of better behaviour from LAPACK 3.1 (which is what I
suppose you are talking about, that is R compiled with its internal
LAPACK).
But in at least one case that means finding a complex set of eigenvalues
where previously a real one was found.
On Tue, 10 Apr 2007, Paul Gilber
I am having some trouble with a case where eigen in R-beta gives a
different largest value than in previous versions of R. Other values
seem to be the same. Before I spend too much time, is anyone aware of a
problem (symmetric = FALSE, only.values = TRUE).
Paul Gilbert
==
15 matches
Mail list logo
