OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
getuloc.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine getuloc (sav, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, ulx1, ulx2, ulx3, ulx4, ulx5, ulx6, ulx7, ulx8, uly1, uly2, uly3, uly4, uly5, uly6, uly7, uly8, ulz1, ulz2, ulz3, ulz4, ulz5, ulz6, ulz7, ulz8, r, nel)

Function/Subroutine Documentation

◆ getuloc()

subroutine getuloc ( double precision, dimension(nel,21) sav,
double precision, dimension(mvsiz) xd1,
double precision, dimension(mvsiz) xd2,
double precision, dimension(mvsiz) xd3,
double precision, dimension(mvsiz) xd4,
double precision, dimension(mvsiz) xd5,
double precision, dimension(mvsiz) xd6,
double precision, dimension(mvsiz) xd7,
double precision, dimension(mvsiz) xd8,
double precision, dimension(mvsiz) yd1,
double precision, dimension(mvsiz) yd2,
double precision, dimension(mvsiz) yd3,
double precision, dimension(mvsiz) yd4,
double precision, dimension(mvsiz) yd5,
double precision, dimension(mvsiz) yd6,
double precision, dimension(mvsiz) yd7,
double precision, dimension(mvsiz) yd8,
double precision, dimension(mvsiz) zd1,
double precision, dimension(mvsiz) zd2,
double precision, dimension(mvsiz) zd3,
double precision, dimension(mvsiz) zd4,
double precision, dimension(mvsiz) zd5,
double precision, dimension(mvsiz) zd6,
double precision, dimension(mvsiz) zd7,
double precision, dimension(mvsiz) zd8,
double precision, dimension(mvsiz) ulx1,
double precision, dimension(mvsiz) ulx2,
double precision, dimension(mvsiz) ulx3,
double precision, dimension(mvsiz) ulx4,
double precision, dimension(mvsiz) ulx5,
double precision, dimension(mvsiz) ulx6,
double precision, dimension(mvsiz) ulx7,
double precision, dimension(mvsiz) ulx8,
double precision, dimension(mvsiz) uly1,
double precision, dimension(mvsiz) uly2,
double precision, dimension(mvsiz) uly3,
double precision, dimension(mvsiz) uly4,
double precision, dimension(mvsiz) uly5,
double precision, dimension(mvsiz) uly6,
double precision, dimension(mvsiz) uly7,
double precision, dimension(mvsiz) uly8,
double precision, dimension(mvsiz) ulz1,
double precision, dimension(mvsiz) ulz2,
double precision, dimension(mvsiz) ulz3,
double precision, dimension(mvsiz) ulz4,
double precision, dimension(mvsiz) ulz5,
double precision, dimension(mvsiz) ulz6,
double precision, dimension(mvsiz) ulz7,
double precision, dimension(mvsiz) ulz8,
double precision, dimension(3,3,mvsiz) r,
integer nel )

Definition at line 29 of file getuloc.F.

43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C G l o b a l P a r a m e t e r s
49C-----------------------------------------------
50#include "mvsiz_p.inc"
51C-----------------------------------------------
52C C o m m o n B l o c k s
53C-----------------------------------------------
54C-----------------------------------------------
55C D u m m y A r g u m e n t s
56C-----------------------------------------------
57 INTEGER NEL
58C REAL
59! my_real
60! . OFFG(*)
61
62 DOUBLE PRECISION
63 . XD1(MVSIZ), XD2(MVSIZ), XD3(MVSIZ), XD4(MVSIZ),
64 . XD5(MVSIZ), XD6(MVSIZ), XD7(MVSIZ), XD8(MVSIZ),
65 . YD1(MVSIZ), YD2(MVSIZ), YD3(MVSIZ), YD4(MVSIZ),
66 . YD5(MVSIZ), YD6(MVSIZ), YD7(MVSIZ), YD8(MVSIZ),
67 . ZD1(MVSIZ), ZD2(MVSIZ), ZD3(MVSIZ), ZD4(MVSIZ),
68 . ZD5(MVSIZ), ZD6(MVSIZ), ZD7(MVSIZ), ZD8(MVSIZ),
69 . SAV(NEL,21), R(3,3,MVSIZ),
70 . ULX1(MVSIZ), ULX2(MVSIZ), ULX3(MVSIZ), ULX4(MVSIZ),
71 . ULX5(MVSIZ), ULX6(MVSIZ), ULX7(MVSIZ), ULX8(MVSIZ),
72 . ULY1(MVSIZ), ULY2(MVSIZ), ULY3(MVSIZ), ULY4(MVSIZ),
73 . ULY5(MVSIZ), ULY6(MVSIZ), ULY7(MVSIZ), ULY8(MVSIZ),
74 . ULZ1(MVSIZ), ULZ2(MVSIZ), ULZ3(MVSIZ), ULZ4(MVSIZ),
75 . ULZ5(MVSIZ), ULZ6(MVSIZ), ULZ7(MVSIZ), ULZ8(MVSIZ)
76C-----------------------------------------------
77C L o c a l V a r i a b l e s
78C-----------------------------------------------
79 INTEGER I
80C=======================================================================
81 DOUBLE PRECISION U,V,W,UL,VL,WL
82
83 DO i=1,nel
84 ulx1(i) = zero ! sb : On pourrait eviter de stocker ces 3 tables (toujours zero)
85 uly1(i) = zero
86 ulz1(i) = zero
87
88 u = (xd2(i)-xd1(i))
89 v = (yd2(i)-yd1(i))
90 w = (zd2(i)-zd1(i))
91 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
92 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
93 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
94 ulx2(i) = ul-sav(i,1)
95 uly2(i) = vl-sav(i,2)
96 ulz2(i) = wl-sav(i,3)
97
98 u = (xd3(i)-xd1(i))
99 v = (yd3(i)-yd1(i))
100 w = (zd3(i)-zd1(i))
101 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
102 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
103 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
104 ulx3(i) = ul-sav(i,4)
105 uly3(i) = vl-sav(i,5)
106 ulz3(i) = wl-sav(i,6)
107
108 u = (xd4(i)-xd1(i))
109 v = (yd4(i)-yd1(i))
110 w = (zd4(i)-zd1(i))
111 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
112 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
113 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
114 ulx4(i) = ul-sav(i,7)
115 uly4(i) = vl-sav(i,8)
116 ulz4(i) = wl-sav(i,9)
117
118 u = (xd5(i)-xd1(i))
119 v = (yd5(i)-yd1(i))
120 w = (zd5(i)-zd1(i))
121 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
122 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
123 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
124 ulx5(i) = ul-sav(i,10)
125 uly5(i) = vl-sav(i,11)
126 ulz5(i) = wl-sav(i,12)
127
128 u = (xd6(i)-xd1(i))
129 v = (yd6(i)-yd1(i))
130 w = (zd6(i)-zd1(i))
131 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
132 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
133 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
134 ulx6(i) = ul-sav(i,13)
135 uly6(i) = vl-sav(i,14)
136 ulz6(i) = wl-sav(i,15)
137
138 u = (xd7(i)-xd1(i))
139 v = (yd7(i)-yd1(i))
140 w = (zd7(i)-zd1(i))
141 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
142 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
143 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
144 ulx7(i) = ul-sav(i,16)
145 uly7(i) = vl-sav(i,17)
146 ulz7(i) = wl-sav(i,18)
147
148 u = (xd8(i)-xd1(i))
149 v = (yd8(i)-yd1(i))
150 w = (zd8(i)-zd1(i))
151 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
152 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
153 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
154 ulx8(i) = ul-sav(i,19)
155 uly8(i) = vl-sav(i,20)
156 ulz8(i) = wl-sav(i,21)
157
158 ENDDO
159 RETURN