657 	double sum = 0., sum2 = 0., xtemp;  in FnMeanSigma_short()  local
665 		        xtemp = (double) *value;  in FnMeanSigma_short()
666 		        sum += xtemp;  in FnMeanSigma_short()
667 		        sum2 += (xtemp * xtemp);  in FnMeanSigma_short()
673 		        xtemp = (double) *value;  in FnMeanSigma_short()
674 		        sum += xtemp;  in FnMeanSigma_short()
675 		        sum2 += (xtemp * xtemp);  in FnMeanSigma_short()
681 		xtemp = sum / ngood;  in FnMeanSigma_short()
682 		if (mean)     *mean = xtemp;  in FnMeanSigma_short()
683 		if (sigma)    *sigma = sqrt((sum2 / ngood) - (xtemp * xtemp));  in FnMeanSigma_short()
715 	double sum = 0., sum2 = 0., xtemp;  in FnMeanSigma_int()  local
723 		        xtemp = (double) *value;  in FnMeanSigma_int()
724 		        sum += xtemp;  in FnMeanSigma_int()
725 		        sum2 += (xtemp * xtemp);  in FnMeanSigma_int()
731 		        xtemp = (double) *value;  in FnMeanSigma_int()
732 		        sum += xtemp;  in FnMeanSigma_int()
733 		        sum2 += (xtemp * xtemp);  in FnMeanSigma_int()
739 		xtemp = sum / ngood;  in FnMeanSigma_int()
740 		if (mean)     *mean = xtemp;  in FnMeanSigma_int()
741 		if (sigma)    *sigma = sqrt((sum2 / ngood) - (xtemp * xtemp));  in FnMeanSigma_int()
773 	double sum = 0., sum2 = 0., xtemp;  in FnMeanSigma_float()  local
781 		        xtemp = (double) *value;  in FnMeanSigma_float()
782 		        sum += xtemp;  in FnMeanSigma_float()
783 		        sum2 += (xtemp * xtemp);  in FnMeanSigma_float()
789 		        xtemp = (double) *value;  in FnMeanSigma_float()
790 		        sum += xtemp;  in FnMeanSigma_float()
791 		        sum2 += (xtemp * xtemp);  in FnMeanSigma_float()
797 		xtemp = sum / ngood;  in FnMeanSigma_float()
798 		if (mean)     *mean = xtemp;  in FnMeanSigma_float()
799 		if (sigma)    *sigma = sqrt((sum2 / ngood) - (xtemp * xtemp));  in FnMeanSigma_float()
831 	double sum = 0., sum2 = 0., xtemp;  in FnMeanSigma_double()  local
839 		        xtemp = *value;  in FnMeanSigma_double()
840 		        sum += xtemp;  in FnMeanSigma_double()
841 		        sum2 += (xtemp * xtemp);  in FnMeanSigma_double()
847 		        xtemp = *value;  in FnMeanSigma_double()
848 		        sum += xtemp;  in FnMeanSigma_double()
849 		        sum2 += (xtemp * xtemp);  in FnMeanSigma_double()
855 		xtemp = sum / ngood;  in FnMeanSigma_double()
856 		if (mean)     *mean = xtemp;  in FnMeanSigma_double()
857 		if (sigma)    *sigma = sqrt((sum2 / ngood) - (xtemp * xtemp));  in FnMeanSigma_double()