Télécharger eqhom3.eso

Retour à la liste

Numérotation des lignes :

  1. C EQHOM3 SOURCE CHAT 05/01/12 23:43:56 5004
  2. C SOURCE MARTINEZ
  3. SUBROUTINE EQHOM3 (X,A)
  4. C
  5. C resolution d'une equation sans second membre dans un espace
  6. C a 3 dimensions
  7. C
  8. C
  9. C
  10. C
  11. IMPLICIT INTEGER(I-N)
  12. IMPLICIT REAL*8 (A-H,O-Z)
  13. REAL*8 A(3,3),X(3)
  14. C
  15. A11=A(1,1)
  16. A12=A(1,2)
  17. A13=A(1,3)
  18. A21=A(2,1)
  19. A22=A(2,2)
  20. A23=A(2,3)
  21. A31=A(3,1)
  22. A32=A(3,2)
  23. A33=A(3,3)
  24. C
  25. C------- ON RESOUD ApX=Yp, Ap ET Yp CALC. A PARTIR DE A
  26. C
  27. YP1=-1.D0*(A11*A13+A21*A23+A31*A33)
  28. YP2=-1.D0*(A12*A13+A22*A23+A32*A33)
  29. C
  30. AP11=A11*A11+A21*A21+A31*A31
  31. AP12=A12*A11+A22*A21+A32*A31
  32. AP21=A12*A11+A22*A21+A32*A31
  33. AP22=A12*A12+A22*A22+A32*A32
  34. C
  35. IF (ABS(AP11*AP22-AP21*AP12).GT.1.D-12) THEN
  36. X(1)=(YP1*AP22-YP2*AP12)/(AP11*AP22-AP21*AP12)
  37. X(2)=(YP1*AP21-YP2*AP11)/(AP12*AP21-AP11*AP22)
  38. X(3)=1.D0
  39. C
  40. ELSE
  41. YP1=-1.D0*(A11*A12+A21*A22+A31*A32)
  42. YP2=-1.D0*(A12*A13+A22*A23+A32*A33)
  43. C
  44. AP11=A11*A11+A21*A21+A31*A31
  45. AP12=A13*A11+A23*A21+A33*A31
  46. AP21=A13*A11+A23*A21+A33*A31
  47. AP22=A13*A13+A23*A23+A33*A33
  48. C
  49. IF (ABS(AP11*AP22-AP21*AP12).GT.1.D-12) THEN
  50. X(1)=(YP1*AP22-YP2*AP12)/(AP11*AP22-AP21*AP12)
  51. X(2)=1.D0
  52. X(3)=(YP1*AP21-YP2*AP11)/(AP12*AP21-AP11*AP22)
  53. C
  54. ELSE
  55. YP1=-1.D0*(A11*A13+A21*A23+A31*A33)
  56. YP2=-1.D0*(A12*A11+A22*A21+A32*A31)
  57. C
  58. AP11=A13*A13+A23*A23+A33*A33
  59. AP12=A12*A13+A22*A23+A32*A33
  60. AP21=A12*A13+A22*A23+A32*A33
  61. AP22=A12*A12+A22*A22+A32*A32
  62. C
  63. IF (ABS(AP11*AP22-AP21*AP12).GT.1.D-12) THEN
  64. X(1)=1.D0
  65. X(2)=(YP1*AP21-YP2*AP11)/(AP12*AP21-AP11*AP22)
  66. X(3)=(YP1*AP22-YP2*AP12)/(AP11*AP22-AP21*AP12)
  67. ELSE
  68. X(1)=0.D0
  69. X(2)=0.D0
  70. X(3)=0.D0
  71. ENDIF
  72. ENDIF
  73. ENDIF
  74. C
  75. C
  76. RETURN
  77. END
  78.  
  79.  
  80.  

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