header <napatool.hdr>

fs    1.0e+6

directive   WINDOW   NONE

ivar  npts  POWEROF2(16)

dvar  freq coherent(1230.0, npts)
dvar  ph   rand_uniform(0.0, _2pi_)

#* make a noisy distorted signal
node  t sine      0.0 1.000000 freq ph
node  d poly      0.0001 0.01 1.000 0.000  t
node  n noise     0.0 0.000001

node  s sum       d   n

#* analyze

#* tool lsp   "lsp.out" s 1.0  100000.0 npts  // bugged ????
   tool fft   "fft.out" s 1.0  100000.0 npts
   tool tsnr "tsnr.out" s 1.0  100000.0 npts

terminate 1 <= TOOL_INDEX