Note
Go to the end to download the full example code.
Reading Integral Indices
This example demonstrates how to read/get TSUNAMI-IP integral indices (like \(c_k\) and \(E\)) computed from a set of experiments and applications.
Reading Integral Indices From a TSUNAMI-IP .out
File
The first way to obtain integral indices is to first create a TSUNAMI-IP input file (see the SCALE manual), run it with SCALE, and then read the integral indices from the
resulting .out
file. There’s a built in function to do this in tsunami_ip_utils
, provided the TSUNAMI-IP input file
follows the expected format
:
tsnami_ip_utils.readers.read_integral_indices()
. This function just takes a path to the .out
file and returns a
dictionary of integral index matrices. Here’s an example of how to use it:
from tsunami_ip_utils.readers import read_integral_indices
from paths import EXAMPLES
# Path to the .out file
out_file = EXAMPLES / 'data' / 'tsunami_ip_hmf.out'
# Read the integral indices
integral_indices = read_integral_indices(out_file)
# Print the integral indices
for key, value in integral_indices.items():
print(f'{key}:\n{value}\n')
c_k:
[[1.0+/-0.0024 0.9986+/-0.0024 0.9941+/-0.0025 0.9892+/-0.0025
0.9729+/-0.0026 0.9701+/-0.0026 0.9582+/-0.0026 0.9032+/-0.0009
0.896+/-0.0008 0.8902+/-0.0008 0.8883+/-0.0007 0.8353+/-0.0008]
[0.9986+/-0.0024 1.0+/-0.0025 0.9976+/-0.0025 0.9939+/-0.0026
0.9798+/-0.0026 0.9772+/-0.0026 0.9663+/-0.0026 0.9036+/-0.0009
0.898+/-0.0008 0.8939+/-0.0008 0.8924+/-0.0007 0.8347+/-0.0008]
[0.9941+/-0.0025 0.9976+/-0.0025 1.0+/-0.0026 0.9991+/-0.0026
0.9909+/-0.0025 0.9891+/-0.0026 0.9811+/-0.0026 0.9192+/-0.0008
0.9145+/-0.0008 0.9115+/-0.0008 0.9102+/-0.0007 0.8491+/-0.0008]
[0.9892+/-0.0025 0.9939+/-0.0026 0.9991+/-0.0026 1.0+/-0.0025
0.9956+/-0.0025 0.9943+/-0.0025 0.9882+/-0.0025 0.927+/-0.0008
0.9229+/-0.0008 0.9204+/-0.0007 0.9192+/-0.0007 0.8564+/-0.0008]
[0.9729+/-0.0026 0.9798+/-0.0026 0.9909+/-0.0025 0.9956+/-0.0025
1.0+/-0.0024 0.9998+/-0.0024 0.9981+/-0.0024 0.9428+/-0.0008
0.9397+/-0.0008 0.9386+/-0.0007 0.9374+/-0.0007 0.871+/-0.0008]
[0.9701+/-0.0026 0.9772+/-0.0026 0.9891+/-0.0026 0.9943+/-0.0025
0.9998+/-0.0024 1.0+/-0.0024 0.9988+/-0.0024 0.9437+/-0.0008
0.9406+/-0.0008 0.9396+/-0.0007 0.9385+/-0.0007 0.8717+/-0.0008]
[0.9582+/-0.0026 0.9663+/-0.0026 0.9811+/-0.0026 0.9882+/-0.0025
0.9981+/-0.0024 0.9988+/-0.0024 1.0+/-0.0024 0.9482+/-0.0008
0.9455+/-0.0008 0.945+/-0.0007 0.9439+/-0.0007 0.8761+/-0.0008]
[0.9032+/-0.001 0.9036+/-0.001 0.9192+/-0.001 0.927+/-0.001
0.9428+/-0.0009 0.9437+/-0.0009 0.9482+/-0.0009 1.0+/-0.001
0.9988+/-0.0009 0.996+/-0.0009 0.9945+/-0.0009 0.9226+/-0.0008]
[0.896+/-0.001 0.898+/-0.001 0.9145+/-0.0009 0.9229+/-0.0009
0.9397+/-0.0009 0.9406+/-0.0009 0.9455+/-0.0009 0.9988+/-0.0009
1.0+/-0.0009 0.999+/-0.0009 0.9982+/-0.0008 0.92+/-0.0008]
[0.8902+/-0.0009 0.8939+/-0.0009 0.9115+/-0.0009 0.9204+/-0.0009
0.9386+/-0.0008 0.9396+/-0.0008 0.945+/-0.0008 0.996+/-0.0009
0.999+/-0.0009 1.0+/-0.0008 0.9998+/-0.0008 0.9169+/-0.0008]
[0.8883+/-0.0009 0.8924+/-0.0009 0.9102+/-0.0009 0.9192+/-0.0009
0.9374+/-0.0008 0.9385+/-0.0008 0.9439+/-0.0008 0.9945+/-0.0009
0.9982+/-0.0008 0.9998+/-0.0008 1.0+/-0.0008 0.9157+/-0.0008]
[0.8353+/-0.001 0.8347+/-0.001 0.8491+/-0.001 0.8564+/-0.0009
0.871+/-0.0009 0.8717+/-0.0009 0.8761+/-0.0009 0.9226+/-0.0008
0.92+/-0.0007 0.9169+/-0.0007 0.9157+/-0.0007 1.0+/-0.0011]]
E_total:
[[1.0+/-0.0024 0.9986+/-0.0024 0.9965+/-0.0024 0.9943+/-0.0024
0.9916+/-0.0024 0.9911+/-0.0024 0.9896+/-0.0024 0.9825+/-0.0023
0.9792+/-0.0023 0.9768+/-0.0023 0.9758+/-0.0023 0.9789+/-0.0023]
[0.9986+/-0.0024 1.0+/-0.0024 0.9992+/-0.0024 0.9981+/-0.0024
0.9964+/-0.0024 0.996+/-0.0024 0.9949+/-0.0025 0.9767+/-0.0023
0.9745+/-0.0023 0.9726+/-0.0023 0.972+/-0.0023 0.9721+/-0.0023]
[0.9965+/-0.0024 0.9992+/-0.0024 1.0+/-0.0024 0.9995+/-0.0024
0.9985+/-0.0024 0.9983+/-0.0024 0.9974+/-0.0025 0.9729+/-0.0023
0.9717+/-0.0023 0.9703+/-0.0023 0.9699+/-0.0023 0.9675+/-0.0023]
[0.9943+/-0.0024 0.9981+/-0.0024 0.9995+/-0.0024 1.0+/-0.0024
0.9994+/-0.0025 0.9992+/-0.0025 0.9986+/-0.0025 0.9694+/-0.0023
0.9689+/-0.0023 0.9678+/-0.0023 0.9676+/-0.0023 0.9634+/-0.0023]
[0.9916+/-0.0024 0.9964+/-0.0024 0.9985+/-0.0024 0.9994+/-0.0025
1.0+/-0.0025 0.9997+/-0.0025 0.9996+/-0.0025 0.967+/-0.0023
0.9671+/-0.0023 0.9664+/-0.0023 0.9664+/-0.0023 0.9604+/-0.0023]
[0.9911+/-0.0024 0.996+/-0.0024 0.9983+/-0.0024 0.9992+/-0.0025
0.9997+/-0.0025 1.0+/-0.0025 0.9996+/-0.0025 0.9666+/-0.0023
0.9668+/-0.0023 0.9662+/-0.0023 0.9662+/-0.0023 0.9599+/-0.0023]
[0.9896+/-0.0024 0.9949+/-0.0025 0.9974+/-0.0025 0.9986+/-0.0025
0.9996+/-0.0025 0.9996+/-0.0025 1.0+/-0.0026 0.9655+/-0.0024
0.9661+/-0.0024 0.9657+/-0.0023 0.9658+/-0.0023 0.9585+/-0.0023]
[0.9825+/-0.0023 0.9767+/-0.0023 0.9729+/-0.0023 0.9694+/-0.0023
0.967+/-0.0023 0.9666+/-0.0023 0.9655+/-0.0024 1.0+/-0.0023
0.9988+/-0.0023 0.9971+/-0.0023 0.9961+/-0.0023 0.9871+/-0.0023]
[0.9792+/-0.0023 0.9745+/-0.0023 0.9717+/-0.0023 0.9689+/-0.0023
0.9671+/-0.0023 0.9668+/-0.0023 0.9661+/-0.0024 0.9988+/-0.0023
1.0+/-0.0023 0.9992+/-0.0023 0.9986+/-0.0023 0.9836+/-0.0023]
[0.9768+/-0.0023 0.9726+/-0.0023 0.9703+/-0.0023 0.9678+/-0.0023
0.9664+/-0.0023 0.9662+/-0.0023 0.9657+/-0.0023 0.9971+/-0.0023
0.9992+/-0.0023 1.0+/-0.0023 0.9996+/-0.0023 0.9814+/-0.0023]
[0.9758+/-0.0023 0.972+/-0.0023 0.9699+/-0.0023 0.9676+/-0.0023
0.9664+/-0.0023 0.9662+/-0.0023 0.9658+/-0.0023 0.9961+/-0.0023
0.9986+/-0.0023 0.9996+/-0.0023 1.0+/-0.0023 0.9805+/-0.0023]
[0.9789+/-0.0023 0.9721+/-0.0023 0.9675+/-0.0023 0.9634+/-0.0023
0.9604+/-0.0023 0.9599+/-0.0023 0.9585+/-0.0023 0.9871+/-0.0023
0.9836+/-0.0023 0.9814+/-0.0023 0.9805+/-0.0023 1.0+/-0.0023]]
E_fission:
[[1.0+/-0.0118 0.9991+/-0.0119 0.9973+/-0.012 0.9957+/-0.012
0.9935+/-0.0122 0.9931+/-0.0121 0.9918+/-0.0123 0.9944+/-0.0116
0.9931+/-0.0117 0.9914+/-0.0116 0.9904+/-0.0117 0.994+/-0.0116]
[0.9991+/-0.0119 1.0+/-0.0121 0.9994+/-0.0121 0.9986+/-0.0122
0.9973+/-0.0124 0.9971+/-0.0123 0.9961+/-0.0125 0.9908+/-0.0116
0.9905+/-0.0117 0.9892+/-0.0117 0.9885+/-0.0118 0.9898+/-0.0117]
[0.9973+/-0.012 0.9994+/-0.0121 1.0+/-0.0122 0.9997+/-0.0123
0.999+/-0.0125 0.9989+/-0.0124 0.9983+/-0.0126 0.9876+/-0.0116
0.988+/-0.0118 0.9872+/-0.0118 0.9868+/-0.0118 0.9858+/-0.0117]
[0.9957+/-0.012 0.9986+/-0.0122 0.9997+/-0.0123 1.0+/-0.0124
0.9997+/-0.0126 0.9996+/-0.0125 0.9992+/-0.0127 0.985+/-0.0117
0.986+/-0.0118 0.9855+/-0.0118 0.9852+/-0.0119 0.9827+/-0.0117]
[0.9935+/-0.0122 0.9973+/-0.0124 0.999+/-0.0125 0.9997+/-0.0126
1.0+/-0.0128 0.9999+/-0.0127 0.9998+/-0.0129 0.982+/-0.0118
0.9835+/-0.012 0.9833+/-0.012 0.9833+/-0.012 0.9791+/-0.0119]
[0.9931+/-0.0121 0.9971+/-0.0123 0.9989+/-0.0124 0.9996+/-0.0125
0.9999+/-0.0127 1.0+/-0.0127 0.9999+/-0.0129 0.9813+/-0.0118
0.983+/-0.0119 0.9829+/-0.0119 0.9828+/-0.012 0.9785+/-0.0118]
[0.9918+/-0.0123 0.9961+/-0.0125 0.9983+/-0.0126 0.9992+/-0.0127
0.9998+/-0.0129 0.9999+/-0.0129 1.0+/-0.013 0.9797+/-0.012
0.9817+/-0.0121 0.9817+/-0.0121 0.9817+/-0.0122 0.9765+/-0.012]
[0.9944+/-0.0116 0.9908+/-0.0116 0.9876+/-0.0116 0.985+/-0.0117
0.982+/-0.0118 0.9813+/-0.0118 0.9797+/-0.012 1.0+/-0.0117
0.9993+/-0.0118 0.9981+/-0.0118 0.9974+/-0.0118 0.9988+/-0.0116]
[0.9931+/-0.0117 0.9905+/-0.0117 0.988+/-0.0118 0.986+/-0.0118
0.9835+/-0.012 0.983+/-0.0119 0.9817+/-0.0121 0.9993+/-0.0118
1.0+/-0.0119 0.9996+/-0.0119 0.9992+/-0.012 0.9971+/-0.0117]
[0.9914+/-0.0116 0.9892+/-0.0117 0.9872+/-0.0118 0.9855+/-0.0118
0.9833+/-0.012 0.9829+/-0.0119 0.9817+/-0.0121 0.9981+/-0.0118
0.9996+/-0.0119 1.0+/-0.012 0.9999+/-0.012 0.9952+/-0.0117]
[0.9904+/-0.0117 0.9885+/-0.0118 0.9868+/-0.0118 0.9852+/-0.0119
0.9833+/-0.012 0.9828+/-0.012 0.9817+/-0.0122 0.9974+/-0.0118
0.9992+/-0.012 0.9999+/-0.012 1.0+/-0.0121 0.9941+/-0.0118]
[0.994+/-0.0116 0.9898+/-0.0117 0.9858+/-0.0117 0.9827+/-0.0117
0.9791+/-0.0119 0.9785+/-0.0118 0.9765+/-0.012 0.9988+/-0.0116
0.9971+/-0.0117 0.9952+/-0.0117 0.9941+/-0.0118 1.0+/-0.0119]]
E_capture:
[[1.0+/-1.5451 0.9865+/-1.3235 0.9541+/-1.1437 0.9192+/-1.0314
0.8608+/-0.9172 0.8444+/-0.8909 0.82+/-0.86 0.9642+/-1.3939
0.9406+/-1.242 0.8989+/-1.1248 0.8763+/-1.0677 0.9488+/-1.4066]
[0.9865+/-1.3235 1.0+/-1.1314 0.99+/-0.9749 0.9708+/-0.8788
0.9313+/-0.7832 0.919+/-0.7613 0.9004+/-0.737 0.9178+/-1.1853
0.8973+/-1.0544 0.8588+/-0.952 0.8375+/-0.9027 0.9012+/-1.1939]
[0.9541+/-1.1437 0.99+/-0.9749 1.0+/-0.8382 0.9948+/-0.7549
0.9728+/-0.6741 0.9647+/-0.6557 0.9517+/-0.6364 0.8591+/-1.0202
0.8415+/-0.9061 0.8064+/-0.8159 0.7867+/-0.7728 0.8416+/-1.0258]
[0.9192+/-1.0314 0.9708+/-0.8788 0.9948+/-0.7549 1.0+/-0.6803
0.9911+/-0.608 0.9862+/-0.5915 0.9776+/-0.5747 0.8062+/-0.9179
0.7907+/-0.8152 0.7585+/-0.7333 0.7403+/-0.6944 0.789+/-0.9221]
[0.8608+/-0.9172 0.9313+/-0.7832 0.9728+/-0.6741 0.9911+/-0.608
1.0+/-0.5439 0.9994+/-0.529 0.9968+/-0.5138 0.7277+/-0.817
0.7152+/-0.7266 0.6872+/-0.6546 0.671+/-0.6202 0.7111+/-0.8204]
[0.8444+/-0.8909 0.919+/-0.7613 0.9647+/-0.6557 0.9862+/-0.5915
0.9994+/-0.529 1.0+/-0.5147 0.9989+/-0.4997 0.7071+/-0.7939
0.6953+/-0.7064 0.6683+/-0.6367 0.6526+/-0.6033 0.6906+/-0.7973]
[0.82+/-0.86 0.9004+/-0.737 0.9517+/-0.6364 0.9776+/-0.5747
0.9968+/-0.5138 0.9989+/-0.4997 1.0+/-0.4849 0.6767+/-0.7677
0.6659+/-0.6844 0.6405+/-0.6181 0.6255+/-0.5863 0.6608+/-0.7711]
[0.9642+/-1.3939 0.9178+/-1.1853 0.8591+/-1.0202 0.8062+/-0.9179
0.7277+/-0.817 0.7071+/-0.7939 0.6767+/-0.7677 1.0+/-1.3006
0.9921+/-1.1596 0.965+/-1.0489 0.9473+/-0.9948 0.9681+/-1.2998]
[0.9406+/-1.242 0.8973+/-1.0544 0.8415+/-0.9061 0.7907+/-0.8152
0.7152+/-0.7266 0.6953+/-0.7064 0.6659+/-0.6844 0.9921+/-1.1596
1.0+/-1.0347 0.9897+/-0.9348 0.9786+/-0.8862 0.9481+/-1.1556]
[0.8989+/-1.1248 0.8588+/-0.952 0.8064+/-0.8159 0.7585+/-0.7333
0.6872+/-0.6546 0.6683+/-0.6367 0.6405+/-0.6181 0.965+/-1.0489
0.9897+/-0.9348 1.0+/-0.8435 0.9978+/-0.799 0.9104+/-1.0439]
[0.8763+/-1.0677 0.8375+/-0.9027 0.7867+/-0.7728 0.7403+/-0.6944
0.671+/-0.6202 0.6526+/-0.6033 0.6255+/-0.5863 0.9473+/-0.9948
0.9786+/-0.8862 0.9978+/-0.799 1.0+/-0.757 0.8891+/-0.9899]
[0.9488+/-1.4066 0.9012+/-1.1939 0.8416+/-1.0258 0.789+/-0.9221
0.7111+/-0.8204 0.6906+/-0.7973 0.6608+/-0.7711 0.9681+/-1.2998
0.9481+/-1.1556 0.9104+/-1.0439 0.8891+/-0.9899 1.0+/-1.3328]]
E_scatter:
[[1.0+/-0.496 0.9869+/-0.4988 0.9784+/-0.5328 0.9694+/-0.565
0.958+/-0.6741 0.9534+/-0.7057 0.9495+/-0.78 0.2886+/-0.6303
0.2551+/-0.5931 0.2396+/-0.6574 0.2461+/-0.6978 0.2374+/-0.5398]
[0.9869+/-0.4988 1.0+/-0.4947 0.9947+/-0.5223 0.9899+/-0.5489
0.982+/-0.6464 0.9791+/-0.6748 0.9755+/-0.7424 0.22+/-0.6316
0.1944+/-0.5943 0.1826+/-0.6516 0.1872+/-0.689 0.1812+/-0.5362]
[0.9784+/-0.5328 0.9947+/-0.5223 1.0+/-0.5465 0.9942+/-0.5697
0.9868+/-0.663 0.9871+/-0.6904 0.9827+/-0.7561 0.2029+/-0.6742
0.1821+/-0.6344 0.1685+/-0.6897 0.1735+/-0.7271 0.166+/-0.5684]
[0.9694+/-0.565 0.9899+/-0.5489 0.9942+/-0.5697 1.0+/-0.59
0.9917+/-0.6794 0.9917+/-0.7058 0.9872+/-0.7697 0.1851+/-0.7145
0.166+/-0.6723 0.1545+/-0.7261 0.1587+/-0.7636 0.152+/-0.5991]
[0.958+/-0.6741 0.982+/-0.6464 0.9868+/-0.663 0.9917+/-0.6794
1.0+/-0.7701 0.9921+/-0.7967 0.9926+/-0.8628 0.183+/-0.8533
0.164+/-0.8027 0.1524+/-0.8593 0.1575+/-0.9006 0.1507+/-0.7107]
[0.9534+/-0.7057 0.9791+/-0.6748 0.9871+/-0.6904 0.9917+/-0.7058
0.9921+/-0.7967 1.0+/-0.8239 0.992+/-0.8905 0.1767+/-0.8935
0.1589+/-0.8406 0.1467+/-0.898 0.1508+/-0.9405 0.1436+/-0.7432]
[0.9495+/-0.78 0.9755+/-0.7424 0.9827+/-0.7561 0.9872+/-0.7697
0.9926+/-0.8628 0.992+/-0.8905 1.0+/-0.96 0.1823+/-0.9886 0.1628+/-0.93
0.1524+/-0.9904 0.1557+/-1.0358 0.1479+/-0.8201]
[0.2886+/-0.6303 0.22+/-0.6316 0.2029+/-0.6742 0.1851+/-0.7145
0.183+/-0.8533 0.1767+/-0.8935 0.1823+/-0.9886 1.0+/-0.8296
0.9866+/-0.7826 0.9807+/-0.8688 0.9731+/-0.9224 0.2956+/-0.7044]
[0.2551+/-0.5931 0.1944+/-0.5943 0.1821+/-0.6344 0.166+/-0.6723
0.164+/-0.8027 0.1589+/-0.8406 0.1628+/-0.93 0.9866+/-0.7826
1.0+/-0.7398 0.9847+/-0.8214 0.9791+/-0.8721 0.2582+/-0.6624]
[0.2396+/-0.6574 0.1826+/-0.6516 0.1685+/-0.6897 0.1545+/-0.7261
0.1524+/-0.8593 0.1467+/-0.898 0.1524+/-0.9904 0.9807+/-0.8688
0.9847+/-0.8214 1.0+/-0.9062 0.9828+/-0.9596 0.2471+/-0.731]
[0.2461+/-0.6978 0.1872+/-0.689 0.1735+/-0.7271 0.1587+/-0.7636
0.1575+/-0.9006 0.1508+/-0.9405 0.1557+/-1.0358 0.9731+/-0.9224
0.9791+/-0.8721 0.9828+/-0.9596 1.0+/-1.016 0.2526+/-0.7747]
[0.2374+/-0.5398 0.1812+/-0.5362 0.166+/-0.5684 0.152+/-0.5991
0.1507+/-0.7107 0.1436+/-0.7432 0.1479+/-0.8201 0.2956+/-0.7044
0.2582+/-0.6624 0.2471+/-0.731 0.2526+/-0.7747 1.0+/-0.6116]]
From the output, we can see that the integral indices are stored in a dictionary where the keys are the names of the integral
indices and the values are the matrices of the integral indices. The integral indices are stored as numpy arrays, and have shapes
(num_applications, num_experiments)
.
Getting Integral Indices From a Set of Application and Experiment SDFs
Another way to get the integral indices from a given set of applications and experiments that requires less manual effort is to use
the tsunami_ip_utils.integral_indices.get_integral_indices()
function. This function takes a list of application SDFs and
experiment SDFs, and the covariance library (which the user would have to manually input when creating their TSUNAMI-IP input
file in the previous section anyways). This function just creates a TSUNAMI-IP input file, runs it, then extracts the results, simply
abstracting away the manual effort required in the previous method.
from tsunami_ip_utils.integral_indices import get_integral_indices
application_sdfs = [ EXAMPLES / 'data' / 'example_sdfs' / 'HMF' / f'HEU-MET-FAST-003-00{i}.sdf' for i in range(1, 13) ]
experiment_sdfs = application_sdfs
coverx_library = '252groupcov7.1'
integral_indices = get_integral_indices(application_sdfs, experiment_sdfs, coverx_library=coverx_library)
# Print the integral indices
for key, value in integral_indices.items():
print(f'{key}:\n{value}\n')
c_k:
[[1.0+/-0.0024 0.9986+/-0.0024 0.9941+/-0.0025 0.9892+/-0.0025
0.9729+/-0.0026 0.9701+/-0.0026 0.9582+/-0.0026 0.9032+/-0.0009
0.896+/-0.0008 0.8902+/-0.0008 0.8883+/-0.0007 0.8353+/-0.0008]
[0.9986+/-0.0024 1.0+/-0.0025 0.9976+/-0.0025 0.9939+/-0.0026
0.9798+/-0.0026 0.9772+/-0.0026 0.9663+/-0.0026 0.9036+/-0.0009
0.898+/-0.0008 0.8939+/-0.0008 0.8924+/-0.0007 0.8347+/-0.0008]
[0.9941+/-0.0025 0.9976+/-0.0025 1.0+/-0.0026 0.9991+/-0.0026
0.9909+/-0.0025 0.9891+/-0.0026 0.9811+/-0.0026 0.9192+/-0.0008
0.9145+/-0.0008 0.9115+/-0.0008 0.9102+/-0.0007 0.8491+/-0.0008]
[0.9892+/-0.0025 0.9939+/-0.0026 0.9991+/-0.0026 1.0+/-0.0025
0.9956+/-0.0025 0.9943+/-0.0025 0.9882+/-0.0025 0.927+/-0.0008
0.9229+/-0.0008 0.9204+/-0.0007 0.9192+/-0.0007 0.8564+/-0.0008]
[0.9729+/-0.0026 0.9798+/-0.0026 0.9909+/-0.0025 0.9956+/-0.0025
1.0+/-0.0024 0.9998+/-0.0024 0.9981+/-0.0024 0.9428+/-0.0008
0.9397+/-0.0008 0.9386+/-0.0007 0.9374+/-0.0007 0.871+/-0.0008]
[0.9701+/-0.0026 0.9772+/-0.0026 0.9891+/-0.0026 0.9943+/-0.0025
0.9998+/-0.0024 1.0+/-0.0024 0.9988+/-0.0024 0.9437+/-0.0008
0.9406+/-0.0008 0.9396+/-0.0007 0.9385+/-0.0007 0.8717+/-0.0008]
[0.9582+/-0.0026 0.9663+/-0.0026 0.9811+/-0.0026 0.9882+/-0.0025
0.9981+/-0.0024 0.9988+/-0.0024 1.0+/-0.0024 0.9482+/-0.0008
0.9455+/-0.0008 0.945+/-0.0007 0.9439+/-0.0007 0.8761+/-0.0008]
[0.9032+/-0.001 0.9036+/-0.001 0.9192+/-0.001 0.927+/-0.001
0.9428+/-0.0009 0.9437+/-0.0009 0.9482+/-0.0009 1.0+/-0.001
0.9988+/-0.0009 0.996+/-0.0009 0.9945+/-0.0009 0.9226+/-0.0008]
[0.896+/-0.001 0.898+/-0.001 0.9145+/-0.0009 0.9229+/-0.0009
0.9397+/-0.0009 0.9406+/-0.0009 0.9455+/-0.0009 0.9988+/-0.0009
1.0+/-0.0009 0.999+/-0.0009 0.9982+/-0.0008 0.92+/-0.0008]
[0.8902+/-0.0009 0.8939+/-0.0009 0.9115+/-0.0009 0.9204+/-0.0009
0.9386+/-0.0008 0.9396+/-0.0008 0.945+/-0.0008 0.996+/-0.0009
0.999+/-0.0009 1.0+/-0.0008 0.9998+/-0.0008 0.9169+/-0.0008]
[0.8883+/-0.0009 0.8924+/-0.0009 0.9102+/-0.0009 0.9192+/-0.0009
0.9374+/-0.0008 0.9385+/-0.0008 0.9439+/-0.0008 0.9945+/-0.0009
0.9982+/-0.0008 0.9998+/-0.0008 1.0+/-0.0008 0.9157+/-0.0008]
[0.8353+/-0.001 0.8347+/-0.001 0.8491+/-0.001 0.8564+/-0.0009
0.871+/-0.0009 0.8717+/-0.0009 0.8761+/-0.0009 0.9226+/-0.0008
0.92+/-0.0007 0.9169+/-0.0007 0.9157+/-0.0007 1.0+/-0.0011]]
E_total:
[[1.0+/-0.0024 0.9986+/-0.0024 0.9965+/-0.0024 0.9943+/-0.0024
0.9916+/-0.0024 0.9911+/-0.0024 0.9896+/-0.0024 0.9825+/-0.0023
0.9792+/-0.0023 0.9768+/-0.0023 0.9758+/-0.0023 0.9789+/-0.0023]
[0.9986+/-0.0024 1.0+/-0.0024 0.9992+/-0.0024 0.9981+/-0.0024
0.9964+/-0.0024 0.996+/-0.0024 0.9949+/-0.0025 0.9767+/-0.0023
0.9745+/-0.0023 0.9726+/-0.0023 0.972+/-0.0023 0.9721+/-0.0023]
[0.9965+/-0.0024 0.9992+/-0.0024 1.0+/-0.0024 0.9995+/-0.0024
0.9985+/-0.0024 0.9983+/-0.0024 0.9974+/-0.0025 0.9729+/-0.0023
0.9717+/-0.0023 0.9703+/-0.0023 0.9699+/-0.0023 0.9675+/-0.0023]
[0.9943+/-0.0024 0.9981+/-0.0024 0.9995+/-0.0024 1.0+/-0.0024
0.9994+/-0.0025 0.9992+/-0.0025 0.9986+/-0.0025 0.9694+/-0.0023
0.9689+/-0.0023 0.9678+/-0.0023 0.9676+/-0.0023 0.9634+/-0.0023]
[0.9916+/-0.0024 0.9964+/-0.0024 0.9985+/-0.0024 0.9994+/-0.0025
1.0+/-0.0025 0.9997+/-0.0025 0.9996+/-0.0025 0.967+/-0.0023
0.9671+/-0.0023 0.9664+/-0.0023 0.9664+/-0.0023 0.9604+/-0.0023]
[0.9911+/-0.0024 0.996+/-0.0024 0.9983+/-0.0024 0.9992+/-0.0025
0.9997+/-0.0025 1.0+/-0.0025 0.9996+/-0.0025 0.9666+/-0.0023
0.9668+/-0.0023 0.9662+/-0.0023 0.9662+/-0.0023 0.9599+/-0.0023]
[0.9896+/-0.0024 0.9949+/-0.0025 0.9974+/-0.0025 0.9986+/-0.0025
0.9996+/-0.0025 0.9996+/-0.0025 1.0+/-0.0026 0.9655+/-0.0024
0.9661+/-0.0024 0.9657+/-0.0023 0.9658+/-0.0023 0.9585+/-0.0023]
[0.9825+/-0.0023 0.9767+/-0.0023 0.9729+/-0.0023 0.9694+/-0.0023
0.967+/-0.0023 0.9666+/-0.0023 0.9655+/-0.0024 1.0+/-0.0023
0.9988+/-0.0023 0.9971+/-0.0023 0.9961+/-0.0023 0.9871+/-0.0023]
[0.9792+/-0.0023 0.9745+/-0.0023 0.9717+/-0.0023 0.9689+/-0.0023
0.9671+/-0.0023 0.9668+/-0.0023 0.9661+/-0.0024 0.9988+/-0.0023
1.0+/-0.0023 0.9992+/-0.0023 0.9986+/-0.0023 0.9836+/-0.0023]
[0.9768+/-0.0023 0.9726+/-0.0023 0.9703+/-0.0023 0.9678+/-0.0023
0.9664+/-0.0023 0.9662+/-0.0023 0.9657+/-0.0023 0.9971+/-0.0023
0.9992+/-0.0023 1.0+/-0.0023 0.9996+/-0.0023 0.9814+/-0.0023]
[0.9758+/-0.0023 0.972+/-0.0023 0.9699+/-0.0023 0.9676+/-0.0023
0.9664+/-0.0023 0.9662+/-0.0023 0.9658+/-0.0023 0.9961+/-0.0023
0.9986+/-0.0023 0.9996+/-0.0023 1.0+/-0.0023 0.9805+/-0.0023]
[0.9789+/-0.0023 0.9721+/-0.0023 0.9675+/-0.0023 0.9634+/-0.0023
0.9604+/-0.0023 0.9599+/-0.0023 0.9585+/-0.0023 0.9871+/-0.0023
0.9836+/-0.0023 0.9814+/-0.0023 0.9805+/-0.0023 1.0+/-0.0023]]
E_fission:
[[1.0+/-0.0118 0.9991+/-0.0119 0.9973+/-0.012 0.9957+/-0.012
0.9935+/-0.0122 0.9931+/-0.0121 0.9918+/-0.0123 0.9944+/-0.0116
0.9931+/-0.0117 0.9914+/-0.0116 0.9904+/-0.0117 0.994+/-0.0116]
[0.9991+/-0.0119 1.0+/-0.0121 0.9994+/-0.0121 0.9986+/-0.0122
0.9973+/-0.0124 0.9971+/-0.0123 0.9961+/-0.0125 0.9908+/-0.0116
0.9905+/-0.0117 0.9892+/-0.0117 0.9885+/-0.0118 0.9898+/-0.0117]
[0.9973+/-0.012 0.9994+/-0.0121 1.0+/-0.0122 0.9997+/-0.0123
0.999+/-0.0125 0.9989+/-0.0124 0.9983+/-0.0126 0.9876+/-0.0116
0.988+/-0.0118 0.9872+/-0.0118 0.9868+/-0.0118 0.9858+/-0.0117]
[0.9957+/-0.012 0.9986+/-0.0122 0.9997+/-0.0123 1.0+/-0.0124
0.9997+/-0.0126 0.9996+/-0.0125 0.9992+/-0.0127 0.985+/-0.0117
0.986+/-0.0118 0.9855+/-0.0118 0.9852+/-0.0119 0.9827+/-0.0117]
[0.9935+/-0.0122 0.9973+/-0.0124 0.999+/-0.0125 0.9997+/-0.0126
1.0+/-0.0128 0.9999+/-0.0127 0.9998+/-0.0129 0.982+/-0.0118
0.9835+/-0.012 0.9833+/-0.012 0.9833+/-0.012 0.9791+/-0.0119]
[0.9931+/-0.0121 0.9971+/-0.0123 0.9989+/-0.0124 0.9996+/-0.0125
0.9999+/-0.0127 1.0+/-0.0127 0.9999+/-0.0129 0.9813+/-0.0118
0.983+/-0.0119 0.9829+/-0.0119 0.9828+/-0.012 0.9785+/-0.0118]
[0.9918+/-0.0123 0.9961+/-0.0125 0.9983+/-0.0126 0.9992+/-0.0127
0.9998+/-0.0129 0.9999+/-0.0129 1.0+/-0.013 0.9797+/-0.012
0.9817+/-0.0121 0.9817+/-0.0121 0.9817+/-0.0122 0.9765+/-0.012]
[0.9944+/-0.0116 0.9908+/-0.0116 0.9876+/-0.0116 0.985+/-0.0117
0.982+/-0.0118 0.9813+/-0.0118 0.9797+/-0.012 1.0+/-0.0117
0.9993+/-0.0118 0.9981+/-0.0118 0.9974+/-0.0118 0.9988+/-0.0116]
[0.9931+/-0.0117 0.9905+/-0.0117 0.988+/-0.0118 0.986+/-0.0118
0.9835+/-0.012 0.983+/-0.0119 0.9817+/-0.0121 0.9993+/-0.0118
1.0+/-0.0119 0.9996+/-0.0119 0.9992+/-0.012 0.9971+/-0.0117]
[0.9914+/-0.0116 0.9892+/-0.0117 0.9872+/-0.0118 0.9855+/-0.0118
0.9833+/-0.012 0.9829+/-0.0119 0.9817+/-0.0121 0.9981+/-0.0118
0.9996+/-0.0119 1.0+/-0.012 0.9999+/-0.012 0.9952+/-0.0117]
[0.9904+/-0.0117 0.9885+/-0.0118 0.9868+/-0.0118 0.9852+/-0.0119
0.9833+/-0.012 0.9828+/-0.012 0.9817+/-0.0122 0.9974+/-0.0118
0.9992+/-0.012 0.9999+/-0.012 1.0+/-0.0121 0.9941+/-0.0118]
[0.994+/-0.0116 0.9898+/-0.0117 0.9858+/-0.0117 0.9827+/-0.0117
0.9791+/-0.0119 0.9785+/-0.0118 0.9765+/-0.012 0.9988+/-0.0116
0.9971+/-0.0117 0.9952+/-0.0117 0.9941+/-0.0118 1.0+/-0.0119]]
E_capture:
[[1.0+/-1.5451 0.9865+/-1.3235 0.9541+/-1.1437 0.9192+/-1.0314
0.8608+/-0.9172 0.8444+/-0.8909 0.82+/-0.86 0.9642+/-1.3939
0.9406+/-1.242 0.8989+/-1.1248 0.8763+/-1.0677 0.9488+/-1.4066]
[0.9865+/-1.3235 1.0+/-1.1314 0.99+/-0.9749 0.9708+/-0.8788
0.9313+/-0.7832 0.919+/-0.7613 0.9004+/-0.737 0.9178+/-1.1853
0.8973+/-1.0544 0.8588+/-0.952 0.8375+/-0.9027 0.9012+/-1.1939]
[0.9541+/-1.1437 0.99+/-0.9749 1.0+/-0.8382 0.9948+/-0.7549
0.9728+/-0.6741 0.9647+/-0.6557 0.9517+/-0.6364 0.8591+/-1.0202
0.8415+/-0.9061 0.8064+/-0.8159 0.7867+/-0.7728 0.8416+/-1.0258]
[0.9192+/-1.0314 0.9708+/-0.8788 0.9948+/-0.7549 1.0+/-0.6803
0.9911+/-0.608 0.9862+/-0.5915 0.9776+/-0.5747 0.8062+/-0.9179
0.7907+/-0.8152 0.7585+/-0.7333 0.7403+/-0.6944 0.789+/-0.9221]
[0.8608+/-0.9172 0.9313+/-0.7832 0.9728+/-0.6741 0.9911+/-0.608
1.0+/-0.5439 0.9994+/-0.529 0.9968+/-0.5138 0.7277+/-0.817
0.7152+/-0.7266 0.6872+/-0.6546 0.671+/-0.6202 0.7111+/-0.8204]
[0.8444+/-0.8909 0.919+/-0.7613 0.9647+/-0.6557 0.9862+/-0.5915
0.9994+/-0.529 1.0+/-0.5147 0.9989+/-0.4997 0.7071+/-0.7939
0.6953+/-0.7064 0.6683+/-0.6367 0.6526+/-0.6033 0.6906+/-0.7973]
[0.82+/-0.86 0.9004+/-0.737 0.9517+/-0.6364 0.9776+/-0.5747
0.9968+/-0.5138 0.9989+/-0.4997 1.0+/-0.4849 0.6767+/-0.7677
0.6659+/-0.6844 0.6405+/-0.6181 0.6255+/-0.5863 0.6608+/-0.7711]
[0.9642+/-1.3939 0.9178+/-1.1853 0.8591+/-1.0202 0.8062+/-0.9179
0.7277+/-0.817 0.7071+/-0.7939 0.6767+/-0.7677 1.0+/-1.3006
0.9921+/-1.1596 0.965+/-1.0489 0.9473+/-0.9948 0.9681+/-1.2998]
[0.9406+/-1.242 0.8973+/-1.0544 0.8415+/-0.9061 0.7907+/-0.8152
0.7152+/-0.7266 0.6953+/-0.7064 0.6659+/-0.6844 0.9921+/-1.1596
1.0+/-1.0347 0.9897+/-0.9348 0.9786+/-0.8862 0.9481+/-1.1556]
[0.8989+/-1.1248 0.8588+/-0.952 0.8064+/-0.8159 0.7585+/-0.7333
0.6872+/-0.6546 0.6683+/-0.6367 0.6405+/-0.6181 0.965+/-1.0489
0.9897+/-0.9348 1.0+/-0.8435 0.9978+/-0.799 0.9104+/-1.0439]
[0.8763+/-1.0677 0.8375+/-0.9027 0.7867+/-0.7728 0.7403+/-0.6944
0.671+/-0.6202 0.6526+/-0.6033 0.6255+/-0.5863 0.9473+/-0.9948
0.9786+/-0.8862 0.9978+/-0.799 1.0+/-0.757 0.8891+/-0.9899]
[0.9488+/-1.4066 0.9012+/-1.1939 0.8416+/-1.0258 0.789+/-0.9221
0.7111+/-0.8204 0.6906+/-0.7973 0.6608+/-0.7711 0.9681+/-1.2998
0.9481+/-1.1556 0.9104+/-1.0439 0.8891+/-0.9899 1.0+/-1.3328]]
E_scatter:
[[1.0+/-0.496 0.9869+/-0.4988 0.9784+/-0.5328 0.9694+/-0.565
0.958+/-0.6741 0.9534+/-0.7057 0.9495+/-0.78 0.2886+/-0.6303
0.2551+/-0.5931 0.2396+/-0.6574 0.2461+/-0.6978 0.2374+/-0.5398]
[0.9869+/-0.4988 1.0+/-0.4947 0.9947+/-0.5223 0.9899+/-0.5489
0.982+/-0.6464 0.9791+/-0.6748 0.9755+/-0.7424 0.22+/-0.6316
0.1944+/-0.5943 0.1826+/-0.6516 0.1872+/-0.689 0.1812+/-0.5362]
[0.9784+/-0.5328 0.9947+/-0.5223 1.0+/-0.5465 0.9942+/-0.5697
0.9868+/-0.663 0.9871+/-0.6904 0.9827+/-0.7561 0.2029+/-0.6742
0.1821+/-0.6344 0.1685+/-0.6897 0.1735+/-0.7271 0.166+/-0.5684]
[0.9694+/-0.565 0.9899+/-0.5489 0.9942+/-0.5697 1.0+/-0.59
0.9917+/-0.6794 0.9917+/-0.7058 0.9872+/-0.7697 0.1851+/-0.7145
0.166+/-0.6723 0.1545+/-0.7261 0.1587+/-0.7636 0.152+/-0.5991]
[0.958+/-0.6741 0.982+/-0.6464 0.9868+/-0.663 0.9917+/-0.6794
1.0+/-0.7701 0.9921+/-0.7967 0.9926+/-0.8628 0.183+/-0.8533
0.164+/-0.8027 0.1524+/-0.8593 0.1575+/-0.9006 0.1507+/-0.7107]
[0.9534+/-0.7057 0.9791+/-0.6748 0.9871+/-0.6904 0.9917+/-0.7058
0.9921+/-0.7967 1.0+/-0.8239 0.992+/-0.8905 0.1767+/-0.8935
0.1589+/-0.8406 0.1467+/-0.898 0.1508+/-0.9405 0.1436+/-0.7432]
[0.9495+/-0.78 0.9755+/-0.7424 0.9827+/-0.7561 0.9872+/-0.7697
0.9926+/-0.8628 0.992+/-0.8905 1.0+/-0.96 0.1823+/-0.9886 0.1628+/-0.93
0.1524+/-0.9904 0.1557+/-1.0358 0.1479+/-0.8201]
[0.2886+/-0.6303 0.22+/-0.6316 0.2029+/-0.6742 0.1851+/-0.7145
0.183+/-0.8533 0.1767+/-0.8935 0.1823+/-0.9886 1.0+/-0.8296
0.9866+/-0.7826 0.9807+/-0.8688 0.9731+/-0.9224 0.2956+/-0.7044]
[0.2551+/-0.5931 0.1944+/-0.5943 0.1821+/-0.6344 0.166+/-0.6723
0.164+/-0.8027 0.1589+/-0.8406 0.1628+/-0.93 0.9866+/-0.7826
1.0+/-0.7398 0.9847+/-0.8214 0.9791+/-0.8721 0.2582+/-0.6624]
[0.2396+/-0.6574 0.1826+/-0.6516 0.1685+/-0.6897 0.1545+/-0.7261
0.1524+/-0.8593 0.1467+/-0.898 0.1524+/-0.9904 0.9807+/-0.8688
0.9847+/-0.8214 1.0+/-0.9062 0.9828+/-0.9596 0.2471+/-0.731]
[0.2461+/-0.6978 0.1872+/-0.689 0.1735+/-0.7271 0.1587+/-0.7636
0.1575+/-0.9006 0.1508+/-0.9405 0.1557+/-1.0358 0.9731+/-0.9224
0.9791+/-0.8721 0.9828+/-0.9596 1.0+/-1.016 0.2526+/-0.7747]
[0.2374+/-0.5398 0.1812+/-0.5362 0.166+/-0.5684 0.152+/-0.5991
0.1507+/-0.7107 0.1436+/-0.7432 0.1479+/-0.8201 0.2956+/-0.7044
0.2582+/-0.6624 0.2471+/-0.731 0.2526+/-0.7747 1.0+/-0.6116]]
Spaces in Annotated Names
In each SDF file, the first line contains an annotated name for the SDF file, that is used to identify the SDF in the TSUNAMI-IP
output tables. If these names contain spaces (as some of those provided in the MCT directory do), they can cause the
tsnami_ip_utils.readers.read_integral_indices()
function to fail. For example, the following code will raise an error:
try:
integral_indices = read_integral_indices(EXAMPLES / 'data' / 'tsunami_ip_mct.out')
print(integral_indices)
except Exception as e:
print(e)
No integral values found in the output file. Please ensure that the output/TSUNAMI-IP input files are formatted correctly.
This is because of the spaces in the following SDF annotated names (that were used to generate the TSUNAMI-IP output file):
mixed-comp-therm-002-001s (experiment pnl-30)
mixed-comp-therm-002-002s (experiment pnl-31)
mixed-comp-therm-002-003s (experiment pnl-32)
mixed-comp-therm-002-004s (experiment pnl-33)
mixed-comp-therm-002-005s (experiment pnl-34)
mixed-comp-therm-002-006s (experiment pnl-35)
So if you have SDFs with names that contain spaces, you should either use the
tsunami_ip_utils.integral_indices.get_integral_indices()
, which will handle these names by replacing whitespace with '-'
,
or you can manually replace the whitespaces in the SDF annotated names using the tsunami_ip_utils.utils.modify_sdf_names()
function. Below is an example showing the correct output using the tsunami_ip_utils.integral_indices.get_integral_indices()
function:
from tsunami_ip_utils.integral_indices import get_integral_indices
application_sdfs = [ EXAMPLES / 'data' / 'example_sdfs' / 'MCT' / f'MIX-COMP-THERM-002-00{i}S.sdf' for i in range(1, 7) ]
experiment_sdfs = application_sdfs
coverx_library = '252groupcov7.1'
integral_indices = get_integral_indices(application_sdfs, experiment_sdfs, coverx_library)
# Print the integral indices
for key, value in integral_indices.items():
print(f'{key}:\n{value}\n')
c_k:
[[1.0+/-0.0005 0.9897+/-0.0005 0.9694+/-0.0004 0.9381+/-0.0014
0.9392+/-0.0003 0.8983+/-0.0003]
[0.9897+/-0.0005 1.0+/-0.0005 0.9405+/-0.0004 0.9573+/-0.0016
0.9123+/-0.0004 0.916+/-0.0003]
[0.9694+/-0.0004 0.9405+/-0.0004 1.0+/-0.0003 0.9374+/-0.0008
0.9917+/-0.0003 0.9215+/-0.0003]
[0.9381+/-0.0003 0.9573+/-0.0003 0.9374+/-0.0003 1.0+/-0.0003
0.9434+/-0.0003 0.9919+/-0.0002]
[0.9392+/-0.0003 0.9123+/-0.0004 0.9917+/-0.0003 0.9434+/-0.0007
1.0+/-0.0002 0.9431+/-0.0002]
[0.8983+/-0.0003 0.916+/-0.0003 0.9215+/-0.0003 0.9919+/-0.0009
0.9431+/-0.0002 1.0+/-0.0002]]
E_total:
[[1.0+/-0.0007 0.991+/-0.0005 0.974+/-0.0006 0.9603+/-0.0005
0.9541+/-0.0006 0.942+/-0.0005]
[0.991+/-0.0005 1.0+/-0.0004 0.9473+/-0.0004 0.967+/-0.0004
0.9217+/-0.0004 0.9408+/-0.0004]
[0.974+/-0.0006 0.9473+/-0.0004 1.0+/-0.0005 0.9634+/-0.0004
0.9966+/-0.0005 0.9668+/-0.0004]
[0.9603+/-0.0005 0.967+/-0.0004 0.9634+/-0.0004 1.0+/-0.0003
0.9544+/-0.0004 0.9948+/-0.0003]
[0.9541+/-0.0006 0.9217+/-0.0004 0.9966+/-0.0005 0.9544+/-0.0004
1.0+/-0.0005 0.9658+/-0.0004]
[0.942+/-0.0005 0.9408+/-0.0004 0.9668+/-0.0004 0.9948+/-0.0003
0.9658+/-0.0004 1.0+/-0.0004]]
E_fission:
[[1.0+/-0.0074 0.9984+/-0.0066 0.9716+/-0.0072 0.9829+/-0.0056
0.9525+/-0.0067 0.9643+/-0.0055]
[0.9984+/-0.0066 1.0+/-0.0059 0.9581+/-0.0063 0.9736+/-0.005
0.9356+/-0.0058 0.9504+/-0.0049]
[0.9716+/-0.0072 0.9581+/-0.0063 1.0+/-0.0073 0.997+/-0.0056
0.9975+/-0.007 0.999+/-0.0056]
[0.9829+/-0.0056 0.9736+/-0.005 0.997+/-0.0056 1.0+/-0.0045
0.9901+/-0.0053 0.9961+/-0.0045]
[0.9525+/-0.0067 0.9356+/-0.0058 0.9975+/-0.007 0.9901+/-0.0053
1.0+/-0.0067 0.9984+/-0.0053]
[0.9643+/-0.0055 0.9504+/-0.0049 0.999+/-0.0056 0.9961+/-0.0045
0.9984+/-0.0053 1.0+/-0.0045]]
E_capture:
[[1.0+/-0.0062 0.9286+/-0.006 0.9463+/-0.0058 0.7276+/-0.005
0.8919+/-0.0053 0.7141+/-0.0047]
[0.9286+/-0.006 1.0+/-0.006 0.8122+/-0.0055 0.8725+/-0.0049
0.7326+/-0.005 0.8197+/-0.0045]
[0.9463+/-0.0058 0.8122+/-0.0055 1.0+/-0.0057 0.7085+/-0.0047
0.9889+/-0.0053 0.7464+/-0.0045]
[0.7276+/-0.005 0.8725+/-0.0049 0.7085+/-0.0047 1.0+/-0.0042
0.6732+/-0.0044 0.9861+/-0.004]
[0.8919+/-0.0053 0.7326+/-0.005 0.9889+/-0.0053 0.6732+/-0.0044
1.0+/-0.0049 0.7341+/-0.0042]
[0.7141+/-0.0047 0.8197+/-0.0045 0.7464+/-0.0045 0.9861+/-0.004
0.7341+/-0.0042 1.0+/-0.0038]]
E_scatter:
[[1.0+/-0.169 0.983+/-0.1599 0.9461+/-0.1543 0.9375+/-0.2897
0.9113+/-0.1747 0.8837+/-0.4451]
[0.983+/-0.1599 1.0+/-0.1448 0.8962+/-0.1359 0.9311+/-0.2822
0.8507+/-0.1562 0.8563+/-0.4593]
[0.9461+/-0.1543 0.8962+/-0.1359 1.0+/-0.142 0.918+/-0.2754
0.9943+/-0.1666 0.9056+/-0.4582]
[0.9375+/-0.2897 0.9311+/-0.2822 0.918+/-0.2754 1.0+/-0.5368
0.9008+/-0.3139 0.981+/-0.8229]
[0.9113+/-0.1747 0.8507+/-0.1562 0.9943+/-0.1666 0.9008+/-0.3139
1.0+/-0.1952 0.9044+/-0.5153]
[0.8837+/-0.4451 0.8563+/-0.4593 0.9056+/-0.4582 0.981+/-0.8229
0.9044+/-0.5153 1.0+/-1.146]]
Total running time of the script: (0 minutes 51.992 seconds)