libblas3gf crash caused by AMD prefetchw instruction on Intel CPU

[Davy Zhu]

Hi, I run Python3.2/Numpy/Scipy on Ubuntu 12.04, I call spsolve(sparse solver), and libblas crashed python. The reason is libblas invoke an AMD only instruction prefetchw.

Backtrace:
Program terminated with signal 4, Illegal instruction.
#0 0x00007f6f079fe304 in ATL_dupKBmm9_1_1_b1 () from /usr/lib/libblas.so.3gf
(gdb) bt
#0 0x00007f6f079fe304 in ATL_dupKBmm9_1_1_b1 () from /usr/lib/libblas.so.3gf
#1 0x00007f6f0799abbc in ATL_dmmIJK2 () from /usr/lib/libblas.so.3gf
#2 0x00007f6f0799b331 in ATL_dmmIJK () from /usr/lib/libblas.so.3gf
#3 0x00007f6f07988c6b in ATL_dgemm () from /usr/lib/libblas.so.3gf
#4 0x00007f6f079a375d in ATL_dptgemm0 () from /usr/lib/libblas.so.3gf
#5 0x00007f6f0b553e9a in start_thread (arg=0x7f6f03132700)
    at pthread_create.c:308
#6 0x00007f6f0a116cbd in clone ()
    at ../sysdeps/unix/sysv/linux/x86_64/clone.S:112
#7 0x0000000000000000 in ?? ()

(gdb) x/i $pc
=> 0x7f6f079fe304 <ATL_dupKBmm9_1_1_b1+164>: prefetchw (%rsi)

My enviroment: Intel P4 cpu, x86_64, Ubuntu 12.04

[Davy Zhu (zhushenli)]

Pickle file to recreate bug.

[Davy Zhu (zhushenli)]

Pickle file to recreate bug.

[Pauli Virtanen (pauli-virtanen)]

@Davy Zhu: list `dpkg --list|grep -i libatlas` and check that `readlink -f /usr/lib/libblas.so.3gf` points to an Ubuntu-provided BLAS.

[Pauli Virtanen (pauli-virtanen)]

The `libatlas3gf-base_3.8.4-3build1_amd64.deb` package in Precise does contain 3DNow! instructions:

$ objdump -d libblas.so.3gf | grep -c prefetchw
348

[Davy Zhu (zhushenli)]

$ dpkg --list|grep -i libatlas
ii libatlas-base-dev 3.8.4-3build1
ii libatlas-dev 3.8.4-3build1
ii libatlas3gf-base 3.8.4-3build1

$ readlink -f /usr/lib/libblas.so.3gf
/usr/lib/atlas-base/atlas/libblas.so.3gf.0