r_el – shape [N_batch x n_el x 3]
R_ion – shape [N_ion x 3]
shape [N_batch x n_el x N_ion x 3] dist: shape [N_batch x n_el x N_ion]
diff