[Pkg-octave-devel] Bug#613070: Bug#613070: octave3.2: Error in big matrix multiplication

[Thomas Weber]

Hi, 
On Sat, Feb 12, 2011 at 08:40:15PM +0100, José Luis García Pallero wrote:
> Hello,
> I think the problem is not in octave but is in ATLAS. The simple program:
> 
> #include<cblas.h>
> #define N 52
> int main()
> {
>     double A[N*N],B[N*N],C[N*N];
>     cblas_dgemm(CblasColMajor,CblasNoTrans,CblasNoTrans,N,N,N,1.0,A,N,B,N,1.0,C,N);
>     return 0;
> }
> 
> shows the error: with #define N 51 or less it rins OK, but with
> #define N 52 or more the program crash with the message "Illegal
> instruction". The version of ATLAS is 3.8.2-29 and in Debian GNU/Linux
> runs without problems but in kFreeBSD crashes

I think this bug has been fixed with the newer versions of Atlas we have
in unstable and testing. José, can you test this on your system? Please
note that you have to manually switch the LAPACK library on your system
to Atlas as described in /usr/share/doc/libatlas3-base/README.Debian.

Thanks
	Thomas