Esope source files

   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 :