use interpos after avoiding NaNs, use interpos_nan to do this automatically deployed since now ne_rho does not work