Would anyone be interested in a quaternion version of this for nx4 arrays with nx3 as a special case where the scalar part == 0? Looking at the the cross product implementation, it shouldn't be to hard to duplicate this for quaternions. What should such a product be called? Something like qprod? Chuck