[Numpy-discussion] numarray bug: dot product between 2x2 and 3x2x3 on Mac different from PC

[Christopher Barker]

Sebastian Haase wrote:
> Hi!
> This is what I got:

> Can someone confirm this ?

This is what I get on Linux (ix86, Fedora core4, python 2.4.3)

 >>> import numarray as na
 >>> na.__version__
'1.5.0'

 >>> bbb=na.zeros((2,3,3), na.Float32)
 >>> bbb[0,:,:]=1
 >>> bbb[1,:,:]=2
 >>> bbb
array([[[ 1.,  1.,  1.],
         [ 1.,  1.,  1.],
         [ 1.,  1.,  1.]],

        [[ 2.,  2.,  2.],
         [ 2.,  2.,  2.],
         [ 2.,  2.,  2.]]], type=Float32)
 >>> ccc=na.transpose(bbb,(1,0,2))
 >>> ccc
array([[[ 1.,  1.,  1.],
         [ 2.,  2.,  2.]],

        [[ 1.,  1.,  1.],
         [ 2.,  2.,  2.]],

        [[ 1.,  1.,  1.],
         [ 2.,  2.,  2.]]], type=Float32)
 >>> d=na.array([[1,0],[0,1]])
 >>> d
array([[1, 0],
        [0, 1]])
 >>> na.dot(d,ccc)
array([[[ 1.,  1.,  1.],
         [ 2.,  2.,  2.],
         [ 1.,  1.,  1.]],

        [[ 2.,  2.,  2.],
         [ 1.,  1.,  1.],
         [ 2.,  2.,  2.]]], type=Float32)
 >>>

-- 
Christopher Barker, Ph.D.
Oceanographer
                                     		
NOAA/OR&R/HAZMAT         (206) 526-6959   voice
7600 Sand Point Way NE   (206) 526-6329   fax
Seattle, WA  98115       (206) 526-6317   main reception

Chris.Barker at noaa.gov