Title: An efficient and high order accurate solution technique for variable coefficient elliptic partial differential equations
Abstract: For many applications in science and engineering, the ability to efficiently and accurately approximate solutions to elliptic PDEs dictates
what physical phenomena can be simulated numerically. In this talk, we present a high-order accurate discretization technique for variable coefficient PDEs with smooth coefficients. For two dimensional problems, a direct solver that is inspired by the multifrontal method is scales linearly or nearly linear with respect to the number of unknowns.
Unlike the application of multifrontal methods to classic discretization techniques, the constant prefactors do not grow with the order of the discretization. The discretization is robust even for problems with highly oscillatory solutions. For example, a problem 100 wavelengths in size can be solved to 9 digits of accuracy with 3.7 million unknowns on a desktop computer. The precomputation of the direct solver takes 6 minutes on a desktop computer. Then applying the computed solver takes 3 seconds. A parallel implementation of the solution technique reduces the precomputation time to roughly 30 seconds and halves the time it takes to apply the solver. Recent work making the solution technique viable for three dimensional problems and some applications will be presented.
If you would like to attend, please email crampersad@flatironinstitute.org for the Zoom details.