Esope source files

   1 : C BDBSTS    SOURCE    CHAT      05/01/12    21:38:00     5004           
   2 :       SUBROUTINE BDBSTS(BGENE,DJAC,DHOOK,LRE,NSTRS,REL)
   3 : C=======================================================================
   4 : C
   5 : C
   6 : C      INTEGRE LE PRODUIT BDB
   7 : C      "EXTENSION" A PARTIR DE BDBST AU CAS NON-SYMETRIQUE
   8 : C      WARNING INB<-->INC
   9 : C
  10 : C      VERSION SEPTEMBRE 92 PIERRE PEGON
  11 : C      ROUTINE FORTRAN PUR
  12 : C
  13 : C=======================================================================
  14 : C     BGENE(NSTRS,LRE)=MATRICE B
  15 : C     DJAC =JACOBIEN
  16 : C     DHOOK(NSTRS,NSTRS)=MATRICE DE HOOKE
  17 : C     LRE=TAILLE DE LA MATRICE DE RIGIDITE
  18 : C     REL=MATRICE DE RIGIDITE
  19 : C           LES KIJ SONT AJOUTES SUR REL
  20 : C=======================================================================
  21 :       IMPLICIT INTEGER(I-N)
  22 :       IMPLICIT REAL*8 (A-H,O-Z)
  23 :       DIMENSION DHOOK(NSTRS,1),REL(LRE,1)
  24 :       DIMENSION WORK1(10),WORK2(10),BGENE(NSTRS,1)
  25 : C
  26 :       DO   301 IX=1,LRE
  27 : C
  28 : C     BOUCLE SUR LES DDL
  29 : C
  30 :         DO   401 INA=1,NSTRS
  31 : C
  32 : C     CALCUL DE JACOBIAN*POIDS*B   POUR LE DDL IX
  33 : C
  34 :           WORK1(INA)=DJAC*BGENE(INA,IX)
  35 :   401   CONTINUE
  36 :         DO   501 INB=1,NSTRS
  37 :           CC=0.D0
  38 :           DO   601 INC=1,NSTRS
  39 :             CC=CC+WORK1(INC)*DHOOK(INC,INB)
  40 :   601     CONTINUE
  41 :           WORK2(INB)=CC
  42 :   501   CONTINUE
  43 :         DO   701 IY=1,LRE
  44 : C
  45 : C     CALCUL DE BDB POUR LES DDL IX ET IY
  46 : C
  47 :           CC=0.D0
  48 :           DO   801 IND=1,NSTRS
  49 :             CC=CC+BGENE(IND,IY)*WORK2(IND)
  50 :   801     CONTINUE
  51 :           REL(IX,IY)=REL(IX,IY)+CC
  52 :   701   CONTINUE
  53 :   301 CONTINUE
  54 : C
  55 :       RETURN
  56 :       END
  57 :  
  58 :