function z=iqforma(A,u), z=sum(u.*(A\u),1);
% iqforma - Quadratic form (array version) iqforma(A,x)=diag(x'*inv(A)*x)
%
% iqforma :: [[N,N]], [[N,T]] -> [[1,T]].
