Télécharger r1f3kb.eso

Retour à la liste

Numérotation des lignes :

  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