add diffusion.py which calcualte relative pair wise diffusion