def calculate_observables_rep(data_out, Tr):
pred_samp = (data_out > np.random.random(data_out.shape)).astype(np.int)
obs = np.zeros([3, 7])
obs[0] = get_observables(data_out, Tr)
obs[1] = get_observables(np.round(data_out), Tr)
obs[2] = get_observables(pred_samp, Tr)
return obs
def calculate_observables(data_or, data_in, data_out, T):
pred_samp = (data_out > np.random.random(data_out.shape)).astype(np.int)
obs = np.zeros([5, 7])
obs[0] = get_observables(data_or, T)
obs[1] = get_observables(data_in, T)
obs[2] = get_observables(data_out, T)
obs[3] = get_observables(np.round(data_out), T)
obs[4] = get_observables(pred_samp, T)
return obs