Hey Gerard,
Thanks for the suggestion! It took me a while to figure out how to get
this to work. Two things were important: I needed to use the
matrixmultiply() function, and the order of the two matrices being
multiplied is critcial. Here's how I got the example to work.
>>> from Numeric import *
>>> array1 = array([[0.000000, 0.500000, 1.000000],
... [0.500000, 0.000000, 0.500000],
... [1.000000, 0.500000, 0.000000]])
>>> array1
array([[ 0. , 0.5, 1. ],
[ 0.5, 0. , 0.5],
[ 1. , 0.5, 0. ]])
>>> premul = array([[1, 0, 0], [0, 1, 0], [0, 1, 0], [0, 0, 1]])
>>> premul
array([[1, 0, 0],
[0, 1, 0],
[0, 1, 0],
[0, 0, 1]])
>>> res1 = matrixmultiply(premul, array1)
>>> res1
array([[ 0. , 0.5, 1. ],
[ 0.5, 0. , 0.5],
[ 0.5, 0. , 0.5],
[ 1. , 0.5, 0. ]])
>>> postmul = array([[1, 0, 0, 0], [0, 1, 1, 0], [0, 0, 0, 1]])
>>> postmul
array([[1, 0, 0, 0],
[0, 1, 1, 0],
[0, 0, 0, 1]])
>>> array2 = matrixmultiply(res1, postmul)
>>> array2
array([[ 0. , 0.5, 0.5, 1. ],
[ 0.5, 0. , 0. , 0.5],
[ 0.5, 0. , 0. , 0.5],
[ 1. , 0.5, 0.5, 0. ]])
Now I have to figure out how to expand this concept to more complex
cases, and how to generate appropriate pre- and post-multiply
matrices... Hmm...
Thanks so much for getting me started,
Chris
--
http://mail.python.org/mailman/listinfo/python-list
