There should be a multiplication, imho.
int32_t vertical_thrust = stabilization_cmd[COMMAND_THRUST] * get_vertical_thrust_coeff();
But it is still not fully correct, because we multiplying here both feedback and feedforward commands by the thrust coeff, while initially we divide only feedforward command into thrust coefficient.