# this is <Makefile>
# ----------------------------------------------------------------------------
# $Id: Makefile,v 1.37 2006-06-29 07:11:44 tforb Exp $
# 
# Copyright (c) 2005 by Thomas Forbriger (BFO Schiltach) 
# 
# time series waveform filters
#
# Makefile fuer tools /src/ts/wf
#
# 22/01/2001      introduced smoos and did some "tidy up"
# 24/10/2001      introduced siggen
# 24/05/2001      introduced tsfilt
# 
# REVISIONS and CHANGES 
#    11/07/2005   V1.0   Thomas Forbriger
# 
# ============================================================================
#

F2CLIB=-lf2c -lm -L${SERVERLIBDIR}
PGPLOTLIB=-lf2cpgplot52 -lX11 -L/usr/X11/lib
CC=gcc
CFLAGS=-O2
F2CFLAGS=-f -u

FLAGS=
CXXFLAGS=-fhonor-std -Wall $(FLAGS)
CXXFLAGS=-Wall $(FLAGS)
LDFLAGS=-L$(LOCLIBDIR)
CPPFLAGS=-I$(LOCINCLUDEDIR) $(FLAGS)

.f.o:
	f2c $(F2CFLAGS) $<
	gcc $(CFLAGS) $(<:.f=.c) -c -I${SERVERINCLUDEDIR} -I${LOCINCLUDEDIR}
	@rm  $(<:.f=.c)

all:

.PHONY: clean
clean: ; 
	-find . -name \*.bak | xargs --no-run-if-empty /bin/rm -v
	-/bin/rm -vf flist *.o *.xxx.* *.ps *~

flist: $(wildcard *.f *.inc Makefile *.cc)
	echo $^ | tr ' ' '\n' | sort > $@

.PHONY: edit
edit: flist; vim $<

tsfilt: tsfilt.o
	$(CC) $(CFLAGS) $< -o $@ -lseife -ltf \
        $(F2CLIB) -L$(LOCLIBDIR) -static
	newprog $@

stufi: stufi.o
	$(CC) $(CFLAGS) $< -o $@ -lseife -ltf -lsff \
        $(F2CLIB) -L$(LOCLIBDIR) -static
	newprog $@

rotate coro xyz2uvw detect: %: %.o
	$(CC) $(CFLAGS) $< -o $@ -lsff $(F2CLIB) -L$(LOCLIBDIR) -ltf
	newprog $@

susei evelo tesiff teswf: %: %.o
	$(CC) $(CFLAGS) $< -o $@ -ltf -lsff $(F2CLIB) -L$(LOCLIBDIR)
	newprog $@

phasedsignals hamres siggen smoos dise: %: %.o
	$(CC) $(CFLAGS) $< -o $@ -ltf -lsffu -ltime -lsff $(F2CLIB) \
           -L$(LOCLIBDIR)
	newprog $@

noisymize teseco resaseda gatherdiff autocorr cross tidofi fredofi sigfit: \
  %: %.o
	$(CXX) -o $@ $^ -I$(LOCINCLUDEDIR) -lsffxx -ldatreadxx -llinearxx \
                  -ltsxx -ltfxx -lsffxx -lgsexx -ltime++ -laff \
                  -llapack -lblas -lcseife -lsffxx -lg2c -lm\
                  -L$(LOCLIBDIR) $(CXXFLAGS) $(FLAGS)
	newprog $@

phasetest%.sff: phasedsignals
	phasedsignals $@ -v -o -p $(patsubst phasetest%.sff,%,$@) -n 16 \
        -t 10. -b 3. -s 0. -O 2 -N 6 -h 0.5
P%.sff: %.sff; evelo $< $@ -P

%.s.ps: %.sff; splot $< 0. $@/ps 1.e-2

RVAL=0.
%.ps: %.sff; stuplo -d $@/ps -R $(RVAL) $<
%.psp: %.ps; gv $<; /bin/rm -fv $<

#======================================================================
# test section

#----------------------------------------------------------------------
# test teseco
t1_sig.xxx.sff:
	/bin/rm -fv $@
	tesiff $@
t1_sig.xxx.ps: t1_sig.xxx.sff
	stuplo -d $@/cps -t -a -s x -i -Y 0.6 -V -c fTdtF \
        -X "Zeit (UT)" -l 3,2,1 -h 1.,1.,1.4,1.3 $<
t1_teseco1.xxx.sff: t1_sig.xxx.sff
	teseco -v -o $@ $< t:3-5 $< t:3,5 $< t:4
t1_teseco2.xxx.sff: t1_sig.xxx.sff
	teseco -v -o $@ $< t:3-5 $< f:-3. t:5 $< f:0.7 t:4 $< t:3 f:0.7
t1_teseco1.xxx.ps t1_teseco2.xxx.ps: \
  t1_teseco%.xxx.ps: t1_sig.xxx.sff t1_teseco%.xxx.sff
	stuplo -d $@/cps -t -a -s x -i -Y 0.6 -V -c fTdtF \
        -X "Zeit (UT)" -l 3,2,1 -h 1.,1.,1.4,1.3 -C \
        $(word 1,$^) t:3-5 nc: $(word 2,$^)

#----------------------------------------------------------------------
# test tidofi against stufi
sig.xxx.sff: ;siggen 6 $@  -Tm 0.2 -Te 1.

tid.t.xxx.fil tid.s.xxx.fil: ; echo -e "DBG  0\ntid  270\nend" > $@
lpb.t.xxx.fil lpb.s.xxx.fil: ; echo -e "DBG  0\nlpb  0.05,5\nend" > $@
hpb.t.xxx.fil hpb.s.xxx.fil: ; echo -e "DBG  0\nhpb  0.05,5\nend" > $@
lp2.t.xxx.fil lp2.s.xxx.fil: ; echo -e "DBG  0\nlp2  0.2,0.707\nend" > $@
hp2.t.xxx.fil hp2.s.xxx.fil: ; echo -e "DBG  0\nhp2  0.2,0.707\nend" > $@
bp2.t.xxx.fil bp2.s.xxx.fil: ; echo -e "DBG  0\nbp2  0.2,0.707\nend" > $@
he2.t.xxx.fil he2.s.xxx.fil: ; echo -e "DBG  0\nhe2  0.2,0.707,1.,0.1\nend" > $@
le2.t.xxx.fil le2.s.xxx.fil: ; echo -e "DBG 0\nle2  0.2,0.707,.01,0.1\nend" > $@
he1.t.xxx.fil he1.s.xxx.fil: ; echo -e "DBG  0\nhe1  0.2,1.,0.1\nend" > $@
le1.t.xxx.fil le1.s.xxx.fil: ; echo -e "DBG 0\nle1  0.2,.01,0.1\nend" > $@
lp1.t.xxx.fil lp1.s.xxx.fil: ; echo -e "DBG  0\nlp1  0.2\nend" > $@
hp1.t.xxx.fil hp1.s.xxx.fil: ; echo -e "DBG  0\nhp1  0.2\nend" > $@
fac.t.xxx.fil fac.s.xxx.fil: ; echo -e "DBG  0\nfac  -10.\nend" > $@
add.t.xxx.fil add.s.xxx.fil: ; echo -e "DBG  0\nadd  -1.e2\nend" > $@
add1.t.xxx.fil add1.s.xxx.fil: ; echo -e "DBG  0\nadd  -1.e-1\nend" > $@
add2.t.xxx.fil add2.s.xxx.fil: ; echo -e "DBG  0\nadd  +1.e-1\nend" > $@
rev.t.xxx.fil rev.s.xxx.fil: ; echo -e "DBG  0\nrev  200\nend" > $@
# 0.2/(2*pi)=0.031830989
dif.t.xxx.fil dif.s.xxx.fil: ; echo -e "DBG  0\ndif  0.031830989\nend" > $@
int.t.xxx.fil int.s.xxx.fil: ; echo -e "DBG  0\nint  0.031830989\nend" > $@
avg.t.xxx.fil avg.s.xxx.fil: ; echo -e "DBG  0\navg  30\nend" > $@
tre.t.xxx.fil tre.s.xxx.fil: ; echo -e "DBG  0\ntre  100\nend" > $@
han.t.xxx.fil: ; echo -e "DBG  0\nhan  0\nend" > $@
han.s.xxx.fil: ; echo -e "DBG  0\nsis  1,2000000\nend" > $@
first.t.xxx.fil: ; echo -e "DBG  0\nfirst  0\nend" > $@
first.s.xxx.fil: ; echo -e "DBG  0\navg 1\nend" > $@
fbl.t.xxx.fil: ; echo -e "DBG  0\nfbl 0,00,600\nend" > $@
fbl.s.xxx.fil: ; echo -e "DBG  0\nfbl 200,20,600\nend" > $@

%.xxx.rules: Makefile
	( echo "%.$(patsubst %.xxx.rules,%,$@).s.xxx.sff:" \
             " %.xxx.sff $(patsubst %.xxx.rules,%,$@).s.xxx.fil"; \
        echo -e '\tstufi $$(word 2,$$^) -v -o $$< ' \
                '>$$(patsubst %.sff,%.out,$$@) 2>&1'; \
        echo -e '\t/bin/mv -fv $$<.sfi $$@'; \
	  echo "%.$(patsubst %.xxx.rules,%,$@).t.xxx.sff:" \
             " %.xxx.sff $(patsubst %.xxx.rules,%,$@).t.xxx.fil"; \
        echo -e '\ttidofi -cf $$(word 2,$$^)' \
                '-v -o $$@ $$< ' \
                '>$$(patsubst %.sff,%.out,$$@) 2>&1'; \
        echo "%.$(patsubst %.xxx.rules,%,$@).t.xxx.fil:" \
             " $(patsubst %.xxx.rules,%,$@).t.xxx.fil" \
             '; /bin/ln -sfv $$< $$@'; \
        echo "%.$(patsubst %.xxx.rules,%,$@).s.xxx.fil:" \
             " $(patsubst %.xxx.rules,%,$@).s.xxx.fil" \
             '; /bin/ln -sfv $$< $$@'; \
        ) > $@

FILTERS=lp2 hp2 lp1 hp1 fac add rev dif int avg tre hpb lpb han bp2 he2 le2 \
  he1 le1 tid first fbl add1 add2

ifdef filters
include $(addsuffix .xxx.rules,$(FILTERS))
endif

REF=sig.xxx.sff
cmp.%.ps: %.s.xxx.sff %.t.xxx.sff $(REF)
	stuplo -d $@/cps -C -t -s x -i -c fT -l 3,2,1 -h 1.,1.2,1.4,1.2 \
        -Y 0.8 -V -X "Zeit (s)" \
        $(word 1,$^) nc: $(word 2,$^) nc: $(word 3,$^) $^
control.%.xxx.out: %.s.xxx.sff %.t.xxx.sff %.s.xxx.fil %.t.xxx.fil
	tail -999 $(patsubst %.sff,%.out,$^) > $@
control.%.ps: control.%.xxx.out
	cat $< \
	  | a2ps -A virtual -o $@ -1 -r -l 100 \
          --center-title="Control files for $(patsubst control.%,cmp.%,$@)"

report.%.ps: cmp.%.rw.ps control.%.ps
	a2ps -A virtual -o $@ -1 $^ 

#----------------------------------------------------------------------
# test equivalence of lp2, hp2 and bp2
sig1.xxx.sff:
	siggen 5 $@ -v -o -f1 0.01 -f2 10. -n 50 -T 10. -d 0.001
eq2.ps: sig1.bp2.t.xxx.sff sig1.lp2.t.dif.t.xxx.sff sig1.hp2.t.int.t.xxx.sff \
        sig1.lp2.t.dif.t.add1.t.xxx.sff sig1.hp2.t.int.t.add2.t.xxx.sff
	stuplo -d $@/cps -C -t -s x -i -c fT -l 3,2,1 -h 1.,1.2,1.4,1.2 \
        -Y 0.8 -V -X "Zeit (s)" \
        $(word 5,$^) $(wordlist 1,5,$^)  \
        nc: $(word 4,$^) $(wordlist 1,4,$^) \
        nc: $(word 3,$^) $(wordlist 1,3,$^) \
        nc: $(word 2,$^) $(wordlist 1,2,$^) \
        nc: $(word 1,$^) $(wordlist  1,1,$^)

#----------------------------------------------------------------------
# test resampling algorithm
rsa.xxx.sff:
	siggen 5 $@ -v -o -f1 1.e-5 -f2 20.e-3 -n 10 -T 1.e3 -d 10.
rsa.xxx.rsa.sff: rsa.xxx.sff
	resaseda -v -o -tf 2001/2/9/0/1/1.23 -dt 0.05 $@ $<
rsa.ps: rsa.xxx.sff rsa.xxx.rsa.sff
	stuplo -d $@/cps -t -a \
        -c dtfT -s x -i -l 3,2,1 -h 1.3,1.3,1.3,1.3 -N -C \
        $(word 1,$^) nc: $(word 2,$^)

#----------------------------------------------------------------------

%.psp: %.ps; gv $<; /bin/rm -fv $<
%.pdf: %.ps; ps2pdf $<
%.rw.ps: %.ps; gs -sDEVICE=pswrite -dNOPAUSE -dBATCH -sOutputFile=$@ $< 
%.fix.ps: %.ps; fixps $< > $@
%.2up.ps: %.ps; psnup -r -l -2 -Pa4 -pa4 $< > $@
%.rs.ps: %.ps; psresize -Pa4 -pa5 $< > $@
%.pngp: %.png; xv $<; /bin/rm -fv $<

# ----- END OF Makefile -----