Sources Esope | Cast3M

ffqua8
C FFQUA8    SOURCE    CHAT      05/01/12    23:59:13     5004
      SUBROUTINE FFQUA8(XA)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C************************************************************************
C
C
C
C************************************************************************
      DIMENSION FN(8),COOR(3,1),XA(3,9)
      DATA COOR/0.5D0,0.5D0,0.D0 /
C
      DO 1 I=1,1
      X=COOR(1,I)
      Y=COOR(2,I)
      Z=COOR(3,I)
 
      FN(1)=-2.D0*(1.D0-X)*(1.D0-Y)*(1.D0-Z)*(X+Y+Z-0.5D0)
      FN(2)= 4.D0*X*(1.D0-X)*(1.D0-Y)*(1.D0-Z)
      FN(3)=-2.D0*X*(1.D0-Y)*(1.D0-Z)*(Z+Y-X+0.5D0)
      FN(4)= 4.D0*X*Y*(1.D0-Y)*(1.D0-Z)
      FN(5)=-2.D0*X*Y*(1.D0-Z)*(1.5D0-X-Y+Z)
      FN(6)= 4.D0*X*Y*(1.D0-X)*(1.D0-Z)
      FN(7)=-2.D0*Y*(1.D0-X)*(1.D0-Z)*(X-Y+Z+0.5D0)
      FN(8)= 4.D0*Y*(1.D0-Y)*(1.D0-X)*(1.D0-Z)
 
      DO 2 N=1,3
      XA(N,I+8)=
     &         FN(1)*XA(N,1)+FN(2)*XA(N,2)+FN(3)*XA(N,3)+FN(4)*XA(N,4)
     &+        FN(5)*XA(N,5)+FN(6)*XA(N,6)+FN(7)*XA(N,7)+FN(8)*XA(N,8)
 2    CONTINUE
 1    CONTINUE
 
 
 
 
1002  format(10(1x,1pe11.4))
      RETURN
      END