A variable MIP neighborhood descent for the multi-attribute inventory routing problem