Esope source files

   1 : C HPRIMEX   SOURCE    AS218369  10/01/29    21:15:01     6623           
   2 :       SUBROUTINE HPRIMEX(XE,NBNN,IDIM,NBSH,SHPTOT,IGAU,SHPWRK,DJAC)                   
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C   ENTREES                                                                     
   6 : C     XE(3,NBNN)= COORDONNEES DES NOEUDS                                        
   7 : C     NBNN      = NOMBRE DES NOEUDS       
   8 : C     NBSH      = NOMBRE de FONCTIONS de FORME                                      
   9 : C     IDIM      = DIMENSION DU PROBLEME                                         
  10 : C     SHPTOT(6,NBNO,IGAU)= FONCTIONS DE FORME REDUITES (VIENT DE MINTE          
  11 : C     IGAU      =NUMERO DU POINT DE GAUUS                                       
  12 : C  SORTIES                                                                      
  13 : C     SHPWRK(6,NBNO)= H  H,X   H,Y   H,Z                                        
  14 : C     DJAC      =JACOBIEN                                                       
  15 : C                                                                               
  16 : C=======================================================================        
  17 :       IMPLICIT INTEGER(I-N)
  18 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  19 :       DIMENSION XE(3,*),SHPTOT(6,NBSH,*),SHPWRK(6,*)                            
  20 :       IDIM1=IDIM+1                                                              
  21 :       DO 100 IA=1,NBNN                                                          
  22 :       DO 100 IB=1,IDIM1                                                         
  23 :       SHPWRK(IB,IA)=SHPTOT(IB,IA,IGAU)                                          
  24 :   100 CONTINUE                                                                  
  24 :       CALL JACOBI(XE,SHPWRK,IDIM,NBNN,DJAC)                                     
  25 :       RETURN                                                                    
  26 :       END                                                                       
  27 :  
  28 :  
  29 :