function wrapper_modifyX
input Real X_real[:, size(X_real, 1)] "Eigenvector matrix, real part";
input Real X_imag[size(X_real, 1), size(X_real, 2)] "Eigenvector matrix, imaginary part";
input Integer n "Order of X";
input Real S_real[size(X_real, 1), :] "Eigenvector bases, real part";
input Real S_imag[size(S_real, 1), size(S_real, 2)] "Eigenvector bases, imaginary part";
input Integer m "Rank of the system input matrix B; S_real and S_imag must have n*m columns";
input Integer ncp "number of complex pairs";
input Integer steps "Number of iterations";
output Real Xm_real[size(X_real, 1), size(X_real, 2)] = X_real;
output Real Xm_imag[size(X_imag, 1), size(X_imag, 2)] = X_imag;
end wrapper_modifyX;