uu.seUppsala University Publications

References$(function(){PrimeFaces.cw("TieredMenu","widget_formSmash_upper_j_idt145",{id:"formSmash:upper:j_idt145",widgetVar:"widget_formSmash_upper_j_idt145",autoDisplay:true,overlay:true,my:"left top",at:"left bottom",trigger:"formSmash:upper:referencesLink",triggerEvent:"click"});}); $(function(){PrimeFaces.cw("OverlayPanel","widget_formSmash_upper_j_idt146_j_idt148",{id:"formSmash:upper:j_idt146:j_idt148",widgetVar:"widget_formSmash_upper_j_idt146_j_idt148",target:"formSmash:upper:j_idt146:permLink",showEffect:"blind",hideEffect:"fade",my:"right top",at:"right bottom",showCloseIcon:true});});

Adaptive Solvers for High-Dimensional PDE Problems on Clusters of Multicore ProcessorsPrimeFaces.cw("AccordionPanel","widget_formSmash_some",{id:"formSmash:some",widgetVar:"widget_formSmash_some",multiple:true}); PrimeFaces.cw("AccordionPanel","widget_formSmash_all",{id:"formSmash:all",widgetVar:"widget_formSmash_all",multiple:true});
function selectAll()
{
var panelSome = $(PrimeFaces.escapeClientId("formSmash:some"));
var panelAll = $(PrimeFaces.escapeClientId("formSmash:all"));
panelAll.toggle();
toggleList(panelSome.get(0).childNodes, panelAll);
toggleList(panelAll.get(0).childNodes, panelAll);
}
/*Toggling the list of authorPanel nodes according to the toggling of the closeable second panel */
function toggleList(childList, panel)
{
var panelWasOpen = (panel.get(0).style.display == 'none');
// console.log('panel was open ' + panelWasOpen);
for (var c = 0; c < childList.length; c++) {
if (childList[c].classList.contains('authorPanel')) {
clickNode(panelWasOpen, childList[c]);
}
}
}
/*nodes have styleClass ui-corner-top if they are expanded and ui-corner-all if they are collapsed */
function clickNode(collapse, child)
{
if (collapse && child.classList.contains('ui-corner-top')) {
// console.log('collapse');
child.click();
}
if (!collapse && child.classList.contains('ui-corner-all')) {
// console.log('expand');
child.click();
}
}
PrimeFaces.cw("AccordionPanel","widget_formSmash_responsibleOrgs",{id:"formSmash:responsibleOrgs",widgetVar:"widget_formSmash_responsibleOrgs",multiple:true}); 2014 (English)Doctoral thesis, comprehensive summary (Other academic)
##### Abstract [en]

##### Place, publisher, year, edition, pages

Uppsala: Acta Universitatis Upsaliensis, 2014. , 34 p.
##### Series

Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology, ISSN 1651-6214 ; 1199
##### Keyword [en]

adaptive mesh refinement, anisotropic refinement, exponential integrators, Lanczos' algorithm, hybrid parallelization, time-dependent Schrödinger equation
##### National Category

Computational Mathematics
##### Research subject

Scientific Computing
##### Identifiers

URN: urn:nbn:se:uu:diva-234984ISBN: 978-91-554-9095-9OAI: oai:DiVA.org:uu-234984DiVA: diva2:758606
##### Public defence

2014-12-12, Room 2446, Polacksbacken, Lägerhyddsvägen 2, Uppsala, 10:15 (English)
##### Opponent

PrimeFaces.cw("AccordionPanel","widget_formSmash_j_idt376",{id:"formSmash:j_idt376",widgetVar:"widget_formSmash_j_idt376",multiple:true});
##### Supervisors

PrimeFaces.cw("AccordionPanel","widget_formSmash_j_idt382",{id:"formSmash:j_idt382",widgetVar:"widget_formSmash_j_idt382",multiple:true});
#####

PrimeFaces.cw("AccordionPanel","widget_formSmash_j_idt388",{id:"formSmash:j_idt388",widgetVar:"widget_formSmash_j_idt388",multiple:true});
##### Projects

eSSENCE
Available from: 2014-11-21 Created: 2014-10-27 Last updated: 2015-02-03Bibliographically approved
##### List of papers

Accurate numerical solution of time-dependent, high-dimensional partial differential equations (PDEs) usually requires efficient numerical techniques and massive-scale parallel computing. In this thesis, we implement and evaluate discretization schemes suited for PDEs of higher dimensionality, focusing on high order of accuracy and low computational cost.

Spatial discretization is particularly challenging in higher dimensions. The memory requirements for uniform grids quickly grow out of reach even on large-scale parallel computers. We utilize high-order discretization schemes and implement adaptive mesh refinement on structured hyperrectangular domains in order to reduce the required number of grid points and computational work. We allow for anisotropic (non-uniform) refinement by recursive bisection and show how to construct, manage and load balance such grids efficiently. In our numerical examples, we use finite difference schemes to discretize the PDEs. In the adaptive case we show how a stable discretization can be constructed using SBP-SAT operators. However, our adaptive mesh framework is general and other methods of discretization are viable.

For integration in time, we implement exponential integrators based on the Lanczos/Arnoldi iterative schemes for eigenvalue approximations. Using adaptive time stepping and a truncated Magnus expansion, we attain high levels of accuracy in the solution at low computational cost. We further investigate alternative implementations of the Lanczos algorithm with reduced communication costs.

As an example application problem, we have considered the time-dependent Schrödinger equation (TDSE). We present solvers and results for the solution of the TDSE on equidistant as well as adaptively refined Cartesian grids.

1. An implementation framework for solving high-dimensional PDEs on massively parallel computers$(function(){PrimeFaces.cw("OverlayPanel","overlay359735",{id:"formSmash:j_idt424:0:j_idt428",widgetVar:"overlay359735",target:"formSmash:j_idt424:0:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

2. Communication-efficient algorithms for numerical quantum dynamics$(function(){PrimeFaces.cw("OverlayPanel","overlay375913",{id:"formSmash:j_idt424:1:j_idt428",widgetVar:"overlay375913",target:"formSmash:j_idt424:1:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

3. Numerical evaluation of the Communication-Avoiding Lanczos algorithm$(function(){PrimeFaces.cw("OverlayPanel","overlay505799",{id:"formSmash:j_idt424:2:j_idt428",widgetVar:"overlay505799",target:"formSmash:j_idt424:2:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

4. Stable difference methods for block-oriented adaptive grids$(function(){PrimeFaces.cw("OverlayPanel","overlay758564",{id:"formSmash:j_idt424:3:j_idt428",widgetVar:"overlay758564",target:"formSmash:j_idt424:3:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

5. Data structures and algorithms for high-dimensional structured adaptive mesh refinement$(function(){PrimeFaces.cw("OverlayPanel","overlay758570",{id:"formSmash:j_idt424:4:j_idt428",widgetVar:"overlay758570",target:"formSmash:j_idt424:4:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

6. Parallel data structures and algorithms for high-dimensional structured adaptive mesh refinement$(function(){PrimeFaces.cw("OverlayPanel","overlay758573",{id:"formSmash:j_idt424:5:j_idt428",widgetVar:"overlay758573",target:"formSmash:j_idt424:5:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

References$(function(){PrimeFaces.cw("TieredMenu","widget_formSmash_lower_j_idt1101",{id:"formSmash:lower:j_idt1101",widgetVar:"widget_formSmash_lower_j_idt1101",autoDisplay:true,overlay:true,my:"left top",at:"left bottom",trigger:"formSmash:lower:referencesLink",triggerEvent:"click"});}); $(function(){PrimeFaces.cw("OverlayPanel","widget_formSmash_lower_j_idt1102_j_idt1104",{id:"formSmash:lower:j_idt1102:j_idt1104",widgetVar:"widget_formSmash_lower_j_idt1102_j_idt1104",target:"formSmash:lower:j_idt1102:permLink",showEffect:"blind",hideEffect:"fade",my:"right top",at:"right bottom",showCloseIcon:true});});