Standard spectral methods are capable of providing very accurate approximations to well-behaved smooth functions with significantly less degrees of freedom when compared with finite difference or finite element methods (cf. [6,7,11]). However, if a function exhibits localized behaviors such as sharp interfaces, very thin internal or boundary layers, using a standard Gauss-type grid usually fail...
