Add tr_bregman which can be used to compute one component, not finished yet