samer@4: function Z=vecop1(F,X,Y)
samer@4: % vecop1 - apply binary function to different sized arrays
samer@4: %
samer@4: % vecop1 ::
samer@4: %    ([[D]],[[D]]->[[D]])  ~'some function requiring equal size args',
samer@4: %    [[DX]]                ~'first arg of size DX',
samer@4: %    [[DY]]                ~'second arg of size DY'
samer@4: % -> [[DX]]                ~'result of size DX' :- all(DX>=DY).
samer@4: %
samer@4: % this is a slightly simpler version of vecop that assumes that the
samer@4: % first argument is at least as big as the second in every dimension.
samer@4: Z=feval(F,X,repmat_to(Y,size(X)));
samer@4: