Before revision:  The method is in some sense a generalization of the Optimized Gradient Method of Kim and Fessler (2016), and asymptotically corresponds to the Triple Momentum Method (Van Scoy et al., 2017), in the presence of strong convexity. Furthermore, the method is numerically stable to arbitrarily large condition numbers and admits a conceptually very simple proof, which involves a Lyapunov argument and a sum of two inequalities. Finally, we provide a numerical recipe for obtaining the algorithmic parameters of the method , using semidefinite programming, and illustrate that it can be used for developing other methods as well.
After revision:  In addition, we provide a numerical recipe for obtaining the algorithmic parameters of the method , using semidefinite programming, and illustrate that it can be used for developing other methods as well.