1 : C HPRIME SOURCE CHAT 05/01/13 00:26:20 5004 2 : SUBROUTINE HPRIME(XE,NBNN,IDIM,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 IDIM = DIMENSION DU PROBLEME 9 : C SHPTOT(6,NBNO,IGAU)= FONCTIONS DE FORME REDUITES (VIENT DE MINTE 10 : C IGAU =NUMERO DU POINT DE GAUUS 11 : C SORTIES 12 : C SHPWRK(6,NBNO)= H H,X H,Y H,Z 13 : C DJAC =JACOBIEN 14 : C 15 : C======================================================================= 16 : IMPLICIT INTEGER(I-N) 17 : IMPLICIT REAL*8(A-H,O-Z) 18 : DIMENSION XE(3,*),SHPTOT(6,NBNN,*),SHPWRK(6,*) 19 : IDIM1=IDIM+1 20 : DO 100 IA=1,NBNN 21 : DO 100 IB=1,IDIM1 22 : SHPWRK(IB,IA)=SHPTOT(IB,IA,IGAU) 23 : 100 CONTINUE 23 : CALL JACOBI(XE,SHPWRK,IDIM,NBNN,DJAC) 24 : RETURN 25 : END 26 :
© Cast3M 2003 - All rights reserved.
Disclaimer