This technical note is dedicated to model order reduction of linear time-invariant systems. The main contribution of this technical note is the derivation of two scalable stability-preserving model reduction algorithms. Both algorithms constitute a development of a recently proposed model reduction method. The algorithms perform a curve fitting procedure using frequency response samples of a model and semidefinite programming methods. Computation of these samples can be done efficiently even for large scale models. Both algorithms are obtained from a reformulation of the model reduction problem. One proposes a semidefinite relaxation, while the other is an iterative semidefinite approach. The relaxation approach is similar to Hankel model reduction, which is a well-known and established method in the control literature. Due to this resemblance, the accuracy of approximation is also similar to the one of Hankel model reduction. An appealing quality of the proposed algorithms is the ability to easily perform extensions, e. g., introduce frequency-weighting, positive-real and bounded-real constraints.