fmax
extern double fmax(double x, double y);
double stddev = sqrt(fmax(0, tsumsq / num - avg * avg));
return fmax(x, y);