A variable MIP neighborhood descent for the multi-attribute inventory routing problem. (December 2020)