_RL mult_wdrift
& mult_wdrift,
& mult_sdrift, mult_tdrift, mult_wdrift,
mult_wdrift = 0. _d 0
& + mult_wdrift * objf_wdrift(bi,bj)