from rsfproj import *

matlab = WhereIs('matlab')
if not matlab:
    raise Exception('Cannot find Matlab')

matfun1 = 'GenPhaseDiag1'

for mode in ['REG','IRREG','JIT1','JIT2']:
    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="Ratio k/n" label2="Nb. of erroneous entries"
           min2=0 max2=400
           ''')

End()