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):

sdf_paths = [ EXAMPLES / 'data' / 'example_sdfs' / 'MCT' / f'MIX-COMP-THERM-002-00{i}S.sdf' for i in range(1, 7) ]

for sdf_path in sdf_paths:
    with open(sdf_path, 'r') as f:
        print(f.readline())
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)

Gallery generated by Sphinx-Gallery