A feasible method for optimization with orthogonality constraints
On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming