Numérotation des lignes :

C CUBKSI    SOURCE    CHAT      07/10/22    21:15:34     5921      SUBROUTINE CUBKSI(NPINT,XXG,BKSI)      implicit real*8(A-H,O-Z)      implicit integer (I-N) *      SUBROUTINE CUB_KSI(NPINT,XXG,BKSI)**     ------------------------------------------------------------------**         DERIVEE P/R VARIABLES CANONIQUE AUX POINTS D INTEGRATION*          POUR LES ELEMENTS CUBES A 8 NOEUDS*                                                    H. BUNG     03-93*     ------------------------------------------------------------------***   ENTREES*     NPINT   : NBRE DE  PTS D INTEGRATION*     XG      : COOR. CANONIQUES DES  PTS D INTEGRATION**   SORTIES :*     BKSI(3,8,NPINT) : LES DERIVEES**      IMPLICIT NONE**---    VARIABLES GLOBALES**      INTEGER NPINT*      REAL *8 XG(3,5),XXG(5),GS(8,4)*      REAL *8 BKSI(3,8,5)       dimension XG(3,5),XXG(*),GS(8,4),BKSI(3,8,*) **---    VARIABLES LOCALES**      REAL *8 AUX*      REAL *8 UN,UNS8*      INTEGER IP,I       dimension XP(3),XM(3)       UN=1.D0      UNS8=0.125D0*C MODIF ANTOINE:C ON REMPLIT LE TABLEAU XGC      DO IP=1,NPINT      XG(1,IP)=0.d0      XG(2,IP)=0.d0      XG(3,IP)=XXG(IP)      END DO**      DO IP=1,NPINT*-        DO  I=1,3          XP(I) = UN + XG(I,IP)          XM(I) = UN - XG(I,IP)        END DO*-        AUX= -UNS8*XM(3)        BKSI(1,1,IP)= AUX*XM(2)        BKSI(1,2,IP)=-BKSI(1,1,IP)        BKSI(1,3,IP)=-AUX*XP(2)        BKSI(1,4,IP)=-BKSI(1,3,IP)        BKSI(2,1,IP)= AUX*XM(1)        BKSI(2,2,IP)= AUX*XP(1)        BKSI(2,3,IP)=-BKSI(2,2,IP)        BKSI(2,4,IP)=-BKSI(2,1,IP)        AUX= -UNS8*XP(3)        BKSI(1,5,IP)= AUX*XM(2)        BKSI(1,6,IP)=-BKSI(1,5,IP)        BKSI(1,7,IP)=-AUX*XP(2)        BKSI(1,8,IP)=-BKSI(1,7,IP)        BKSI(2,5,IP)= AUX*XM(1)        BKSI(2,6,IP)= AUX*XP(1)        BKSI(2,7,IP)=-BKSI(2,6,IP)        BKSI(2,8,IP)=-BKSI(2,5,IP)        AUX= -UNS8*XM(2)        BKSI(3,1,IP)= AUX*XM(1)        BKSI(3,2,IP)= AUX*XP(1)        BKSI(3,5,IP)=-BKSI(3,1,IP)        BKSI(3,6,IP)=-BKSI(3,2,IP)        AUX= -UNS8*XP(2)        BKSI(3,3,IP)= AUX*XP(1)        BKSI(3,4,IP)= AUX*XM(1)        BKSI(3,7,IP)=-BKSI(3,3,IP)        BKSI(3,8,IP)=-BKSI(3,4,IP)*      END DO*      END

© Cast3M 2003 - Tous droits réservés.
Mentions légales