from rsfproj import *

matlab = WhereIs('matlab')
matfun1 = 'GenPhaseDiag1'

for mode in ['REG','IRREG']:
    Flow('phasediag%(mode)s backup%(mode)s.mat'%vars() ,matfun1+'.m',
         matlab + '''
         -nosplash -nojvm -r "%(matfun1)s('${TARGETS[0]}','${TARGETS[1]}','%(mode)s');quit"
         '''%vars(),stdin=0,stdout=-1)
    Result('phasediag%(mode)s'%vars(),
           '''
           sfput o1=.05 d1=.0380 |
           sfgraph title='' label1='k/n' label2='Nb. of false detection'
           min2=0 max2=400
           ''')

End()