# matrixfunction

28 Feb 2009 (Updated )

summary

matrixfunn(hFunction,A,B,Dim)
```function Y = matrixfunn(hFunction,A,B,Dim)
% matrixfunn
%   Y = matrixfunn(hFunction,A,B,Dim)
%
%   Invokes hFunction on the columns of matrices A and B.
%   hFunction is a function handle that accepts two vectors.
%   hFunction is invoked aCols*(2*bCols-aCols+1)/2 times.
%   Y is a cell matrix of size([aCols bCols]).
%   If Dim is 2, hFunction operates on the rows of matrices A and B.
%
%   Example
%       A = [1 1 1; 2 2 2; 3 3 3]
%       Y = matrixfunn(@(X,Y)sum(abs(X-Y)),A,A);
%       Y = [0]    [0]    [0]
%            []     []    [0]
%            []     []    [0]
%
%       Y = matrixfunn(@(X,Y)sum(abs(X-Y)),A,A,2);
%       Y = [0]    [-3]    [-6]
%            []      []    [-3]
%            []      []    [ 0]
%
%       matrixfun matrixfunrr matrixfunrc matrixfuncc matrixfuncr matrixfunt
%       matrixfunn matrixfunrrn matrixfunrcn matrixfunccn matrixfuncrn matrixfuntn
%       matrixfuni matrixfunir matrixfunic
%

%% Author Information
%   Smart Systems Integration Group
%	Cork Institute of Technology
%

%% Assign defaults
if nargin<4 || isempty(Dim), Dim = 1; end

%%
nA = size(A.',Dim);     % Matrix size
nB = size(B.',Dim);     % Matrix size
Y = cell([nA nB]);    	% Predefine loop variable
if Dim==2
for i = 1:nA                             	% Loop through each row of A
for j = i:nB                            % Loop through each row of B
Y{i,j} = hFunction(A(i,:),B(j,:));  % Invoke function
end
end
elseif Dim==1
for i = 1:nA                             	% Loop through each column of A
for j = i:nB                            % Loop through each column of B
Y{i,j} = hFunction(A(:,i),B(:,j));  % Invoke function
end
end
end
end

```