Télécharger ffqua8.eso

Retour à la liste

Numérotation des lignes :

  1. C FFQUA8 SOURCE CHAT 05/01/12 23:59:13 5004
  2. SUBROUTINE FFQUA8(XA)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. C************************************************************************
  6. C
  7. C
  8. C
  9. C************************************************************************
  10. DIMENSION FN(8),COOR(3,1),XA(3,9)
  11. DATA COOR/0.5D0,0.5D0,0.D0 /
  12. C
  13. DO 1 I=1,1
  14. X=COOR(1,I)
  15. Y=COOR(2,I)
  16. Z=COOR(3,I)
  17.  
  18. FN(1)=-2.D0*(1.D0-X)*(1.D0-Y)*(1.D0-Z)*(X+Y+Z-0.5D0)
  19. FN(2)= 4.D0*X*(1.D0-X)*(1.D0-Y)*(1.D0-Z)
  20. FN(3)=-2.D0*X*(1.D0-Y)*(1.D0-Z)*(Z+Y-X+0.5D0)
  21. FN(4)= 4.D0*X*Y*(1.D0-Y)*(1.D0-Z)
  22. FN(5)=-2.D0*X*Y*(1.D0-Z)*(1.5D0-X-Y+Z)
  23. FN(6)= 4.D0*X*Y*(1.D0-X)*(1.D0-Z)
  24. FN(7)=-2.D0*Y*(1.D0-X)*(1.D0-Z)*(X-Y+Z+0.5D0)
  25. FN(8)= 4.D0*Y*(1.D0-Y)*(1.D0-X)*(1.D0-Z)
  26.  
  27. DO 2 N=1,3
  28. XA(N,I+8)=
  29. & FN(1)*XA(N,1)+FN(2)*XA(N,2)+FN(3)*XA(N,3)+FN(4)*XA(N,4)
  30. &+ FN(5)*XA(N,5)+FN(6)*XA(N,6)+FN(7)*XA(N,7)+FN(8)*XA(N,8)
  31. 2 CONTINUE
  32. 1 CONTINUE
  33.  
  34.  
  35.  
  36.  
  37. 1002 format(10(1x,1pe11.4))
  38. RETURN
  39. END
  40.  
  41.  
  42.  

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