29
30
31
32#include "implicit_f.inc"
33
34
35
36#include "mvsiz_p.inc"
37
38
39
40 INTEGER :: NEL
41C
43 . sig(mvsiz,6), sign(mvsiz,6), dir(mvsiz,2)
44
45
46
47 INTEGER I
48
50 . x2,y2,xy,xys2
51
52 DO i=1,nel
53 x2=dir(i,1)*dir(i,1)
54 y2=dir(i,2)*dir(i,2)
55 xy=dir(i,1)*dir(i,2)
56 xys2=two*xy*sig(i,4)
57 sign(i,1) = x2*sig(i,1)+y2*sig(i,2)-xys2
58 sign(i,2) = y2*sig(i,1)+x2*sig(i,2)+xys2
59 sign(i,3) = sig(i,3)
60 sign(i,4) =(sig(i,1)-sig(i,2))*xy+(x2-y2)*sig(i,4)
61 sign(i,5) = dir(i,2)*sig(i,6)+dir(i,1)*sig(i,5)
62 sign(i,6) = dir(i,1)*sig(i,6)-dir(i,2)*sig(i,5)
63 ENDDO
64
65 RETURN