
##S = ones((7,4,3),Float32)
##XC = zeros((7,7),Float32)
##YC = zeros((4,4),Float32)
##
##print '\nS forward'
##print S.shape
##print S
##
##print '\nS lattice reversed'
##St = transpose(S,(1,0,2))
##print St.shape
##print St
##
##MX = matrixmultiply(XC,St)
##
##print '\nMX'
##print MX.shape
##print MX
##
##MY = matrixmultiply(YC,S)
##MY = transpose(MY,(1,0,2))
##
##print '\nMY'
##print MY.shape
##print MY
