1 : C VALPRP SOURCE CHAT 05/01/13 04:01:13 5004 2 : SUBROUTINE VALPRP(S11,S22,S33,S12,S23,S13,S1,S2,S3) 3 : C======================================================================= 4 : C 5 : C CALCUL DES VALEURS PROPRES 6 : C 7 : C======================================================================= 8 : C 9 : C CREATION : F.CORMERY 10 : C E.N.S.M.A - LMPM 11 : C DEC 1992 12 : C 13 : C======================================================================= 14 : IMPLICIT INTEGER(I-N) 15 : IMPLICIT REAL*8 (A-H,O-Z) 16 : C---------------- 17 : SN11=S11 18 : SN22=S22 19 : SN33=S33 20 : SN12=S12 21 : SN23=S23 22 : SN13=S13 23 : C------- CALCUL DES 3 INVARIANTS A,B,C 24 : A=SN11+SN22+SN33 25 : B=SN22*SN33-SN23*SN23+SN33*SN11-SN13*SN13+SN11*SN22-SN12*SN12 26 : C=SN11*(SN22*SN33-SN23*SN23)-SN12*(SN12*SN33-SN13*SN23)+ 27 : 1 SN13*(SN12*SN23-SN13*SN22) 28 : A=-A 29 : C=-C 30 : C------- CALCUL DES CONTRAINTES PRINCIPALES 30 : CALL EQUAT3(A,B,C,SN1,SN2,SN3) 31 : C------- 32 : S1=SN1 33 : S2=SN2 34 : S3=SN3 35 : C------ FIN 36 : 10 RETURN 37 : END 38 : 39 : 40 :
© Cast3M 2003 - All rights reserved.
Disclaimer