% Evaluate different function handle for each block in a multiblock.Grid
% Function handles may optionaly take a time argument
% f -- cell array of function handles
%       f{i} = f_i(t,x,y,...)
% t -- optional time point. If not specified, it is assumed that the functions take only spatial arguments.
function gf = evalOn(g, f, t)
    assertType(g, 'multiblock.Grid');
    assertType(f, 'cell');

    default_arg('t', []);

    grids = g.grids;
    nBlocks = length(grids);
    gf = cell(nBlocks, 1);

    if isempty(t)
        for i = 1:nBlocks
            if isa(f{i}, 'function_handle')
                gf{i} = grid.evalOn(grids{i}, f{i});
            else
                gf{i} = f{i};
            end
        end
    else
        for i = 1:nBlocks
            if isa(f{i}, 'function_handle')
                gf{i} = grid.evalOn(grids{i}, @(varargin)f{i}(t,varargin{:}));
            else
                gf{i} = f{i};
            end
        end
    end

    gf = blockmatrix.toMatrix(gf);
end