ENH: making computation of Bi Bij faster by exploiting that R_i is diagonal