numarray dotblas problem on OSX
Hi, I have a colleague running Mac OS 10.3, running numarray-1.3.1 (from fink) who has managed to bomb on this little code example:
import numarray as na import numarray.random_array as ra a = ra.random(shape=(257,256)) b = ra.random(shape=(1,256)) na.innerproduct(a, b)
He gets a blas error: ldc must be >= MAX(N,1): ldc=256 N=257Parameter 14 to routine cblas_dgemm was incorrect Mac OS BLAS parameter error in cblas_dgemm, parameter #0, (unavailable), is 0 Simon. -- Simon Burton, B.Sc. Licensed PO Box 8066 ANU Canberra 2601 Australia Ph. 61 02 6249 6940 http://arrowtheory.com
Simon Burton wrote:
Hi,
I have a colleague running Mac OS 10.3, running numarray-1.3.1 (from fink) who has managed to bomb on this little code example:
import numarray as na import numarray.random_array as ra a = ra.random(shape=(257,256)) b = ra.random(shape=(1,256)) na.innerproduct(a, b)
He gets a blas error:
ldc must be >= MAX(N,1): ldc=256 N=257Parameter 14 to routine cblas_dgemm was incorrect Mac OS BLAS parameter error in cblas_dgemm, parameter #0, (unavailable), is 0
On OS X 10.3, numarray 1.3.0, self-compiled for the Apple-installed Python with vecLib as the BLAS, I don't get an error. I don't get a result that's sensible to me, either; I get a (257,1)-shape array with only the first and last entries non-zero. Your colleague might want to reconsider whether he wants innerproduct() or dot(), with the appropriate change of shape for b. -- Robert Kern rkern@ucsd.edu "In the fields of hell where the grass grows high Are the graves of dreams allowed to die." -- Richard Harter
Robert Kern wrote:
Simon Burton wrote:
Hi,
I have a colleague running Mac OS 10.3, running numarray-1.3.1 (from fink) who has managed to bomb on this little code example:
import numarray as na import numarray.random_array as ra a = ra.random(shape=(257,256)) b = ra.random(shape=(1,256)) na.innerproduct(a, b)
He gets a blas error:
ldc must be >= MAX(N,1): ldc=256 N=257Parameter 14 to routine cblas_dgemm was incorrect Mac OS BLAS parameter error in cblas_dgemm, parameter #0, (unavailable), is 0
On OS X 10.3, numarray 1.3.0, self-compiled for the Apple-installed Python with vecLib as the BLAS, I don't get an error.
I don't get a result that's sensible to me, either; I get a (257,1)-shape array with only the first and last entries non-zero.
Oh yes, and apparently a segfault on exit, too. -- Robert Kern rkern@ucsd.edu "In the fields of hell where the grass grows high Are the graves of dreams allowed to die." -- Richard Harter
participants (2)
-
Robert Kern -
Simon Burton