genop.m

Nonlinear dynamical factor analysis Matlab package

function y=genop(f,x1,x2)
% GENOP - Generalized operation
%   
%   C = GENOP(F,A,B) Call function F with exapanded matrices Y and X.
%   The dimensions of the two operands are compared and singleton
%   dimensions in one are copied to match the size of the other.
%   Returns a matrix having dimension lengths equal to
%   MAX(SIZE(A),SIZE(B))
%
% See also GENOPS

% Copyright (C) 2003 Aki Vehtari
%
% This software is distributed under the GNU General Public 
% Licence (version 2 or later); please refer to the file 
% Licence.txt, included with the software, for details.

s1=size(x1);
s2=size(x2);
ls1=numel(s1);
ls2=numel(s2);
l=max(ls1,ls2);
d=ls1-ls2;
if d<0
  s1(ls1+1:ls1+d)=1;
elseif d>0
  s2(ls2+1:ls2+d)=1;
end
if any(s1>1 & s2>1 & s1~=s2)
  error('Array dimensions are not appropriate.');
end
r1=ones(1,l);
r2=r1;
r1(s1==1)=s2(s1==1);
r2(s2==1)=s1(s2==1);
y=feval(f,repmat(x1,r1),repmat(x2,r2));