_RL mult_ies
& mult_ies,
& mult_ctdt, mult_ctds, mult_ies,
mult_ies = 0. _d 0
& + mult_ies * objf_ies(bi,bj)