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