view general/numerical/matrix/minner.m @ 29:61921dceded1

More documentation on type system.

function C=minner(A,B)
% minner - inner product for arbitrary dimension arrays
%
% ttimes :: [[I,J]], [[J,K]] -> [[I,K]].
%
% ie just like mtimes, but I and K can be row vectors

Adom=size1(A);
Bdom=size1(B);
A=reshape(A,[],Bdom(1));
B=reshape(B,Bdom(1),[]),
C=reshape(A*B,tosize([Adom(1:end-1) Bdom(2:end)]));