If you could help by elaborating a bit more on how did you calculate the denominators with `denominators = sum_of_rows — exp_sim_by_tau_diag`, that would be very helpful.