MINICERN/mathlib/gen/d/synt.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:19  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE SYNT (MP,N,IBCX,POT)
  11       COMMON/FW1/W1(774)/FW2/W2(100)
  12       DIMENSION POT(MP,N)
  13       M=MP-1
  14       MD=MP+M
  15       AX=(1./LOG10(2.))+.01
  16       XM=M
  17       AM=AX*LOG10(XM)
  18       IM=INT(AM)
  19       IBC=IBCX+3
  20       IF(IBCX.EQ.4) GO TO 400
  21 C
  22 C   THE RFT ROUTINE IS A MODIFIED
  23 C   VERSION OF HOCKNEY FAST FOURIER
  24 C   TRANSFORM PACKAGE AND CORRESPONDS TO
  25 C   ISELIN LIBRARY PROGRAM D 700
  26 C
  27       DO 1 J=1,N
  28  1    CALL RFT(IM,POT(1,J),1,POT(1,J),1,IBC)
  29       RETURN
  30 C
  31 C   THE MFT ROUTINE IS USED FOR
  32 C   THE FOURTH TYPE OF BOUNDARY
  33 C   CONDITION ALONG X
  34 C
  35  400  DO 2 J=1,N
  36       DO 3 I=1,MP
  37  3    W1(I)=POT(I,J)
  38 C
  39       CALL MFT(MP,MD,W1,W1(MP+1),W1(4*MP+1),2)
  40 C
  41       DO 4 I=1,MP
  42  4    POT(I,J)=W1(I+MP)
  43  2    CONTINUE
  44       WRITE(6,1008)(POT(I,1),I=1,MP)
  45  1008 FORMAT(7E11.4)
  46       RETURN
  47       END