mcarrayMax.m

function maxOut = mcarrayMax(d,maxtype)
% 
% Max value across all mocap structs in an array.
% If the number of frames in the structs are different, the longer structs
% will be cut.
%
% maxOut = mcarrayMax(d)
% maxOut = mcarrayMax(d,maxtype)
%
% d: array of mocap structs
% maxtype: 'normal' (Default) choose max value across all data streams
%          'absmax' choose most extreme value (disregard sign)
%
% maxout: a single mocap struct
%
% By Kristian Nymoen, RITMO/University of Oslo, 2019
%

    

if nargin < 2
    maxtype = 'normal';
end


    dframes = min([d.nFrames]);
    dmarkers = max([d.nMarkers]);

    dl = length(d);

    q = nan(dframes,dmarkers*3,dl);

    for i = 1:dl

        q(1:dframes,1:1:d(i).nMarkers*3,i) = d(i).data(1:dframes,:);

    end

    maxOut = d(1);
    maxOut.nFrames = dframes;
    maxOut.nMarkers = dmarkers;
    maxOut.filename = 'arraymax';
    
    if strcmpi(maxtype,'absmax')

        
        [~,y] = nanmax(abs(q),[],3);
    
        x= (1:dframes*dmarkers*3)' + (y(:)-1)*(dframes*dmarkers*3);
        
        maxOut.data = reshape(q(x(:)),dframes,dmarkers*3);
        
    else
        
        maxOut.data = nanmax(q,[],3);
        
    end
    
       
    

end