#pragma TextEncoding = "UTF-8"
#pragma rtGlobals=3		// Use modern global access method and strict wave access.


Function average(wave0)
	Wave wave0
	Variable numrows = dimsize(wave0, 0)
	Wavestats/PCST wave0
	Wave M_WaveStats
	Make/D/O/N = (numrows) results_avg = M_WaveStats[%avg][p]
End

Function stdev(wave0)
	Wave wave0
	Variable numrows = dimsize(wave0, 0)
	Wavestats/PCST wave0
	Wave M_WaveStats
	Make/D/O/N = (numrows) results_stdev = M_WaveStats[%sdev][p]
End

Function count(wave0)
	Wave wave0
	Variable numrows = dimsize(wave0, 0)
	Wavestats/PCST wave0
	Wave M_WaveStats
	Make/D/O/N = (numrows) results_count = M_WaveStats[%numpoints][p]
End

Function sumalong(wave0)
	Wave wave0
	Variable numrows = dimsize(wave0, 0)
	Wavestats/PCST wave0
	Wave M_WaveStats
	Make/D/O/N = (numrows) results_sum = M_WaveStats[%avg][p]
End