Esope source files

   1 : C CALUHG    SOURCE    CHAT      05/01/12    21:48:06     5004           
   2 :       SUBROUTINE CALUHG(U,H,NG)                                                 
   3 : C************************************************************************       
   4 : C     CALCUL DES POIDS ET ABCISSES DES POINTS DE GAUSS SUR LE SEGMENT           
   5 : C     (-1,1) DE 1 A 6 RACINES                                                   
   6 : C                                                                               
   7 : C     NG   NOMBRE DE RACINES                                                    
   8 : C     U    ABCISSES                                                             
   9 : C     H    POIDS                                                                
  10 : C************************************************************************       
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  13 :       REAL*8 U(NG),H(NG),USQ                                                    
  14 : C***                                                                            
  15 :                                                                                 
  16 :                                                                                 
  17 :                                                                                 
  18 :                                                                                 
  19 :       GO TO (10,20,30,40,50,60),NG                                              
  19 :       CALL ARRET(0)                                                             
  20 : C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>        
  21 :  10   CONTINUE                                                                  
  22 :       U(1)=0.D0                                                                 
  23 :       H(1)=2.D0                                                                 
  24 :       RETURN                                                                    
  25 : C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>        
  26 :  20   CONTINUE                                                                  
  27 :       U(1)=-SQRT(3.D0)/3.D0                                                    
  28 :       U(2)=-U(1)                                                                
  29 :       H(1)=1.D0                                                                 
  30 :       H(2)=1.D0                                                                 
  31 :       RETURN                                                                    
  32 : C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>        
  33 :  30   CONTINUE                                                                  
  34 :       U(1)=-SQRT(0.6D0)                                                        
  35 :       U(2)=0.D0                                                                 
  36 :       U(3)=-U(1)                                                                
  37 :       H(1)=5.D0/9.D0                                                            
  38 :       H(3)=H(1)                                                                 
  39 :       H(2)=8.D0/9.D0                                                            
  40 :       RETURN                                                                    
  41 : C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>        
  42 :  40   CONTINUE                                                                  
  43 :       USQ=SQRT(4.8D0)                                                          
  44 :       U(1)=(3.D0+USQ)/7.D0                                                      
  45 :       U(1)=-SQRT(U(1))                                                         
  46 :       U(4)=-U(1)                                                                
  47 :       U(2)=(3.D0-USQ)/7.D0                                                      
  48 :       U(2)=-SQRT(U(2))                                                         
  49 :       U(3)=-U(2)                                                                
  50 :       H(1)=(7.2D0-USQ)/14.4D0                                                   
  51 :       H(2)=(7.2D0+USQ)/14.4D0                                                   
  52 :       H(3)=H(2)                                                                 
  53 :       H(4)=H(1)                                                                 
  54 :       RETURN                                                                    
  55 : C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>        
  56 :  50   CONTINUE                                                                  
  57 :       USQ=SQRT(1120.D0)                                                        
  58 :       U(1)=(70.D0+USQ)/126.D0                                                   
  59 :       U(1)=-SQRT(U(1))                                                         
  60 :       U(5)=-U(1)                                                                
  61 :       U(3)=0.D0                                                                 
  62 :       U(2)=(70.D0-USQ)/126.D0                                                   
  63 :       U(2)=-SQRT(U(2))                                                         
  64 :       U(4)=-U(2)                                                                
  65 :       H(1)=21.D0/5.D0*(28.D0+5.D0*USQ)/(1120.D0+70.D0*USQ)                      
  66 :       H(2)=21.D0/5.D0*(28.D0-5.D0*USQ)/(1120.D0-70.D0*USQ)                      
  67 :       H(3)=128.D0/225.D0                                                        
  68 :       H(4)=H(2)                                                                 
  69 :       H(5)=H(1)                                                                 
  70 :       RETURN                                                                    
  71 : C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>         
  72 :  60   CONTINUE                                                                  
  73 :       U(1)=-0.932469514203152D0
  74 :       U(2)=-0.661209386466265D0
  75 :       U(3)=-0.238619186083197D0
  76 :       U(4)=-U(3)                                                                
  77 :       U(5)=-U(2)                                                                
  78 :       U(6)=-U(1)                                                                
  79 :       H(1)=0.171324492379170D0
  80 :       H(2)=0.360761573048139D0
  81 :       H(3)=0.467913934572691D0
  82 :       H(4)=H(3)                                                                 
  83 :       H(5)=H(2)                                                                 
  84 :       H(6)=H(1)                                                                 
  85 :       RETURN                                                                    
  86 :       END                                                                       
  87 :  
  88 :