Esope source files

   1 : C CALG1     SOURCE    CHAT      05/01/12    21:45:59     5004           
   2 :       SUBROUTINE CALG1(A,B,NG,H,U,X,PG)                                         
   3 : C************************************************************************       
   4 : C   CALCUL DES COORDONNEES ET DES POIDS DES POINTS DE GAUSS SUR LE PAVE         
   5 : C   (A,B)  LES CONNAISSANT SUR LE SEGMENT (-1,1)                                
   6 : C                                                                               
   7 : C   NG     NOMBRE DE POINTS DE GAUSS SUR LE SEGMENT (-1,1)                      
   8 : C   U(NG)  ABCISSES DES POINTS DE GAUSS                                         
   9 : C   H(NG)  POIDS DES POINTS DE GAUSS                                            
  10 : C                                                                               
  11 : C   X      COORDONNEES DES POINTS DE GAUSS SUR LE PAVE (A,B)                    
  12 : C   PG     POIDS DES POINTS DE GAUSS SUR LE PAVE (A,B)                          
  13 : C                                                                               
  14 : C************************************************************************       
  15 :       IMPLICIT INTEGER(I-N)
  16 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  17 :       REAL*8 A,B,H(NG),U(NG),X(NG)                                              
  18 :       REAL*8 PG(NG)                                                             
  19 :                                                                                 
  20 :                                                                                 
  21 :       BMA=(B-A)/2.D0                                                            
  22 :       BPA=(B+A)/2.D0                                                            
  23 : C                                                                               
  24 :       DO 1 LX=1,NG                                                              
  25 :       X(LX)=BMA*U(LX)+BPA                                                       
  26 :       PG(LX)=BMA*H(LX)                                                          
  27 :  1    CONTINUE                                                                  
  28 : 
  29 : C     WRITE(6,*)' SUB CALG1 : POIDS ET POSITIONS DES PTS DE GAUSS '
  30 : C     WRITE(6,*)' SUB CALG1 : H,U,X,Y ET PG '
  31 : C     WRITE(6,1002)H
  32 : C     WRITE(6,1002)U
  33 : C     WRITE(6,1002)X
  34 : C     WRITE(6,1002)Y
  35 : C     WRITE(6,1002)PG
  36 :  1002 FORMAT(10(1X,1PD11.4))                                                    
  37 :       RETURN                                                                    
  38 :       END                                                                       
  39 :  
  40 :