cosmo average samples hdr

function ds_avg = cosmo_average_samples(ds, varargin)