up [pdf]
from rsf.proj import *

# ===========================================================================
# Comparion of mode separation using low-rank decomp. of different error level

SConscript('../bptti2007/SConstruct')
# SConscript('../bptti2007-10e-3/SConstruct')

eps6p = '../bptti2007/ElasticSepP.rsf'
eps3p = '../bptti2007-10e-3/ElasticSepP.rsf'
eps6s = '../bptti2007/ElasticSepSV.rsf'
eps3s = '../bptti2007-10e-3/ElasticSepSV.rsf'

Flow('ElasticSepP6',eps6p,
     '''
     cp  
         ''')
Flow('ElasticSepP3',eps3p,
     '''
     cp
         ''')
Flow('ElasticSepSV6',eps6s,
     '''
     cp  
         ''')
Flow('ElasticSepSV3',eps3s,
     '''
     cp
         ''')

Flow('ElasticSepP6vs3Dif','ElasticSepP6 ElasticSepP3',
     '''
         add scale=1,-1 ${SOURCES[1]} 
         ''')
Flow('ElasticSepSV6vs3Dif','ElasticSepSV6 ElasticSepSV3',
     '''
         add scale=1,-1 ${SOURCES[1]} 
         ''')

name='''ElasticSepP3 ElasticSepSV3 ElasticSepP6vs3Dif ElasticSepSV6vs3Dif'''

for qq in Split(name):
        Result(qq,
        '''
        grey color=j scalebar=y allpos=n barreverse=y bartype=h dbarnum=5 barlabel='' barunit='' screenratio=0.87 wanttitle=n axisfat=5 axiscol=7 labelsz=10
        ''')

End()

sfsegyread
sfput
sfmath
sfgrey
sfwindow
sfsmooth
sftti2delrsep
sftti2delrdecomp2p
sfcp
sfadd