Array-matrix product (R implementation)

This R implementation is used for operations not covered by the specialized Cpp optimizations (like 4D Z).

Usage

amprod(A, M, k)