1       subroutine e3StsLhs( xl,  lStsVec )  argument
12       real*8  lDir(npro,nshl,3), lStsVec(npro,nshl,nResDims),  local
18          lStsVec(:,i,1) = lDir(:,i,1) * lDir(:,i,1)
19          lStsVec(:,i,2) = lDir(:,i,2) * lDir(:,i,2)
20          lStsVec(:,i,3) = lDir(:,i,3) * lDir(:,i,3)
22          lStsVec(:,i,4) = lDir(:,i,1) * lDir(:,i,2)
23          lStsVec(:,i,5) = lDir(:,i,2) * lDir(:,i,3)
24          lStsVec(:,i,6) = lDir(:,i,3) * lDir(:,i,1)
26          lStsVec(:,i,7) = 0.0
27          lStsVec(:,i,8) = 0.0
28          lStsVec(:,i,9) = 0.0
29          lStsVec(:,i,10) = 0.0
30          lStsVec(:,i,11) = 0.0
37       subroutine e3StsRes( xl, rl, lStsVec )  argument
47       real*8  lDir(npro,nshl,3), lStsVec(npro,nshl,nResDims)  local
52          lStsVec(:,i,1) = lDir(:,i,1) * rl(:,i,4)
53          lStsVec(:,i,2) = lDir(:,i,2) * rl(:,i,4)
54          lStsVec(:,i,3) = lDir(:,i,3) * rl(:,i,4)
56          lStsVec(:,i,4) = lDir(:,i,1) * rl(:,i,1)
57          lStsVec(:,i,5) = lDir(:,i,2) * rl(:,i,2)
58          lStsVec(:,i,6) = lDir(:,i,3) * rl(:,i,3)
60          lStsVec(:,i,7) = lDir(:,i,1) * rl(:,i,2)
62          lStsVec(:,i,8) = lDir(:,i,2) * rl(:,i,3)
64          lStsVec(:,i,9) = lDir(:,i,3) * rl(:,i,1)
66          lStsVec(:,i,10) = 0
67          lStsVec(:,i,11) = 0