The linear algebra system is solved by a combination of the preconditioned conjugate gradient method (PCG) with the fast Fourier transform (FFT), where T.