Add pair_sph_rhosum_multiphase