LAUSR

Abstract Spectral/pseudospectral integration preconditioning matrices are useful tools for solving differential equations involving pure differential operators d m / d x m . In this study we construct well-conditioned inverse pseudospectral matrices for the basic differential operator and the mixed differential operator d d r a ( r ) d d r on Gauss-Radau-Legendre points based on stepwise integrations and low-rank updates. The inverse matrices can be used either as a solution operator or an effective preconditioner for variable coefficient differential equations of first and second order in polar coordinates. Numerical experiments were conducted and we observed the performance of the inverse operator as expected.