Package test_suite :: Package unit_tests :: Package _lib :: Package _dispersion :: Module test_ns_mmq_3site
[hide private]
[frames] | no frames]

Source Code for Module test_suite.unit_tests._lib._dispersion.test_ns_mmq_3site

   1  ############################################################################### 
   2  #                                                                             # 
   3  # Copyright (C) 2019 Edward d'Auvergne                                        # 
   4  #                                                                             # 
   5  # This file is part of the program relax (http://www.nmr-relax.com).          # 
   6  #                                                                             # 
   7  # This program is free software: you can redistribute it and/or modify        # 
   8  # it under the terms of the GNU General Public License as published by        # 
   9  # the Free Software Foundation, either version 3 of the License, or           # 
  10  # (at your option) any later version.                                         # 
  11  #                                                                             # 
  12  # This program is distributed in the hope that it will be useful,             # 
  13  # but WITHOUT ANY WARRANTY; without even the implied warranty of              # 
  14  # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the               # 
  15  # GNU General Public License for more details.                                # 
  16  #                                                                             # 
  17  # You should have received a copy of the GNU General Public License           # 
  18  # along with this program.  If not, see <http://www.gnu.org/licenses/>.       # 
  19  #                                                                             # 
  20  ############################################################################### 
  21   
  22  # Python module imports. 
  23  from numpy import array, int16 
  24  from unittest import TestCase 
  25   
  26  # relax module imports. 
  27  from lib.dispersion.ns_mmq_3site import r2eff_ns_mmq_3site_mq, r2eff_ns_mmq_3site_sq_dq_zq 
  28   
  29   
30 -class Test_ns_mmq_3site(TestCase):
31 """Unit tests for the lib.dispersion.ns_mmq_3site relax module.""" 32
33 - def setUp(self):
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
1115 - def test_ns_mmq_3site_mq(self):
1116 """Test the r2eff_ns_mmq_3site_mq() function for no exchange when pA = 1.0 and pB = 0.0.""" 1117 1118 # Parameter reset. 1119 self.pA = 1.0 1120 self.pB = 0.0 1121 1122 # Calculate and check the R2eff values. 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
1126 - def test_ns_mmq_3site_sq_dq_zq(self):
1127 """Test the r2eff_ns_mmq_3site_sq_dq_zq() function for no exchange when pA = 1.0 and pB = 0.0.""" 1128 1129 # Parameter reset. 1130 self.pA = 1.0 1131 self.pB = 0.0 1132 1133 # Calculate and check the R2eff values. 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