#include "implicit_f.inc"
#include "mvsiz_p.inc"
Go to the source code of this file.
|
| subroutine | qrota2 (sig, s1, s2, s3, s4, s5, s6, wyz, nel, jcvt) |
◆ qrota2()
| subroutine qrota2 |
( |
| sig, |
|
|
| s1, |
|
|
| s2, |
|
|
| s3, |
|
|
| s4, |
|
|
| s5, |
|
|
| s6, |
|
|
| wyz, |
|
|
integer | nel, |
|
|
integer, intent(in) | jcvt ) |
Definition at line 30 of file qrota2.F.
34
35
36
37#include "implicit_f.inc"
38
39
40
41#include "mvsiz_p.inc"
42
43
44
45 INTEGER, INTENT(IN) :: JCVT
46 INTEGER NEL
47 my_real sig(nel,6),s1(*), s2(*), s3(*), s4(*), s5(*), s6(*), wyz(*)
48
49
50
51 INTEGER I
53
54 DO i=1,nel
55 s1(i)=sig(i,1)
56 s2(i)=sig(i,2)
57 s3(i)=sig(i,3)
58 s4(i)=sig(i,4)
59 s5(i)=sig(i,5)
60 s6(i)=sig(i,6)
61 ENDDO
62
63 IF(jcvt/=0) RETURN
64
65 DO i=1,nel
66 q1(i)=two*s4(i)*wyz(i)
67 ENDDO
68
69 DO i=1,nel
70 sig(i,1)=s1(i)+q1(i)
71 sig(i,2)=s2(i)-q1(i)
72 sig(i,4)=s4(i)+wyz(i)*(s2(i)-s1(i))
73 ENDDO
74
75 RETURN