Esope source files

   1 : C FFQUA8    SOURCE    CHAT      05/01/12    23:59:13     5004           
   2 :       SUBROUTINE FFQUA8(XA)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
   5 : C************************************************************************       
   6 : C                                                                               
   7 : C                                                                               
   8 : C                                                                               
   9 : C************************************************************************       
  10 :       DIMENSION FN(8),COOR(3,1),XA(3,9)
  11 :       DATA COOR/0.5D0,0.5D0,0.D0 /
  12 : C                                                                               
  13 :       DO 1 I=1,1
  14 :       X=COOR(1,I)
  15 :       Y=COOR(2,I)
  16 :       Z=COOR(3,I)
  17 : 
  18 :       FN(1)=-2.D0*(1.D0-X)*(1.D0-Y)*(1.D0-Z)*(X+Y+Z-0.5D0)
  19 :       FN(2)= 4.D0*X*(1.D0-X)*(1.D0-Y)*(1.D0-Z)
  20 :       FN(3)=-2.D0*X*(1.D0-Y)*(1.D0-Z)*(Z+Y-X+0.5D0)
  21 :       FN(4)= 4.D0*X*Y*(1.D0-Y)*(1.D0-Z)
  22 :       FN(5)=-2.D0*X*Y*(1.D0-Z)*(1.5D0-X-Y+Z)
  23 :       FN(6)= 4.D0*X*Y*(1.D0-X)*(1.D0-Z)
  24 :       FN(7)=-2.D0*Y*(1.D0-X)*(1.D0-Z)*(X-Y+Z+0.5D0)
  25 :       FN(8)= 4.D0*Y*(1.D0-Y)*(1.D0-X)*(1.D0-Z)
  26 : 
  27 :       DO 2 N=1,3
  28 :       XA(N,I+8)=
  29 :      &         FN(1)*XA(N,1)+FN(2)*XA(N,2)+FN(3)*XA(N,3)+FN(4)*XA(N,4)
  30 :      &+        FN(5)*XA(N,5)+FN(6)*XA(N,6)+FN(7)*XA(N,7)+FN(8)*XA(N,8)
  31 :  2    CONTINUE
  32 :  1    CONTINUE
  33 : 
  34 : 
  35 : 
  36 : 
  37 : 1002  format(10(1x,1pe11.4))
  38 :       RETURN                                                                    
  39 :       END                                                                       
  40 :  
  41 :