finite difference scheme is exploited for the model equations. The resulting coefficient matrix
is a nonsymmetric ill-conditioned Toeplitz-like matrix. Symmetrized strategies are proposed
for the nonsymmetric system so that the conjugate gradient method can be utilized to derive
the numerical solutions. Moreover, preconditioners based on the sine transform are
designed to speed up the convergence rate of the proposed methods. Theoretically, we …