[R-SIG-Mac] eigen() with complex matrices

[Simon Urbanek]

On Feb 22, 2007, at 7:31 AM, Tony Chiang wrote:

> Hi Robin,
>
> I have 2.5 (pretty recent) installed my my mac-mini and I don't get  
> the error. I do not run R-2.4.1 so I cannot see if I can reproduce  
> your error. I build R from source rather than using the universal  
> binary though, so try that and see if that helps. Can anyone with  
> R-2.4.1 reproduce this error?
>

Yes, I can reproduce it with R 2.4.1 on G5 when using vecLib BLAS. It  
doesn't crash when using R's BLAS, so that is a work-around. Given  
that the crash occurs in vecLib and works with with R's BLAS, I'm  
inclined to think that it's a bug in vecLib, but it will require some  
more digging to see if that's really the case...

*Please* when reporting problems include exact information on the  
hardware used - there are simply too many different Mac platforms now...

Thanks,
Simon



>> A <- matrix(c(1,0.2+0.1i,0.2-0.1i,1),2,2)
>> A
>          [,1]     [,2]
> [1,] 1.0+0.0i 0.2-0.1i
> [2,] 0.2+0.1i 1.0+0.0i
>> eigen(A)
> $values
> [1] 1.2236068 0.7763932
>
> $vectors
>                       [,1]                  [,2]
> [1,] -0.7071068+0.0000000i -0.7071068+0.0000000i
> [2,] -0.6324555-0.3162278i  0.6324555+0.3162278i
>
>> sessionInfo()
> R version 2.5.0 Under development (unstable) (2007-02-08 r40671)
> i386-apple-darwin8.8.1
>
> locale:
> C
>
> attached base packages:
> [1] "stats"     "graphics"  "grDevices" "utils"     "datasets"   
> "methods"
> [7] "base"
>>
>
> Cheers,
> --Tony
>
> On 2/22/07, Robin Hankin <r.hankin at noc.soton.ac.uk> wrote:
>>
>> Hello
>>
>> I just now get a bus error with R-2.4.1[Universal binary downloaded
>> from CRAN just now], and also with (recently compiled) R-2.5.0.
>> I have Xcode 2.4.1.
>>
>>
>>
>> Under R-2.4.1:
>>
>>
>> octopus:~/scratch% R --vanilla
>>
>> R version 2.4.1 (2006-12-18)
>> Copyright (C) 2006 The R Foundation for Statistical Computing
>> ISBN 3-900051-07-0
>>
>> [snip]
>>
>> 'help.start()' for an HTML browser interface to help.
>> Type 'q()' to quit R.
>>
>>
>>>
>>> A <- matrix(c(1,0.2+0.1i,0.2-0.1i,1),2,2)
>>> eigen(A)
>>
>> *** caught bus error ***
>> address 0x1, cause 'invalid alignment'
>>
>> Traceback:
>> 1: .Call("La_rs_cmplx", x, only.values, PACKAGE = "base")
>> 2: eigen(A)
>>
>> Possible actions:
>> 1: abort (with core dump)
>> 2: normal R exit
>> 3: exit R without saving workspace
>> 4: exit R saving workspace
>> Selection:
>>
>>
>>
>> the same thing happens with R-2.5.0
>>
>> Anyone?
>>
>>
>>
>>
>>
>>
>>
>> --
>> Robin Hankin
>> Uncertainty Analyst
>> National Oceanography Centre, Southampton
>> European Way, Southampton SO14 3ZH, UK
>>   tel  023-8059-7743
>>
>> _______________________________________________
>> R-SIG-Mac mailing list
>> R-SIG-Mac at stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/r-sig-mac
>>
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-SIG-Mac mailing list
> R-SIG-Mac at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-mac
>
>