An efficient implementation of two-component relativistic density functional theory with torque-free auxiliary variables