1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23 from numpy import array, int16
24 from unittest import TestCase
25
26
27 from lib.dispersion.ns_mmq_3site import r2eff_ns_mmq_3site_mq, r2eff_ns_mmq_3site_sq_dq_zq
28
29
31 """Unit tests for the lib.dispersion.ns_mmq_3site relax module."""
32
34 """Set up for all unit tests.
35
36 The variable values have been taken from the Relax_disp.test_ns_mmq_3site system test.
37 """
38
39 self.M0 = array([ 0.85, 0.05, 0.1 ])
40 self.F_vector = array([ 1., 0., 0.])
41 self.R20A = array([[[[ 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
42 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
43 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
44 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
45 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
46 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
47 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
48 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
49 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5]],
50
51 [[ 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
52 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
53 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
54 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
55 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
56 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
57 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
58 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
59 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5]],
60
61 [[ 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
62 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
63 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
64 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
65 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
66 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
67 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
68 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
69 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5]],
70
71 [[ 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
72 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
73 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
74 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
75 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
76 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
77 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
78 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
79 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5]]]])
80 self.R20B = array([[[[ 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
81 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
82 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
83 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
84 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
85 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
86 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
87 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
88 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5]],
89
90 [[ 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
91 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
92 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
93 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
94 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
95 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
96 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
97 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
98 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5]],
99
100 [[ 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
101 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
102 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
103 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
104 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
105 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
106 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
107 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
108 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5]],
109
110 [[ 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
111 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
112 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
113 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
114 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
115 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
116 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
117 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
118 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5]]]])
119 self.R20C = array([[[[ 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
120 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
121 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
122 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
123 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
124 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
125 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
126 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5,
127 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5, 7.5]],
128
129 [[ 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
130 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
131 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
132 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
133 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
134 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
135 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
136 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5,
137 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5, 8.5]],
138
139 [[ 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
140 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
141 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
142 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
143 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
144 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
145 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
146 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5,
147 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5, 11.5]],
148
149 [[ 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
150 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
151 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
152 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
153 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
154 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
155 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
156 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5,
157 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5, 13.5]]]])
158 self.pA = 0.85
159 self.pB = 0.05
160 self.dw_AB = array([[[[ 1256.637061435917076, 1256.637061435917076,
161 1256.637061435917076, 1256.637061435917076,
162 1256.637061435917076, 1256.637061435917076,
163 1256.637061435917076, 1256.637061435917076,
164 1256.637061435917076, 1256.637061435917076,
165 1256.637061435917076, 1256.637061435917076,
166 1256.637061435917076, 1256.637061435917076,
167 1256.637061435917076, 1256.637061435917076,
168 1256.637061435917076, 1256.637061435917076,
169 1256.637061435917076, 1256.637061435917076,
170 1256.637061435917076, 1256.637061435917076,
171 1256.637061435917076, 1256.637061435917076,
172 1256.637061435917076, 1256.637061435917076,
173 1256.637061435917076, 1256.637061435917076,
174 1256.637061435917076, 1256.637061435917076,
175 1256.637061435917076, 1256.637061435917076,
176 1256.637061435917076, 1256.637061435917076,
177 1256.637061435917076, 1256.637061435917076,
178 1256.637061435917076, 1256.637061435917076,
179 1256.637061435917076, 1256.637061435917076,
180 1256.637061435917076, 1256.637061435917076,
181 1256.637061435917076, 1256.637061435917076,
182 1256.637061435917076, 1256.637061435917076,
183 1256.637061435917076, 1256.637061435917076,
184 1256.637061435917076, 1256.637061435917076,
185 1256.637061435917076, 1256.637061435917076,
186 1256.637061435917076, 1256.637061435917076,
187 1256.637061435917076, 1256.637061435917076,
188 1256.637061435917076, 1256.637061435917076,
189 1256.637061435917076, 1256.637061435917076,
190 1256.637061435917076, 1256.637061435917076,
191 1256.637061435917076, 1256.637061435917076,
192 1256.637061435917076, 1256.637061435917076,
193 1256.637061435917076, 1256.637061435917076,
194 1256.637061435917076, 1256.637061435917076,
195 1256.637061435917076, 1256.637061435917076,
196 1256.637061435917076, 1256.637061435917076,
197 1256.637061435917076, 1256.637061435917076,
198 1256.637061435917076, 1256.637061435917076,
199 1256.637061435917076, 1256.637061435917076]],
200
201 [[ 1884.955592153875841, 1884.955592153875841,
202 1884.955592153875841, 1884.955592153875841,
203 1884.955592153875841, 1884.955592153875841,
204 1884.955592153875841, 1884.955592153875841,
205 1884.955592153875841, 1884.955592153875841,
206 1884.955592153875841, 1884.955592153875841,
207 1884.955592153875841, 1884.955592153875841,
208 1884.955592153875841, 1884.955592153875841,
209 1884.955592153875841, 1884.955592153875841,
210 1884.955592153875841, 1884.955592153875841,
211 1884.955592153875841, 1884.955592153875841,
212 1884.955592153875841, 1884.955592153875841,
213 1884.955592153875841, 1884.955592153875841,
214 1884.955592153875841, 1884.955592153875841,
215 1884.955592153875841, 1884.955592153875841,
216 1884.955592153875841, 1884.955592153875841,
217 1884.955592153875841, 1884.955592153875841,
218 1884.955592153875841, 1884.955592153875841,
219 1884.955592153875841, 1884.955592153875841,
220 1884.955592153875841, 1884.955592153875841,
221 1884.955592153875841, 1884.955592153875841,
222 1884.955592153875841, 1884.955592153875841,
223 1884.955592153875841, 1884.955592153875841,
224 1884.955592153875841, 1884.955592153875841,
225 1884.955592153875841, 1884.955592153875841,
226 1884.955592153875841, 1884.955592153875841,
227 1884.955592153875841, 1884.955592153875841,
228 1884.955592153875841, 1884.955592153875841,
229 1884.955592153875841, 1884.955592153875841,
230 1884.955592153875841, 1884.955592153875841,
231 1884.955592153875841, 1884.955592153875841,
232 1884.955592153875841, 1884.955592153875841,
233 1884.955592153875841, 1884.955592153875841,
234 1884.955592153875841, 1884.955592153875841,
235 1884.955592153875841, 1884.955592153875841,
236 1884.955592153875841, 1884.955592153875841,
237 1884.955592153875841, 1884.955592153875841,
238 1884.955592153875841, 1884.955592153875841,
239 1884.955592153875841, 1884.955592153875841,
240 1884.955592153875841, 1884.955592153875841]],
241
242 [[ 2513.274122871834152, 2513.274122871834152,
243 2513.274122871834152, 2513.274122871834152,
244 2513.274122871834152, 2513.274122871834152,
245 2513.274122871834152, 2513.274122871834152,
246 2513.274122871834152, 2513.274122871834152,
247 2513.274122871834152, 2513.274122871834152,
248 2513.274122871834152, 2513.274122871834152,
249 2513.274122871834152, 2513.274122871834152,
250 2513.274122871834152, 2513.274122871834152,
251 2513.274122871834152, 2513.274122871834152,
252 2513.274122871834152, 2513.274122871834152,
253 2513.274122871834152, 2513.274122871834152,
254 2513.274122871834152, 2513.274122871834152,
255 2513.274122871834152, 2513.274122871834152,
256 2513.274122871834152, 2513.274122871834152,
257 2513.274122871834152, 2513.274122871834152,
258 2513.274122871834152, 2513.274122871834152,
259 2513.274122871834152, 2513.274122871834152,
260 2513.274122871834152, 2513.274122871834152,
261 2513.274122871834152, 2513.274122871834152,
262 2513.274122871834152, 2513.274122871834152,
263 2513.274122871834152, 2513.274122871834152,
264 2513.274122871834152, 2513.274122871834152,
265 2513.274122871834152, 2513.274122871834152,
266 2513.274122871834152, 2513.274122871834152,
267 2513.274122871834152, 2513.274122871834152,
268 2513.274122871834152, 2513.274122871834152,
269 2513.274122871834152, 2513.274122871834152,
270 2513.274122871834152, 2513.274122871834152,
271 2513.274122871834152, 2513.274122871834152,
272 2513.274122871834152, 2513.274122871834152,
273 2513.274122871834152, 2513.274122871834152,
274 2513.274122871834152, 2513.274122871834152,
275 2513.274122871834152, 2513.274122871834152,
276 2513.274122871834152, 2513.274122871834152,
277 2513.274122871834152, 2513.274122871834152,
278 2513.274122871834152, 2513.274122871834152,
279 2513.274122871834152, 2513.274122871834152,
280 2513.274122871834152, 2513.274122871834152,
281 2513.274122871834152, 2513.274122871834152]],
282
283 [[ 3141.592653589792917, 3141.592653589792917,
284 3141.592653589792917, 3141.592653589792917,
285 3141.592653589792917, 3141.592653589792917,
286 3141.592653589792917, 3141.592653589792917,
287 3141.592653589792917, 3141.592653589792917,
288 3141.592653589792917, 3141.592653589792917,
289 3141.592653589792917, 3141.592653589792917,
290 3141.592653589792917, 3141.592653589792917,
291 3141.592653589792917, 3141.592653589792917,
292 3141.592653589792917, 3141.592653589792917,
293 3141.592653589792917, 3141.592653589792917,
294 3141.592653589792917, 3141.592653589792917,
295 3141.592653589792917, 3141.592653589792917,
296 3141.592653589792917, 3141.592653589792917,
297 3141.592653589792917, 3141.592653589792917,
298 3141.592653589792917, 3141.592653589792917,
299 3141.592653589792917, 3141.592653589792917,
300 3141.592653589792917, 3141.592653589792917,
301 3141.592653589792917, 3141.592653589792917,
302 3141.592653589792917, 3141.592653589792917,
303 3141.592653589792917, 3141.592653589792917,
304 3141.592653589792917, 3141.592653589792917,
305 3141.592653589792917, 3141.592653589792917,
306 3141.592653589792917, 3141.592653589792917,
307 3141.592653589792917, 3141.592653589792917,
308 3141.592653589792917, 3141.592653589792917,
309 3141.592653589792917, 3141.592653589792917,
310 3141.592653589792917, 3141.592653589792917,
311 3141.592653589792917, 3141.592653589792917,
312 3141.592653589792917, 3141.592653589792917,
313 3141.592653589792917, 3141.592653589792917,
314 3141.592653589792917, 3141.592653589792917,
315 3141.592653589792917, 3141.592653589792917,
316 3141.592653589792917, 3141.592653589792917,
317 3141.592653589792917, 3141.592653589792917,
318 3141.592653589792917, 3141.592653589792917,
319 3141.592653589792917, 3141.592653589792917,
320 3141.592653589792917, 3141.592653589792917,
321 3141.592653589792917, 3141.592653589792917,
322 3141.592653589792917, 3141.592653589792917]]]])
323 self.dw_BC = array([[[[ -5026.548245743668303, -5026.548245743668303,
324 -5026.548245743668303, -5026.548245743668303,
325 -5026.548245743668303, -5026.548245743668303,
326 -5026.548245743668303, -5026.548245743668303,
327 -5026.548245743668303, -5026.548245743668303,
328 -5026.548245743668303, -5026.548245743668303,
329 -5026.548245743668303, -5026.548245743668303,
330 -5026.548245743668303, -5026.548245743668303,
331 -5026.548245743668303, -5026.548245743668303,
332 -5026.548245743668303, -5026.548245743668303,
333 -5026.548245743668303, -5026.548245743668303,
334 -5026.548245743668303, -5026.548245743668303,
335 -5026.548245743668303, -5026.548245743668303,
336 -5026.548245743668303, -5026.548245743668303,
337 -5026.548245743668303, -5026.548245743668303,
338 -5026.548245743668303, -5026.548245743668303,
339 -5026.548245743668303, -5026.548245743668303,
340 -5026.548245743668303, -5026.548245743668303,
341 -5026.548245743668303, -5026.548245743668303,
342 -5026.548245743668303, -5026.548245743668303,
343 -5026.548245743668303, -5026.548245743668303,
344 -5026.548245743668303, -5026.548245743668303,
345 -5026.548245743668303, -5026.548245743668303,
346 -5026.548245743668303, -5026.548245743668303,
347 -5026.548245743668303, -5026.548245743668303,
348 -5026.548245743668303, -5026.548245743668303,
349 -5026.548245743668303, -5026.548245743668303,
350 -5026.548245743668303, -5026.548245743668303,
351 -5026.548245743668303, -5026.548245743668303,
352 -5026.548245743668303, -5026.548245743668303,
353 -5026.548245743668303, -5026.548245743668303,
354 -5026.548245743668303, -5026.548245743668303,
355 -5026.548245743668303, -5026.548245743668303,
356 -5026.548245743668303, -5026.548245743668303,
357 -5026.548245743668303, -5026.548245743668303,
358 -5026.548245743668303, -5026.548245743668303,
359 -5026.548245743668303, -5026.548245743668303,
360 -5026.548245743668303, -5026.548245743668303,
361 -5026.548245743668303, -5026.548245743668303,
362 -5026.548245743668303, -5026.548245743668303]],
363
364 [[ -7539.822368615503365, -7539.822368615503365,
365 -7539.822368615503365, -7539.822368615503365,
366 -7539.822368615503365, -7539.822368615503365,
367 -7539.822368615503365, -7539.822368615503365,
368 -7539.822368615503365, -7539.822368615503365,
369 -7539.822368615503365, -7539.822368615503365,
370 -7539.822368615503365, -7539.822368615503365,
371 -7539.822368615503365, -7539.822368615503365,
372 -7539.822368615503365, -7539.822368615503365,
373 -7539.822368615503365, -7539.822368615503365,
374 -7539.822368615503365, -7539.822368615503365,
375 -7539.822368615503365, -7539.822368615503365,
376 -7539.822368615503365, -7539.822368615503365,
377 -7539.822368615503365, -7539.822368615503365,
378 -7539.822368615503365, -7539.822368615503365,
379 -7539.822368615503365, -7539.822368615503365,
380 -7539.822368615503365, -7539.822368615503365,
381 -7539.822368615503365, -7539.822368615503365,
382 -7539.822368615503365, -7539.822368615503365,
383 -7539.822368615503365, -7539.822368615503365,
384 -7539.822368615503365, -7539.822368615503365,
385 -7539.822368615503365, -7539.822368615503365,
386 -7539.822368615503365, -7539.822368615503365,
387 -7539.822368615503365, -7539.822368615503365,
388 -7539.822368615503365, -7539.822368615503365,
389 -7539.822368615503365, -7539.822368615503365,
390 -7539.822368615503365, -7539.822368615503365,
391 -7539.822368615503365, -7539.822368615503365,
392 -7539.822368615503365, -7539.822368615503365,
393 -7539.822368615503365, -7539.822368615503365,
394 -7539.822368615503365, -7539.822368615503365,
395 -7539.822368615503365, -7539.822368615503365,
396 -7539.822368615503365, -7539.822368615503365,
397 -7539.822368615503365, -7539.822368615503365,
398 -7539.822368615503365, -7539.822368615503365,
399 -7539.822368615503365, -7539.822368615503365,
400 -7539.822368615503365, -7539.822368615503365,
401 -7539.822368615503365, -7539.822368615503365,
402 -7539.822368615503365, -7539.822368615503365,
403 -7539.822368615503365, -7539.822368615503365]],
404
405 [[-10053.096491487336607, -10053.096491487336607,
406 -10053.096491487336607, -10053.096491487336607,
407 -10053.096491487336607, -10053.096491487336607,
408 -10053.096491487336607, -10053.096491487336607,
409 -10053.096491487336607, -10053.096491487336607,
410 -10053.096491487336607, -10053.096491487336607,
411 -10053.096491487336607, -10053.096491487336607,
412 -10053.096491487336607, -10053.096491487336607,
413 -10053.096491487336607, -10053.096491487336607,
414 -10053.096491487336607, -10053.096491487336607,
415 -10053.096491487336607, -10053.096491487336607,
416 -10053.096491487336607, -10053.096491487336607,
417 -10053.096491487336607, -10053.096491487336607,
418 -10053.096491487336607, -10053.096491487336607,
419 -10053.096491487336607, -10053.096491487336607,
420 -10053.096491487336607, -10053.096491487336607,
421 -10053.096491487336607, -10053.096491487336607,
422 -10053.096491487336607, -10053.096491487336607,
423 -10053.096491487336607, -10053.096491487336607,
424 -10053.096491487336607, -10053.096491487336607,
425 -10053.096491487336607, -10053.096491487336607,
426 -10053.096491487336607, -10053.096491487336607,
427 -10053.096491487336607, -10053.096491487336607,
428 -10053.096491487336607, -10053.096491487336607,
429 -10053.096491487336607, -10053.096491487336607,
430 -10053.096491487336607, -10053.096491487336607,
431 -10053.096491487336607, -10053.096491487336607,
432 -10053.096491487336607, -10053.096491487336607,
433 -10053.096491487336607, -10053.096491487336607,
434 -10053.096491487336607, -10053.096491487336607,
435 -10053.096491487336607, -10053.096491487336607,
436 -10053.096491487336607, -10053.096491487336607,
437 -10053.096491487336607, -10053.096491487336607,
438 -10053.096491487336607, -10053.096491487336607,
439 -10053.096491487336607, -10053.096491487336607,
440 -10053.096491487336607, -10053.096491487336607,
441 -10053.096491487336607, -10053.096491487336607,
442 -10053.096491487336607, -10053.096491487336607,
443 -10053.096491487336607, -10053.096491487336607,
444 -10053.096491487336607, -10053.096491487336607]],
445
446 [[-12566.370614359171668, -12566.370614359171668,
447 -12566.370614359171668, -12566.370614359171668,
448 -12566.370614359171668, -12566.370614359171668,
449 -12566.370614359171668, -12566.370614359171668,
450 -12566.370614359171668, -12566.370614359171668,
451 -12566.370614359171668, -12566.370614359171668,
452 -12566.370614359171668, -12566.370614359171668,
453 -12566.370614359171668, -12566.370614359171668,
454 -12566.370614359171668, -12566.370614359171668,
455 -12566.370614359171668, -12566.370614359171668,
456 -12566.370614359171668, -12566.370614359171668,
457 -12566.370614359171668, -12566.370614359171668,
458 -12566.370614359171668, -12566.370614359171668,
459 -12566.370614359171668, -12566.370614359171668,
460 -12566.370614359171668, -12566.370614359171668,
461 -12566.370614359171668, -12566.370614359171668,
462 -12566.370614359171668, -12566.370614359171668,
463 -12566.370614359171668, -12566.370614359171668,
464 -12566.370614359171668, -12566.370614359171668,
465 -12566.370614359171668, -12566.370614359171668,
466 -12566.370614359171668, -12566.370614359171668,
467 -12566.370614359171668, -12566.370614359171668,
468 -12566.370614359171668, -12566.370614359171668,
469 -12566.370614359171668, -12566.370614359171668,
470 -12566.370614359171668, -12566.370614359171668,
471 -12566.370614359171668, -12566.370614359171668,
472 -12566.370614359171668, -12566.370614359171668,
473 -12566.370614359171668, -12566.370614359171668,
474 -12566.370614359171668, -12566.370614359171668,
475 -12566.370614359171668, -12566.370614359171668,
476 -12566.370614359171668, -12566.370614359171668,
477 -12566.370614359171668, -12566.370614359171668,
478 -12566.370614359171668, -12566.370614359171668,
479 -12566.370614359171668, -12566.370614359171668,
480 -12566.370614359171668, -12566.370614359171668,
481 -12566.370614359171668, -12566.370614359171668,
482 -12566.370614359171668, -12566.370614359171668,
483 -12566.370614359171668, -12566.370614359171668,
484 -12566.370614359171668, -12566.370614359171668,
485 -12566.370614359171668, -12566.370614359171668]]]])
486
487 self.dwH_AB = array([[[[ 764.516785511118997, 764.516785511118997,
488 764.516785511118997, 764.516785511118997,
489 764.516785511118997, 764.516785511118997,
490 764.516785511118997, 764.516785511118997,
491 764.516785511118997, 764.516785511118997,
492 764.516785511118997, 764.516785511118997,
493 764.516785511118997, 764.516785511118997,
494 764.516785511118997, 764.516785511118997,
495 764.516785511118997, 764.516785511118997,
496 764.516785511118997, 764.516785511118997,
497 764.516785511118997, 764.516785511118997,
498 764.516785511118997, 764.516785511118997,
499 764.516785511118997, 764.516785511118997,
500 764.516785511118997, 764.516785511118997,
501 764.516785511118997, 764.516785511118997,
502 764.516785511118997, 764.516785511118997,
503 764.516785511118997, 764.516785511118997,
504 764.516785511118997, 764.516785511118997,
505 764.516785511118997, 764.516785511118997,
506 764.516785511118997, 764.516785511118997,
507 764.516785511118997, 764.516785511118997,
508 764.516785511118997, 764.516785511118997,
509 764.516785511118997, 764.516785511118997,
510 764.516785511118997, 764.516785511118997,
511 764.516785511118997, 764.516785511118997,
512 764.516785511118997, 764.516785511118997,
513 764.516785511118997, 764.516785511118997,
514 764.516785511118997, 764.516785511118997,
515 764.516785511118997, 764.516785511118997,
516 764.516785511118997, 764.516785511118997,
517 764.516785511118997, 764.516785511118997,
518 764.516785511118997, 764.516785511118997,
519 764.516785511118997, 764.516785511118997,
520 764.516785511118997, 764.516785511118997,
521 764.516785511118997, 764.516785511118997,
522 764.516785511118997, 764.516785511118997,
523 764.516785511118997, 764.516785511118997,
524 764.516785511118997, 764.516785511118997,
525 764.516785511118997, 764.516785511118997,
526 764.516785511118997, 764.516785511118997]],
527
528 [[ 1146.775178266678722, 1146.775178266678722,
529 1146.775178266678722, 1146.775178266678722,
530 1146.775178266678722, 1146.775178266678722,
531 1146.775178266678722, 1146.775178266678722,
532 1146.775178266678722, 1146.775178266678722,
533 1146.775178266678722, 1146.775178266678722,
534 1146.775178266678722, 1146.775178266678722,
535 1146.775178266678722, 1146.775178266678722,
536 1146.775178266678722, 1146.775178266678722,
537 1146.775178266678722, 1146.775178266678722,
538 1146.775178266678722, 1146.775178266678722,
539 1146.775178266678722, 1146.775178266678722,
540 1146.775178266678722, 1146.775178266678722,
541 1146.775178266678722, 1146.775178266678722,
542 1146.775178266678722, 1146.775178266678722,
543 1146.775178266678722, 1146.775178266678722,
544 1146.775178266678722, 1146.775178266678722,
545 1146.775178266678722, 1146.775178266678722,
546 1146.775178266678722, 1146.775178266678722,
547 1146.775178266678722, 1146.775178266678722,
548 1146.775178266678722, 1146.775178266678722,
549 1146.775178266678722, 1146.775178266678722,
550 1146.775178266678722, 1146.775178266678722,
551 1146.775178266678722, 1146.775178266678722,
552 1146.775178266678722, 1146.775178266678722,
553 1146.775178266678722, 1146.775178266678722,
554 1146.775178266678722, 1146.775178266678722,
555 1146.775178266678722, 1146.775178266678722,
556 1146.775178266678722, 1146.775178266678722,
557 1146.775178266678722, 1146.775178266678722,
558 1146.775178266678722, 1146.775178266678722,
559 1146.775178266678722, 1146.775178266678722,
560 1146.775178266678722, 1146.775178266678722,
561 1146.775178266678722, 1146.775178266678722,
562 1146.775178266678722, 1146.775178266678722,
563 1146.775178266678722, 1146.775178266678722,
564 1146.775178266678722, 1146.775178266678722,
565 1146.775178266678722, 1146.775178266678722,
566 1146.775178266678722, 1146.775178266678722,
567 1146.775178266678722, 1146.775178266678722]],
568
569 [[ 1529.033571022237993, 1529.033571022237993,
570 1529.033571022237993, 1529.033571022237993,
571 1529.033571022237993, 1529.033571022237993,
572 1529.033571022237993, 1529.033571022237993,
573 1529.033571022237993, 1529.033571022237993,
574 1529.033571022237993, 1529.033571022237993,
575 1529.033571022237993, 1529.033571022237993,
576 1529.033571022237993, 1529.033571022237993,
577 1529.033571022237993, 1529.033571022237993,
578 1529.033571022237993, 1529.033571022237993,
579 1529.033571022237993, 1529.033571022237993,
580 1529.033571022237993, 1529.033571022237993,
581 1529.033571022237993, 1529.033571022237993,
582 1529.033571022237993, 1529.033571022237993,
583 1529.033571022237993, 1529.033571022237993,
584 1529.033571022237993, 1529.033571022237993,
585 1529.033571022237993, 1529.033571022237993,
586 1529.033571022237993, 1529.033571022237993,
587 1529.033571022237993, 1529.033571022237993,
588 1529.033571022237993, 1529.033571022237993,
589 1529.033571022237993, 1529.033571022237993,
590 1529.033571022237993, 1529.033571022237993,
591 1529.033571022237993, 1529.033571022237993,
592 1529.033571022237993, 1529.033571022237993,
593 1529.033571022237993, 1529.033571022237993,
594 1529.033571022237993, 1529.033571022237993,
595 1529.033571022237993, 1529.033571022237993,
596 1529.033571022237993, 1529.033571022237993,
597 1529.033571022237993, 1529.033571022237993,
598 1529.033571022237993, 1529.033571022237993,
599 1529.033571022237993, 1529.033571022237993,
600 1529.033571022237993, 1529.033571022237993,
601 1529.033571022237993, 1529.033571022237993,
602 1529.033571022237993, 1529.033571022237993,
603 1529.033571022237993, 1529.033571022237993,
604 1529.033571022237993, 1529.033571022237993,
605 1529.033571022237993, 1529.033571022237993,
606 1529.033571022237993, 1529.033571022237993,
607 1529.033571022237993, 1529.033571022237993,
608 1529.033571022237993, 1529.033571022237993]],
609
610 [[ 1911.291963777797719, 1911.291963777797719,
611 1911.291963777797719, 1911.291963777797719,
612 1911.291963777797719, 1911.291963777797719,
613 1911.291963777797719, 1911.291963777797719,
614 1911.291963777797719, 1911.291963777797719,
615 1911.291963777797719, 1911.291963777797719,
616 1911.291963777797719, 1911.291963777797719,
617 1911.291963777797719, 1911.291963777797719,
618 1911.291963777797719, 1911.291963777797719,
619 1911.291963777797719, 1911.291963777797719,
620 1911.291963777797719, 1911.291963777797719,
621 1911.291963777797719, 1911.291963777797719,
622 1911.291963777797719, 1911.291963777797719,
623 1911.291963777797719, 1911.291963777797719,
624 1911.291963777797719, 1911.291963777797719,
625 1911.291963777797719, 1911.291963777797719,
626 1911.291963777797719, 1911.291963777797719,
627 1911.291963777797719, 1911.291963777797719,
628 1911.291963777797719, 1911.291963777797719,
629 1911.291963777797719, 1911.291963777797719,
630 1911.291963777797719, 1911.291963777797719,
631 1911.291963777797719, 1911.291963777797719,
632 1911.291963777797719, 1911.291963777797719,
633 1911.291963777797719, 1911.291963777797719,
634 1911.291963777797719, 1911.291963777797719,
635 1911.291963777797719, 1911.291963777797719,
636 1911.291963777797719, 1911.291963777797719,
637 1911.291963777797719, 1911.291963777797719,
638 1911.291963777797719, 1911.291963777797719,
639 1911.291963777797719, 1911.291963777797719,
640 1911.291963777797719, 1911.291963777797719,
641 1911.291963777797719, 1911.291963777797719,
642 1911.291963777797719, 1911.291963777797719,
643 1911.291963777797719, 1911.291963777797719,
644 1911.291963777797719, 1911.291963777797719,
645 1911.291963777797719, 1911.291963777797719,
646 1911.291963777797719, 1911.291963777797719,
647 1911.291963777797719, 1911.291963777797719,
648 1911.291963777797719, 1911.291963777797719,
649 1911.291963777797719, 1911.291963777797719]]]])
650
651 self.dwH_BC = array([[[[-2803.228213540769957, -2803.228213540769957,
652 -2803.228213540769957, -2803.228213540769957,
653 -2803.228213540769957, -2803.228213540769957,
654 -2803.228213540769957, -2803.228213540769957,
655 -2803.228213540769957, -2803.228213540769957,
656 -2803.228213540769957, -2803.228213540769957,
657 -2803.228213540769957, -2803.228213540769957,
658 -2803.228213540769957, -2803.228213540769957,
659 -2803.228213540769957, -2803.228213540769957,
660 -2803.228213540769957, -2803.228213540769957,
661 -2803.228213540769957, -2803.228213540769957,
662 -2803.228213540769957, -2803.228213540769957,
663 -2803.228213540769957, -2803.228213540769957,
664 -2803.228213540769957, -2803.228213540769957,
665 -2803.228213540769957, -2803.228213540769957,
666 -2803.228213540769957, -2803.228213540769957,
667 -2803.228213540769957, -2803.228213540769957,
668 -2803.228213540769957, -2803.228213540769957,
669 -2803.228213540769957, -2803.228213540769957,
670 -2803.228213540769957, -2803.228213540769957,
671 -2803.228213540769957, -2803.228213540769957,
672 -2803.228213540769957, -2803.228213540769957,
673 -2803.228213540769957, -2803.228213540769957,
674 -2803.228213540769957, -2803.228213540769957,
675 -2803.228213540769957, -2803.228213540769957,
676 -2803.228213540769957, -2803.228213540769957,
677 -2803.228213540769957, -2803.228213540769957,
678 -2803.228213540769957, -2803.228213540769957,
679 -2803.228213540769957, -2803.228213540769957,
680 -2803.228213540769957, -2803.228213540769957,
681 -2803.228213540769957, -2803.228213540769957,
682 -2803.228213540769957, -2803.228213540769957,
683 -2803.228213540769957, -2803.228213540769957,
684 -2803.228213540769957, -2803.228213540769957,
685 -2803.228213540769957, -2803.228213540769957,
686 -2803.228213540769957, -2803.228213540769957,
687 -2803.228213540769957, -2803.228213540769957,
688 -2803.228213540769957, -2803.228213540769957,
689 -2803.228213540769957, -2803.228213540769957,
690 -2803.228213540769957, -2803.228213540769957]],
691
692 [[-4204.842320311155163, -4204.842320311155163,
693 -4204.842320311155163, -4204.842320311155163,
694 -4204.842320311155163, -4204.842320311155163,
695 -4204.842320311155163, -4204.842320311155163,
696 -4204.842320311155163, -4204.842320311155163,
697 -4204.842320311155163, -4204.842320311155163,
698 -4204.842320311155163, -4204.842320311155163,
699 -4204.842320311155163, -4204.842320311155163,
700 -4204.842320311155163, -4204.842320311155163,
701 -4204.842320311155163, -4204.842320311155163,
702 -4204.842320311155163, -4204.842320311155163,
703 -4204.842320311155163, -4204.842320311155163,
704 -4204.842320311155163, -4204.842320311155163,
705 -4204.842320311155163, -4204.842320311155163,
706 -4204.842320311155163, -4204.842320311155163,
707 -4204.842320311155163, -4204.842320311155163,
708 -4204.842320311155163, -4204.842320311155163,
709 -4204.842320311155163, -4204.842320311155163,
710 -4204.842320311155163, -4204.842320311155163,
711 -4204.842320311155163, -4204.842320311155163,
712 -4204.842320311155163, -4204.842320311155163,
713 -4204.842320311155163, -4204.842320311155163,
714 -4204.842320311155163, -4204.842320311155163,
715 -4204.842320311155163, -4204.842320311155163,
716 -4204.842320311155163, -4204.842320311155163,
717 -4204.842320311155163, -4204.842320311155163,
718 -4204.842320311155163, -4204.842320311155163,
719 -4204.842320311155163, -4204.842320311155163,
720 -4204.842320311155163, -4204.842320311155163,
721 -4204.842320311155163, -4204.842320311155163,
722 -4204.842320311155163, -4204.842320311155163,
723 -4204.842320311155163, -4204.842320311155163,
724 -4204.842320311155163, -4204.842320311155163,
725 -4204.842320311155163, -4204.842320311155163,
726 -4204.842320311155163, -4204.842320311155163,
727 -4204.842320311155163, -4204.842320311155163,
728 -4204.842320311155163, -4204.842320311155163,
729 -4204.842320311155163, -4204.842320311155163,
730 -4204.842320311155163, -4204.842320311155163,
731 -4204.842320311155163, -4204.842320311155163]],
732
733 [[-5606.456427081539914, -5606.456427081539914,
734 -5606.456427081539914, -5606.456427081539914,
735 -5606.456427081539914, -5606.456427081539914,
736 -5606.456427081539914, -5606.456427081539914,
737 -5606.456427081539914, -5606.456427081539914,
738 -5606.456427081539914, -5606.456427081539914,
739 -5606.456427081539914, -5606.456427081539914,
740 -5606.456427081539914, -5606.456427081539914,
741 -5606.456427081539914, -5606.456427081539914,
742 -5606.456427081539914, -5606.456427081539914,
743 -5606.456427081539914, -5606.456427081539914,
744 -5606.456427081539914, -5606.456427081539914,
745 -5606.456427081539914, -5606.456427081539914,
746 -5606.456427081539914, -5606.456427081539914,
747 -5606.456427081539914, -5606.456427081539914,
748 -5606.456427081539914, -5606.456427081539914,
749 -5606.456427081539914, -5606.456427081539914,
750 -5606.456427081539914, -5606.456427081539914,
751 -5606.456427081539914, -5606.456427081539914,
752 -5606.456427081539914, -5606.456427081539914,
753 -5606.456427081539914, -5606.456427081539914,
754 -5606.456427081539914, -5606.456427081539914,
755 -5606.456427081539914, -5606.456427081539914,
756 -5606.456427081539914, -5606.456427081539914,
757 -5606.456427081539914, -5606.456427081539914,
758 -5606.456427081539914, -5606.456427081539914,
759 -5606.456427081539914, -5606.456427081539914,
760 -5606.456427081539914, -5606.456427081539914,
761 -5606.456427081539914, -5606.456427081539914,
762 -5606.456427081539914, -5606.456427081539914,
763 -5606.456427081539914, -5606.456427081539914,
764 -5606.456427081539914, -5606.456427081539914,
765 -5606.456427081539914, -5606.456427081539914,
766 -5606.456427081539914, -5606.456427081539914,
767 -5606.456427081539914, -5606.456427081539914,
768 -5606.456427081539914, -5606.456427081539914,
769 -5606.456427081539914, -5606.456427081539914,
770 -5606.456427081539914, -5606.456427081539914,
771 -5606.456427081539914, -5606.456427081539914,
772 -5606.456427081539914, -5606.456427081539914]],
773
774 [[-7008.070533851924665, -7008.070533851924665,
775 -7008.070533851924665, -7008.070533851924665,
776 -7008.070533851924665, -7008.070533851924665,
777 -7008.070533851924665, -7008.070533851924665,
778 -7008.070533851924665, -7008.070533851924665,
779 -7008.070533851924665, -7008.070533851924665,
780 -7008.070533851924665, -7008.070533851924665,
781 -7008.070533851924665, -7008.070533851924665,
782 -7008.070533851924665, -7008.070533851924665,
783 -7008.070533851924665, -7008.070533851924665,
784 -7008.070533851924665, -7008.070533851924665,
785 -7008.070533851924665, -7008.070533851924665,
786 -7008.070533851924665, -7008.070533851924665,
787 -7008.070533851924665, -7008.070533851924665,
788 -7008.070533851924665, -7008.070533851924665,
789 -7008.070533851924665, -7008.070533851924665,
790 -7008.070533851924665, -7008.070533851924665,
791 -7008.070533851924665, -7008.070533851924665,
792 -7008.070533851924665, -7008.070533851924665,
793 -7008.070533851924665, -7008.070533851924665,
794 -7008.070533851924665, -7008.070533851924665,
795 -7008.070533851924665, -7008.070533851924665,
796 -7008.070533851924665, -7008.070533851924665,
797 -7008.070533851924665, -7008.070533851924665,
798 -7008.070533851924665, -7008.070533851924665,
799 -7008.070533851924665, -7008.070533851924665,
800 -7008.070533851924665, -7008.070533851924665,
801 -7008.070533851924665, -7008.070533851924665,
802 -7008.070533851924665, -7008.070533851924665,
803 -7008.070533851924665, -7008.070533851924665,
804 -7008.070533851924665, -7008.070533851924665,
805 -7008.070533851924665, -7008.070533851924665,
806 -7008.070533851924665, -7008.070533851924665,
807 -7008.070533851924665, -7008.070533851924665,
808 -7008.070533851924665, -7008.070533851924665,
809 -7008.070533851924665, -7008.070533851924665,
810 -7008.070533851924665, -7008.070533851924665,
811 -7008.070533851924665, -7008.070533851924665,
812 -7008.070533851924665, -7008.070533851924665,
813 -7008.070533851924665, -7008.070533851924665]]]])
814 self.kex_AB = 500.0
815 self.kex_BC = 2000.0
816 self.kex_AC = 1000.0
817 self.inv_tcpmg = array([[[[ 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
818 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
819 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
820 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
821 50., 50., 50., 50., 50., 50., 0., 0., 0., 0., 0.,
822 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
823 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
824 0., 0., 0.]],
825
826 [[ 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
827 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
828 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
829 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
830 50., 50., 50., 50., 50., 50., 0., 0., 0., 0., 0.,
831 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
832 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
833 0., 0., 0.]],
834
835 [[ 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
836 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
837 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
838 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
839 50., 50., 50., 50., 50., 50., 0., 0., 0., 0., 0.,
840 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
841 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
842 0., 0., 0.]],
843
844 [[ 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
845 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
846 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
847 50., 50., 50., 50., 50., 50., 50., 50., 50., 50., 50.,
848 50., 50., 50., 50., 50., 50., 0., 0., 0., 0., 0.,
849 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
850 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
851 0., 0., 0.]]]])
852 self.tcp = array([[[[ 0.005 , 0.0025 , 0.001666666666667,
853 0.00125 , 0.001 , 0.000833333333333,
854 0.000714285714286, 0.000625 , 0.000555555555556,
855 0.0005 , 0.000454545454545, 0.000416666666667,
856 0.000384615384615, 0.000357142857143, 0.000333333333333,
857 0.0003125 , 0.000294117647059, 0.000277777777778,
858 0.000263157894737, 0.00025 , 0.000238095238095,
859 0.000227272727273, 0.000217391304348, 0.000208333333333,
860 0.0002 , 0.000192307692308, 0.000185185185185,
861 0.000178571428571, 0.000172413793103, 0.000166666666667,
862 0.000161290322581, 0.00015625 , 0.000151515151515,
863 0.000147058823529, 0.000142857142857, 0.000138888888889,
864 0.000135135135135, 0.000131578947368, 0.000128205128205,
865 0.000125 , 0.000121951219512, 0.000119047619048,
866 0.000116279069767, 0.000113636363636, 0.000111111111111,
867 0.000108695652174, 0.000106382978723, 0.000104166666667,
868 0.000102040816327, 0.0001 , 0. ,
869 0. , 0. , 0. ,
870 0. , 0. , 0. ,
871 0. , 0. , 0. ,
872 0. , 0. , 0. ,
873 0. , 0. , 0. ,
874 0. , 0. , 0. ,
875 0. , 0. , 0. ,
876 0. , 0. , 0. ,
877 0. , 0. , 0. ,
878 0. , 0. ]],
879
880 [[ 0.005 , 0.0025 , 0.001666666666667,
881 0.00125 , 0.001 , 0.000833333333333,
882 0.000714285714286, 0.000625 , 0.000555555555556,
883 0.0005 , 0.000454545454545, 0.000416666666667,
884 0.000384615384615, 0.000357142857143, 0.000333333333333,
885 0.0003125 , 0.000294117647059, 0.000277777777778,
886 0.000263157894737, 0.00025 , 0.000238095238095,
887 0.000227272727273, 0.000217391304348, 0.000208333333333,
888 0.0002 , 0.000192307692308, 0.000185185185185,
889 0.000178571428571, 0.000172413793103, 0.000166666666667,
890 0.000161290322581, 0.00015625 , 0.000151515151515,
891 0.000147058823529, 0.000142857142857, 0.000138888888889,
892 0.000135135135135, 0.000131578947368, 0.000128205128205,
893 0.000125 , 0.000121951219512, 0.000119047619048,
894 0.000116279069767, 0.000113636363636, 0.000111111111111,
895 0.000108695652174, 0.000106382978723, 0.000104166666667,
896 0.000102040816327, 0.0001 , 0. ,
897 0. , 0. , 0. ,
898 0. , 0. , 0. ,
899 0. , 0. , 0. ,
900 0. , 0. , 0. ,
901 0. , 0. , 0. ,
902 0. , 0. , 0. ,
903 0. , 0. , 0. ,
904 0. , 0. , 0. ,
905 0. , 0. , 0. ,
906 0. , 0. ]],
907
908 [[ 0.005 , 0.0025 , 0.001666666666667,
909 0.00125 , 0.001 , 0.000833333333333,
910 0.000714285714286, 0.000625 , 0.000555555555556,
911 0.0005 , 0.000454545454545, 0.000416666666667,
912 0.000384615384615, 0.000357142857143, 0.000333333333333,
913 0.0003125 , 0.000294117647059, 0.000277777777778,
914 0.000263157894737, 0.00025 , 0.000238095238095,
915 0.000227272727273, 0.000217391304348, 0.000208333333333,
916 0.0002 , 0.000192307692308, 0.000185185185185,
917 0.000178571428571, 0.000172413793103, 0.000166666666667,
918 0.000161290322581, 0.00015625 , 0.000151515151515,
919 0.000147058823529, 0.000142857142857, 0.000138888888889,
920 0.000135135135135, 0.000131578947368, 0.000128205128205,
921 0.000125 , 0.000121951219512, 0.000119047619048,
922 0.000116279069767, 0.000113636363636, 0.000111111111111,
923 0.000108695652174, 0.000106382978723, 0.000104166666667,
924 0.000102040816327, 0.0001 , 0. ,
925 0. , 0. , 0. ,
926 0. , 0. , 0. ,
927 0. , 0. , 0. ,
928 0. , 0. , 0. ,
929 0. , 0. , 0. ,
930 0. , 0. , 0. ,
931 0. , 0. , 0. ,
932 0. , 0. , 0. ,
933 0. , 0. , 0. ,
934 0. , 0. ]],
935
936 [[ 0.005 , 0.0025 , 0.001666666666667,
937 0.00125 , 0.001 , 0.000833333333333,
938 0.000714285714286, 0.000625 , 0.000555555555556,
939 0.0005 , 0.000454545454545, 0.000416666666667,
940 0.000384615384615, 0.000357142857143, 0.000333333333333,
941 0.0003125 , 0.000294117647059, 0.000277777777778,
942 0.000263157894737, 0.00025 , 0.000238095238095,
943 0.000227272727273, 0.000217391304348, 0.000208333333333,
944 0.0002 , 0.000192307692308, 0.000185185185185,
945 0.000178571428571, 0.000172413793103, 0.000166666666667,
946 0.000161290322581, 0.00015625 , 0.000151515151515,
947 0.000147058823529, 0.000142857142857, 0.000138888888889,
948 0.000135135135135, 0.000131578947368, 0.000128205128205,
949 0.000125 , 0.000121951219512, 0.000119047619048,
950 0.000116279069767, 0.000113636363636, 0.000111111111111,
951 0.000108695652174, 0.000106382978723, 0.000104166666667,
952 0.000102040816327, 0.0001 , 0. ,
953 0. , 0. , 0. ,
954 0. , 0. , 0. ,
955 0. , 0. , 0. ,
956 0. , 0. , 0. ,
957 0. , 0. , 0. ,
958 0. , 0. , 0. ,
959 0. , 0. , 0. ,
960 0. , 0. , 0. ,
961 0. , 0. , 0. ,
962 0. , 0. ]]]])
963 self.back_calc = array([[[[ 105.465130755614368, 104.042707539434943, 103.236522151168145,
964 101.241760501741368, 102.177693763105637, 101.271226767091861,
965 99.655513326683703, 100.173149424447033, 102.350262471078565,
966 104.736498244558589, 106.309167799583904, 107.115597853339921,
967 107.213812890830923, 106.92699609238494 , 106.353039273978524,
968 105.663097189073696, 104.893939951800405, 104.12086480872405 ,
969 103.351551547094829, 102.617870273970823, 101.917169400320432,
970 101.262359496147241, 100.647905487734462, 100.078407319932339,
971 99.548131031951044, 99.058049306637614, 98.603089443554339,
972 98.182698645073231, 97.792649079381135, 97.431816773388078,
973 97.096766520526188, 96.786233745696151, 96.497455919599133,
974 96.229214870279094, 95.979295195234513, 95.746598061970218,
975 95.529344635040516, 95.326575475401128, 95.136857348597232,
976 94.959366680883434, 94.792943290840469, 94.636886426261086,
977 94.490251937114081, 94.352445999707101, 94.22269594627295 ,
978 94.100499178311438, 93.985219749085772, 93.876432033876611,
979 93.773609641473669, 93.676391526141146, 0. ,
980 0. , 0. , 0. ,
981 0. , 0. , 0. ,
982 0. , 0. , 0. ,
983 0. , 0. , 0. ,
984 0. , 0. , 0. ,
985 0. , 0. , 0. ,
986 0. , 0. , 0. ,
987 0. , 0. , 0. ,
988 0. , 0. , 0. ,
989 0. , 0. ]],
990
991 [[ 120.287133737199881, 119.518159476124779, 118.524338753085502,
992 118.067520837287617, 116.924450261129678, 114.631983039295065,
993 116.353417238519754, 116.655796886645959, 114.13165367565648 ,
994 111.567800462356587, 111.115724985737813, 112.837427738883363,
995 115.506176747110857, 118.055631512485377, 119.933100915347609,
996 121.183201536122681, 121.861210479651504, 122.177323205824635,
997 122.204942008882966, 122.07332867972454 , 121.815593652192604,
998 121.499263795711016, 121.135420910556803, 120.758929148120558,
999 120.371572411978335, 119.99163217063699 , 119.617391031834373,
1000 119.258577182809717, 118.912397680843512, 118.584054342201341,
1001 118.270650454227848, 117.974936036959903, 117.69426929374842 ,
1002 117.43004912442639 , 117.179991823721679, 116.944744718073594,
1003 116.722382442578336, 116.513134461940808, 116.315396455204308,
1004 116.12916883456505 , 115.953121190633567, 115.787132854908208,
1005 115.630101591537411, 115.481847655036191, 115.3414566516675 ,
1006 115.208725226122127, 115.08289281265742 , 114.963752426326309,
1007 114.850669286714606, 114.743443663330467, 0. ,
1008 0. , 0. , 0. ,
1009 0. , 0. , 0. ,
1010 0. , 0. , 0. ,
1011 0. , 0. , 0. ,
1012 0. , 0. , 0. ,
1013 0. , 0. , 0. ,
1014 0. , 0. , 0. ,
1015 0. , 0. , 0. ,
1016 0. , 0. , 0. ,
1017 0. , 0. ]],
1018
1019 [[ 129.92115804111495 , 129.500108932148464, 128.949277054807254,
1020 128.383816107598591, 128.303930781717241, 127.995351116526606,
1021 125.89637750932134 , 124.438507694837952, 126.447208093921944,
1022 127.647124035218567, 126.264385528177598, 123.256370173959823,
1023 120.453742994435657, 119.437129654549196, 120.507883560068279,
1024 122.834280518382428, 125.394552902335832, 127.610049881724706,
1025 129.263675393752578, 130.430563222285713, 131.175097923692022,
1026 131.632964048434786, 131.860563026044844, 131.944279721263882,
1027 131.911807835944899, 131.811928038273379, 131.656138264887545,
1028 131.471901979718751, 131.263117736595291, 131.045580500290811,
1029 130.819830994752522, 130.595170705884129, 130.370747043133548,
1030 130.152148133579146, 129.938023638940223, 129.731780077454886,
1031 129.531970647265808, 129.340689451500594, 129.156564870780699,
1032 128.98088580619762 , 128.812420512875946, 128.651956113681308,
1033 128.498416205920762, 128.352271579976644, 128.212594916086772,
1034 128.07965609099503 , 127.952662021771829, 127.831754602366047,
1035 127.716257757614997, 127.606232051795487, 0. ,
1036 0. , 0. , 0. ,
1037 0. , 0. , 0. ,
1038 0. , 0. , 0. ,
1039 0. , 0. , 0. ,
1040 0. , 0. , 0. ,
1041 0. , 0. , 0. ,
1042 0. , 0. , 0. ,
1043 0. , 0. , 0. ,
1044 0. , 0. , 0. ,
1045 0. , 0. ]],
1046
1047 [[ 135.631068708819186, 135.336515276208303, 135.116760936415631,
1048 134.811908231880381, 134.317982194993533, 134.503459261963087,
1049 133.999005946859 , 133.810386986666686, 131.195208527296614,
1050 130.791255682473064, 132.752846219057972, 134.219412845316469,
1051 133.91657850121851 , 131.892571309240424, 128.842801014884088,
1052 126.064993180705869, 124.753418794368415, 125.318622395484525,
1053 127.237582101023094, 129.627434949689274, 131.857145915200732,
1054 133.679722786026815, 135.045284048098182, 136.03731512062339 ,
1055 136.717332165036652, 137.177266641818591, 137.460318566100796,
1056 137.624845888398539, 137.692936568914121, 137.698940894509462,
1057 137.653127412775802, 137.575848217939438, 137.471480323167981,
1058 137.352378459386614, 137.220069884526737, 137.082250146055628,
1059 136.93909581388462 , 136.795510875902238, 136.65105892343189 ,
1060 136.508924525045018, 136.368422639016302, 136.231652286447456,
1061 136.097858779535784, 135.968441018811234, 135.842660314740556,
1062 135.721455794620084, 135.604142975805132, 135.491354449276258,
1063 135.382474063238163, 135.277927474700078, 0. ,
1064 0. , 0. , 0. ,
1065 0. , 0. , 0. ,
1066 0. , 0. , 0. ,
1067 0. , 0. , 0. ,
1068 0. , 0. , 0. ,
1069 0. , 0. , 0. ,
1070 0. , 0. , 0. ,
1071 0. , 0. , 0. ,
1072 0. , 0. , 0. ,
1073 0. , 0. ]]]])
1074 self.num_points = array([[[50],
1075 [50],
1076 [50],
1077 [50]]], dtype=int16)
1078 self.power = array([[[[ 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11.,
1079 12., 13., 14., 15., 16., 17., 18., 19., 20., 21., 22.,
1080 23., 24., 25., 26., 27., 28., 29., 30., 31., 32., 33.,
1081 34., 35., 36., 37., 38., 39., 40., 41., 42., 43., 44.,
1082 45., 46., 47., 48., 49., 50., 0., 0., 0., 0., 0.,
1083 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
1084 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
1085 0., 0., 0.]],
1086
1087 [[ 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11.,
1088 12., 13., 14., 15., 16., 17., 18., 19., 20., 21., 22.,
1089 23., 24., 25., 26., 27., 28., 29., 30., 31., 32., 33.,
1090 34., 35., 36., 37., 38., 39., 40., 41., 42., 43., 44.,
1091 45., 46., 47., 48., 49., 50., 0., 0., 0., 0., 0.,
1092 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
1093 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
1094 0., 0., 0.]],
1095
1096 [[ 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11.,
1097 12., 13., 14., 15., 16., 17., 18., 19., 20., 21., 22.,
1098 23., 24., 25., 26., 27., 28., 29., 30., 31., 32., 33.,
1099 34., 35., 36., 37., 38., 39., 40., 41., 42., 43., 44.,
1100 45., 46., 47., 48., 49., 50., 0., 0., 0., 0., 0.,
1101 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
1102 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
1103 0., 0., 0.]],
1104
1105 [[ 1., 2., 3., 4., 5., 6., 7., 8., 9., 10., 11.,
1106 12., 13., 14., 15., 16., 17., 18., 19., 20., 21., 22.,
1107 23., 24., 25., 26., 27., 28., 29., 30., 31., 32., 33.,
1108 34., 35., 36., 37., 38., 39., 40., 41., 42., 43., 44.,
1109 45., 46., 47., 48., 49., 50., 0., 0., 0., 0., 0.,
1110 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
1111 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
1112 0., 0., 0.]]]])
1113
1114
1116 """Test the r2eff_ns_mmq_3site_mq() function for no exchange when pA = 1.0 and pB = 0.0."""
1117
1118
1119 self.pA = 1.0
1120 self.pB = 0.0
1121
1122
1123 r2eff_ns_mmq_3site_mq(M0=self.M0, F_vector=self.F_vector, R20A=self.R20A, R20B=self.R20B, R20C=self.R20C, pA=self.pA, pB=self.pB, dw_AB=self.dw_AB, dw_BC=self.dw_BC, dwH_AB=self.dwH_AB, dwH_BC=self.dwH_BC, kex_AB=self.kex_AB, kex_BC=self.kex_BC, kex_AC=self.kex_AC, inv_tcpmg=self.inv_tcpmg, tcp=self.tcp, back_calc=self.back_calc, num_points=self.num_points, power=self.power)
1124
1125
1127 """Test the r2eff_ns_mmq_3site_sq_dq_zq() function for no exchange when pA = 1.0 and pB = 0.0."""
1128
1129
1130 self.pA = 1.0
1131 self.pB = 0.0
1132
1133
1134 r2eff_ns_mmq_3site_sq_dq_zq(M0=self.M0, F_vector=self.F_vector, R20A=self.R20A, R20B=self.R20B, R20C=self.R20C, pA=self.pA, pB=self.pB, dw_AB=self.dw_AB, dw_BC=self.dw_BC, dwH_AB=self.dwH_AB, dwH_BC=self.dwH_BC, kex_AB=self.kex_AB, kex_BC=self.kex_BC, kex_AC=self.kex_AC, inv_tcpmg=self.inv_tcpmg, tcp=self.tcp, back_calc=self.back_calc, num_points=self.num_points, power=self.power)
1135