Télécharger la plaquette Cast3M


fr en

Télécharger r1f3kb.eso

Retour à la liste

Numérotation des lignes : 

r1f3kb

  1. C R1F3KB    SOURCE    BP208322  18/10/08    21:15:09     9952           
  2. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  3. C
  4. C   FFTPACK 5.1 
  5. C
  6. C   Authors:  Paul N. Swarztrauber and Richard A. Valent
  7. C
  8. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  9.  
  10.       SUBROUTINE R1F3KB (IDO,L1,CC,IN1,CH,IN2,WA1,WA2)
  11.       IMPLICIT INTEGER(I-N)
  12.       IMPLICIT REAL*8(A-H,O-Z)
  13.       REAL*8     CC(IN1,IDO,3,L1)    ,CH(IN2,IDO,L1,3),
  14.      1           WA1(IDO)   ,WA2(IDO)
  15. C
  16.       ARG=8.D0*ATAN(1.0D0)/3.D0
  17.       TAUR=COS(ARG)
  18.       TAUI=SIN(ARG)
  19.       DO 101 K=1,L1
  20.          CH(1,1,K,1) = CC(1,1,1,K)+2.*CC(1,IDO,2,K)
  21.          CH(1,1,K,2) = CC(1,1,1,K)+(2.*TAUR)*CC(1,IDO,2,K)
  22.      1   -(2.*TAUI)*CC(1,1,3,K)
  23.          CH(1,1,K,3) = CC(1,1,1,K)+(2.*TAUR)*CC(1,IDO,2,K)
  24.      1   +2.*TAUI*CC(1,1,3,K)
  25.   101 CONTINUE
  26.       IF (IDO .EQ. 1) RETURN
  27.       IDP2 = IDO+2
  28.       DO 103 K=1,L1
  29.          DO 102 I=3,IDO,2
  30.             IC = IDP2-I
  31.         CH(1,I-1,K,1) = CC(1,I-1,1,K)+(CC(1,I-1,3,K)+CC(1,IC-1,2,K))
  32.         CH(1,I,K,1) = CC(1,I,1,K)+(CC(1,I,3,K)-CC(1,IC,2,K))
  33.         CH(1,I-1,K,2) = WA1(I-2)*
  34.      1 ((CC(1,I-1,1,K)+TAUR*(CC(1,I-1,3,K)+CC(1,IC-1,2,K)))-
  35.      * (TAUI*(CC(1,I,3,K)+CC(1,IC,2,K))))
  36.      2                   -WA1(I-1)*
  37.      3 ((CC(1,I,1,K)+TAUR*(CC(1,I,3,K)-CC(1,IC,2,K)))+
  38.      * (TAUI*(CC(1,I-1,3,K)-CC(1,IC-1,2,K))))
  39.             CH(1,I,K,2) = WA1(I-2)*
  40.      4 ((CC(1,I,1,K)+TAUR*(CC(1,I,3,K)-CC(1,IC,2,K)))+
  41.      8 (TAUI*(CC(1,I-1,3,K)-CC(1,IC-1,2,K))))
  42.      5                  +WA1(I-1)*
  43.      6 ((CC(1,I-1,1,K)+TAUR*(CC(1,I-1,3,K)+CC(1,IC-1,2,K)))-
  44.      8 (TAUI*(CC(1,I,3,K)+CC(1,IC,2,K))))
  45.               CH(1,I-1,K,3) = WA2(I-2)*
  46.      7 ((CC(1,I-1,1,K)+TAUR*(CC(1,I-1,3,K)+CC(1,IC-1,2,K)))+
  47.      8 (TAUI*(CC(1,I,3,K)+CC(1,IC,2,K))))
  48.      8   -WA2(I-1)*
  49.      9 ((CC(1,I,1,K)+TAUR*(CC(1,I,3,K)-CC(1,IC,2,K)))-
  50.      8 (TAUI*(CC(1,I-1,3,K)-CC(1,IC-1,2,K))))
  51.             CH(1,I,K,3) = WA2(I-2)*
  52.      1 ((CC(1,I,1,K)+TAUR*(CC(1,I,3,K)-CC(1,IC,2,K)))-
  53.      8 (TAUI*(CC(1,I-1,3,K)-CC(1,IC-1,2,K))))
  54.      2                 +WA2(I-1)*
  55.      3 ((CC(1,I-1,1,K)+TAUR*(CC(1,I-1,3,K)+CC(1,IC-1,2,K)))+
  56.      8 (TAUI*(CC(1,I,3,K)+CC(1,IC,2,K))))
  57.   102    CONTINUE
  58.   103 CONTINUE
  59.       RETURN
  60.       END
  61.  
  62.  
  63.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales