elementwise product of 3d arrays `a` and `b`,"np.einsum('ijk,ikl->ijl', A, B)"