Télécharger ksym.eso

Retour à la liste

Numérotation des lignes :

ksym
  1. C KSYM SOURCE CHAT 05/01/13 01:09:27 5004
  2. SUBROUTINE KSYM (IES,IES1,A1,U1,A2,U2,UA,UN)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. C----------------------------------------------------------------------
  6. C CALCUL DES SYMETRIQUES PAR RAPPORT A UN AXE(2D) OU UN PLAN (3D)
  7. C A1 : SOMMETS---> A2
  8. C U1 : NORMALE--->U2
  9. C L'AXE OU LE PLAN EST DEFINIPAR UN POINT ET UN VECTEUR NORMAL
  10. C IES1 = IES+1
  11. C LE NOMBRE DE SOMMETS EN 2D OU 3D EST EGAL A IES
  12. C----------------------------------------------------------------------
  13. C
  14. DIMENSION A1(IES,IES),A2(IES,IES)
  15. DIMENSION U1(IES1),U2(IES1),UA(IES),UN(IES)
  16. C
  17. C VECTEUR NORMAL
  18. C
  19. VN = 0.
  20. DO 1 K=1,IES
  21. VN = VN +U1(K)*UN(K)
  22. 1 CONTINUE
  23. DO 2 K=1,IES
  24. U2(K) = U1(K) - 2*VN*UN(K)
  25. 2 CONTINUE
  26. C
  27. C SOMMETS
  28. C
  29. DO 30 NS = 1,IES
  30. VN = 0.
  31. DO 10 K=1,IES
  32. VN = VN +(A1(K,NS)-UA(K))*UN(K)
  33. 10 CONTINUE
  34. DO 20 K=1,IES
  35. A2(K,NS) = UA(K) + (A1(K,NS)-UA(K)) - 2*VN*UN(K)
  36. 20 CONTINUE
  37.  
  38. 30 CONTINUE
  39. C
  40. U2(IES1)=0.
  41. DO 40 K=1,IES
  42. U2(IES1) = U2(IES1)- U2(K)*A2(K,1)
  43. 40 CONTINUE
  44. C
  45. C
  46. RETURN
  47. END
  48.  
  49.  

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