An accelerated Poisson solver based on multidomain spectral discretization

dc.citation.journalTitleBIT Numerical Mathematicsen_US
dc.contributor.authorBabb, Tracyen_US
dc.contributor.authorGillman, Adriannaen_US
dc.contributor.authorHao, Sijiaen_US
dc.contributor.authorMartinsson, Per-Gunnaren_US
dc.date.accessioned2018-09-27T17:33:25Zen_US
dc.date.available2018-09-27T17:33:25Zen_US
dc.date.issued2018en_US
dc.description.abstractThis paper presents a numerical method for variable coefficient elliptic PDEs with mostly smooth solutions on two dimensional domains. The method works best for domains that can readily be mapped onto a rectangle, or a collection of nonoverlapping rectangles. The PDE is discretized via a multi-domain spectral collocation method of high local order (order 30 and higher have been tested and work well). Local mesh refinement results in highly accurate solutions even in the presence of local irregular behavior due to corner singularities, localized loads, etc. The system of linear equations attained upon discretization is solved using a direct (as opposed to iterative) solver with O(N1.5)O(N1.5) complexity for the factorization stage and O(NlogN)O(Nlogā”N) complexity for the solve. The scheme is ideally suited for executing the elliptic solve required when parabolic problems are discretized via time-implicit techniques. In situations where the geometry remains unchanged between time-steps, very fast execution speeds are obtained since the solution operator for each implicit solve can be pre-computed.en_US
dc.identifier.citationBabb, Tracy, Gillman, Adrianna, Hao, Sijia, et al.. "An accelerated Poisson solver based on multidomain spectral discretization." <i>BIT Numerical Mathematics,</i> (2018) Springer: https://doi.org/10.1007/s10543-018-0714-0.en_US
dc.identifier.digitalBabb2018en_US
dc.identifier.doihttps://doi.org/10.1007/s10543-018-0714-0en_US
dc.identifier.urihttps://hdl.handle.net/1911/102722en_US
dc.language.isoengen_US
dc.publisherSpringeren_US
dc.rightsThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.en_US
dc.rights.urihttp://creativecommons.org/licenses/by/4.0/en_US
dc.titleAn accelerated Poisson solver based on multidomain spectral discretizationen_US
dc.typeJournal articleen_US
dc.type.dcmiTexten_US
dc.type.publicationpublisher versionen_US
Files
Original bundle
Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
Babb2018.pdf
Size:
1.44 MB
Format:
Adobe Portable Document Format