The least squares mean of positive Hilbert–Schmidt operators