V3FIT
fftfax.f90
1  SUBROUTINE fftfax_g(n, ifax, trigs)
2  USE stel_kinds
3  IMPLICIT NONE
4 !-----------------------------------------------
5 ! D u m m y A r g u m e n t s
6 !-----------------------------------------------
7  INTEGER n
8  INTEGER, DIMENSION(13) :: ifax
9  REAL(rprec), DIMENSION(*) :: trigs
10 #if !defined(CRAY) || defined(LONESTAR) || defined(MCURIE)
11 !-----------------------------------------------
12 ! L o c a l V a r i a b l e s
13 !-----------------------------------------------
14  INTEGER :: mode=3, i
15 !-----------------------------------------------
16 !
17 ! mode 3 is used for REAL/half-complex transforms. it is possible
18 ! to do complex/complex transforms with other values of mode, but
19 ! documentation of the details were not available when this routine
20 ! was written.
21 !
22  CALL fax (ifax, n, mode)
23  i = ifax(1)
24  IF (ifax(i+1)>5 .or. n<=4) ifax(1) = -99
25  IF (ifax(1) <= 0) stop 'IFAX(1) <= 0 in fftfax'
26  CALL fftrig_g (trigs, n, mode)
27 #else
28  CALL fftfax (n, ifax, trigs)
29 #endif
30  END SUBROUTINE fftfax_g