function n_part = count(ipr)
    % returns the number of particles in a reconstruction
    n_part  = size(ipr.re_part, 2);
end