1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46
|
function resu = bsxfun(func, A, B)
%# function resu = bsxfun(func, A, B)
%# Implements a wrapper around internal bsxfun
%% Copyright (C) 2009-2017 Pascal Dupuis <cdemills@gmail.com>
%%
%% This file is part of the dataframe package for Octave.
%%
%% This package is free software; you can redistribute it and/or
%% modify it under the terms of the GNU General Public
%% License as published by the Free Software Foundation;
%% either version 3, or (at your option) any later version.
%%
%% This package is distributed in the hope that it will be useful,
%% but WITHOUT ANY WARRANTY; without even the implied
%% warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
%% PURPOSE. See the GNU General Public License for more
%% details.
%%
%% You should have received a copy of the GNU General Public
%% License along with this package; see the file COPYING. If not,
%% see <http://www.gnu.org/licenses/>.
try
[A, B, resu] = df_basecomp (A, B, true, @bsxfun);
for indi = (1:max (A.x_cnt(2), B.x_cnt(2)))
indA = min (indi, A.x_cnt(2));
indB = min (indi, B.x_cnt(2));
Au = A.x_data{indA}(:, A.x_rep{indA});
Bu = B.x_data{indB}(:, B.x_rep{indB});
resu.x_data{indi} = bsxfun(func, Au, Bu);
resu.x_rep{indi} = 1:size(resu.x_data{indi}, 2);
end
resu = df_thirddim (resu);
catch
disp (lasterr ());
error ('bsxfun: non-compatible dimensions')
end
end
|