Efficient vector and parallel manipulation of tensor products

We present efficient vector and parallel methods formanipulating tensor products of matrices. We consider bothcomputing the matrix-vector product (A<sub>1</sub><font face="Symbol">Ä</font><font face="Symbol">¼</font><fontface="Symbol">Ä</font>A<sub>K</sub>)x and solving the systemof linear equations (A<sub>1</sub><font face="Symbol">Ä</font><font face="Symbol">¼</font><fontface="Symbol">Ä</font>A<sub>K</sub>)x=b. The methods describedare independent of <i>K</i>. We accompany this article with acompanion algorithm which describes an implementation of a completeset of tensor product routines based on LAPACK and the Level 2 and3 Basic Linear Algebra Subprograms (BLAS) which providevectorization and parallelization.