# Author: G. Hennenfent
#         Seismic Laboratory for Imaging and Modeling
#         Department of Earch & Ocean Sciences
#         The University of British Columbia
#         
# Date  : February, 07

from rsfproj import *

n1 = 2**13
k11 = n1/3
k12 = 2*n1/3
k13 = 3*n1/4

d1 = 1./(n1-1)

graphx = '''
sfput o1=-.5 d1=%(d1)f |
sfgraph label1="Frequency" label2="Amplitude" unit1="" unit2=""
title=""
'''%vars()

# generate sparse DCT vector
Flow('x',None,
     'sfspike n1=%(n1)i n2=1 k1=%(k11)i,%(k12)i,%(k13)i mag=50,10,30 nsp=3'%vars() )

Result('model',
       'sfgraph title="" label1="Sample" label2="Amplitude" unit1="" unit2="" ')
Flow('spec specREG specIRREG model','x script.py',
     './${SOURCES[1]} specREG=${TARGETS[1]} specIRREG=${TARGETS[2]} model=${TARGETS[3]}')
Result('spec',graphx)
Result('specREG',graphx)
Result('specIRREG',graphx)


End()