From dffb9858e80f052bc4152b5d1602df9c9c5e06ea Mon Sep 17 00:00:00 2001 From: Adam Nelson Date: Wed, 8 Jun 2016 19:38:44 -0400 Subject: [PATCH] Fixed per @paulromano comments --- .../pythonapi/examples/mgxs-part-iv.ipynb | 172 +-- docs/source/usersguide/mgxs_library.rst | 4 +- openmc/mgxs/library.py | 42 +- src/mgxs_header.F90 | 14 +- src/scattdata_header.F90 | 20 +- tests/test_mg_basic/results_true.dat | 2 +- tests/test_mg_max_order/results_true.dat | 2 +- tests/test_mg_nuclide/results_true.dat | 2 +- tests/test_mg_tallies/results_true.dat | 1306 ++++++++--------- .../results_true.dat | 2 +- 10 files changed, 771 insertions(+), 795 deletions(-) diff --git a/docs/source/pythonapi/examples/mgxs-part-iv.ipynb b/docs/source/pythonapi/examples/mgxs-part-iv.ipynb index 5459c3dfe..e5c80c192 100644 --- a/docs/source/pythonapi/examples/mgxs-part-iv.ipynb +++ b/docs/source/pythonapi/examples/mgxs-part-iv.ipynb @@ -432,7 +432,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD6AgMAAAD1grKuAAAABGdBTUEAALGPC/xhBQAAACBjSFJN\nAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEX////pgJFyEhJNv8RV\nUZDeAAAAAWJLR0QAiAUdSAAAAAd0SU1FB+AGBxMuGM4M8GgAAAWFSURBVGje7Zs7cttADIZ9CSvX\ncrP0iCxUqbBc8Ag6xR6BhV2EvYvwFD4CCx1ABT1jMdgndpegRQnOrCbjpPlGESISC4A/gd27e8H5\n83CX3b4+iKJrRHkS4vkghMPBonRYWGwtfgD2YN+dRDUOoh6lACw0Noi9w2fESuEoAR/uVuMolX03\n9oXGT7F3eFL2iEfhUX1f4cPdL/ishs+68ai+udE4xPhexbjX2FfjGNoPj/DPNX4Tsd+EODr8FvsV\ndf1Hd9P2VvCi4+s/aXvrf+upAD+1/9GV1mkOH5X9vV6THtfvACslcaUCbESL61drBPtdI8SrFMWr\nELsXCkuFDYW75gbiP7d9Cf7bAYI/aCwUShrBvh30+lWQkzVgZ/HD4OixNCgcQpJ3BxU/Ln91elKo\nM5VEE38QtJ+Yv6cQ9xjKNYayyl8TypP8DfJnQ2H/b/N3ye9P83cT33SQv/sQh9gV7zZ/0dNj5HQa\nC5vVzv9+/WFN2w8KVaZ2BwL1+pv4g0x1QRfjq0dB4Q3kT277oP6VNL6gKxNU9a8zK+WLbi/Wwpdi\nhbboKqyxFOulHMj6v4W/AXbmUeAxrv9J/CqEBXaRKsXaodD4nsYvkT/G6H1D4SR/iPy1Roj9JsQ5\ne18/7EUHv1+Fvx/Xj5V9Ugb5K8TW4TZEEdcvoz/up0VTe9qsVIppKVX6a7D6y9ZvwEKjrtQxPtv6\nfXII9vCxKOGaIeAIfEF8IvAG8ie3vRK9rRQl+PPpSctbhfpTUCpviH+kxsZgpT91+snoX1l49KK3\niUQvICRy5aUw6l8leoVwoo3Uv1rKreF/UFLY6d9QP4L9Wf2r7EP9GOSfcsjZ56f60kz+XmVPXv+R\nuP49ff0T/53Rv6n/7m2lvXT9Wqd/VUz8hvh5M/ED6ILmt4mfHYZSaePnTWpsf/SvqV9O6dLYYClL\nEetnoH/LBLFoBvrX189uTv8++kot5vTvQD4/9jP690g9P/4z/bvo/XVG/xYoZZx+8fr3MxAtsf7t\nUOkG2JqsTtCIpgCt/qX1226KqZS7gfzJbe+c9jLrtIZ8lXD+s4umlW6AKIVrlML2/cXjgPFjlJqI\nRC+Fj0bVJe+vSh56pSdR6YkQ1ygF10Wqf0FeLta/iKn9Mv1L24ti2e+7W4n1b3T/W+L+t9H9T/Sv\nVboUmqJJon1/hZq8LnzRDlDrX1u0xRT1+6vEpomMmyYkqi95vIH8yW1PN+122KkLcNLKi/WTF01z\n/cNASrWE/l3ev6T17zX909z9X27/euK/Rf3zWP+Waf9eEv37KkWJ+rfDl6ZglNDa+cEBhwYDvkoN\nP/rX69814NaI3imq0l7OYDy/qSdDGwr7r+Y3VbzoKZr6XX2lfxfOb87qXzr+b1j/Xlp/nP6dn98M\ncdH7cn7zjPObKsYWS3Eb9w8n85smHtqQuPuZ30T2dlIT6F9xFl+n8xslegL9a4c2KRr9W4rp/GYq\numiM9Nec/j2v/yj9u1h//hv9e93vc++f63/u+rPjL3f+5Lbn1j9m/eXWf+7zh/v8+2b9e/Hzn6s/\nuPqHrb8g71n6L3f+5Lbnvn8w33+4718/+5d47//c/gO7/5E7/nPbc/tv3P4fs//I7X9y+6/fqH+v\n6j9z+9/c/ju3/8+eP+TOn9z23PkXc/7Gnf9x5483q38Xzn+582fu/Js9fy8kb/6fO39y23P3n3S8\n/S/c/Tfc/T83uX/pgv1XE/9duP+Lu/+Mvf8td/znti8kb/8ld/9nx9t/Sjw/Ltr/yt1/+337f6/b\nf0zoB3nJ/ucVc/81d/83e/957vzJbc89/8A8f8E9/5HE78XnT/4H/cs5f8Q9/8Q9f8U+/5U7f3Lb\nc88fdrzzjyvm+cuf/Uu887/c88fs88954/8vO4SjPC+2QRIAAAAldEVYdGRhdGU6Y3JlYXRlADIw\nMTYtMDYtMDdUMTk6NDY6MjQtMDQ6MDAgoUlWAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDE2LTA2LTA3\nVDE5OjQ2OjI0LTA0OjAwUfzx6gAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAD6AgMAAAD1grKuAAAABGdBTUEAALGPC/xhBQAAACBjSFJN\nAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEX////pgJFyEhJNv8RV\nUZDeAAAAAWJLR0QAiAUdSAAAAAd0SU1FB+AGCBMhJtCW4VoAAAWFSURBVGje7Zs7cttADIZ9CSvX\ncrP0iCxUqbBc8Ag6xR6BhV2EvYvwFD4CCx1ABT1jMdgndpegRQnOrCbjpPlGESISC4A/gd27e8H5\n83CX3b4+iKJrRHkS4vkghMPBonRYWGwtfgD2YN+dRDUOoh6lACw0Noi9w2fESuEoAR/uVuMolX03\n9oXGT7F3eFL2iEfhUX1f4cPdL/ishs+68ai+udE4xPhexbjX2FfjGNoPj/DPNX4Tsd+EODr8FvsV\ndf1Hd9P2VvCi4+s/aXvrf+upAD+1/9GV1mkOH5X9vV6THtfvACslcaUCbESL61drBPtdI8SrFMWr\nELsXCkuFDYW75gbiP7d9Cf7bAYI/aCwUShrBvh30+lWQkzVgZ/HD4OixNCgcQpJ3BxU/Ln91elKo\nM5VEE38QtJ+Yv6cQ9xjKNYayyl8TypP8DfJnQ2H/b/N3ye9P83cT33SQv/sQh9gV7zZ/0dNj5HQa\nC5vVzv9+/WFN2w8KVaZ2BwL1+pv4g0x1QRfjq0dB4Q3kT277oP6VNL6gKxNU9a8zK+WLbi/Wwpdi\nhbboKqyxFOulHMj6v4W/AXbmUeAxrv9J/CqEBXaRKsXaodD4nsYvkT/G6H1D4SR/iPy1Roj9JsQ5\ne18/7EUHv1+Fvx/Xj5V9Ugb5K8TW4TZEEdcvoz/up0VTe9qsVIppKVX6a7D6y9ZvwEKjrtQxPtv6\nfXII9vCxKOGaIeAIfEF8IvAG8ie3vRK9rRQl+PPpSctbhfpTUCpviH+kxsZgpT91+snoX1l49KK3\niUQvICRy5aUw6l8leoVwoo3Uv1rKreF/UFLY6d9QP4L9Wf2r7EP9GOSfcsjZ56f60kz+XmVPXv+R\nuP49ff0T/53Rv6n/7m2lvXT9Wqd/VUz8hvh5M/ED6ILmt4mfHYZSaePnTWpsf/SvqV9O6dLYYClL\nEetnoH/LBLFoBvrX189uTv8++kot5vTvQD4/9jP690g9P/4z/bvo/XVG/xYoZZx+8fr3MxAtsf7t\nUOkG2JqsTtCIpgCt/qX1226KqZS7gfzJbe+c9jLrtIZ8lXD+s4umlW6AKIVrlML2/cXjgPFjlJqI\nRC+Fj0bVJe+vSh56pSdR6YkQ1ygF10Wqf0FeLta/iKn9Mv1L24ti2e+7W4n1b3T/W+L+t9H9T/Sv\nVboUmqJJon1/hZq8LnzRDlDrX1u0xRT1+6vEpomMmyYkqi95vIH8yW1PN+122KkLcNLKi/WTF01z\n/cNASrWE/l3ev6T17zX909z9X27/euK/Rf3zWP+Waf9eEv37KkWJ+rfDl6ZglNDa+cEBhwYDvkoN\nP/rX69814NaI3imq0l7OYDy/qSdDGwr7r+Y3VbzoKZr6XX2lfxfOb87qXzr+b1j/Xlp/nP6dn98M\ncdH7cn7zjPObKsYWS3Eb9w8n85smHtqQuPuZ30T2dlIT6F9xFl+n8xslegL9a4c2KRr9W4rp/GYq\numiM9Nec/j2v/yj9u1h//hv9e93vc++f63/u+rPjL3f+5Lbn1j9m/eXWf+7zh/v8+2b9e/Hzn6s/\nuPqHrb8g71n6L3f+5Lbnvn8w33+4718/+5d47//c/gO7/5E7/nPbc/tv3P4fs//I7X9y+6/fqH+v\n6j9z+9/c/ju3/8+eP+TOn9z23PkXc/7Gnf9x5483q38Xzn+582fu/Js9fy8kb/6fO39y23P3n3S8\n/S/c/Tfc/T83uX/pgv1XE/9duP+Lu/+Mvf8td/znti8kb/8ld/9nx9t/Sjw/Ltr/yt1/+337f6/b\nf0zoB3nJ/ucVc/81d/83e/957vzJbc89/8A8f8E9/5HE78XnT/4H/cs5f8Q9/8Q9f8U+/5U7f3Lb\nc88fdrzzjyvm+cuf/Uu887/c88fs88954/8vO4SjPC+2QRIAAAAldEVYdGRhdGU6Y3JlYXRlADIw\nMTYtMDYtMDhUMTk6MzM6MzgtMDQ6MDBqU5pBAAAAJXRFWHRkYXRlOm1vZGlmeQAyMDE2LTA2LTA4\nVDE5OjMzOjM4LTA0OjAwGw4i/QAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] @@ -722,8 +722,8 @@ " Copyright: 2011-2016 Massachusetts Institute of Technology\n", " License: http://openmc.readthedocs.io/en/latest/license.html\n", " Version: 0.7.1\n", - " Git SHA1: b41ec2178cc78c6c5968b47ebf09d926cdc44d00\n", - " Date/Time: 2016-06-07 19:46:24\n", + " Git SHA1: 826d5a43d85eaec1b6c7b4ce22e1a8f5e9336a4f\n", + " Date/Time: 2016-06-08 19:33:38\n", " OpenMP Threads: 4\n", "\n", " ===========================================================================\n", @@ -810,20 +810,20 @@ "\n", " =======================> TIMING STATISTICS <=======================\n", "\n", - " Total time for initialization = 1.4790E+00 seconds\n", - " Reading cross sections = 1.1810E+00 seconds\n", - " Total time in simulation = 1.7439E+01 seconds\n", - " Time in transport only = 1.7334E+01 seconds\n", - " Time in inactive batches = 2.1960E+00 seconds\n", - " Time in active batches = 1.5243E+01 seconds\n", - " Time synchronizing fission bank = 6.0000E-03 seconds\n", - " Sampling source sites = 5.0000E-03 seconds\n", + " Total time for initialization = 1.4220E+00 seconds\n", + " Reading cross sections = 1.1320E+00 seconds\n", + " Total time in simulation = 1.6571E+01 seconds\n", + " Time in transport only = 1.6501E+01 seconds\n", + " Time in inactive batches = 2.1010E+00 seconds\n", + " Time in active batches = 1.4470E+01 seconds\n", + " Time synchronizing fission bank = 5.0000E-03 seconds\n", + " Sampling source sites = 4.0000E-03 seconds\n", " SEND/RECV source sites = 1.0000E-03 seconds\n", " Time accumulating tallies = 0.0000E+00 seconds\n", " Total time for finalization = 0.0000E+00 seconds\n", - " Total time elapsed = 1.8926E+01 seconds\n", - " Calculation Rate (inactive) = 22768.7 neutrons/second\n", - " Calculation Rate (active) = 13120.8 neutrons/second\n", + " Total time elapsed = 1.8002E+01 seconds\n", + " Calculation Rate (inactive) = 23798.2 neutrons/second\n", + " Calculation Rate (active) = 13821.7 neutrons/second\n", "\n", " ============================> RESULTS <============================\n", "\n", @@ -1098,8 +1098,8 @@ " Copyright: 2011-2016 Massachusetts Institute of Technology\n", " License: http://openmc.readthedocs.io/en/latest/license.html\n", " Version: 0.7.1\n", - " Git SHA1: b41ec2178cc78c6c5968b47ebf09d926cdc44d00\n", - " Date/Time: 2016-06-07 19:46:43\n", + " Git SHA1: 826d5a43d85eaec1b6c7b4ce22e1a8f5e9336a4f\n", + " Date/Time: 2016-06-08 19:33:56\n", " OpenMP Threads: 4\n", "\n", " ===========================================================================\n", @@ -1124,56 +1124,56 @@ "\n", " Bat./Gen. k Average k \n", " ========= ======== ==================== \n", - " 1/1 1.02438 \n", - " 2/1 1.04495 \n", - " 3/1 1.03438 \n", - " 4/1 0.99945 \n", - " 5/1 1.00818 \n", - " 6/1 1.00174 \n", - " 7/1 1.03430 \n", - " 8/1 1.02233 \n", - " 9/1 1.01200 \n", - " 10/1 1.02727 \n", - " 11/1 1.04051 \n", - " 12/1 1.02665 1.03358 +/- 0.00693\n", - " 13/1 1.01925 1.02880 +/- 0.00623\n", - " 14/1 1.03520 1.03040 +/- 0.00469\n", - " 15/1 1.03185 1.03069 +/- 0.00364\n", - " 16/1 1.03729 1.03179 +/- 0.00317\n", - " 17/1 1.00718 1.02827 +/- 0.00442\n", - " 18/1 1.00733 1.02566 +/- 0.00464\n", - " 19/1 1.02123 1.02516 +/- 0.00412\n", - " 20/1 1.01673 1.02432 +/- 0.00378\n", - " 21/1 1.03445 1.02524 +/- 0.00354\n", - " 22/1 1.03063 1.02569 +/- 0.00326\n", - " 23/1 1.00806 1.02433 +/- 0.00329\n", - " 24/1 1.03707 1.02524 +/- 0.00318\n", - " 25/1 1.01907 1.02483 +/- 0.00299\n", - " 26/1 1.01970 1.02451 +/- 0.00282\n", - " 27/1 1.00208 1.02319 +/- 0.00296\n", - " 28/1 1.02836 1.02348 +/- 0.00280\n", - " 29/1 1.02679 1.02365 +/- 0.00266\n", - " 30/1 1.03997 1.02447 +/- 0.00265\n", - " 31/1 0.99825 1.02322 +/- 0.00281\n", - " 32/1 1.01865 1.02301 +/- 0.00269\n", - " 33/1 1.00846 1.02238 +/- 0.00265\n", - " 34/1 1.01864 1.02222 +/- 0.00254\n", - " 35/1 1.02137 1.02219 +/- 0.00243\n", - " 36/1 1.02781 1.02241 +/- 0.00235\n", - " 37/1 1.02395 1.02246 +/- 0.00226\n", - " 38/1 1.00883 1.02198 +/- 0.00223\n", - " 39/1 0.99776 1.02114 +/- 0.00231\n", - " 40/1 1.01546 1.02095 +/- 0.00224\n", - " 41/1 1.04425 1.02170 +/- 0.00229\n", - " 42/1 1.02481 1.02180 +/- 0.00222\n", - " 43/1 1.00925 1.02142 +/- 0.00219\n", - " 44/1 1.04389 1.02208 +/- 0.00222\n", - " 45/1 1.04451 1.02272 +/- 0.00225\n", - " 46/1 1.02519 1.02279 +/- 0.00219\n", - " 47/1 1.03656 1.02316 +/- 0.00216\n", - " 48/1 1.04310 1.02369 +/- 0.00217\n", - " 49/1 1.03258 1.02392 +/- 0.00212\n", - " 50/1 1.01093 1.02359 +/- 0.00210\n", + " 1/1 0.99367 \n", + " 2/1 1.03173 \n", + " 3/1 1.01999 \n", + " 4/1 1.01421 \n", + " 5/1 1.03980 \n", + " 6/1 1.04540 \n", + " 7/1 1.04199 \n", + " 8/1 1.02680 \n", + " 9/1 1.01267 \n", + " 10/1 1.03420 \n", + " 11/1 1.05773 \n", + " 12/1 1.03475 1.04624 +/- 0.01149\n", + " 13/1 1.03632 1.04293 +/- 0.00741\n", + " 14/1 0.99297 1.03044 +/- 0.01355\n", + " 15/1 1.02413 1.02918 +/- 0.01057\n", + " 16/1 1.02359 1.02825 +/- 0.00868\n", + " 17/1 0.99913 1.02409 +/- 0.00843\n", + " 18/1 1.01493 1.02294 +/- 0.00739\n", + " 19/1 1.03010 1.02374 +/- 0.00657\n", + " 20/1 1.04890 1.02626 +/- 0.00639\n", + " 21/1 1.01267 1.02502 +/- 0.00591\n", + " 22/1 1.02637 1.02513 +/- 0.00540\n", + " 23/1 1.01374 1.02426 +/- 0.00504\n", + " 24/1 1.06661 1.02728 +/- 0.00556\n", + " 25/1 1.03212 1.02760 +/- 0.00519\n", + " 26/1 1.05433 1.02927 +/- 0.00513\n", + " 27/1 0.99891 1.02749 +/- 0.00514\n", + " 28/1 1.00616 1.02630 +/- 0.00499\n", + " 29/1 1.04583 1.02733 +/- 0.00483\n", + " 30/1 1.01512 1.02672 +/- 0.00462\n", + " 31/1 0.98104 1.02455 +/- 0.00491\n", + " 32/1 1.04202 1.02534 +/- 0.00474\n", + " 33/1 1.00779 1.02458 +/- 0.00460\n", + " 34/1 1.02450 1.02457 +/- 0.00440\n", + " 35/1 0.98882 1.02314 +/- 0.00446\n", + " 36/1 1.01541 1.02285 +/- 0.00429\n", + " 37/1 1.02050 1.02276 +/- 0.00413\n", + " 38/1 1.03573 1.02322 +/- 0.00401\n", + " 39/1 1.03649 1.02368 +/- 0.00389\n", + " 40/1 1.01434 1.02337 +/- 0.00378\n", + " 41/1 1.02345 1.02337 +/- 0.00365\n", + " 42/1 1.01900 1.02323 +/- 0.00354\n", + " 43/1 1.01450 1.02297 +/- 0.00344\n", + " 44/1 1.03127 1.02321 +/- 0.00335\n", + " 45/1 1.01598 1.02301 +/- 0.00326\n", + " 46/1 1.00851 1.02260 +/- 0.00319\n", + " 47/1 1.03406 1.02291 +/- 0.00312\n", + " 48/1 1.02373 1.02294 +/- 0.00303\n", + " 49/1 1.04066 1.02339 +/- 0.00299\n", + " 50/1 1.02011 1.02331 +/- 0.00292\n", " Creating state point statepoint.50.h5...\n", "\n", " ===========================================================================\n", @@ -1183,27 +1183,27 @@ "\n", " =======================> TIMING STATISTICS <=======================\n", "\n", - " Total time for initialization = 3.7000E-02 seconds\n", + " Total time for initialization = 3.1000E-02 seconds\n", " Reading cross sections = 5.0000E-03 seconds\n", - " Total time in simulation = 1.1738E+01 seconds\n", - " Time in transport only = 1.1701E+01 seconds\n", - " Time in inactive batches = 1.3240E+00 seconds\n", - " Time in active batches = 1.0414E+01 seconds\n", - " Time synchronizing fission bank = 8.0000E-03 seconds\n", - " Sampling source sites = 5.0000E-03 seconds\n", - " SEND/RECV source sites = 3.0000E-03 seconds\n", + " Total time in simulation = 1.1867E+01 seconds\n", + " Time in transport only = 1.1830E+01 seconds\n", + " Time in inactive batches = 1.2670E+00 seconds\n", + " Time in active batches = 1.0600E+01 seconds\n", + " Time synchronizing fission bank = 7.0000E-03 seconds\n", + " Sampling source sites = 7.0000E-03 seconds\n", + " SEND/RECV source sites = 0.0000E+00 seconds\n", " Time accumulating tallies = 0.0000E+00 seconds\n", " Total time for finalization = 0.0000E+00 seconds\n", - " Total time elapsed = 1.1785E+01 seconds\n", - " Calculation Rate (inactive) = 37764.4 neutrons/second\n", - " Calculation Rate (active) = 19204.9 neutrons/second\n", + " Total time elapsed = 1.1907E+01 seconds\n", + " Calculation Rate (inactive) = 39463.3 neutrons/second\n", + " Calculation Rate (active) = 18867.9 neutrons/second\n", "\n", " ============================> RESULTS <============================\n", "\n", - " k-effective (Collision) = 1.02444 +/- 0.00180\n", - " k-effective (Track-length) = 1.02359 +/- 0.00210\n", - " k-effective (Absorption) = 1.02530 +/- 0.00175\n", - " Combined k-effective = 1.02471 +/- 0.00130\n", + " k-effective (Collision) = 1.02638 +/- 0.00260\n", + " k-effective (Track-length) = 1.02331 +/- 0.00292\n", + " k-effective (Absorption) = 1.02579 +/- 0.00132\n", + " Combined k-effective = 1.02558 +/- 0.00136\n", " Leakage Fraction = 0.00000 +/- 0.00000\n", "\n" ] @@ -1285,8 +1285,8 @@ "output_type": "stream", "text": [ "Continuous-Energy keff = 1.024295\n", - "Multi-Group keff = 1.024707\n", - "bias [pcm]: -41.2\n" + "Multi-Group keff = 1.025577\n", + "bias [pcm]: -128.2\n" ] } ], @@ -1383,7 +1383,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 40, @@ -1392,9 +1392,9 @@ }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAADDCAYAAACS2+oqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXmYFNXV/z8HUFQEBBQUFVxwQYwoKoJJzLhh8OcWgwpq\n3LefGolJ3KLR0RjXvO4ajeIaFeOu0Sj64piYuCMYjahEQBEBBRFccGHO+0fVYDNM96lhpumZ8vt5\nnnmmu+637j1dderUrVt16pq7I4QQovXTptIGCCGEaB4U0IUQIicooAshRE5QQBdCiJyggC6EEDlB\nAV0IIXJCiwvoZvaamW1XaTu+y5jZo2b2syas/0czO705bfouYma1ZrZeifLcHivywaXE3cM/YH/g\nRWA+8D7wCPD9LOsG9d4EnNPUeir5l/6GL4F56d984JVK25XB7rOArwpsngf8utJ2NcLmOcAzwKBG\nrP8UcNgysHMKsADoWm/5eKAW6JWxnoXAegV+1qhjBVgOOBOYmO7j99Jjd+dK78sG9qd8sBn+wh66\nmf0SuAQ4F+gO9AKuAfaI1v0OcaG7d0r/Orr7Fs3dgJm1be46gdEFNndy9z+UoY3mZrS7dwJWBWqA\nuytrToM4MBkYUbfAzDYFVkjLsmJNtONeYHfgQKALsC5wObBrg42Vx8ci5IPNSXA26URy5ty7hGZ5\n4DKSnvs04FJgubTsRyS9gl8CM1PNIWnZkSRnugUkZ7sH0+WTgR0KzoZ3Abekmn8DAwrariXtwaTf\nF+vFpG28DXwEPACskS7vna7bpqEzJ7A+yY6aC8wC7izx+4v2nAraOQiYmtb1m4JyA04FJgEfAqOB\nVeqte1i6bk26/CCSHuCHwBl12wvoAXwGdCmof8u0zbZFehq3Rr2IUtsi3dcz07LxwCaN2Q8F+/Bo\n4C1gNnBV0Du6teB7X5JebLf0+yrAw6mds9PPPdOyc4FvgM9TX7oiXb4xMCbVvwHsU1D/rsDrqf49\n4JcZe2GTgd8ALxQsuxg4LbW3V0O9NeBg4B/1/ZsMx0oDNuyU+sMaGWw9GZgAfEEyDNs3te1jkmNu\n94Z8o4TNPwf+m+6Hi7LuT/lg030w6qEPBtqnG6AYZwADgc2A/unnMwrKVwc6Aj2BI4Crzayzu18P\n3E6ywzu5+55F6t8duAPonG6cqwvKivZ2zGwH4DxgGLAG8C5JwAzXBX4HPO7uqwBrAVeW0Gbh+8AG\nJAfZmWa2Ubp8JMmVzg9Jts/HJFc/hWxHssN3MbO+JL9/BMlv6pyuh7vPJDkI9i1Y9wAS51/YBNsb\n3BZmNgT4AdAnLduPxCEXI8N+APh/JCefzYF907pLYmbLkwST2STbDZJgdCOwNsmV5Oek/uLuZwD/\nAI5P/e0EM1uJ5ED6M0lvawRwTbqdAW4AjvSkN7YpMDayq4DngI5mtpGZtSHZL38m7nUv4ZeNOFYK\n2RF43t0/yKAdDgwlCUZtgIeAx4DVgBOA281sg0bYvBcwIP3b08wOy2BDKeSDGX0wCujdgI/cvbaE\nZn/gbHef7e6zgbOBwpsZXwG/c/eF7v434FNgowbqKcYz7v64J6er20hOHHWUOjj2B0a5+wR3/5qk\ndzTYzHplaPNroLeZrenuX7n7vwL9SWY2x8w+Tv/fVFDmQHVaz6skPaH+adlRwOnu/kFq4znAsDQA\n1K17lrt/4e5fkjjkQ+7+rLt/QzI+WsitpNs+rWMEyTYrxn717F69Edvia5IT9SZmZu7+ZnpSqU+W\n/XC+u8939/dITkqbRzaTHCiHA8Pq/NPd57j7/e7+pbt/BpxPckIsxm7AZHe/1RPGkwxTDEvLvwL6\nmVlHd/8kLW8Mt5Ec8DuTjGNPb+T6TWFVYEbdFzPrku7nuWb2RT3t5e4+PfWxQUAHd7/Q3b9x96eA\nv1IwfJSBC9LtNY3k6r3UuvLBZvTBKKDPBlYtCDAN0ZPkjFfH1HTZojrqnRA+B1YO2i1kRsHnz4EV\nAnsK7Zpa9yXduLOBNTOsexLJtnnBzP5tZocCmNlpZjbfzOaZWWFP+mJ37+ruXdL/h9arr9DJCn9/\nb+D+1JHnAP8hcdIeBfpp9X7TewW/6QsW75E8CPQ1s3WAIcBcd3+pxO+8q57dMxrQNLgt0gP9KpLe\nxwwzu9bMGtqvWfZDse1T1GaS+zmvAVvVFZjZimZ2nZlNMbO5wNPAKmZW7MTfGxhUt/3N7GOSg79u\n+/+UpOc21cyeMrNBJexqiD+n9R1CcrItG6lf1vnmWiTbeI26cnf/2N27kPRCl6+3elEfS5lKtuOm\nofrqx4P6yAeb0QejwPgsybjdXiU076dGFRqYtSfSmBtEDfE5sFLB98Kz+/RCu8ysA8kVxzSSsUWK\nrevus9z9KHdfEziG5BJoPXc/37+9eXNsE22H5EQ4NHXkOqfuUO8yuXAbfUByyVn3m1ZMf1Od3V8C\nfyG5CXYgpXvnmSi2LdKyq9x9K6AfyVXXSQ1UUWo/NMWuOak91WZW5/y/Ihna2jq9BK/rGdUdTPX9\n7T2SexOF27+Tux+ftvGyu+9FMvTwIMm2bYyN75KMUQ8F7mtA8hnF/XeJ6oK2Ohb45jTgf4Gtzayh\nYFo/uBTWPZ1kuKCQXiTHeVabC9fvRROvTOSD2X2wZEB393kkNwGuNrM907NPOzMbamYXpLLRwBlm\ntqqZrQr8luyBZCbJTZ/GUOiMrwD7m1kbM/sxyU3YOu4ADjWzzcysPckY2nPu/p67f0TioAem6x5G\ncuMlacBsmJnVnb3nktw0Wdpx6FLDQtcB59Vd+pnZamZW+PRQ/XXvAXY3s0FmthzJ8FZ9biPpEe5O\n0kNsEsW2hZltZWYDzawdyc20BTS8jYruh6ba5u5vkoz1npIu6pjaMs/MugLV9Vap729/BTY0swNT\nv14u/V0bp5/3N7NOntyDmE9yQ6uxHEZy47L+MAckN/H2To+rPiSX78Vo1LHi7k+QDB08kO6n5dJ9\nNZjSJ4fngc/M7OR0m1SRDAvc2QibTzKzVcxsbZL7RPXHqxuFfDC7D4ZDF+5+KclTKmeQ3Ll9FziW\nb2+Ungu8BNSND78E/L5UlQWfR5GMD80xs/saKI/W/wXJTcWPScbp7i+weyzJyeU+kuC9LsnNnzqO\nJLm7/xHJnep/FpRtDTxvZvPS33mCu0+lOCenl7rz0sveWUXsrf/9cpKz7hgz+wT4F8lN5QbXdff/\nkDxBcBdJr+MTkn3yZYHmXyQOPy7tIS4Nhe0W2xadgOtJnsWdTLIdl3jkLMN+KLV9svAH4Mi0M3EZ\nSe/xI5Jt+Wg97eXAPmY228wuc/dPSYamhpNsz+nABXw7JPEzYHJ66XwUyU3mLCz6De4+2d3HNVRG\n8oTG1yTDijex5Am4qcfK3iQB488kx8g7JMfJLkXaIB1j3oPk6YqPSIY0fubub2e0GRKffhkYR/Ig\nw42BnQ0hH0xolA+ae1NHPUSlSC8d55Lc5Z9asPx/gdvdfWkOJCGWGjOrJfHHdypty3eRFpf6L0pj\nZrull7sdgP8BXq0XzLcGtiDpxQshvkMooLc+9iS5LJtGMu6/6NLRzG4meaZ1ZHonX4hljS75K4iG\nXIQQIieohy6EEDmhXTkqTR8hvIzkhDHK3S9sQKNLA1FW3L2pL7daAvm2aAkU8+1mH3KxJIvzLZJ3\nSUwnee3ucHefWE/nnFLQ9jPV8IPqxSu7JUODl2bQZEla/ncGzUv1ttUD1bBX9eLLDo1zFS6tPTfU\nnPjWtSXLr9jwyLCOE8Zev+TCW6rh4OpFX9fa8e0lNfUYv+hNBcW5gSNCzVe0DzXX+dFLLJtXfQWd\nqk9Y9P39F0q9ViRlkDV7QG+Ub/NJwZLzSbLNU57sGLZ16477hJpZdA81j/rQUPPktCVfnFp9CVT/\n8tvvm631fFjPqy9sE2qGD7wp1Nw2r36i9ZI83HnJ33VX9dvsV/2tb1Qv8Qj4kgz2Z0PN9TeODDV0\nyxBHf/J5vQW/B+q/sv3iklUMGbI+Y8YcVNS3yzHkMhB4292nps+0jia5kSdEa0e+LVo05Qjoa7L4\nuyCm0bj3QAjRUpFvixZNOcbQG7oUaPh65Jnqbz+3X6UMppSZjasqbUHj6V9VaQsaTfuq+FKel2tg\nXE25Tcnu25xf8LlTOWwpK1WDK21B4+lX1bXSJjSSH2bUTUn/YNKkCSWV5Qjo00heyFPHWhR7OU/9\nMfPWRmsM6JtXVdqCRpMpoG9ZlfzVMaqh19w0mey+XThm3gppjQF906pusahFkXU62HXSP+jTZ33e\neaf49BTlGHJ5EehjZr0teQH8cJIX5gvR2pFvixZNs/fQ3X2hmR1PkrFY92jXG83djhDLGvm2aOmU\n5Tl0d3+Mxs1KJESrQL4tWjIVS/03M2edoO2dYtvsd/Wf7VwSv2qlULNweDzh+Qqrzw01HVeZF2q2\naRc/09s36PhdPit+NnZk98tDTZU9FWreZsNQM2GxmQEb5lnigdm5dAk1a2SYJvPVNoPLkliUBTPz\nK2qLv9p8PvFz6L0sfvPxYKKZEcn0ZpX1PmhokqDFeaznj0JNR499fwErhpq9P29oLpDFGdBhXKjZ\nquRkXQkXvXxWqOk3IK5n4uYDQg2/zbAzgkfwh+wIYx5qs0yfQxdCCFEBFNCFECInKKALIUROUEAX\nQoicoIAuhBA5QQFdCCFyggK6EELkBAV0IYTICWXJFM1M6Xe5w5dxFUd3Lz0RBECP834daqxjnIPy\n5clxQkhVhmSPnl7kfU4F/GGJF98vzkbdDwrr2JTXQs1gjxM0nmOLUHPiJfF+uO9X8WQLF3FyqNnN\n/hpqXg0V5aXByUVShu8YT/JwlR8Xat6fFU/04Q9lyK06ojaUHJvhDQf/vXrTUGPHLQw1c2fFSX4j\n1rkx1FzMGaGGu88MJdtuGR/TEw/KkFj0foZ9sXFQvnbpYvXQhRAiJyigCyFETlBAF0KInKCALoQQ\nOUEBXQghcoICuhBC5AQFdCGEyAmVneDixaDtlWPbvlktfma1Tdf42Vf+FJ/bPjhqlVBzol0aag7x\nm0PNLmP/XrL8qB3jySsu9PiZ7q4D44f9vXMowZ6It7GPirfxTkc8HGrmEE8GPMG2regEF0Nq7y9a\nfu+XPw3rWOnyDJO7nBxv8wkZJle63Q8MNRdNjSeCsHgeGfzaDLtkRvzbd7omnsr1yYf2jNsaH0s4\nM97OazA51Hw4s3uoObDHbSXLN2VtTm6zuya4EEKIvKOALoQQOUEBXQghcoICuhBC5AQFdCGEyAkK\n6EIIkRMU0IUQIicooAshRE6o6AQXL27Zr2T51u+/FNbRZk6chLB1l3+Empfi/AvW+NvcUDO6/6Fx\nRWNjCaeVLl5u2tdhFftyV6jZ4aU4+eL0p0IJtefFCV4v/6b0/gb4E0dl0BwdaiaEivIy14onoV3Z\n/ufh+qeuEieOffVJvM3Hdd4v1Fx0Y5w0xMRY8tm5cdLQSq/Hx+sZo+OJKZ7sv0ds0E9iyTvVq4ea\nxz0+pt//+NZQ025UnKB06+rHlCwfsmbp9dVDF0KInKCALoQQOUEBXQghcoICuhBC5AQFdCGEyAkK\n6EIIkRMU0IUQIicooAshRE4oS2KRmU0BPgFqga/dfWBDuuE2umQ943punqGx2lAykEtCzRdbxkkR\nK6ycYeajl+Nz5J0HxMk8Iw4sPuMNwKd+fVjHE7P2CjVWG28/fyD+Tf88bYtQ8wPiRDGejtva5kfP\nx/WUiay+fRf7Fq1jPPG24qh4v7zK90LNIY/GyWUcHrd1hcfJXD+/+Ia4rbvits69JEM/c0yGmY96\nxMfrV7ZuqGlHXE+WWdHO/c2vQs0Z9waxqkPp4nJlitYCVe7+cZnqF6JSyLdFi6VcQy5WxrqFqCTy\nbdFiKZdjOvC4mb1oZkeWqQ0hKoF8W7RYyjXksq27zzCz1YAnzOwNd3+mTG0JsSyRb4sWS1kCurvP\nSP9/aGb3AwOBJZx+TvU1iz6vWLU1K1ZtXQ5zxHeA12pm83rNnLK3k9W3L6v+bNHnQVXLMahq+bLb\nJnLK6zXJHzBpWd8UNbOVgDbu/qmZdQCGAGc3pO1afWxzNy++o2xa1Y1Nq7ot+n732ZOavY3G+PYv\nqoMjT4is9KtK/oA+PeCdmxt0OaA8PfQewP1m5mn9t7v7mDK0I8SyRr4tWjTNHtDdfTKQ4QFyIVoX\n8m3R0jH3eAaRsjRs5jf4iJKabswO6/mbDw013W1mqNnf7ww1L/lWoeZnb94barx9KGFCn9KJEzW1\n8QMW7X1BqDnmpNtiY4bEPrL6zpNDzYzn14vbah+31eatODmF4W1w9wzZJ82PmfmXJZ5SX+7iDJUM\nzZI4E2+r0X3imX2GX/Rw3NaguK3J2/UINd0XfBhqOlwQ799HqrcPNVULno7bei2DL/0tlpz+29+G\nmlcy9AWe+qSqZPlO7drxSMfORX1bz9MKIUROUEAXQoicoIAuhBA5QQFdCCFyggK6EELkBAV0IYTI\nCQroQgiRExTQhRAiJ1Q0schGl57l45t+cSLr4E3HhpoHiGfu+TlXhpqfck+o6ejzQ81Onz0Vatqf\nWLr8/uvjhKq1eC/UrOJzQ83LDAg1P1zy/VRL2rN2/PIsPyCUsPsFfwk1j7TZt6KJRfO/KN5XGtI+\nflvAFrwSaq56+aRQc/dWu4WaT7xzqBnC46Gm1wMfhZpoxh0gU/76hB02CDXjfMtQ05c3Qs2jxMfa\nlowLNSPm3RFqNuz0dsnywXTmujabKLFICCHyjgK6EELkBAV0IYTICQroQgiRExTQhRAiJyigCyFE\nTlBAF0KInKCALoQQOaGiiUXDa0eV1Ozqj4b1HGh3x41dEZ+3Xjjhe6FmIBPitm6J27r7oDjZYx97\nqLTgibidrwbGeTXLdy6d3AXANnFb/nDclnXP0NbkuK2J6/UONZvY1IomFh1Re3nR8hO5NKyjL++E\nmuP5n1Bz5fMnhxrbJt4vw7k51PRte3ioOWth3JaNjX2g9t8Z/G1k3NZTtm2o6eTzQs2WvBZqhnF7\nqLnv76Uz64Z0gTH9i8/GpR66EELkBAV0IYTICQroQgiRExTQhRAiJyigCyFETlBAF0KInKCALoQQ\nOUEBXQghckKGuUHKx1tWetaR79naYR0jam8ONXf+JLZlHp1CjZ/dNtRMOmut2B5GhJphfUq39YdJ\nx4V1bMRboeb7GaaQeeT5fULN56wYatpyUKhZd93vh5rp9Aw1MDWDpny8Tr+iZX3fnBKuf/lGx4Sa\niRYnqB2yzR9Dzc3jY7/uv8XIUNNu4bGhZgqrh5o3d9gu1EzeYZ1Qc/Td8e/quk+fULPF3yeGmvW3\nixOLdsswgxhPBglTgbnqoQshRE5QQBdCiJyggC6EEDlBAV0IIXKCAroQQuQEBXQhhMgJCuhCCJET\nFNCFECInLPWMRWY2CtgNmOnum6XLugB3Ab2BKcC+7v5JkfV9j9o7Srbx8N/3C+34plecG3XWOqeE\nmv3sL6Hmv75+qJll3UPNeh7PRrPjdc+WLK8dFFbByf3PCTW/++zMULNCnFfB6C33CDXDhwazMAHv\nPBYnnvyU+0LNBNt2qWcsag7fZmht0foPeeSa0IYbDzs+1PimoYR7fxknHw1776+hpu3kePaf537U\nP9T8jaGhpn+GmcF+zR9CzRjfOdTMtNjfRvidoWZzeyXUHMItoeZx36VkeV96M7LNsLLMWHQTUL/1\nU4En3X0jYCxwWhPqF6JSyLdFq2SpA7q7PwN8XG/xnrDoNHQLsNfS1i9EpZBvi9ZKc4+hd3f3mQDu\nPgNYrZnrF6JSyLdFi0c3RYUQIic099sWZ5pZD3efaWarA7NKiSdW37Po86pVm7Bq1SbNbI74rvBp\nzTg+rRlXziYa5du8Xf3t565V0K2qjKaJPPN+zSSm1/wXgCl0LqltakC39K+Oh4BDgAuBg4EHS628\ncfWwJjYvRMLKVQNYuWrAou8zz76xqVU2ybfZoLqp7QsBwJpVfVizKnlvbl9689g5dxfVLvWQi5nd\nAfwL2NDM3jWzQ4ELgJ3N7E1gp/S7EK0K+bZorSx1D93d9y9StNPS1ilES0C+LVorS51Y1OSGzfzJ\n2sElNdvf/Vxczz5xwoPfHl+I/OqA34eaSzI8ejzDuoSaBd4+1KzDjJLl9sP4N9UeHefV2IHx9rN7\nMrS1oHna4sO4rbHdS/sNwE727FInFjUVM/MrFh5etPy4B+LhIPtJhv1yVIb9sleG/bJrhv0yK27L\nH8rQ1hEZjteD45mGuv3p/VAzp32cNHStHRJqdvCnQs2GGWbIuogTYk1t6STI7WnPve26lyWxSAgh\nRAtCAV0IIXKCAroQQuQEBXQhhMgJCuhCCJETFNCFECInKKALIUROUEAXQoic0Nwv52oUOz/6TMny\nk/eJZ9w57/I4CeG2kfuEmvXtv6FmqvcINe2I7ZlkG4Sap314yfKd/xEnMK32Sf1Xei/JW/QJNT2H\nxYlQXQ74MtTU7hpvm686hBKqqY5FS8xPsWzZvk3xZJQr9zoiXP/nq8Xbyk6N7WjzTZw4WHtg3NaY\nP28XanZZ5e+h5hI/LtQcd01s8+wPeoaaD3uvHGqO+fDTUHNB95GhprMfGmpeZvtQM2f8miXL53cq\nvb566EIIkRMU0IUQIicooAshRE5QQBdCiJyggC6EEDlBAV0IIXKCAroQQuQEBXQhhMgJFZ2xiFNr\nS2p+cP4TYT1dauPkmQdt31BzpF0VatZgeqg5xq8LNT0nxzbzcuniOcNWCKuY0nZBqBmwbmyKj84w\nE81W8Uw0t7JfqOnP+FBzg8WJOVfbKRWdsaj2reLl4/r0Devo882kUNOpXZzM9ULbuM/W7/NYM32F\nOJHHMsSSPrfHMw09c8CAUPM0VaHmWK4ONcsv+DrUdFjhm1DT5t24nt694+TFb7x0ruePWIk72qyl\nGYuEECLvKKALIUROUEAXQoicoIAuhBA5QQFdCCFyggK6EELkBAV0IYTICQroQgiREyo6YxHblM77\neObsncMqbM84mcEfjGdkOe7MeBahzZkYav4/3UPNH9/4VahhWOlEna4/js/FXVbLkFczKU4Isifi\nttpOiPfDCf23CDWH/md0qDmm36WhpuLcWbzo/t/uFa5+3thzQ83YnQeFmu32ypAU1j5OnLnKLwo1\nl915WqjhgNLJhABj7ZRQM9ifDTVd+CK2J4Nv37J76dnDAPboFe/Th+6J6+HC0sUfBbtcPXQhhMgJ\nCuhCCJETFNCFECInKKALIUROUEAXQoicoIAuhBA5QQFdCCFyggK6EELkhKWescjMRgG7ATPdfbN0\n2VnAkcCsVPYbd3+syPrO0NJJBgc9cm1ox32f/TTU/KJDnIjyb/teqNmh9qm4rVfiGYseH7BdqOnH\n6yXLb+HgsI6fcH+o6fXZe6Hmsg4jQ83pYy4JNf13eS7UdPXZoWZC7eahZk67tZd6xqLm8O07avco\nWv9btlFowxY+LtTcRZyocsNnR4aaVzvEvv8njgo1N407NtSMHHBBqHnXeoWaF9k61Jzp54SaGzks\n1PRnQqj5S208K1pn+yTUTO7dr2T5kO1gzB1WlhmLbgJ2aWD5Je4+IP1r0OGFaOHIt0WrZKkDurs/\nAzQ0OWZF5nEUormQb4vWSjnG0I8zs/FmdoOZdS5D/UJUCvm2aNE098u5rgHOcXc3s3OBS4DDi6rf\nrv72c9cq6FbVzOaI7wpf1zzL10/HL2xqAo3y7Xuqv32R2yZVq7JJ1arltE3kmQU18GUNAJNeLS1t\n1oDu7h8WfL0eeLjkChtUN2fz4jvMclWDWa5q8KLvC37XvG9kbKxvD6veuFnbF99hVqhK/oA+m8E7\nr51dVNrUIRejYFzRzFYvKNsbeK2J9QtRKeTbotWx1D10M7sDqAK6mdm7wFnA9ma2OVALTAGObgYb\nhVimyLdFa2WpA7q779/A4puaYIsQLQL5tmitVHTGopGPlE4y2MUeD+v4osNKoeZRdg014y78Qah5\ncOKIUMM38ZNtQ/b8R1zPvUHC111x0sS9tbeFmis7nBBqTj8zThri61gyYeg2sWj9LLMsPRlrKkyV\nPV20bArrhutPtXVCzcr+aahpG09GxKD/xIkzZKhniwH/DDXTrWeoWdnnh5r3/rRhqBlz1JBQ82zb\nrULN4QtHhZprnzkx1MxZsGao4YCgfAPgjuLFSv0XQoicoIAuhBA5QQFdCCFyggK6EELkhBYT0KfV\nvFNpExrPBzWVtqDRzK95pdImNJ7PayptQZP4V02GO4otjRdqKm1B43mzptIWNI4JNc1epQJ6U5hR\nU2kLGs38mvGVNqHxfFFTaQuaRKsM6C/WVNqCxvNWTaUtaByv1jR7lS0moAshhGgaFX0OfS2+zabu\nxMqLfU+Wxc+arkvXULMcK8fG9IglfLH41+mToWfvepqFGepZJYMmelR5wBphFRs38Lu/YvnFlneg\nb1jPgPjR4UzPKTMgg2atJRdNnwg9C1+N0qljWM24eH6IstKO/os+t2EK7Vhn0ffVG/qR9ViRz0NN\nL7qHGmsbb/QBKyy5bHo76Fm4PINfb0SHUNOLbqFmeRaEmgENvOts+krQs2B5ltgwYEDcp+1WsO+K\n1pMhxLD84l+nt4eeneppgkfV+3SDMSXKl3rGoqZiZpVpWHxnWNoZi5qKfFuUm2K+XbGALoQQonnR\nGLoQQuQEBXQhhMgJLSKgm9mPzWyimb1lZqdU2p4smNkUM5tgZq+Y2QuVtqchzGyUmc00s1cLlnUx\nszFm9qaZPd6SplIrYu9ZZjbNzMalfz+upI2NpbX5tvy6PCwr3654QDezNsBVJLOs9wNGmFlrmO6l\nFqhy9y3cfWCljSlCQ7PXnwo86e4bAWOB05a5VcVpyF6AS9x9QPr32LI2amlppb4tvy4Py8S3Kx7Q\ngYHA2+4+1d2/BkYDe1bYpiwYLWP7FaXI7PV7Arekn28B9lqmRpWgiL1QMHNQK6M1+rb8ugwsK99u\nCTtuTeC9gu/TCJ/GbBE48LiZvWhmR1bamEbQ3d1nArj7DGC1CtuThePMbLyZ3dDSLqUDWqNvy6+X\nLc3q2y3VFwiAAAABSklEQVQhoDd0hmoNz1Ju6+5bAbuS7JR4hgyxNFwDrO/umwMzgAyzbbQYWqNv\ny6+XHc3u2y0hoE8DehV8XwuYXiFbMpP2Aupmg7+f5PK6NTDTzHrAoomPZ1XYnpK4+4f+bbLE9cDW\nlbSnkbQ635ZfLzvK4dstIaC/CPQxs95mtjwwHHiowjaVxMxWMrOV088dgCG03FngF5u9nmTbHpJ+\nPhh4cFkbFLCYvenBWcfetNzt3BCtyrfl12Wn7L5d0Xe5ALj7QjM7nuQVBW2AUe7+RoXNiugB3J+m\neLcDbnf3Uq9YqAhFZq+/ALjbzA4D3gX2qZyFi1PE3u3NbHOSpy+mAEdXzMBG0gp9W35dJpaVbyv1\nXwghckJLGHIRQgjRDCigCyFETlBAF0KInKCALoQQOUEBXQghcoICuhBC5AQFdCGEyAkK6EIIkRP+\nD7ZNrAVaR8hwAAAAAElFTkSuQmCC\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAADDCAYAAACS2+oqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXm8H9P5x99PZEGIrEIaCRFrVAgiVLm1RPlRikSC2tqg\nraJaSqlcS2tpG0vRKqEootTaKlFcrRJbSEutFZGIBEkk9kTu8/tj5ibf3Nzv95m7fPO9d/J5v173\ndb8z5zPnPHPmmWfOnJkzx9wdIYQQbZ92lTZACCFEy6CALoQQOUEBXQghcoICuhBC5AQFdCGEyAkK\n6EIIkRNaXUA3sxfMbOdK27EyY2b3mdm3mrH9b83sjJa0aWXEzGrNbECJ9NyeK/LBJuLu4R9wCPA0\n8CHwNvBX4CtZtg3yvQ44p7n5VPIv3YfPgQXp34fAc5W2K4PdY4GFBTYvAH5cabsaYfNc4DFgWCO2\nfwQ4egXY+SbwGdC93vrngVqgX8Z8FgMDCvysUecK0AE4C3g5PcbT03N3j0ofywaOp3ywBf7CFrqZ\nnQyMA84D1gb6AVcC34i2XYm40N27pH9ruvvWLV2Ama3S0nkCEwps7uLuvypDGS3NBHfvAvQEaoDb\nKmtOgzgwFRhdt8LMtgBWTdOyYs2048/AvsBhQDdgA+BSYO8GCyuPj0XIB1uS4GrSheTKeUAJTUfg\nEpKW+wzgYqBDmrYLSavgZGB2qjkyTRtDcqX7jORqd3e6fiqwa8HV8Fbg+lTzH2BIQdm1pC2YdHmZ\nVkxaxmvA+8BdwLrp+v7ptu0aunICG5IcqA+Ad4FbSux/0ZZTQTmHA9PSvH5akG7AacDrwHvABKBr\nvW2PTretSdcfTtICfA84s66+gN7Ax0C3gvy3SctcpUhL44aoFVGqLtJjPTtNex7YvDHHoeAYHgu8\nCswBLg9aRzcULG9G0ortkS53Be5N7ZyT/u6Tpp0HfAF8kvrSZen6TYGJqf4lYERB/nsDL6b66cDJ\nGVthU4GfAk8VrPslcHpqb7+GWmvAEcA/6/s3Gc6VBmzYPfWHdTPYeiowBfiUpBt2s9S2eSTn3L4N\n+UYJm38A/C89DhdlPZ7yweb7YNRC3wHolFZAMc4EhgJbAoPT32cWpK8DrAn0Ab4DXGFma7n71cBN\nJAe8i7vvVyT/fYGbgbXSyrmiIK1oa8fMdgV+ARwErAu8RRIww22Bc4EH3L0r0Bf4TQltFr4CbERy\nkp1lZpuk608kudP5Kkn9zCO5+ylkZ5IDvqeZbUay/6NJ9mmtdDvcfTbJSTCyYNtDSZx/cTNsb7Au\nzGw4sBMwME07mMQhlyHDcQD4P5KLz1bAyDTvkphZR5JgMoek3iAJRtcC65HcSX5C6i/ufibwT+D4\n1N9OMLPVSU6kP5K0tkYDV6b1DHANMMaT1tgWwMORXQVMAtY0s03MrB3Jcfkjcat7Ob9sxLlSyG7A\nk+7+TgbtKGAvkmDUDrgHuB/oBZwA3GRmGzXC5v2BIenffmZ2dAYbSiEfzOiDUUDvAbzv7rUlNIcA\nZ7v7HHefA5wNFD7MWAic6+6L3f1vwEfAJg3kU4zH3P0BTy5XN5JcOOoodXIcAox39ynuvoikdbSD\nmfXLUOYioL+ZfcndF7r744H+FDOba2bz0v/XFaQ5UJ3m82+SltDgNO0Y4Ax3fye18RzgoDQA1G07\n1t0/dffPSRzyHnd/wt2/IOkfLeQG0rpP8xhNUmfFOLie3es0oi4WkVyoNzczc/dX0otKfbIch/Pd\n/UN3n05yUdoqspnkRPk2cFCdf7r7XHe/090/d/ePgfNJLojF2AeY6u43eMLzJN0UB6XpC4FBZram\nu89P0xvDjSQn/B4k/dgzG7l9c+gJzKpbMLNu6XH+wMw+rae91N1npj42DOjs7he6+xfu/gjwFwq6\njzJwQVpfM0ju3kttKx9sQR+MAvocoGdBgGmIPiRXvDqmpeuW5FHvgvAJsEZQbiGzCn5/Aqwa2FNo\n17S6hbRy5wBfyrDtKSR185SZ/cfMjgIws9PN7EMzW2BmhS3pX7p7d3fvlv4/ql5+hU5WuP/9gTtT\nR54L/JfESXsX6GfU26fpBfv0Kcu2SO4GNjOz9YHhwAfu/kyJ/by1nt2zGtA0WBfpiX45Setjlpn9\nzswaOq5ZjkOx+ilqM8nznBeAbesSzGw1M7vKzN40sw+AR4GuZlbswt8fGFZX/2Y2j+Tkr6v/A0la\nbtPM7BEzG1bCrob4Y5rfkSQX27KR+mWdb/YlqeN169LdfZ67dyNphXast3lRH0uZRrbzpqH86seD\n+sgHW9AHo8D4BEm/3f4lNG+nRhUamLUl0pgHRA3xCbB6wXLh1X1moV1m1pnkjmMGSd8ixbZ193fd\n/Rh3/xJwHMkt0AB3P9+XPrz5XjNth+RCuFfqyHVO3bnebXJhHb1DcstZt0+rpftUZ/fnwJ9IHoId\nRunWeSaK1UWadrm7bwsMIrnrOqWBLEodh+bYNTe1p9rM6pz/RyRdW9ult+B1LaO6k6m+v00neTZR\nWP9d3P34tIxn3X1/kq6Hu0nqtjE2vkXSR70XcEcDko8p7r/LZReUtWaBb84AHgK2M7OGgmn94FKY\n90yS7oJC+pGc51ltLty+H828M5EPZvfBkgHd3ReQPAS4wsz2S68+7c1sLzO7IJVNAM40s55m1hP4\nGdkDyWyShz6NodAZnwMOMbN2ZvZ1koewddwMHGVmW5pZJ5I+tEnuPt3d3ydx0MPSbY8mefCSFGB2\nkJnVXb0/IHlo0tR+6FLdQlcBv6i79TOzXmZW+PZQ/W1vB/Y1s2Fm1oGke6s+N5K0CPclaSE2i2J1\nYWbbmtlQM2tP8jDtMxquo6LHobm2ufsrJH29P0lXrZnassDMugPV9Tap729/ATY2s8NSv+6Q7tem\n6e9DzKyLJ88gPiR5oNVYjiZ5cFm/mwOSh3gHpOfVQJLb92I06lxx9wdJug7uSo9Th/RY7UDpi8OT\nwMdmdmpaJ1Uk3QK3NMLmU8ysq5mtR/KcqH5/daOQD2b3wbDrwt0vJnlL5UySJ7dvAd9j6YPS84Bn\ngLr+4WeAn5fKsuD3eJL+oblmdkcD6dH2J5E8VJxH0k93Z4HdD5NcXO4gCd4bkDz8qWMMydP990me\nVP+rIG074EkzW5Du5wnuPo3inJre6i5Ib3vfLWJv/eVLSa66E81sPvA4yUPlBrd19/+SvEFwK0mr\nYz7JMfm8QPM4icNPTluITaGw3GJ10QW4muRd3Kkk9bjcK2cZjkOp+snCr4AxaWPiEpLW4/skdXlf\nPe2lwAgzm2Nml7j7RyRdU6NI6nMmcAFLuyS+BUxNb52PIXnInIUl++DuU919ckNpJG9oLCLpVryO\n5S/AzT1XDiAJGH8kOUfeIDlP9ixSBmkf8zdI3q54n6RL41vu/lpGmyHx6WeBySQvMlwb2NkQ8sGE\nRvmguTe310NUivTW8QOSp/zTCtY/BNzk7k05kYRoMmZWS+KPb1TalpWRVjf0X5TGzPZJb3c7A78G\n/l0vmG8HbE3SihdCrEQooLc99iO5LZtB0u+/5NbRzP5A8k7riemTfCFWNLrlryDqchFCiJygFroQ\nQuSE9uXINH2F8BKSC8Z4d7+wAY1uDURZcffmftxqOeTbojVQzLdbvMvFklGcr5J8S2ImyWd3R7n7\ny/V0zk8Kyn6sGnaqXjaz6zMUeHEGTZZBy//JoHmmXl3dVQ37Vy+77qh4rMLFteeFmh+++ruS6Zdt\nPCbM44SHr15+5fXVcET1ksW+u722vKYezy/5UkFxruE7oWYhnULNVX7scusWVF9Gl+oTliy//VSp\nz4qkDLMWD+iN8e0LknEhADxY/SR7VG+/ZPlyfhCWdZjHQwh+PffHoSbL6f3FQ2suv/K2ahhRvWTx\n+JEXhfn85p7TYntKfRUqpdOvPgg1b3Vf/gsev6r+jB9Xr7pkeY2FH4X5DO40JdR8x68JNac9fmmo\nOX/Hk5ZZfqj6CXar3mGZdWdMLh3QhneBiRsX9+1ydLkMBV5z92npO60TSB7kCdHWkW+LVk05AvqX\nWPZbEDNo3HcghGityLdFq6YcfegN3Qo0fOP3WPXS3526lsGUMrNpVaUtaDyDqyptQaPpVLV9LHq2\nBibXlNuUzL79YPWTS36v2rX+t7DaAJtXVdqCRrNjVVkeCZaNDar6xiKAZ2oS/wZeD3oty1EDM0g+\nyFNHX4p9nKd+n3lboy0G9K2qKm1Bo8kU0LepSv7qGN/QZ26aTWbfLuwzb5MMqqq0BY2mrQX0AVX1\nv4FWhG2rkj9gYBd44/Livl2OLpengYFm1t+SD8CPIvlgvhBtHfm2aNW0+CXN3Reb2fEkIxbrXu16\nqaXLEWJFI98WrZ2y3KO4+/00blYiIdoE8m3Rmqlsp1P0+ah94pdoreqTUONnrB5qFhNPeL7qPvH7\nsWsuil99nhhPV8jJG5f6AjH88N34BfyTdy2dB0AVj4Sa6zky1PyXzUPNE+wQahZa/ABxy6GTQs2/\nQ0V5eY2Ni6Z9m/Hh9s9b/O7/9O7xQ7V1Xp0fai4bGY9pmM9aocYXhhIWZRg3cu9a+4aawyyecuHE\njvG74Xst93Xb5XnA9gw1XXecF2pOeye2p2ZIVcn0IfRiYol0Df0XQoicoIAuhBA5QQFdCCFyggK6\nEELkBAV0IYTICQroQgiRExTQhRAiJyigCyFETqjYnKJm5twWlP15bNtxo8eFmt7t44kAxp6fYS6E\nUxeHkioeCDUbEU8qcTXHB+mHh3lswQuhZgcmh5pJbB1qho2LJwq440d7hZqLODXU7MNfQs1Z9uuy\nzFiUBTNz/lJbNP2k/zs/zGNjfyXUHHfuDaFm/unxQK212n8Wajg7bvu9dlY80GkjeyvUnGLxBDCj\nfEKo2SaD//NevF8+N3ajKzY9OtRM8/6hpr9NK5nej0Hsbyev0AkuhBBCVAAFdCGEyAkK6EIIkRMU\n0IUQIicooAshRE5QQBdCiJyggC6EEDmhshNcrB+krxFncfm8U0JNu8XF3wlewu/ja9s7dA8161j8\nBf8D/c+hxh8+sWT6U7vFH8s/0G8PNQyNJ/bYPp7bAB6M6/ib18R1fMV3vhdq/syBGQz6dQZN+aid\nV3xffW68/YvdNoxFP4vr/KFV4jo/YI/YB2yj2JyNjpwRavyHcVn9t47f6X6ZTUPNNtfEZXFD/I65\n/SMef7LYjws1vz7hzNie3UsnD+8FcHLRdLXQhRAiJyigCyFETlBAF0KInKCALoQQOUEBXQghcoIC\nuhBC5AQFdCGEyAkK6EIIkRMqOrDo6W0GlUzf7u1nwjzazY0nwdiu2z9DzTObhBLW/dsHoWbC4KPi\njB6OJZxeOrnDjEVhFiO5NdTs+sx+oeaMR0IJtb+IB3E8+9PSxxvg9xyTQXNsqImn2ygvLx46oGja\nLRwSbn/elfEkD5O+G088csDfQglMj8+ht7/dLdT0+GheqOl0cVzWtlvH530P5oQauzbD5D13xZLa\nR2Pf3mGXLeKMBmeYb6VTkN6hdLJa6EIIkRMU0IUQIicooAshRE5QQBdCiJyggC6EEDlBAV0IIXKC\nAroQQuQEBXQhhMgJZRlYZGZvAvOBWmCRuw9tSDfKJpTMZ3KfrTIUFs/aMpRxoebTbeKX/lddI565\nhGfja+Qth8aDeUYfdmfJ9I/86jCPB9/dP9RYbVx/fle8T/86PR7kshPxgBEejcvafpcn43zKRFbf\nvtVGFs1jkP8nLuh78XF5KcMApWGT/h2XdVbs11dnmPnorFMzDJw5Py5r2C1xWf5RXJY/HtehTYjL\nspdDCUN3iYeyDfr206HmfXqUTO/GaiXTyzVStBaocvd46JgQbQv5tmi1lKvLxcqYtxCVRL4tWi3l\nckwHHjCzp81sTJnKEKISyLdFq6VcXS47uvssM+sFPGhmL7n7Y2UqS4gViXxbtFrKEtDdfVb6/z0z\nuxMYCizn9HOrr1zye7Wq7VitartymCNWAl6omcOLNXPLXk5W3360eukXPvtX9WP9qv5lt03kk4U1\nk1hYk7wI8EIQsls8oJvZ6kA7d//IzDoDw4GzG9J2r/5eSxcvVlK2qOrBFlVL3xC47ezXW7yMxvj2\nLtVfbfHyxcpJx6phdKwaBsAWrMaL5xR/a68cLfTewJ1m5mn+N7n7xDKUI8SKRr4tWjUtHtDdfSqQ\n4QVyIdoW8m3R2jH3DLN6lKNgM7/GR5fUZJmV5G++V6hZ22aHmkP8llDzjG8bar71yp9DjUezkgBT\nBpYeOFFTG79g0ck/CzXHnXJjbMzw2EfW2WNqqJn1ZPFZfJbQKS6r3avxgBFGtcPdM4x0aXnMzO/1\nXYumD+WpMI/pteuFmiGzXgo1F617fKg5/V+XhprFveJZe+zZUIKfl+GQPBL7wH/X3iDUbHbPm6Hm\n893ist5dvWeoWe+cOFbdPXZ4qDngvtJTTA3vCROHFfdtvU8rhBA5QQFdCCFyggK6EELkBAV0IYTI\nCQroQgiRExTQhRAiJyigCyFETlBAF0KInFDRgUU2ofTsJV8Migey7rDFw6HmLuKZe37Ab0LNgdwe\natb0D0PN7h8/Emo6/bB0+p1XxwOq+jI91HT1D0LNswwJNV9d/vtUy9uzXvzxLD80lLDvBX8KNX9t\nN7KiA4se98FF059qeJKjZTjh9nhGquMOujjU7G93hZoHfM/YHuLBRwP+EQ/gY5NYwrWxxDeLNU/u\nv2Wo2ZD/hZqer38cat4buEaoyTKgbNoxpXds+OYw8WTTwCIhhMg7CuhCCJETFNCFECInKKALIURO\nUEAXQoicoIAuhBA5QQFdCCFyggK6EELkhIoOLBpVO76kZm+/L8znMLstLuyy+Lr11AlfDjVDmRKX\ndX1c1m2H7xNqRtg9pQUPxuUsHBqPq+m4VunBXQBsH5fl98Zl2doZypoal/XygP6hZnObVtGBRSf5\nz4umb+EvhHkczU2h5m2LZ9K5z/cONWO4IdTMpmuo6X1FPKiO72fwgemxD5zed2yoOd9izUiLZ+y6\novb7oaaXLQg1OxIPKJw0vPhMVwDDt4GJF2pgkRBC5B4FdCGEyAkK6EIIkRMU0IUQIicooAshRE5Q\nQBdCiJyggC6EEDlBAV0IIXJCPCVQGXnVNiqZ/mVbL8xjdO0fQs0t34xtWUCXUONnrxJqXh/bN7aH\n0aHmoIGly/rV6/Fgh014NdR8hc6h5q9Pjgg1n7BaqFmFw0PNBht8JdTMpE+ogWkZNOXjkqtOK5p2\n7rGnhNvP8m6hpiPxIJ2niGdHGnNn7Ne9N4gHINo7oYTa++Kynt17UKgZYfGsVeP83VAzxDuGml79\nPgo1PiLer2PGjQw1kw4sPbCIvsCFxZPVQhdCiJyggC6EEDlBAV0IIXKCAroQQuQEBXQhhMgJCuhC\nCJETFNCFECInKKALIUROaPKMRWY2HtgHmO3uW6brugG3Av2BN4GR7j6/yPb+jdqbS5Zx7z8ODu34\nol88Nmrs+j8JNQdnGKjwP98w1Lxra4eaAf5GqNntqidKptcOC7Pg1MHnhJpzPz4r1Kz6clzWhG2+\nEWpG7RXMwgS8cf86oeZA7gg1U2zHJs9Y1BK+zf21RfNfvG3cjrJjMhh6SSyZsl7pwXsAnf2TULPR\nLW+HGv9DbM+ip2PNxHlfCzXtMwyqmufxLEuj3o190l4LJbBBLKn9MHbHnTd5oGT6UHpwcbttyjJj\n0XXAnvXWnQb83d03AR4GTm9G/kJUCvm2aJM0OaC7+2PAvHqr9wOuT39fD+zf1PyFqBTybdFWaek+\n9LXdfTaAu88CerVw/kJUCvm2aPXooagQQuSElv7a4mwz6+3us81sHaDk585err59ye+eVZvTs2rz\nFjZHrCx8VDOZj2oml7OIRvk2N1Yv/b1lFQyuKqNpIs/Mr5nC/JopACxm9ZLa5gZ0S//quAc4kuQD\nj0cAd5faeNPqg5pZvBAJa1QNYY2qIUuWZ599bXOzbJZv863q5pYvBABrVQ1mrarBQPKWy6Rzfl9U\n2+QuFzO7GXgc2NjM3jKzo4ALgD3M7BVg93RZiDaFfFu0VZrcQnf3Q4ok7d7UPIVoDci3RVulyQOL\nml2wmf+9doeSmq/dNinOZ0Q8wMBvim9EfnToz0PNuAyvHs+yeKaZz7xTqFmfWSXT7avxPtUeGw9k\nsMPi+rPbM5T1WcuUxXtxWQ+vXdpvAHa3J5o8sKi5mJljxQcW/f2LeFamXe1fcUH3xnU1Zd94YNFg\nXgk1X3SLy2o/IkN1/z72gTcsnpHqRo9nvxqb4Sbqr+wWal70eAalU+2yULNoflyHF3b9Ucn0AWzM\nYXZsWQYWCSGEaEUooAshRE5QQBdCiJyggC6EEDlBAV0IIXKCAroQQuQEBXQhhMgJCuhCCJETWvrj\nXI1ij/seK5l+6oh4xp1fXLpKqLnxxBGhZkP7X6iZ5r1DTXtie163eLDHoz6qZPoe/4wHMPWaX/+T\n3svzKgNDTZ+D4oFQ3Q79PNTU7h3XzcLOoYRqqmPRcvNTrGCeKp70WLudws371MZT4Gw2PTZj8Cvx\ndDujN/5DqLlmVtz2u73jPqHm4JNjH3h4XLGBuks5a95FoWZBl3hKp6HtO4aahyweIPzP2u1CzS4v\nFx9stoTSExYxfEOAY4umq4UuhBA5QQFdCCFyggK6EELkBAV0IYTICQroQgiRExTQhRAiJyigCyFE\nTlBAF0KInFDRGYs4rfSL9jud/2CYT7faePDM3TYy1Iyxy0PNuswMNcf5VaGmz9TYZp4tnTz3oFXD\nLN5c5bNQMyQev4JPyDAb0bbxTDQ3cHCoGczzoeYa+06oucJ+UtkZix4pcV79KT7n1rvi1VAzjY1D\nzQBeCjVj7exQM8UHh5rRdkuomeF9Q82mvBxqtp7zXKj57MV48J13yuDb28e+fRS/CzXX33BcqOl5\neOnRYlWsyu3temvGIiGEyDsK6EIIkRMU0IUQIicooAshRE5QQBdCiJyggC6EEDlBAV0IIXKCAroQ\nQuSEis5YxPalX+p/7Ow9wixsv3iQht8dz5Ly/bPiWYS2yjDg4busHWp++9KPQg0HlR7M0P3r8bW4\nW68M42pejwdN2INxWatMiY/DCYO3DjVH/XdCqDlu0MWhptJcusuYomkntv99uH2WwVMv+YxQM5Ij\nQs0RHtc5r8c+8OjAoaHmm9wXat6wPqFmXPcfhhp2jmcIuoevh5oDXomL8g7HhJprDx8dao5+vPTg\nrAVdS2+vFroQQuQEBXQhhMgJCuhCCJETFNCFECInKKALIUROUEAXQoicoIAuhBA5QQFdCCFyQpNn\nLDKz8cA+wGx33zJdNxYYA7ybyn7q7vcX2d7Zq/SL/4f/NZ4F5I6PDww1J3WOB6L8x74canatfSQu\n67l4xqIHhuwcagbxYsn06zMMGPkmd4aafh+XniEF4JLOJ4aaMyaOCzWD95wUarr7nFAzpXarUDO3\n/XpNnrGoRXx76ufFC3imY2zElfF52euht0LNe2+sF2q+OiCeGexk4uP7W/tuqBnht4WaM+znoWY4\nE0NN1wyzmV1+xCmhhj/GNtMznhWNOX+LNT32Lpk8vAom3m5lmbHoOmDPBtaPc/ch6V+DDi9EK0e+\nLdokTQ7o7v4Y0NAlsCLzOArRUsi3RVulHH3o3zez583sGjNbqwz5C1Ep5NuiVdPSH+e6EjjH3d3M\nzgPGAd8uqn6teunv7lXQo6qFzRErC4tqnmDRo0+Us4jG+fYl5y79PWxnGLZLOW0TeWZhDSyqAeD1\n0o/WWjagu/t7BYtXA/eW3GCj6pYsXqzEdKjagQ5VOyxZ/uzclv0iY6N9+6SftWj5YiWmY1XyBwwc\nBG+8dHZRaXO7XIyCfkUzW6cg7QDghWbmL0SlkG+LNkeTW+hmdjNQBfQws7eAscDXzGwroBZ4Ezi2\nBWwUYoUi3xZtlSYHdHc/pIHV1zXDFiFaBfJt0VZp8sCiZhds5ifW/qKkZk97IMznOj8q1PyPAaFm\n8kU7hZoMExbBFxnebMswyxJ/DjS3nhNmsU3trqHmN35CqNnxrOdCDYtiCRdl2O8Ns8yy9PcMhQ1v\n8sCi5mJmDlOKC9bdMs7kpxkKylKfu2XIZ34s6XD1glCzVY/nQ810+oaajh4701v/2iTUMD6W8Ie5\nGUTdMmjOjSU9zoo1q5Z22eG7wMSbyzOwSAghRCtCAV0IIXKCAroQQuQEBXQhhMgJrSagz6h5o9Im\nNJ53aiptQaP5sCbDA87Wxic1lbagmTxdaQMaz3s1lbag8TxXU2kLGkc6+rMlUUBvDrNqKm1Bo/mw\nJn4TodXxaU2lLWgmz1TagMbzfk2lLWg8CuitJ6ALIYRoHi39ca5G0Zelo6m7sMYyy8m6jcM8NqB7\nqOnAGrExvWMJny67OHMq9OlfT7M4Qz5dM2g2CNKHrBtmsWkD+72Qjsus78xmYT5D+oQS+CKDZkgG\nTQOvKc98GfpsWrCiy5phNpMnZyirjAwZstqS3zNntqdPn6XL9MyQQa8Mmi0yaOr7Z0N8tPyqmbOg\nT8HwjfarxG2/TTKcZ12JJ/fokKGd2bPz8utmdoA+heuz7PuQVTKIshCfj9T7PufMqdCn/nneqXQW\nA9en5NQeFR1YVJGCxUpDZQcWCVE+ivl2xQK6EEKIlkV96EIIkRMU0IUQIie0ioBuZl83s5fN7FUz\n+0ml7cmCmb1pZlPM7Dkze6rS9jSEmY03s9lm9u+Cdd3MbKKZvWJmD7SmqdSK2DvWzGaY2eT07+uV\ntLGxtDXfll+XhxXl2xUP6GbWDricZJb1QcBoM9u09Fatglqgyt23dvehlTamCA3NXn8a8Hd33wR4\nGDh9hVtVnIbsBRjn7kPSv/tXtFFNpY36tvy6PKwQ3654QAeGAq+5+zR3XwRMAParsE1ZMFpH/RWl\nyOz1+wHXp7+vB/ZfoUaVoIi9UDBzUBujLfq2/LoMrCjfbg0H7kvA9ILlGem61o4DD5jZ02Y2ptLG\nNIK13X02gLvPItsbz5Xm+2b2vJld09pupQPaom/Lr1csLerbrSGgN3SFagvvUu7o7tsCe5MclAwz\nZIgmcCVQ4s+9AAABMUlEQVSwobtvBcwCxlXYnsbQFn1bfr3iaHHfbg0BfQbQr2C5LzCzQrZkJm0F\n1M0GfyfJ7XVbYLaZ9YYlEx+/W2F7SuLu7/nSwRJXA9tV0p5G0uZ8W3694iiHb7eGgP40MNDM+ptZ\nR2AUcE+FbSqJma1uZmukvzsDw2m9s8AvM3s9Sd0emf4+Arh7RRsUsIy96clZxwG03npuiDbl2/Lr\nslN2367ot1wA3H2xmR1P8omCdsB4d3+pwmZF9AbuTId4twducvdSn1ioCEVmr78AuM3MjgbeAkZU\nzsJlKWLv18xsK5K3L94Ejq2YgY2kDfq2/LpMrCjf1tB/IYTICa2hy0UIIUQLoIAuhBA5QQFdCCFy\nggK6EELkBAV0IYTICQroQgiRExTQhRAiJyigCyFETvh/2lCfM9vFM3EAAAAASUVORK5CYII=\n", "text/plain": [ - "" + "" ] }, "metadata": {}, diff --git a/docs/source/usersguide/mgxs_library.rst b/docs/source/usersguide/mgxs_library.rst index daf99d481..32e4c550c 100644 --- a/docs/source/usersguide/mgxs_library.rst +++ b/docs/source/usersguide/mgxs_library.rst @@ -171,7 +171,9 @@ attributes/sub-elements required to describe the meta-data: provided via the ``scatt_type`` element above, is represented and thus used during the scattering process. Specifically, the options are to either convert the Legendre expansion to a tabular representation or leave it as - a set of Legendre coefficients. Converting to a tabular representation will cost memory but can allow for a decrease in runtime compared to leaving as a set of Legendre coefficients. This element has the following + a set of Legendre coefficients. Converting to a tabular representation + will cost memory but can allow for a decrease in runtime compared to + leaving as a set of Legendre coefficients. This element has the following attributes/sub-elements: :enable: diff --git a/openmc/mgxs/library.py b/openmc/mgxs/library.py index 886eebf45..924189264 100644 --- a/openmc/mgxs/library.py +++ b/openmc/mgxs/library.py @@ -778,18 +778,18 @@ class Library(object): Scattering order for this data entry. Default is None, which will set the XSdata object to use the order of the Library. - tabular_legendre : {None, bool} + tabular_legendre : None or bool Flag to denote whether or not the Legendre expansion of the scattering angular distribution is to be converted to a tabular representation by OpenMC. A value of `True` means that it is to be - converted while a value of 'False' means that it will not be. + converted while a value of `False` means that it will not be. Defaults to `None` which leaves the default behavior of OpenMC in place (the distribution is converted to a tabular representation). - tabular_points : {int} - This parameter is not used unless the `tabular_legendre` is set to - `True`. In this case, this parameter sets the number of - equally-spaced points in the domain of [-1,1] to be used in - building the tabular distribution. Default is `33`. + tabular_points : int + This parameter is not used unless the ``tabular_legendre`` + parameter is set to `True`. In this case, this parameter sets the + number of equally-spaced points in the domain of [-1,1] to be used + in building the tabular distribution. Default is `33`. Returns ------- @@ -947,18 +947,18 @@ class Library(object): Cross section set identifier (i.e., '71c') for all data sets (if only str) or for each individual one (if iterable of str). Defaults to '1m'. - tabular_legendre : {None, bool} + tabular_legendre : None or bool Flag to denote whether or not the Legendre expansion of the scattering angular distribution is to be converted to a tabular representation by OpenMC. A value of `True` means that it is to be - converted while a value of 'False' means that it will not be. + converted while a value of `False` means that it will not be. Defaults to `None` which leaves the default behavior of OpenMC in place (the distribution is converted to a tabular representation). - tabular_points : {int} - This parameter is not used unless the `tabular_legendre` is set to - `True`. In this case, this parameter sets the number of - equally-spaced points in the domain of [-1,1] to be used in - building the tabular distribution. Default is `33`. + tabular_points : int + This parameter is not used unless the ``tabular_legendre`` + parameter is set to `True`. In this case, this parameter sets the + number of equally-spaced points in the domain of [-1,1] to be used + in building the tabular distribution. Default is `33`. Returns ------- @@ -1046,18 +1046,18 @@ class Library(object): Cross section set identifier (i.e., '71c') for all data sets (if only str) or for each individual one (if iterable of str). Defaults to '1m'. - tabular_legendre : {None, bool} + tabular_legendre : None or bool Flag to denote whether or not the Legendre expansion of the scattering angular distribution is to be converted to a tabular representation by OpenMC. A value of `True` means that it is to be - converted while a value of 'False' means that it will not be. + converted while a value of `False` means that it will not be. Defaults to `None` which leaves the default behavior of OpenMC in place (the distribution is converted to a tabular representation). - tabular_points : {int} - This parameter is not used unless the `tabular_legendre` is set to - `True`. In this case, this parameter sets the number of - equally-spaced points in the domain of [-1,1] to be used in - building the tabular distribution. Default is `33`. + tabular_points : int + This parameter is not used unless the ``tabular_legendre`` + parameter is set to `True`. In this case, this parameter sets the + number of equally-spaced points in the domain of [-1,1] to be used + in building the tabular distribution. Default is `33`. Returns ------- diff --git a/src/mgxs_header.F90 b/src/mgxs_header.F90 index 519da708c..753f4e149 100644 --- a/src/mgxs_header.F90 +++ b/src/mgxs_header.F90 @@ -277,7 +277,7 @@ module mgxs_header real(8), allocatable :: scatt_coeffs(:, :, :) real(8), allocatable :: input_scatt(:, :, :) real(8), allocatable :: temp_scatt(:, :, :) - real(8) :: dmu, mu, norm, p0, m, mu0 + real(8) :: dmu, mu, norm integer :: order, order_dim, gin, gout, l, arr_len integer :: legendre_mu_points, imu @@ -498,18 +498,6 @@ module mgxs_header (scatt_coeffs(imu - 1, gout, gin) + & scatt_coeffs(imu, gout, gin)) end if - ! Now create CDF from fmu with the analytical integral of a - ! piecewise linear function - if (imu > 1) then - p0 = scatt_coeffs(imu - 1, gout, gin) - mu0 = mu - dmu - m = (scatt_coeffs(imu, gout, gin) - & - scatt_coeffs(imu - 1, gout, gin)) / dmu - norm = norm + & - HALF * m * mu * mu + & - (p0 - m * mu0) * mu + & - (HALF * m * mu0 * mu0 - p0 * mu0) - end if end do ! Now that we have the integral, lets ensure that the distribution ! is normalized such that it preserves the original scattering xs diff --git a/src/scattdata_header.F90 b/src/scattdata_header.F90 index 6fca89190..ee101d63b 100644 --- a/src/scattdata_header.F90 +++ b/src/scattdata_header.F90 @@ -398,32 +398,18 @@ contains end do ! Re-normalize fmu for numerical integration issues and in case - ! the negative fix-up introduced un-normalized data + ! the negative fix-up introduced un-normalized data while + ! accruing the CDF norm = ZERO do imu = 2, order norm = norm + HALF * this % dmu * & (this % fmu(gin) % data(imu - 1, gout) + & this % fmu(gin) % data(imu, gout)) + this % dist(gin) % data(imu, gout) = norm end do if (norm > ZERO) then this % fmu(gin) % data(:, gout) = & this % fmu(gin) % data(:, gout) / norm - end if - - ! Now create CDF from fmu with the analytical integral - this % dist(gin) % data(1, gout) = ZERO - do imu = 2, order - p0 = this % fmu(gin) % data(imu - 1, gout) - mu0 = this % mu(imu - 1) - mu1 = this % mu(imu) - m = (this % fmu(gin) % data(imu, gout) - p0) / (mu1 - mu0) - this % dist(gin) % data(imu, gout) = HALF * m * mu1 * mu1 + & - (p0 - m * mu0) * mu1 + & - (HALF * m * mu0 * mu0 - p0 * mu0) - end do - ! Ensure we normalize to 1 still - norm = this % dist(gin) % data(order, gout) - if (norm > ZERO) then this % dist(gin) % data(:, gout) = & this % dist(gin) % data(:, gout) / norm end if diff --git a/tests/test_mg_basic/results_true.dat b/tests/test_mg_basic/results_true.dat index 1c58db0d7..7cdda73b5 100644 --- a/tests/test_mg_basic/results_true.dat +++ b/tests/test_mg_basic/results_true.dat @@ -1,2 +1,2 @@ k-combined: -1.003952E+00 4.160185E-02 +1.047136E+00 2.765964E-02 diff --git a/tests/test_mg_max_order/results_true.dat b/tests/test_mg_max_order/results_true.dat index 3b9625201..4db04e089 100644 --- a/tests/test_mg_max_order/results_true.dat +++ b/tests/test_mg_max_order/results_true.dat @@ -1,2 +1,2 @@ k-combined: -1.026398E+00 7.824554E-02 +1.102093E+00 4.962190E-02 diff --git a/tests/test_mg_nuclide/results_true.dat b/tests/test_mg_nuclide/results_true.dat index 1c58db0d7..7cdda73b5 100644 --- a/tests/test_mg_nuclide/results_true.dat +++ b/tests/test_mg_nuclide/results_true.dat @@ -1,2 +1,2 @@ k-combined: -1.003952E+00 4.160185E-02 +1.047136E+00 2.765964E-02 diff --git a/tests/test_mg_tallies/results_true.dat b/tests/test_mg_tallies/results_true.dat index 94ec79286..4a654639a 100644 --- a/tests/test_mg_tallies/results_true.dat +++ b/tests/test_mg_tallies/results_true.dat @@ -1,86 +1,86 @@ k-combined: -1.003952E+00 4.160185E-02 +1.047136E+00 2.765964E-02 tally 1: -3.027550E+00 -2.152400E+00 -8.573015E-02 -1.716238E-03 -4.377236E+00 -4.480918E+00 -3.577293E-02 -2.996921E-04 -8.868872E-02 -1.842057E-03 -3.189586E+00 -2.165036E+00 -9.625303E-02 -1.898094E-03 -4.691910E+00 -4.579570E+00 -4.156421E-02 -3.532611E-04 -1.030465E-01 -2.171318E-03 -2.987648E+00 -1.869303E+00 -9.105156E-02 -1.775628E-03 -4.522682E+00 -4.342733E+00 -3.968997E-02 -3.424105E-04 -9.839988E-02 -2.104625E-03 -3.225648E+00 -2.102269E+00 -8.764034E-02 -1.567260E-03 -4.629948E+00 -4.315388E+00 -3.574267E-02 -2.684687E-04 -8.861368E-02 -1.650142E-03 -3.394894E+00 -2.332105E+00 -8.432398E-02 -1.485500E-03 -4.659277E+00 -4.398368E+00 -3.230884E-02 -2.317511E-04 -8.010050E-02 -1.424457E-03 -3.013158E+00 -2.009646E+00 -7.490989E-02 -1.178333E-03 -4.123109E+00 -3.682531E+00 -2.869447E-02 -1.710786E-04 -7.113970E-02 -1.051534E-03 -3.157181E+00 -2.100927E+00 -8.531655E-02 -1.520341E-03 -4.523209E+00 -4.302173E+00 -3.468239E-02 -2.516726E-04 -8.598504E-02 -1.546905E-03 -2.922890E+00 -1.788820E+00 -8.708247E-02 -1.621130E-03 -4.305352E+00 -3.866481E+00 -3.740038E-02 -3.097101E-04 -9.272351E-02 -1.903632E-03 +2.490837E+00 +1.335515E+00 +8.299408E-02 +1.628477E-03 +3.917892E+00 +3.395475E+00 +3.777816E-02 +3.574150E-04 +9.366009E-02 +2.196850E-03 +2.410307E+00 +1.211320E+00 +7.104710E-02 +1.056830E-03 +3.619111E+00 +2.712466E+00 +3.046789E-02 +2.019483E-04 +7.553639E-02 +1.241275E-03 +2.616860E+00 +1.407243E+00 +7.221083E-02 +1.059301E-03 +3.719840E+00 +2.827103E+00 +2.963219E-02 +1.775955E-04 +7.346451E-02 +1.091590E-03 +2.383417E+00 +1.171863E+00 +6.968054E-02 +9.924374E-04 +3.486989E+00 +2.493855E+00 +2.962148E-02 +1.817051E-04 +7.343795E-02 +1.116850E-03 +2.177844E+00 +9.998272E-01 +6.808817E-02 +9.543925E-04 +3.302614E+00 +2.268097E+00 +3.001326E-02 +1.847730E-04 +7.440926E-02 +1.135706E-03 +2.585363E+00 +1.416531E+00 +7.208386E-02 +1.068367E-03 +3.721777E+00 +2.868016E+00 +2.981132E-02 +1.871051E-04 +7.390860E-02 +1.150041E-03 +2.740199E+00 +1.663131E+00 +8.719535E-02 +1.596752E-03 +4.146740E+00 +3.682203E+00 +3.870935E-02 +3.136109E-04 +9.596871E-02 +1.927608E-03 +2.212799E+00 +1.000035E+00 +7.449165E-02 +1.135538E-03 +3.466616E+00 +2.419923E+00 +3.399202E-02 +2.463541E-04 +8.427346E-02 +1.514215E-03 0.000000E+00 0.000000E+00 0.000000E+00 @@ -171,86 +171,86 @@ tally 1: 0.000000E+00 0.000000E+00 0.000000E+00 -3.136746E+00 -2.097844E+00 -9.001259E-02 -1.679379E-03 -4.662194E+00 -4.592725E+00 -3.801404E-02 -2.981421E-04 -9.424490E-02 -1.832530E-03 -3.106448E+00 -2.078391E+00 -9.012356E-02 -1.762294E-03 -4.513037E+00 -4.378913E+00 -3.813004E-02 -3.167893E-04 -9.453248E-02 -1.947145E-03 -2.830121E+00 -1.799064E+00 -8.110239E-02 -1.428875E-03 -4.186480E+00 -3.869678E+00 -3.421633E-02 -2.502127E-04 -8.482957E-02 -1.537931E-03 -3.117676E+00 -1.961910E+00 -8.682325E-02 -1.551480E-03 -4.541798E+00 -4.177594E+00 -3.598671E-02 -2.719848E-04 -8.921872E-02 -1.671754E-03 -2.860792E+00 -1.699774E+00 -7.734553E-02 -1.212035E-03 -4.045685E+00 -3.352230E+00 -3.136656E-02 -2.020318E-04 -7.776439E-02 -1.241788E-03 -2.807124E+00 -1.635117E+00 -8.312523E-02 -1.403341E-03 -4.138428E+00 -3.503069E+00 -3.560217E-02 -2.611276E-04 -8.826537E-02 -1.605020E-03 -2.924592E+00 -1.733173E+00 -8.343992E-02 -1.403493E-03 -4.345943E+00 -3.817213E+00 -3.513702E-02 -2.517119E-04 -8.711216E-02 -1.547146E-03 -2.646486E+00 -1.418290E+00 -7.950085E-02 -1.309718E-03 -3.966626E+00 -3.215631E+00 -3.437595E-02 -2.503895E-04 -8.522531E-02 -1.539019E-03 +2.099094E+00 +8.919271E-01 +7.540489E-02 +1.194829E-03 +3.366752E+00 +2.321133E+00 +3.536709E-02 +2.711517E-04 +8.768256E-02 +1.666633E-03 +2.413065E+00 +1.186499E+00 +7.505466E-02 +1.154624E-03 +3.609948E+00 +2.664048E+00 +3.294274E-02 +2.242161E-04 +8.167207E-02 +1.378144E-03 +2.536174E+00 +1.301621E+00 +7.790327E-02 +1.242427E-03 +3.836384E+00 +2.980367E+00 +3.409705E-02 +2.443384E-04 +8.453386E-02 +1.501825E-03 +2.666159E+00 +1.492268E+00 +7.417343E-02 +1.141120E-03 +3.912196E+00 +3.186439E+00 +3.075257E-02 +1.953163E-04 +7.624217E-02 +1.200511E-03 +2.459342E+00 +1.286534E+00 +6.409323E-02 +8.519553E-04 +3.511461E+00 +2.568245E+00 +2.550853E-02 +1.387171E-04 +6.324108E-02 +8.526242E-04 +3.227353E+00 +2.189454E+00 +8.161989E-02 +1.366996E-03 +4.465292E+00 +4.143037E+00 +3.169508E-02 +2.069605E-04 +7.857885E-02 +1.272082E-03 +3.006650E+00 +1.972749E+00 +8.846181E-02 +1.643195E-03 +4.449165E+00 +4.222523E+00 +3.779771E-02 +2.999515E-04 +9.370858E-02 +1.843651E-03 +2.591189E+00 +1.530310E+00 +7.305643E-02 +1.114250E-03 +3.732882E+00 +3.031120E+00 +3.039294E-02 +1.919600E-04 +7.535057E-02 +1.179882E-03 0.000000E+00 0.000000E+00 0.000000E+00 @@ -341,86 +341,86 @@ tally 1: 0.000000E+00 0.000000E+00 0.000000E+00 -3.421827E+00 -2.488458E+00 -9.626948E-02 -2.065197E-03 -5.040929E+00 -5.451459E+00 -4.018656E-02 -3.779606E-04 -9.963105E-02 -2.323134E-03 -2.869547E+00 -1.698223E+00 -8.662967E-02 -1.551328E-03 -4.227759E+00 -3.702584E+00 -3.740987E-02 -2.903037E-04 -9.274703E-02 -1.784351E-03 -2.588566E+00 -1.357415E+00 -7.706802E-02 -1.216264E-03 -3.821156E+00 -2.967091E+00 -3.309228E-02 -2.275536E-04 -8.204281E-02 -1.398658E-03 -2.777772E+00 -1.584289E+00 -8.469785E-02 -1.456208E-03 -4.195353E+00 -3.574265E+00 -3.693140E-02 -2.790926E-04 -9.156081E-02 -1.715441E-03 -3.118197E+00 -1.989138E+00 -9.265883E-02 -1.753400E-03 -4.574869E+00 -4.260081E+00 -3.969940E-02 -3.241960E-04 -9.842328E-02 -1.992670E-03 -2.591904E+00 -1.370188E+00 -8.148726E-02 -1.342503E-03 -3.988752E+00 -3.228928E+00 -3.608041E-02 -2.655091E-04 -8.945102E-02 -1.631951E-03 -3.123916E+00 -1.985559E+00 -9.297564E-02 -1.754830E-03 -4.650584E+00 -4.387677E+00 -3.995654E-02 -3.270558E-04 -9.906077E-02 -2.010247E-03 -2.507814E+00 -1.313596E+00 -7.534295E-02 -1.160881E-03 -3.761249E+00 -2.915759E+00 -3.258970E-02 -2.199767E-04 -8.079681E-02 -1.352086E-03 +2.667561E+00 +1.528159E+00 +6.981409E-02 +1.036744E-03 +3.696781E+00 +2.912439E+00 +2.768232E-02 +1.628195E-04 +6.863038E-02 +1.000769E-03 +2.579691E+00 +1.399684E+00 +7.252523E-02 +1.085979E-03 +3.762646E+00 +2.966415E+00 +3.021298E-02 +1.872587E-04 +7.490441E-02 +1.150985E-03 +2.326398E+00 +1.111745E+00 +6.380790E-02 +8.422824E-04 +3.373074E+00 +2.338678E+00 +2.619114E-02 +1.463613E-04 +6.493342E-02 +8.996094E-04 +2.689651E+00 +1.483190E+00 +8.007622E-02 +1.328579E-03 +4.015410E+00 +3.294872E+00 +3.443210E-02 +2.525628E-04 +8.536450E-02 +1.552377E-03 +2.712538E+00 +1.496273E+00 +7.040893E-02 +1.033599E-03 +3.801423E+00 +2.967673E+00 +2.785218E-02 +1.687241E-04 +6.905150E-02 +1.037062E-03 +3.101358E+00 +2.022582E+00 +8.316566E-02 +1.469875E-03 +4.433641E+00 +4.167950E+00 +3.363316E-02 +2.428472E-04 +8.338377E-02 +1.492660E-03 +2.669269E+00 +1.482479E+00 +7.250461E-02 +1.097959E-03 +3.815372E+00 +3.039000E+00 +2.952621E-02 +1.837209E-04 +7.320176E-02 +1.129240E-03 +2.719836E+00 +1.670890E+00 +8.176774E-02 +1.472662E-03 +4.030045E+00 +3.607044E+00 +3.530219E-02 +2.721266E-04 +8.752165E-02 +1.672625E-03 0.000000E+00 0.000000E+00 0.000000E+00 @@ -511,86 +511,86 @@ tally 1: 0.000000E+00 0.000000E+00 0.000000E+00 -3.072424E+00 -1.910502E+00 -8.487675E-02 -1.467547E-03 -4.502748E+00 -4.073633E+00 -3.505649E-02 -2.616264E-04 -8.691250E-02 -1.608086E-03 -3.358497E+00 -2.340508E+00 -9.563725E-02 -1.882737E-03 -4.919309E+00 -4.997057E+00 -4.012501E-02 -3.331881E-04 -9.947844E-02 -2.047940E-03 -2.841595E+00 -1.676656E+00 -8.961778E-02 -1.706734E-03 -4.382947E+00 -3.971256E+00 -3.977482E-02 -3.476057E-04 -9.861025E-02 -2.136557E-03 -3.031993E+00 -1.887195E+00 -1.005232E-01 -2.088650E-03 -4.706199E+00 -4.537534E+00 -4.558466E-02 -4.340386E-04 -1.130141E-01 -2.667817E-03 -3.116522E+00 -2.011384E+00 -9.595523E-02 -1.977738E-03 -4.627314E+00 -4.456240E+00 -4.188172E-02 -3.906963E-04 -1.038337E-01 -2.401414E-03 -2.906411E+00 -1.724257E+00 -9.583879E-02 -1.902835E-03 -4.514939E+00 -4.188183E+00 -4.337455E-02 -3.935555E-04 -1.075347E-01 -2.418988E-03 -2.723724E+00 -1.514959E+00 -8.518353E-02 -1.512432E-03 -4.173364E+00 -3.561677E+00 -3.761270E-02 -3.039338E-04 -9.324988E-02 -1.868128E-03 -2.336622E+00 -1.158899E+00 -8.017647E-02 -1.335514E-03 -3.700707E+00 -2.838852E+00 -3.695886E-02 -2.855605E-04 -9.162889E-02 -1.755197E-03 +2.643152E+00 +1.461242E+00 +7.777503E-02 +1.277892E-03 +3.881716E+00 +3.158579E+00 +3.318888E-02 +2.353505E-04 +8.228231E-02 +1.446581E-03 +2.707136E+00 +1.528526E+00 +7.153411E-02 +1.046460E-03 +3.793295E+00 +2.976791E+00 +2.857458E-02 +1.671752E-04 +7.084247E-02 +1.027542E-03 +2.745209E+00 +1.561299E+00 +7.557988E-02 +1.164281E-03 +3.914803E+00 +3.154166E+00 +3.099509E-02 +1.948889E-04 +7.684344E-02 +1.197884E-03 +3.141608E+00 +2.000672E+00 +8.238823E-02 +1.366316E-03 +4.442048E+00 +3.976219E+00 +3.282615E-02 +2.177028E-04 +8.138301E-02 +1.338110E-03 +3.078971E+00 +1.928690E+00 +8.927133E-02 +1.655293E-03 +4.508330E+00 +4.154031E+00 +3.778191E-02 +3.038622E-04 +9.366939E-02 +1.867688E-03 +3.342682E+00 +2.353374E+00 +9.542381E-02 +1.892530E-03 +4.916807E+00 +5.068483E+00 +4.010858E-02 +3.330650E-04 +9.943770E-02 +2.047183E-03 +2.906350E+00 +1.760817E+00 +8.559496E-02 +1.527401E-03 +4.338708E+00 +3.888336E+00 +3.663301E-02 +2.834528E-04 +9.082102E-02 +1.742242E-03 +2.201534E+00 +1.026489E+00 +7.361463E-02 +1.156634E-03 +3.444005E+00 +2.510736E+00 +3.352726E-02 +2.491193E-04 +8.312121E-02 +1.531211E-03 0.000000E+00 0.000000E+00 0.000000E+00 @@ -681,86 +681,86 @@ tally 1: 0.000000E+00 0.000000E+00 0.000000E+00 -3.107110E+00 -2.031843E+00 -8.851130E-02 -1.622664E-03 -4.517810E+00 -4.254353E+00 -3.710637E-02 -2.854685E-04 -9.199458E-02 -1.754631E-03 -3.274545E+00 -2.233343E+00 -1.058426E-01 -2.344479E-03 -4.973131E+00 -5.160486E+00 -4.734962E-02 -4.741951E-04 -1.173898E-01 -2.914639E-03 -2.636225E+00 -1.433972E+00 -7.975913E-02 -1.339963E-03 -3.942728E+00 -3.209671E+00 -3.458062E-02 -2.576972E-04 -8.573273E-02 -1.583935E-03 -2.850580E+00 -1.753714E+00 -8.409223E-02 -1.504889E-03 -4.292353E+00 -3.910272E+00 -3.606743E-02 -2.785145E-04 -8.941884E-02 -1.711888E-03 -3.177076E+00 -2.046322E+00 -9.290111E-02 -1.805479E-03 -4.671510E+00 -4.473692E+00 -3.953032E-02 -3.362191E-04 -9.800407E-02 -2.066570E-03 -2.792020E+00 -1.590839E+00 -8.956724E-02 -1.645807E-03 -4.270468E+00 -3.718893E+00 -3.995877E-02 -3.306439E-04 -9.906631E-02 -2.032302E-03 -2.319545E+00 -1.130506E+00 -7.787530E-02 -1.284059E-03 -3.631235E+00 -2.757748E+00 -3.552455E-02 -2.689956E-04 -8.807293E-02 -1.653381E-03 -2.338383E+00 -1.206709E+00 -7.569363E-02 -1.234797E-03 -3.576288E+00 -2.757162E+00 -3.391726E-02 -2.462226E-04 -8.408811E-02 -1.513407E-03 +2.379478E+00 +1.220072E+00 +7.149541E-02 +1.036720E-03 +3.597267E+00 +2.694760E+00 +3.098378E-02 +1.953397E-04 +7.681538E-02 +1.200655E-03 +3.155079E+00 +2.066218E+00 +9.171950E-02 +1.711836E-03 +4.700403E+00 +4.503318E+00 +3.904100E-02 +3.167906E-04 +9.679096E-02 +1.947152E-03 +2.737388E+00 +1.515779E+00 +7.304714E-02 +1.077128E-03 +3.854064E+00 +2.996781E+00 +2.937470E-02 +1.746013E-04 +7.282614E-02 +1.073187E-03 +3.294568E+00 +2.201280E+00 +8.370459E-02 +1.449531E-03 +4.567777E+00 +4.247816E+00 +3.258414E-02 +2.266289E-04 +8.078302E-02 +1.392974E-03 +3.040055E+00 +1.908677E+00 +8.543146E-02 +1.469528E-03 +4.395363E+00 +3.943052E+00 +3.551022E-02 +2.539139E-04 +8.803740E-02 +1.560681E-03 +2.733720E+00 +1.608810E+00 +7.575232E-02 +1.234497E-03 +3.955386E+00 +3.353987E+00 +3.126476E-02 +2.114507E-04 +7.751200E-02 +1.299681E-03 +2.633801E+00 +1.566751E+00 +7.851898E-02 +1.276803E-03 +3.867974E+00 +3.237952E+00 +3.370851E-02 +2.331450E-04 +8.357057E-02 +1.433025E-03 +2.649820E+00 +1.500516E+00 +6.849773E-02 +9.733415E-04 +3.744876E+00 +2.926356E+00 +2.706771E-02 +1.530408E-04 +6.710661E-02 +9.406650E-04 0.000000E+00 0.000000E+00 0.000000E+00 @@ -851,86 +851,86 @@ tally 1: 0.000000E+00 0.000000E+00 0.000000E+00 -2.615843E+00 -1.511122E+00 -8.062716E-02 -1.419946E-03 -3.984393E+00 -3.462948E+00 -3.538715E-02 -2.775562E-04 -8.773227E-02 -1.705998E-03 -2.882570E+00 -1.761817E+00 -8.002023E-02 -1.361536E-03 -4.138232E+00 -3.590378E+00 -3.301403E-02 -2.428126E-04 -8.184882E-02 -1.492447E-03 -3.223528E+00 -2.083805E+00 -8.733159E-02 -1.544129E-03 -4.604746E+00 -4.254674E+00 -3.551662E-02 -2.594309E-04 -8.805326E-02 -1.594591E-03 -2.927536E+00 -1.741281E+00 -8.346281E-02 -1.434485E-03 -4.324501E+00 -3.806148E+00 -3.508574E-02 -2.599007E-04 -8.698503E-02 -1.597479E-03 -3.057707E+00 -1.944892E+00 -8.714384E-02 -1.580757E-03 -4.450690E+00 -4.103723E+00 -3.652244E-02 -2.807333E-04 -9.054690E-02 -1.725526E-03 -2.719115E+00 -1.510843E+00 -7.867244E-02 -1.258358E-03 -3.983708E+00 -3.239056E+00 -3.327673E-02 -2.253729E-04 -8.250009E-02 -1.385254E-03 -2.396119E+00 -1.153282E+00 -7.566785E-02 -1.182019E-03 -3.610641E+00 -2.634248E+00 -3.348337E-02 -2.377537E-04 -8.301240E-02 -1.461353E-03 -2.154802E+00 -1.067493E+00 -6.324323E-02 -8.945795E-04 -3.142780E+00 -2.238062E+00 -2.691333E-02 -1.628077E-04 -6.672387E-02 -1.000697E-03 +3.128432E+00 +2.123023E+00 +8.683081E-02 +1.593217E-03 +4.568791E+00 +4.439900E+00 +3.595477E-02 +2.733665E-04 +8.913952E-02 +1.680246E-03 +2.955891E+00 +1.796192E+00 +8.526700E-02 +1.464085E-03 +4.355866E+00 +3.839684E+00 +3.605221E-02 +2.663621E-04 +8.938111E-02 +1.637194E-03 +2.731962E+00 +1.617682E+00 +7.685258E-02 +1.239175E-03 +4.057901E+00 +3.505399E+00 +3.210546E-02 +2.150275E-04 +7.959627E-02 +1.321666E-03 +3.058734E+00 +1.911285E+00 +8.780153E-02 +1.572171E-03 +4.531361E+00 +4.177302E+00 +3.703146E-02 +2.828466E-04 +9.180888E-02 +1.738516E-03 +2.854688E+00 +1.652859E+00 +8.501906E-02 +1.459424E-03 +4.215336E+00 +3.600073E+00 +3.650039E-02 +2.691434E-04 +9.049224E-02 +1.654289E-03 +2.532385E+00 +1.341692E+00 +7.908024E-02 +1.289149E-03 +3.816608E+00 +3.020090E+00 +3.480705E-02 +2.494950E-04 +8.629409E-02 +1.533520E-03 +2.946914E+00 +1.857842E+00 +7.586572E-02 +1.217999E-03 +4.089763E+00 +3.522165E+00 +2.980738E-02 +1.901498E-04 +7.389883E-02 +1.168755E-03 +2.788908E+00 +1.667764E+00 +7.477958E-02 +1.148212E-03 +4.029335E+00 +3.403340E+00 +3.033277E-02 +1.914118E-04 +7.520139E-02 +1.176512E-03 0.000000E+00 0.000000E+00 0.000000E+00 @@ -1021,86 +1021,86 @@ tally 1: 0.000000E+00 0.000000E+00 0.000000E+00 -3.073244E+00 -2.018709E+00 -9.073696E-02 -1.753157E-03 -4.516776E+00 -4.329334E+00 -3.876727E-02 -3.208452E-04 -9.611233E-02 -1.972074E-03 -2.671039E+00 -1.446563E+00 -7.823234E-02 -1.250837E-03 -3.946269E+00 -3.138681E+00 -3.334746E-02 -2.319341E-04 -8.267545E-02 -1.425582E-03 -3.349978E+00 -2.319515E+00 -9.950325E-02 -2.061329E-03 -5.001138E+00 -5.168172E+00 -4.276928E-02 -3.828620E-04 -1.060342E-01 -2.353260E-03 -2.978936E+00 -1.799523E+00 -9.017252E-02 -1.638102E-03 -4.406186E+00 -3.924003E+00 -3.904232E-02 -3.071501E-04 -9.679424E-02 -1.887897E-03 -2.702401E+00 -1.471209E+00 -7.923417E-02 -1.296139E-03 -4.023502E+00 -3.272988E+00 -3.383397E-02 -2.416414E-04 -8.388161E-02 -1.485248E-03 -2.759858E+00 -1.531505E+00 -8.419492E-02 -1.488052E-03 -4.150314E+00 -3.494954E+00 -3.667250E-02 -2.957414E-04 -9.091893E-02 -1.817773E-03 -2.601272E+00 -1.380301E+00 -8.755772E-02 -1.625699E-03 -4.019486E+00 -3.299452E+00 -3.989786E-02 -3.528389E-04 -9.891530E-02 -2.168723E-03 -2.510230E+00 -1.275856E+00 -8.725714E-02 -1.601325E-03 -3.973851E+00 -3.204497E+00 -4.039490E-02 -3.543820E-04 -1.001476E-01 -2.178208E-03 +3.320603E+00 +2.324411E+00 +9.169756E-02 +1.769612E-03 +4.811447E+00 +4.885591E+00 +3.780490E-02 +3.050758E-04 +9.372640E-02 +1.875147E-03 +3.157638E+00 +2.088634E+00 +8.917486E-02 +1.627414E-03 +4.530493E+00 +4.262890E+00 +3.711278E-02 +2.791427E-04 +9.201048E-02 +1.715750E-03 +2.826797E+00 +1.654653E+00 +8.798649E-02 +1.562942E-03 +4.238099E+00 +3.671654E+00 +3.864939E-02 +3.007708E-04 +9.582007E-02 +1.848687E-03 +2.639291E+00 +1.470166E+00 +8.380168E-02 +1.469150E-03 +3.995402E+00 +3.337545E+00 +3.720130E-02 +2.907469E-04 +9.222996E-02 +1.787075E-03 +2.563883E+00 +1.375196E+00 +8.901196E-02 +1.610331E-03 +4.042903E+00 +3.361278E+00 +4.116155E-02 +3.428545E-04 +1.020482E-01 +2.107354E-03 +2.651467E+00 +1.472579E+00 +7.547436E-02 +1.201468E-03 +3.877983E+00 +3.154531E+00 +3.166886E-02 +2.150564E-04 +7.851384E-02 +1.321843E-03 +2.884789E+00 +1.801722E+00 +7.517555E-02 +1.168857E-03 +4.085136E+00 +3.522983E+00 +2.985774E-02 +1.859368E-04 +7.402368E-02 +1.142860E-03 +2.810664E+00 +1.609682E+00 +6.917229E-02 +9.875658E-04 +3.854296E+00 +3.029615E+00 +2.634304E-02 +1.483398E-04 +6.531000E-02 +9.117703E-04 0.000000E+00 0.000000E+00 0.000000E+00 @@ -1191,86 +1191,86 @@ tally 1: 0.000000E+00 0.000000E+00 0.000000E+00 -2.355419E+00 -1.227267E+00 -7.298104E-02 -1.204355E-03 -3.584952E+00 -2.839605E+00 -3.206831E-02 -2.384737E-04 -7.950416E-02 -1.465778E-03 -2.885415E+00 -1.707496E+00 -8.937354E-02 -1.638109E-03 -4.328343E+00 -3.842717E+00 -3.915662E-02 -3.169623E-04 -9.707761E-02 -1.948208E-03 -2.817788E+00 -1.659313E+00 -8.795245E-02 -1.652827E-03 -4.286382E+00 -3.831800E+00 -3.878123E-02 -3.288848E-04 -9.614692E-02 -2.021489E-03 -2.987432E+00 -1.822564E+00 -8.500598E-02 -1.500428E-03 -4.380470E+00 -3.921075E+00 -3.565649E-02 -2.685990E-04 -8.840004E-02 -1.650943E-03 -2.685515E+00 -1.443739E+00 -8.434998E-02 -1.460545E-03 -4.166156E+00 -3.492673E+00 -3.738938E-02 -2.956120E-04 -9.269624E-02 -1.816978E-03 -2.686749E+00 -1.496524E+00 -8.914589E-02 -1.663073E-03 -4.156201E+00 -3.578378E+00 -4.042946E-02 -3.475741E-04 -1.002332E-01 -2.136363E-03 -2.760305E+00 -1.579416E+00 -8.914803E-02 -1.630935E-03 -4.214440E+00 -3.628492E+00 -3.989038E-02 -3.331009E-04 -9.889676E-02 -2.047404E-03 -2.681752E+00 -1.472982E+00 -8.701435E-02 -1.526110E-03 -4.175631E+00 -3.518871E+00 -3.914441E-02 -3.118235E-04 -9.704732E-02 -1.916622E-03 +3.709821E+00 +2.847102E+00 +1.037249E-01 +2.264974E-03 +5.329369E+00 +5.878403E+00 +4.295875E-02 +3.964777E-04 +1.065039E-01 +2.436949E-03 +3.237824E+00 +2.159909E+00 +9.708289E-02 +1.898400E-03 +4.765807E+00 +4.617662E+00 +4.179332E-02 +3.503645E-04 +1.036145E-01 +2.153514E-03 +2.503072E+00 +1.321745E+00 +6.623202E-02 +9.120239E-04 +3.541728E+00 +2.623055E+00 +2.654018E-02 +1.469478E-04 +6.579875E-02 +9.032144E-04 +2.360183E+00 +1.196310E+00 +7.336179E-02 +1.123869E-03 +3.552179E+00 +2.637131E+00 +3.222538E-02 +2.218615E-04 +7.989359E-02 +1.363671E-03 +2.353136E+00 +1.260383E+00 +7.514099E-02 +1.159881E-03 +3.611322E+00 +2.808437E+00 +3.349306E-02 +2.280172E-04 +8.303643E-02 +1.401507E-03 +2.404559E+00 +1.241962E+00 +6.952417E-02 +1.036673E-03 +3.551485E+00 +2.687553E+00 +2.947110E-02 +1.929942E-04 +7.306514E-02 +1.186238E-03 +2.740499E+00 +1.588007E+00 +7.592011E-02 +1.232446E-03 +3.959573E+00 +3.327411E+00 +3.134446E-02 +2.121705E-04 +7.770960E-02 +1.304106E-03 +2.647586E+00 +1.416013E+00 +7.421099E-02 +1.112815E-03 +3.826931E+00 +2.947306E+00 +3.083732E-02 +1.923794E-04 +7.645229E-02 +1.182459E-03 0.000000E+00 0.000000E+00 0.000000E+00 @@ -2892,15 +2892,15 @@ tally 1: 0.000000E+00 0.000000E+00 tally 2: -4.417493E+01 -3.967273E+02 -4.418714E+01 -3.969466E+02 -7.878379E+00 -1.248770E+01 -7.878666E+00 -1.248861E+01 -1.271565E+02 -3.251279E+03 -1.271565E+02 -3.251279E+03 +3.995283E+01 +3.210141E+02 +3.996387E+01 +3.211916E+02 +6.783537E+00 +9.416330E+00 +6.783789E+00 +9.417032E+00 +1.239468E+02 +3.134204E+03 +1.239468E+02 +3.134204E+03 diff --git a/tests/test_mgxs_library_ce_to_mg/results_true.dat b/tests/test_mgxs_library_ce_to_mg/results_true.dat index c82da3bce..6710c82be 100644 --- a/tests/test_mgxs_library_ce_to_mg/results_true.dat +++ b/tests/test_mgxs_library_ce_to_mg/results_true.dat @@ -1,2 +1,2 @@ k-combined: -1.151644E+00 3.044607E-02 +1.140804E+00 2.937150E-02