from rsf.proj import *
import math

labelfat=4
titlefat=4

labelsize=6
titlesize=8

nt=1500
dt=0.002
Flow('time',None,'spike n1=%d d1=%g o1=-1.5 | math output=x1'%(nt,dt))


##### Ellipse _________________________________________________________________

label=0
k=0

nearlist=[]
linelist=[]
sinclist=[]

symbolist=['plotfat=1','symbolsz=1 symbol=.','symbol=*']
fatlist=['plotfat=8','plotfat=5','plotfat=1',]

for eps in (0.1,0.5,1.0):


    ei=(nt/2.0)-(eps/dt)
    near='near-%g'%(label)
    line='line-%g'%(label)
    sinc='sinc-%g'%(label)
    label=label+1

    #Nearest-Neighbor
    Flow(near,'time',
         '''
         math output="1.0/%g" |
         window min1=%g max1=%g |
         pad beg1=%d end1=%d
         '''%(eps,-1.0*eps,eps,(ei),(ei)))
    Plot(near,
         '''
         graph plotcol=%d plotfat=8 wanttitle=y wantaxis=y label1="x" label2="Amplitude"
         transp=n yreverse=n  labelfat=%d labelsz=%d titlefat=%d titlesz=%d
         min1=-1.5 max1=1.5 min2=-1.5 max2=15 screenratio=1 %s
         title="Rectangle Function Approximating \F9 d\F2 (t)" font=2
         '''%((label+1),labelfat,labelsize,titlefat,titlesize,fatlist[k]))
    nearlist.append(near)

    #Linear
    Flow(line,'time',
         '''
         math output="1.0/%g-abs(input)/(%g)" |
         window min1=%g max1=%g |
         pad beg1=%d end1=%d
         '''%(eps,eps*eps,-1.0*eps,eps,(ei),(ei)))
    Plot(line,
         '''
         graph plotcol=%d plotfat=3 wanttitle=y wantaxis=y label1="x" label2="Amplitude"
         transp=n yreverse=n  labelfat=%d labelsz=%d titlefat=%d titlesz=%d
         min1=-1.5 max1=1.5 min2=-1.5 max2=15 screenratio=1 %s
         title="Triangle Function Approximating \F9 d\F2 (t)" font=2
         '''%((label+1),labelfat,labelsize,titlefat,titlesize,fatlist[k]))
    linelist.append(line)

    #Sinc
    Flow(sinc,'time',
         '''
         put o1=-6 d1=0.008 |
         math output="(%g*sin(%g*x1))/(%g*x1)" 
         '''%(5*eps,5*eps*math.pi,5*eps*math.pi))
    Plot(sinc,
         '''
         graph plotcol=%d plotfat=3 wanttitle=y wantaxis=y label1="x" label2="Amplitude"
         transp=n yreverse=n  labelfat=%d labelsz=%d titlefat=%d titlesz=%d
         min1=-6.0 max1=6.0 min2=-1.5 max2=6 screenratio=1 %s
         title="Sinc Function Approximating \F9 d\F2 (t)" font=2
         '''%((label+1),labelfat,labelsize,titlefat,titlesize,symbolist[k]))
    sinclist.append(sinc)
    k=k+1

Result('nearest',nearlist,'Overlay')
Result('linear',linelist,'Overlay')
Result('sinc',sinclist,'Overlay')



End()