[Numpy-discussion] numarray dotblas problem on OSX

Robert Kern rkern at ucsd.edu
Thu Apr 28 18:05:30 CDT 2005


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 at ucsd.edu

"In the fields of hell where the grass grows high
  Are the graves of dreams allowed to die."
   -- Richard Harter




More information about the Numpy-discussion mailing list