def potential_numpy(cluster): d = distances_numpy(cluster) dtri = np.triu(d) energy = lj_numpy(dtri[dtri > 1e-6]).sum() return energy #### END: numpy