fixed small bug for noise

DiCarloLab-Delft · Sep 20, 2018 · 3e43e05 · 3e43e05
1 parent 28b3154
commit 3e43e05
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/...uperoperator_simulation_withdistortions_newdevice_singlequbitphases_newcode_fluxnoise2.py b/...uperoperator_simulation_withdistortions_newdevice_singlequbitphases_newcode_fluxnoise2.py
@@ -764,7 +764,7 @@ def simulate_quantities_of_interest_superoperator(tlist, c_ops, noise_parameters
                 if isinstance(c_ops[c],list):
                     c_ops_temp.append(c_ops[c][0]*c_ops[c][1][i])    # c_ops are already in the H_0 basis
                 else:
-                    c_ops_temp.append(S_H * c_ops[c] * S_H.dag())
+                    c_ops_temp.append(c_ops[c])
             liouville_exp_t=(qtp.liouvillian(H,c_ops_temp)*sim_step).expm()
         else:
         	liouville_exp_t=(-1j*H*sim_step).expm()
@@ -1112,7 +1112,7 @@ def average_phases(phases,weights):
 	        qoi_average[5]=average*100
 	        qoi_average[9]=average_partial*100
 
-	        qoi_average[0] = -np.log10(1-qoi_average[5]/100)    # we want log of the average and not average of the log
+	        qoi_average[0] = (-np.log10(1-qoi_average[5]/100))**4    # we want log of the average and not average of the log
 
 	        qoi_plot.append(qoi_average)