From pdalgd at gmail.com  Tue Jun 18 21:49:09 2013
From: pdalgd at gmail.com (peter dalgaard)
Date: Tue, 18 Jun 2013 21:49:09 +0200
Subject: [Rd] eigen(symmetric=TRUE) for complex matrices
In-Reply-To: <35EFC1D4-CA6B-42B6-A15B-4A1FCD0FD687@xs4all.nl>
References: <CAHHjBM4z8d58QNb1HuqT4-k7BKaTo_=zAeuFeDr4pZdjTP2-Tg@mail.gmail.com>
	<5CEBFA58-9D40-42F0-87B8-F006EC3A848B@gmail.com>
	<267A30A5-1865-4D23-8B9A-F04F9BBED1EE@gmail.com>
	<35EFC1D4-CA6B-42B6-A15B-4A1FCD0FD687@xs4all.nl>
Message-ID: <84517510-5BD0-4842-B361-4F336E24E37A@gmail.com>


On Jun 18, 2013, at 21:23 , Berend Hasselman wrote:

> 
> So it seems that the blocked algorithm is the cause of the error and that using the (possibly slow) unblocked algorithm gives the correct result.

Thanks Berend,

The finger seems to be pointing to the internal libRlapack/libRblas. The apparent pattern is that things work when R is linked against external libraries and not when the internal ones are used. So it could be time to start looking for differences between R's lapack module and the original LAPACK code.

-pd 

-- 
Peter Dalgaard, Professor,
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Email: pd.mes at cbs.dk  Priv: PDalgd at gmail.com



[Rd] eigen(symmetric=TRUE) for complex matrices
peter dalgaard