54 research outputs found
Hierarchical interpolative factorization for elliptic operators: differential equations
This paper introduces the hierarchical interpolative factorization for
elliptic partial differential equations (HIF-DE) in two (2D) and three
dimensions (3D). This factorization takes the form of an approximate
generalized LU/LDL decomposition that facilitates the efficient inversion of
the discretized operator. HIF-DE is based on the multifrontal method but uses
skeletonization on the separator fronts to sparsify the dense frontal matrices
and thus reduce the cost. We conjecture that this strategy yields linear
complexity in 2D and quasilinear complexity in 3D. Estimated linear complexity
in 3D can be achieved by skeletonizing the compressed fronts themselves, which
amounts geometrically to a recursive dimensional reduction scheme. Numerical
experiments support our claims and further demonstrate the performance of our
algorithm as a fast direct solver and preconditioner. MATLAB codes are freely
available.Comment: 37 pages, 13 figures, 12 tables; to appear, Comm. Pure Appl. Math.
arXiv admin note: substantial text overlap with arXiv:1307.266
Hierarchical interpolative factorization for elliptic operators: integral equations
This paper introduces the hierarchical interpolative factorization for
integral equations (HIF-IE) associated with elliptic problems in two and three
dimensions. This factorization takes the form of an approximate generalized LU
decomposition that permits the efficient application of the discretized
operator and its inverse. HIF-IE is based on the recursive skeletonization
algorithm but incorporates a novel combination of two key features: (1) a
matrix factorization framework for sparsifying structured dense matrices and
(2) a recursive dimensional reduction strategy to decrease the cost. Thus,
higher-dimensional problems are effectively mapped to one dimension, and we
conjecture that constructing, applying, and inverting the factorization all
have linear or quasilinear complexity. Numerical experiments support this claim
and further demonstrate the performance of our algorithm as a generalized fast
multipole method, direct solver, and preconditioner. HIF-IE is compatible with
geometric adaptivity and can handle both boundary and volume problems. MATLAB
codes are freely available.Comment: 39 pages, 14 figures, 13 tables; to appear, Comm. Pure Appl. Mat
Proceedings for the ICASE Workshop on Heterogeneous Boundary Conditions
Domain Decomposition is a complex problem with many interesting aspects. The choice of decomposition can be made based on many different criteria, and the choice of interface of internal boundary conditions are numerous. The various regions under study may have different dynamical balances, indicating that different physical processes are dominating the flow in these regions. This conference was called in recognition of the need to more clearly define the nature of these complex problems. This proceedings is a collection of the presentations and the discussion groups
Recommended from our members
Final report on the Copper Mountain conference on multigrid methods
The Copper Mountain Conference on Multigrid Methods was held on April 6-11, 1997. It took the same format used in the previous Copper Mountain Conferences on Multigrid Method conferences. Over 87 mathematicians from all over the world attended the meeting. 56 half-hour talks on current research topics were presented. Talks with similar content were organized into sessions. Session topics included: fluids; domain decomposition; iterative methods; basics; adaptive methods; non-linear filtering; CFD; applications; transport; algebraic solvers; supercomputing; and student paper winners
Spectral and High Order Methods for Partial Differential Equations ICOSAHOM 2018
This open access book features a selection of high-quality papers from the presentations at the International Conference on Spectral and High-Order Methods 2018, offering an overview of the depth and breadth of the activities within this important research area. The carefully reviewed papers provide a snapshot of the state of the art, while the extensive bibliography helps initiate new research directions
Trinity: A Unified Treatment of Turbulence, Transport, and Heating in Magnetized Plasmas
To faithfully simulate ITER and other modern fusion devices, one must resolve
electron and ion fluctuation scales in a five-dimensional phase space and time.
Simultaneously, one must account for the interaction of this turbulence with
the slow evolution of the large-scale plasma profiles. Because of the enormous
range of scales involved and the high dimensionality of the problem, resolved
first-principles global simulations are very challenging using conventional
(brute force) techniques. In this thesis, the problem of resolving turbulence
is addressed by developing velocity space resolution diagnostics and an
adaptive collisionality that allow for the confident simulation of velocity
space dynamics using the approximate minimal necessary dissipation. With regard
to the wide range of scales, a new approach has been developed in which
turbulence calculations from multiple gyrokinetic flux tube simulations are
coupled together using transport equations to obtain self-consistent,
steady-state background profiles and corresponding turbulent fluxes and
heating. This approach is embodied in a new code, Trinity, which is capable of
evolving equilibrium profiles for multiple species, including electromagnetic
effects and realistic magnetic geometry, at a fraction of the cost of
conventional global simulations. Furthermore, an advanced model physical
collision operator for gyrokinetics has been derived and implemented, allowing
for the study of collisional turbulent heating, which has not been extensively
studied. To demonstrate the utility of the coupled flux tube approach,
preliminary results from Trinity simulations of the core of an ITER plasma are
presented.Comment: 187 pages, 53 figures, Ph.D. thesis in physics at University of
Maryland, single-space versio
- …