Function pagetensorprod.m, paged tensor product

pagetensorprod generalizes tensorprod for paged operation (analogous to pagemtimes).

The operation pagetensorprod(A,B,dimA,dimB,ndimsA,ndimsB) is equivalent to tensorprod(A,B,dimA,dimB,NumDimensionsA=ndimsA) if ndimsA >= ndims(A) and ndimsB >= ndimsB. Otherwise, the tensor product operation is applied elementwise in dimensions ndimsA+[1,2,...] in A and corresponding dimensions ndimsB+[1,2,...] in B, with singleton expansion.