function obj = partition(obj)
    arguments (Input)
        obj (1, 1) {mustBeA(obj, 'miSim')};
    end
    arguments (Output)
        obj (1, 1) {mustBeA(obj, 'miSim')};
    end

    % Assess sensing performance of each agent at each sample point
    % in the domain
    agentPerformances = cellfun(@(x) reshape(x.sensorModel.sensorPerformance(x.pos, x.pan, x.tilt, [obj.objective.X(:), obj.objective.Y(:), zeros(size(obj.objective.X(:)))]), size(obj.objective.X)), obj.agents, 'UniformOutput', false);
    agentPerformances = cat(3, agentPerformances{:});
    
    % Get highest performance value at each point
    [~, idx] = max(agentPerformances, [], 3);

    % Collect agent indices in the same way
    agentInds = cellfun(@(x) x.index * ones(size(obj.objective.X)), obj.agents, 'UniformOutput', false);
    agentInds = cat(3, agentInds{:});

    % Get highest performing agent's index
    [m,n,~] = size(agentInds);
    [i,j] = ndgrid(1:m, 1:n);
    obj.partitioning = agentInds(sub2ind(size(agentInds), i, j, idx));
end