Diffusion of interstitial carbon atoms in iron is the rate-limiting phenomenon of a number of phase transitions in body-centered (bec) and body-centered tetragonal (bct) phases such as ferrite and martensite. These phases being rarely stress-free and undeformed, the influence of stress/strain on the diffusivity of carbon is essential, although scarcely documented. We developed a model of carbon elastodiffusion in bct-iron. We combined anisotropie linear elasticity theory of point defects, the dilute approximation of regular solutions and the multisite model of random walk into a coherent mean-field theory. The model allows predicting the effects of composition, temperature and mechanical loading on the anisotropy of carbon diffusion. Density functional theory calculations have provided most of the materials parameters. The predictions were successfully tested against kinetic Monte Carlo simulations. Our results show that compression of the crystal increases carbon diffusivity, while tension has the opposite effect. Axial straining is accompanied by a large anisotropy of diffusion. This effect could be exploited to produce stress-controlled diffusion channeling for the engineering of anisotropie micro structures during thermal ageing of martensitic Fe-C alloys.