Sources Esope | Cast3M

xtmxmu
C XTMXMU    SOURCE    CHAT      05/01/13    04:15:09     5004
      FUNCTION XTMXMU(CC,RE,NQ,IANTI)
*
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
-INC CCREEL
      INTEGER IANTI
      REAL*8 CC(*),RE(NQ,*),VB,VC ,XTMXMU
*
*     matrice symetrique
*
      VB=0.D0
      IF (IANTI.EQ.0) THEN
          DO 10 IN=1,NQ
*              IF (ABS(CC(IN)).GT.1.E-10) THEN
              IF (ABS(CC(IN)).GT.XPETIT) THEN
                  VC=0.D0
                  DO 12 JN=1,NQ
                      VC=VC+CC(JN)*RE(JN,IN)
 12               CONTINUE
                  VB=VB+VC*CC(IN)
              ENDIF
 10       CONTINUE
      ELSE IF (IANTI .EQ. 3) THEN
*      matrice diagonale
           DO 15 IN=1,NQ
              VB = VB +  CC(IN)*RE(IN,IN)*CC(IN)
 15        CONTINUE
      ELSE
*     matrice quelconque ou anti symétrique
          DO 20 IN=1,NQ
*              IF (ABS(CC(IN)).GT.1.E-10) THEN
              IF (ABS(CC(IN)).GT.XPETIT) THEN
                  VC=0.D0
                  DO 22 JN = 1,NQ
                      VC = CC(JN)*RE(IN,JN)+ VC
 22               CONTINUE
                  VB = VB + VC*CC(IN)
              ENDIF
 20       CONTINUE
*        end do
      ENDIF
      XTMXMU = VB
*
      END