; testing PKSYM in 2D, unit Hertz and PPM (using SIMUN)
dim 2 specw 5283.0 4230. freq 600.13 600.13 600. offset 500. 500.
chsize 512 1k zero unit h
simun 1 1.0 3461 999 2.1 1.5 0. 0. 
simun 1 1.0 5000 601 2.1 1.5 0. 0. 
sin 0. f12 ft f12 max
minimax 5000 ($max[1]*2)
real f12 peak 0
pksym 1 1.0
set a = (abs($pk2d_f1f[3] - 232.426116943)<0.001)
set b = (abs($pk2d_f2f[3] - 156.606658936)<0.001)
set c = (abs($pk2d_f1f[4] - 251.604873657)<0.001)
set d = (abs($pk2d_f2f[4] + 31.1011352539)<0.001)
report ($a & $b & $c & $d) "PKSYM 2D with option add"
pkselect 1 2 3 0
pksym 0 1.0
set e = (abs($npk2d - 2)<0.001)
set f = (abs($pk2d_f2f[1] - 453)<0.001)
set g = (abs($pk2d_f1f[1] - 114)<0.001)
set h = (abs($pk2d_f2f[2] - 156.606658936)<0.001)
set j = (abs($pk2d_f1f[2] - 232.426116943)<0.001)
report ($e & $f & $g & $h & $j) "PKSYM 2D with option remove"

; testing PKPROJ in 2D, unit Hertz (using SIMUN)
dim 2 specw 5283.0 4230. freq 600.13 600.13 600. offset 500. 500.
chsize 512 1k zero unit h
simun 1 1.0 3461 999 2.1 1.5 0. 0. 
simun 1 1.0 5000 601 2.1 1.5 0. 0. 
sin 0. f12 ft f12 max
minimax 5000 ($max[1]*2)
real f12 peak 0
dim 1 specw 5283.0 freq 600.13 600. offset 500. chsize 256 dim 2
pkproj f2 s 1.0
set c = (abs($pk1d_f[1] - 114)<0.001)
set d = (abs($pk1d_f[2] - 39)<0.001)
report ($c & $d) "PKPROJ 2D along F2"
dim 1 specw 4230.0 freq 600.13 600. offset 500. chsize 512 dim 2
pkproj f1 s 1.0
set a = (abs($pk1d_f[1] - 453.000000000)<0.001)
set b = (abs($pk1d_f[2] - 501.000000000)<0.001)
report ($a & $b) "PKPROJ 2D along F1"

quit
