I appreciate Zu-Qing Qu for the interest shown and for making useful comments on the contents of my paper.
The reduced order matrices and retain the same form as given in Eq. (1) as the iterations progress. The computation results presented by Zu-Qing Qu seem to be incorrect.
With the help of the procedure in Sec. 3 of my paper, the computations for the first three iterations are carried out on the numerical example cited by Zu-Qing Qu. The results of the first three iterations are given in Table 1 hereunder.
In addition to the above, the eigenvalues are also extracted for the full system and the two cases of master selection. The converged eigenvalues are shown in Table 2 below.
Further to the above, I wish to add here that since the formulation in my paper finally falls into the category of unsymmetric matrices—particularly so in the case of the mass matrix in Eq. (2)—the assumption by Zu-Qing Qu that the transformation matrices [R] and [S] are the same even for a symmetric structure is not valid.
