from rsf.proj import *

Flow('spike',None,'spike n1=40 n2=40 n3=40 k1=20 k2=20 k3=20')

for np in (-1,1):
    p = np*0.7

    for case in range(2):
        filt = 'filt%d-%d' % (np,case)
        plane = 'plane%d-%d' % (np,case)

        Flow([filt,'lag'+filt],None,
             '''
             pfactor2 p=%g q=0.5 nt=20 nx=20
             niter=10 lag=${TARGETS[1]} fixed=%d
             ''' % (p,case))
        Flow(plane,['spike',filt],
             '''
             helicon filt=${SOURCES[1]} div=1 adj=0 |
             helicon filt=${SOURCES[1]} div=1 adj=1
             ''')
        Plot(plane,
             '''
             byte gainpanel=all | grey3 frame1=20 frame2=20 frame3=20
             flat=n point1=0.586 point2=0.586 title="Px=%g Py=0.5"
             ''' % p)

Result('eplane','plane1-0 plane-1-0','SideBySideIso')
Result('tplane','plane1-1 plane-1-1','SideBySideIso')

End()