Gradient descent iterative method for energy flow of integrated energy system considering multiple modes of compressors. (1st March 2020)