Sources Esope | Cast3M

kra002
C KRA002    SOURCE    CHAT      05/01/13    01:05:44     5004
C*********************************************************************C
C                                                                     C
C --- KRA002 REPRESENTE LA PARTIE NUMERIQUEMENT INTENSE DU CALCUL     C
C            EFFECTUE DANS CAMDF (2D). ON UTILISE UN PROGRAMME        C
C            FORTRAN PLUTOT QU'ESOPE POUR UNE PLUS GRANDE EFFICACITE. C
C                                                                     C
C            LA BOUCLE SUR LE NOMBRE D'ELEMENTS EST VECTORISABLE      C
C            ET PARALLELISABLE DIRECTEMENT.                           C
C                                                                     C
C            POUR OPTIMISER LES ACCES MEMOIRE, ON BOUCLE SUR          C
C            LE NOMBRE D'ELEMENTS A L'INTERIEUR DE LA BOUCLE SUR LES  C
C            SOMMETS.                                                 C
C                                                                     C
C*********************************************************************C
      SUBROUTINE KRA002(AM,GX,D,B,NUM,IPADL,NP,NEL,NPT)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      DIMENSION AM(NEL,NP,2),
     ;          GX(NPT,2),
     ;          D (NPT,2),
     ;          NUM(NP,1),
     ;          IPADL(1),
     ;          B (NEL)
C
C --- L'AFFECTATION Z=0. N'EST LA QUE POUR EVITER QU'ESOPE N'INTERCALE
C     LA CARTE INTEGEROOA,... ENTRE LA DIRECTIVE DE VECTORISATION ET
C     LA BOUCLE DO 11. |||
      Z=0.
      DO 12 I=1,NP
      DO 11 K=1,NEL
      IU=IPADL(NUM(I,K))
      B(K)=B(K)+AM(K,I,1)*GX(IU,1)/D(IU,1)
     &         +AM(K,I,2)*GX(IU,2)/D(IU,2)
 11   CONTINUE
 12   CONTINUE
 
 1002 FORMAT(10(1X,1PE11.4))
 1001 FORMAT(20(1X,I5))
      RETURN
      END