Diffusion tensor imaging is widely used in brain connectivity research. As more and more studies recruit large numbers of subjects, it is important to design registration methods which are not only theoretically rigorous, but also computationally efficient. However, the requirement of reorienting diffusion tensors complicates and considerably slows down registration procedures, due to the correlated impacts of registration forces at adjacent voxel locations. Based on the diffeomorphic Demons algorithm (Vercauteren , 2009), we propose a fast local trust region algorithm for handling inseparable registration forces for quadratic energy functions. The method guarantees that, at any time and at any voxel location, the velocity is always within its local trust region. This local regularization allows efficient calculation of the transformation update with numeric integration instead of completely solving a large linear system at every iteration. It is able to incorporate exact reorientation and regularization into the velocity optimization, and preserve the linear complexity of the diffeomorphic Demons algorithm. In an experiment with 84 diffusion tensor images involving both pair-wise and group-wise registrations, the proposed algorithm achieves better registration in comparison with other methods solving large linear systems (Yeo , 2009). At the same time, this algorithm reduces the computation time and memory demand tenfold.