This paper presents a primal-dual interior point method based algorithm with sparsity structure, which deals with voltage/reactive power optimization control problems in large scale power systems.