The implementation of dot operations on higher dimensional matrices, is dependent on numpy.einsum, which was implemented in numpy version 1.6. Would it be okay to accept this dependency? Another approach, would be to investigate if numpy,tensordot could be used instead. That should be from numpy version 1.3. I have not tested if this is possible. Best Troels