# This script produces a noise image where each pixel represent the
# standard deviation of that pixel from 10 images.
# n  = number of images.
# s  = sum(d)
# ss = sum(d*d)
# m  = s/n
# sd = sqrt(  (ss - (n*m*m)) / (n-1) )
#---------------------------
# a input buffer
# b = sum(d)
# c = sum(d*d)
# f = tempbuf. ends with variance
# d = std dev
# e = mean
# dv displaytype hd 0
#dv displaytype hd 1
#dv displaytype hd 2
#dv displaytype hd 3
#pre read array
#set itime based on actual signal level to get different data points
#--- sample 1
dvdisplay on
isready
itime 1
#
go
isready
#
m.wait 1
dv copy a to b
m.wait 1
dv c = a * a
m.wait 2
#--- sample 2
isready
itime 1
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 3
isready
itime 1
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 4
isready
itime 1
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 5
isready
itime 1
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 6
isready
itime 1
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 7
isready
itime 1
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 8
isready
itime 1
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 9
isready
itime 1
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 10
isready
itime 1
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#
# calculate
#
# f = m
dv f = b / 10
m.wait 2
dv e = f * 1
# f = (m*m*n)
dv f = f * f
m.wait 2
dv f = f * 10
m.wait 2
# f = ss - (m*m*n)
dv f = c - f
m.wait 2
# f = ss - (m*m*n) / (n-1)
dv f = f / 9
m.wait 2
# f = sqrt( ss - (m*m*n) / (n-1) )
dv sqrt f to d
dv bufinfo b 
dv bufinfo f
#--- sample 1
isready
itime 2
#
go
isready
#
m.wait 1
dv copy a to b
m.wait 1
dv c = a * a
m.wait 2
#--- sample 2
isready
itime 2
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 3
isready
itime 2
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 4
isready
itime 2
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 5
isready
itime 2
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 6
isready
itime 2
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 7
isready
itime 2
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 8
isready
itime 2
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 9
isready
itime 2
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 10
isready
itime 2
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#
# calculate
#
# f = m
dv f = b / 10
m.wait 2
dv e = f * 1
# f = (m*m*n)
dv f = f * f
m.wait 2
dv f = f * 10
m.wait 2
# f = ss - (m*m*n)
dv f = c - f
m.wait 2
# f = ss - (m*m*n) / (n-1)
dv f = f / 9
m.wait 2
# f = sqrt( ss - (m*m*n) / (n-1) )
dv sqrt f to d
dv bufinfo b 
dv bufinfo f
#--- sample 1
isready
itime 3
#
go
isready
#
m.wait 1
dv copy a to b
m.wait 1
dv c = a * a
m.wait 2
#--- sample 2
isready
itime 3
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 3
isready
itime 3
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 4
isready
itime 3
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 5
isready
itime 3
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 6
isready
itime 3
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 7
isready
itime 3
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 8
isready
itime 3
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 9
isready
itime 3
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 10
isready
itime 3
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#
# calculate
#
# f = m
dv f = b / 10
m.wait 2
dv e = f * 1
# f = (m*m*n)
dv f = f * f
m.wait 2
dv f = f * 10
m.wait 2
# f = ss - (m*m*n)
dv f = c - f
m.wait 2
# f = ss - (m*m*n) / (n-1)
dv f = f / 9
m.wait 2
# f = sqrt( ss - (m*m*n) / (n-1) )
dv sqrt f to d
dv bufinfo b 
dv bufinfo f
#--- sample 1
isready
itime 4
#
go
isready
#
m.wait 1
dv copy a to b
m.wait 1
dv c = a * a
m.wait 2
#--- sample 2
isready
itime 4
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 3
isready
itime 4
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 4
isready
itime 4
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 5
isready
itime 4
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 6
isready
itime 4
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 7
isready
itime 4
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 8
isready
itime 4
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 9
isready
itime 4
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 10
isready
itime 4
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#
# calculate
#
# f = m
dv f = b / 10
m.wait 2
dv e = f * 1
# f = (m*m*n)
dv f = f * f
m.wait 2
dv f = f * 10
m.wait 2
# f = ss - (m*m*n)
dv f = c - f
m.wait 2
# f = ss - (m*m*n) / (n-1)
dv f = f / 9
m.wait 2
# f = sqrt( ss - (m*m*n) / (n-1) )
dv sqrt f to d
dv bufinfo b 
dv bufinfo f
#--- sample 1
isready
itime 5
#
go
isready
#
m.wait 1
dv copy a to b
m.wait 1
dv c = a * a
m.wait 2
#--- sample 2
isready
itime 5
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 3
isready
itime 5
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 4
isready
itime 5
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 5
isready
itime 5
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 6
isready
itime 5
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 7
isready
itime 5
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 8
isready
itime 5
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 9
isready
itime 5
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 10
isready
itime 5
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#
# calculate
#
# f = m
dv f = b / 10
m.wait 2
dv e = f * 1
# f = (m*m*n)
dv f = f * f
m.wait 2
dv f = f * 10
m.wait 2
# f = ss - (m*m*n)
dv f = c - f
m.wait 2
# f = ss - (m*m*n) / (n-1)
dv f = f / 9
m.wait 2
# f = sqrt( ss - (m*m*n) / (n-1) )
dv sqrt f to d
dv bufinfo b 
dv bufinfo f
#--- sample 1
isready
itime 6
#
go
isready
#
m.wait 1
dv copy a to b
m.wait 1
dv c = a * a
m.wait 2
#--- sample 2
isready
itime 6
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 3
isready
itime 6
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 4
isready
itime 6
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 5
isready
itime 6
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 6
isready
itime 6
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 7
isready
itime 6
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 8
isready
itime 6
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 9
isready
itime 6
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 10
isready
itime 6
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#
# calculate
#
# f = m
dv f = b / 10
m.wait 2
dv e = f * 1
# f = (m*m*n)
dv f = f * f
m.wait 2
dv f = f * 10
m.wait 2
# f = ss - (m*m*n)
dv f = c - f
m.wait 2
# f = ss - (m*m*n) / (n-1)
dv f = f / 9
m.wait 2
# f = sqrt( ss - (m*m*n) / (n-1) )
dv sqrt f to d
dv bufinfo b 
dv bufinfo f
#--- sample 1
isready
itime 7
#
go
isready
#
m.wait 1
dv copy a to b
m.wait 1
dv c = a * a
m.wait 2
#--- sample 2
isready
itime 7
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 3
isready
itime 7
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 4
isready
itime 7
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 5
isready
itime 7
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 6
isready
itime 7
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 7
isready
itime 7
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 8
isready
itime 7
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 9
isready
itime 7
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 10
isready
itime 7
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#
# calculate
#
# f = m
dv f = b / 10
m.wait 2
dv e = f * 1
# f = (m*m*n)
dv f = f * f
m.wait 2
dv f = f * 10
m.wait 2
# f = ss - (m*m*n)
dv f = c - f
m.wait 2
# f = ss - (m*m*n) / (n-1)
dv f = f / 9
m.wait 2
# f = sqrt( ss - (m*m*n) / (n-1) )
dv sqrt f to d
dv bufinfo b 
dv bufinfo f
#--- sample 1
isready
itime 8
#
go
isready
#
m.wait 1
dv copy a to b
m.wait 1
dv c = a * a
m.wait 2
#--- sample 2
isready
itime 8
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 3
isready
itime 8
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 4
isready
itime 8
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 5
isready
itime 8
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 6
isready
itime 8
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 7
isready
itime 8
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 8
isready
itime 8
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 9
isready
itime 8
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 10
isready
itime 8
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#
# calculate
#
# f = m
dv f = b / 10
m.wait 2
dv e = f * 1
# f = (m*m*n)
dv f = f * f
m.wait 2
dv f = f * 10
m.wait 2
# f = ss - (m*m*n)
dv f = c - f
m.wait 2
# f = ss - (m*m*n) / (n-1)
dv f = f / 9
m.wait 2
# f = sqrt( ss - (m*m*n) / (n-1) )
dv sqrt f to d
dv bufinfo b 
dv bufinfo f
#--- sample 1
isready
itime 9
#
go
isready
#
m.wait 1
dv copy a to b
m.wait 1
dv c = a * a
m.wait 2
#--- sample 2
isready
itime 9
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 3
isready
itime 9
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 4
isready
itime 9
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 5
isready
itime 9
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 6
isready
itime 9
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 7
isready
itime 9
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 8
isready
itime 9
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 9
isready
itime 9
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 10
isready
itime 9
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#
# calculate
#
# f = m
dv f = b / 10
m.wait 2
dv e = f * 1
# f = (m*m*n)
dv f = f * f
m.wait 2
dv f = f * 10
m.wait 2
# f = ss - (m*m*n)
dv f = c - f
m.wait 2
# f = ss - (m*m*n) / (n-1)
dv f = f / 9
m.wait 2
# f = sqrt( ss - (m*m*n) / (n-1) )
dv sqrt f to d
dv bufinfo b 
dv bufinfo f
#--- sample 1
isready
itime 10
#
go
isready
#
m.wait 1
dv copy a to b
m.wait 1
dv c = a * a
m.wait 2
#--- sample 2
isready
itime 10
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 3
isready
itime 10
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 4
isready
itime 10
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 5
isready
itime 10
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 6
isready
itime 10
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 7
isready
itime 10
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 8
isready
itime 10
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 9
isready
itime 10
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#--- sample 10
isready
itime 10
#
go
isready
#
m.wait 1
dv b = b + a
m.wait 1
dv a = a * a
m.wait 1
dv c = c + a
m.wait 2
#
# calculate
#
# f = m
dv f = b / 10
m.wait 2
dv e = f * 1
# f = (m*m*n)
dv f = f * f
m.wait 2
dv f = f * 10
m.wait 2
# f = ss - (m*m*n)
dv f = c - f
m.wait 2
# f = ss - (m*m*n) / (n-1)
dv f = f / 9
m.wait 2
# f = sqrt( ss - (m*m*n) / (n-1) )
dv sqrt f to d
dv bufinfo b 
dv bufinfo f