Doubly stochastic matrix plays an essential role in several areas such as statistics and machine learning. In this paper we consider the optimal approximation of a square set doubly matrices. A structured BFGS method is proposed to solve dual primal problem. The resulting algorithm builds curvature information into diagonal components true Hessian, so that it takes only additional linear cost o...