; benchmark for gifa
verbose 1    ; check what you're doing
; benchmark for gifa
;   run on a 512 x 2k real data set.
;   Processing performed is :
;   - Exponen. Broading in F2
;   - Real FT, phasing, in F2
;   - 5 points Spline base-line correction in F2
;   - cosine apodization in F1, zero-filling in F1
;   - Real FT and phasing in F1.
; The displayed spectrum at the end is 512 by 1k real points
; results are (some figures date from version 1.0 !) :
;  * 3'50" on Dec VS 3100  VMS 5.2   * 1'30"    on SUN Sparc 330
;  * 2'30" on Macintosh Centris 650 - running MachTen
;  * 53"   on Silicon-Graphics 4D/20 * 47"      on S.G. 4D/25
;  * 37"   on SUN SparcStation 2     * 34"      on DEC Station 5000-200
;  * 36.3" on PC Am486DX4/100 Linux
;  * 27"   on S.G. 4D/30             * 31" 	on IBM Risc 6000 320
;  * 20"   on Alliant FX 40 3CE      * 17.5" on SUN Sparc-10/30
;  * 17"   on HP 725/50              * 16"   on Indigo R4000
;  * 15.5" on INDY PC                * 15"   on IBM Risc6000 550
;  * 12.6" on IBM Risc6000 560       * 12.3" on SUN Sparc-20
;  * 11.7" on INDY R4600-132MHz   
;  * 10.7" on SUN SparcCenter 2000E SuperSparc 85 Mhz
;  * 10.5" on INDY R4400-174MHz      * 10"   on HP 715/80
;  * 9.5"  on PC Pentium 133 Linux   * 7.5"  on HP 735/99
;  * 7.4"  on HP PA-8000 160MHz      * 6.8"  on INDY R5000-180MHz
;  * 6.5"  on SGI-O2 R5000-PC        * 5.06" on Pentium Pro 200

; first prepare a pseudo-data set for benchmarking GIFA
; set-up a 512 x 2k real data set

dim 2
chsize 128 128 one sin 0 f12 revf f1 revf f2  ; put something in it
chsize 512 2k itype 0                         ; and zero-fill

sh date
timer 0
timer 1

em 0 1 revf f2 rft f2                                    ; main F2 processing
phase 20 20 f2 real f2               ;throwing Imaginary parts is needed here
bcorr 2 2 f2 10 100 300 600 900 0                       ; for baseline corr.
sin 0 f1 chsize (%*2) % rft f1                              ; main F1 proc
phase 20 20 f1                                          ;done
verbose 0
sh date
exit

