Télécharger la plaquette Cast3M


fr en

Télécharger eqhom3.eso

Retour à la liste

Numérotation des lignes : 

eqhom3

  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