Esope source files

   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 :