Commit a9ca9acf authored by thomas.forbriger's avatar thomas.forbriger

ts/croposp: adjust parameters of test example

parent 4132ac71
......@@ -62,7 +62,7 @@ clean: ;
# set time series parameters
# --------------------------
CPSD_DT=0.01
CPSD_DUR=8640.
CPSD_DUR=8640.065
#CPSD_DUR=864.
# ----------------------------------------------------------------------
......@@ -129,27 +129,27 @@ cpsd1raw.bin:
/bin/bash -c 'export GSL_RNG_SEED=$$(date +"%N"); \
export GSL_RNG_TYPE=$(RNGTYPE); \
siggenx 14 $@ $(CPSD_SIGGENOPT); \
echo seed: $$GSL_RNG_SEED'
echo seed: $$GSL_RNG_SEED > $(patsubst %.bin,%seed.xxx,$@)'
cpsd2raw.bin:
/bin/bash -c 'export GSL_RNG_SEED=$$(date +"%N"); \
export GSL_RNG_TYPE=$(RNGTYPE); \
siggenx 14 $@ $(CPSD_SIGGENOPT); \
echo seed: $$GSL_RNG_SEED'
echo seed: $$GSL_RNG_SEED > $(patsubst %.bin,%seed.xxx,$@)'
cpsd3raw.bin:
/bin/bash -c 'export GSL_RNG_SEED=$$(date +"%N"); \
export GSL_RNG_TYPE=$(RNGTYPE); \
siggenx 14 $@ $(CPSD_SIGGENOPT); \
echo seed: $$GSL_RNG_SEED'
echo seed: $$GSL_RNG_SEED > $(patsubst %.bin,%seed.xxx,$@)'
cpsd4raw.bin:
/bin/bash -c 'export GSL_RNG_SEED=$$(date +"%N"); \
export GSL_RNG_TYPE=$(RNGTYPE); \
siggenx 14 $@ $(CPSD_SIGGENOPT); \
echo seed: $$GSL_RNG_SEED'
echo seed: $$GSL_RNG_SEED > $(patsubst %.bin,%seed.xxx,$@)'
# create time series for common noise signal
cpsdbgraw.bin:
/bin/bash -c 'export GSL_RNG_SEED=$$(date +"%N"); \
siggenx 14 $@ $(CPSD_SIGGENOPT); \
echo seed: $$GSL_RNG_SEED'
echo seed: $$GSL_RNG_SEED > $(patsubst %.bin,%seed.xxx,$@)'
# provide a time series plot to confirm incoherence
rawsig.ps: cpsd1raw.bin cpsd2raw.bin cpsd3raw.bin cpsd4raw.bin cpsdbgraw.bin
......@@ -180,7 +180,7 @@ cpsd2mod.bin: cpsd2raw.bin
printf "fac $(CPSD_AMP2)\nhpb 4.,4\nlpb 2.,4\nend\n" \
| tidofi -type bin -Type bin -v -o -cs $@ $<
cpsd3mod.bin: cpsd3raw.bin
printf "fac $(CPSD_AMP3)\nlp2 100.,3\nend\n" \
printf "fac $(CPSD_AMP3)\nlpb 100.,2\nend\n" \
| tidofi -type bin -Type bin -v -o -cs $@ $<
cpsd4mod.bin: cpsd4raw.bin
printf "fac $(CPSD_AMP4)\nend\n" \
......@@ -192,6 +192,9 @@ cpsdbg2mod.bin: cpsdbg1mod.bin
printf "dif 0.\nend\n" \
| tidofi -type bin -Type bin -v -o -cs $@ $<
%.ascii: %.bin
any2any --verbose --overwrite --itype bin --otype ascii $@ $<
DISCARDS=839.
# discard transient filter response
%X.bin: %.bin
......@@ -435,13 +438,21 @@ SNPSD_%.pdf: psdraw.xxx psd_out.xxx \
REPORTVAR=echo $(1)=$($(1));
REPORTVARS=$(foreach var,$(1),$(call REPORTVAR,$(var)))
CPSDSEEDFILES=$(addprefix cpsd,$(addsuffix rawseed.xxx,$(CPSD_SET) bg))
CPSDASCIIFILES=$(addprefix cpsd,$(addsuffix modhd.xxx,$(CPSD_SET) bg1 bg2))
%hd.xxx: %.ascii
egrep -n '^#' $< > $@
.PHONY: report
report: $(CPSD_BINFILES)
report: $(CPSD_BINFILES) $(CPSDSEEDFILES) $(CPSDASCIIFILES)
sigval -type bin \
-format "File %F: %NS samples, duration: %L" $^
-format "File %F: %NS samples, duration: %L" $(CPSD_BINFILES)
@$(call REPORTVARS,CPSD_NSEGMENTS CPSD_DIVISOR CPSD_OVERLAP)
@$(call REPORTVARS,CPSD_PADFACTOR CPSD_AVGOUT CPSD_NLOG)
@$(call REPORTVARS,CPSD_LOG)
@$(call REPORTVARS,CPSD_DT CPSD_DUR RNGTYPE DISCARDS)
@tail -n 999 $(CPSDSEEDFILES)
@tail -n 999 $(CPSDASCIIFILES)
# ----- END OF Makefile -----
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment