Add a getBoundaryOperator method in multiblock.DiffOp that uses the the diffop get_boundary_operator method rather than its properties. Werpers might not like this.

function benchmark(func, N)
    default_arg('N',100);

    tic
    profile on

    for i = 1:N
        func();
    end

    profile viewer
end