function [avAccels] = computeAvAccels(trajData, IMUaccels)
%COMPUTEAVACCELS Compute average accelerations for a multiposition dataset
%   Detailed explanation goes here

avAccels = zeros(sum(trajData.isPosition), 3);

i_row = 1; i_pos = 1;
for i = 1:length(IMUaccels)
    if ~mod(i, trajData.fs*trajData.avTime) && trajData.isPosition(i_pos)
        avAccels(i_row,:) = mean(IMUaccels((i-((trajData.fs*trajData.avTime) - 2)):i,:), 1); % take average of the samples to reduce effects of noise
        i_row = i_row + 1; i_pos = i_pos + 1;
    elseif ~mod(i, trajData.fs*trajData.avTime) && ~trajData.isPosition(i_pos)
        i_pos = i_pos + 1;
    end
end

end