Esope source files

   1 : C FFPR15    SOURCE    CHAT      05/01/12    23:59:06     5004           
   2 :       SUBROUTINE FFPR15(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(15),COOR(3,3),XA(3,*)
  11 : C                                                                               
  12 :       A=SQRT(2.D0)
  13 :       AS2=A*0.5D0
  14 : C?    CA2=2.D0/A/A
  15 :       CA2=1.D0
  16 : C?    CA4=4.D0/A/A
  17 :       CA4=2.D0
  18 :       A4=4.D0/A
  19 : 
  20 :       COOR(1,1)=AS2
  21 :       COOR(2,1)=0.D0
  22 :       COOR(3,1)=0.5D0
  23 : 
  24 :       COOR(1,2)=AS2
  25 :       COOR(2,2)=AS2
  26 :       COOR(3,2)=0.5D0
  27 : 
  28 :       COOR(1,3)=0.D0
  29 :       COOR(2,3)=AS2
  30 :       COOR(3,3)=0.5D0
  31 : 
  32 :       DO 1 I=1,3
  33 :       X=COOR(1,I)
  34 :       Y=COOR(2,I)
  35 :       Z=COOR(3,I)
  36 : 
  37 :       FN(1)= CA2*(X+Y-A)*(1.D0-Z)*(X+Y+A*Z-AS2)
  38 :       FN(2)=-CA4*X*(X+Y-A)*(1.D0-Z)
  39 :       FN(3)=-CA2*X*(1.D0-Z)*(-X+A*Z+AS2)
  40 :       FN(4)= CA4*X*Y*(1.D0-Z)
  41 :       FN(5)=-CA2*Y*(1.D0-Z)*(-Y+A*Z+AS2)
  42 :       FN(6)=-CA4*Y*(X+Y-A)*(1.D0-Z)
  43 :       FN(7)=-A4*Z*(1.D0-Z)*(X+Y-A)
  44 :       FN(8)= A4*X*Z*(1.D0-Z)
  45 :       FN(9) = A4*Y*Z*(1.D0-Z)
  46 :       FN(10)= CA2*(X+Y-A)*Z*(X+Y+A*(1.D0-Z)-AS2)
  47 :       FN(11)=-CA4*X*(X+Y-A)*Z
  48 :       FN(12)=-CA2*X*Z*(-X+A*(1.D0-Z)+AS2)
  49 :       FN(13)= CA4*X*Y*Z
  50 :       FN(14)=-CA2*Y*Z*(-Y+A*(1.D0-Z)+AS2)
  51 :       FN(15)=-CA4*Y*Z*(X+Y-A)
  52 : 
  53 :       DO 2 N=1,3
  54 :       XA(N,I+15)=
  55 :      &         FN(1)*XA(N,1)+FN(2)*XA(N,2)+FN(3)*XA(N,3)+FN(4)*XA(N,4)
  56 :      &+        FN(5)*XA(N,5)+FN(6)*XA(N,6)+FN(7)*XA(N,7)+FN(8)*XA(N,8)
  57 :      &+  FN(9)*XA(N,9)+FN(10)*XA(N,10)+FN(11)*XA(N,11)+FN(12)*XA(N,12)
  58 :      &+FN(13)*XA(N,13)+FN(14)*XA(N,14)+FN(15)*XA(N,15)
  59 :  2    CONTINUE
  60 :  1    CONTINUE
  61 : 
  62 : 
  63 : 
  64 : 
  65 : 1002  format(10(1x,1pe11.4))
  66 :       RETURN                                                                    
  67 :       END                                                                       
  68 :  
  69 :