discrepancy1 = []
discrepancye1 = []
discrepancye2 = []
discrepancye3 = []
discrepancye4 = []
discrepancye5 = []
discrepancye6 = []
discrepancye7 = []
discrepancye8 = []
discrepancye9 = []
for a, b, c, d, e, f, g, h, i, j, k, s in zip(meHs10, meHs1, meHse1, meHse2,
meHse3, meHse4, meHse5, meHse6,
meHse7, meHse8, meHse9, srgHs):
discrepancy10.append(diff(a, s, 6))
discrepancy1.append(diff(b, s, 6))
discrepancye1.append(diff(c, s, 6))
discrepancye2.append(diff(d, s, 6))
discrepancye3.append(diff(e, s, 6))
discrepancye4.append(diff(f, s, 6))
discrepancye5.append(diff(g, s, 6))
discrepancye6.append(diff(h, s, 6))
discrepancye7.append(diff(i, s, 6))
discrepancye8.append(diff(j, s, 6))
discrepancye9.append(diff(k, s, 6))
#################################################
#
# Graphing
#
# Values taken from srg_pairing. fixed_d = 1.0 # Fixed value of the threshold to be used threshold_value = 0.001 # To hold the final values discrepancy_vary_g = [] # Collects the data for varying the energy spacing for g in g_range: srgHs = srg_main (flowparams, fixed_d, g) meHs = magnus_expansion_main (flowparams, threshold_value, fixed_d, g) discrepancy = [] for i, j in zip (meHs, srgHs): discrepancy.append (diff (i, j, 6)) discrepancy_vary_g.append (discrepancy) # Gathers the data interaction_negative_1 = discrepancy_vary_g[0][3:7] interaction_negative_0_2 = discrepancy_vary_g[4][3:7] interaction_0_2 = discrepancy_vary_g[6][3:7] interaction_1 = discrepancy_vary_g[10][3:7] ################################################# # # GRAPHING # ################################################# # Imports here to prevent error