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

[José Luis García Pallero]

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

-- 
*****************************************
José Luis García Pallero
jgpallero at gmail.com
(o<
/ / \
V_/_
Use Debian GNU/Linux and enjoy!
*****************************************