On other computational platforms such as Matlab or R, the mean or standard deviation of an empty vector is typically computed as NaN. However, in GSL, functions like gsl_stats_mean and gsl_stats_sd return zero instead. May I ask what the rationale is behind this design choice?
Best regards, ZFY