657 research outputs found
An Optimized and Scalable Eigensolver for Sequences of Eigenvalue Problems
In many scientific applications the solution of non-linear differential
equations are obtained through the set-up and solution of a number of
successive eigenproblems. These eigenproblems can be regarded as a sequence
whenever the solution of one problem fosters the initialization of the next. In
addition, in some eigenproblem sequences there is a connection between the
solutions of adjacent eigenproblems. Whenever it is possible to unravel the
existence of such a connection, the eigenproblem sequence is said to be
correlated. When facing with a sequence of correlated eigenproblems the current
strategy amounts to solving each eigenproblem in isolation. We propose a
alternative approach which exploits such correlation through the use of an
eigensolver based on subspace iteration and accelerated with Chebyshev
polynomials (ChFSI). The resulting eigensolver is optimized by minimizing the
number of matrix-vector multiplications and parallelized using the Elemental
library framework. Numerical results show that ChFSI achieves excellent
scalability and is competitive with current dense linear algebra parallel
eigensolvers.Comment: 23 Pages, 6 figures. First revision of an invited submission to
special issue of Concurrency and Computation: Practice and Experienc
An Sn Application of Homotopy Continuation in Neutral Particle Transport
The objective of this dissertation is to investigate the usefulness of homotopy continuation applied in the context of neutral particle transport where traditional methods of acceleration degrade. This occurs in higher dimensional heterogeneous problems [51]. We focus on utilizing homotopy continuation as a means of providing a better initial guess for difficult problems. We investigate various homotopy formulations for two primary diffcult problems: a thick-diffusive fixed internal source, and a k-eigenvalue problem with high dominance ratio. We also investigate the usefulness of homotopy continuation for computationally intensive problems with 30-energy groups. We find that homotopy continuation exhibits usefulness in specific problem formulations. In the thick-diffusive problem it shows benefit when there is a strong internal source in thin materials. In the k-eigenvalue problem, homotopy continuation provides an improvement in convergence speed for fixed point iteration methods in high dominance ratio problems. We also show that one of our imbeddings successfully stabilizes the nonlinear formulation of the k-eigenvalue problem with a high dominance ratio
- …