ctprod function

Compute the contracted tensor product between two multiway arrays.