The current implementation does not allow for variable-coefficient PDOs. Discretizing with the ultraspherical spectral method yields a multiplication operator for variable coefficients which is banded, allowing us to retain an O(p^4) leaf solver in #1.
The current implementation does not allow for variable-coefficient PDOs. Discretizing with the ultraspherical spectral method yields a multiplication operator for variable coefficients which is banded, allowing us to retain an O(p^4) leaf solver in #1.