Karla Morris

Bio: Karla Morris is an academic researcher from Sandia National Laboratories. The author has contributed to research in topics: Fortran & Software. The author has an hindex of 9, co-authored 35 publications receiving 220 citations. Previous affiliations of Karla Morris include City University of New York.

Vortex locking in direct numerical simulations of quantum turbulence.

Abstract: Direct numerical simulations are used to examine the locking of quantized superfluid vortices and normal fluid vorticity in evolving turbulent flows. The superfluid is driven by the normal fluid, which undergoes either a decaying Taylor-Green flow or a linearly forced homogeneous isotropic turbulent flow, although the back reaction of the superfluid on the normal fluid flow is omitted. Using correlation functions and wavelet transforms, we present numerical and visual evidence for vortex locking on length scales above the intervortex spacing.

Extracting UML class diagrams from object-oriented fortran: ForUML

Abstract: Many scientists who implement computational science and engineering software have adopted the object-oriented (OO) Fortran paradigm. One of the challenges faced by OO Fortran developers is the inability to obtain high level software design descriptions of existing applications. Knowledge of the overall software design is not only valuable in the absence of documentation, it can also serve to assist developers with accomplishing different tasks during the software development process, especially maintenance and refactoring. The software engineering community commonly uses reverse engineering techniques to deal with this challenge. A number of reverse engineering-based tools have been proposed, but few of them can be applied to OO Fortran applications. In this paper, we propose a software tool to extract unified modeling language (UML) class diagrams from Fortran code. The UML class diagram facilitates the developers' ability to examine the entities and their relationships in the software system. The extracted diagrams enhance software maintenance and evolution. The experiments carried out to evaluate the proposed tool show its accuracy and a few of the limitations.

Fault Resilient Domain Decomposition Preconditioner for PDEs

Abstract: The move towards extreme-scale computing platforms challenges scientific simulations in many ways Given the recent tendencies in computer architecture development, one needs to reformulate legacy codes in order to cope with large amounts of communication, system faults, and requirements of low-memory usage per core In this work, we develop a novel framework for solving PDEs via domain decomposition that reformulates the solution as a state of knowledge with a probabilistic interpretation Such reformulation allows resiliency with respect to potential faults without having to apply fault detection, avoids unnecessary communication, and is generally well-suited for rigorous uncertainty quantification studies that target improvements of predictive fidelity of scientific models We demonstrate our algorithm for one-dimensional PDE examples where artificial faults have been implemented as bit flips in the binary representation of subdomain solutions

High-Performance design patterns for modern fortran

Abstract: This paper presents ideas for using coordinate-free numerics in modern Fortran to achieve code flexibility in the partial differential equation (PDE) domain. We also show how Fortran, over the last few decades, has changed to become a language well-suited for state-of-the-art software development. Fortran's new coarray distributed data structure, the language's class mechanism, and its side-effect-free, pure procedure capability provide the scaffolding on which we implement HPC software. These features empower compilers to organize parallel computations with efficient communication. We present some programming patterns that support asynchronous evaluation of expressions comprised of parallel operations on distributed data. We implemented these patterns using coarrays and the message passing interface (MPI). We compared the codes' complexity and performance. The MPI code is much more complex and depends on external libraries. The MPI code on Cray hardware using the Cray compiler is 1.5-2 times faster than the coarray code on the same hardware. The Intel compiler implements coarrays atop Intel's MPI library with the result apparently being 2-2.5 times slower than manually coded MPI despite exhibiting nearly linear scaling efficiency. As compilers mature and further improvements to coarrays comes in Fortran 2015, we expect this performance gap to narrow.

Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers.

Abstract: In recent years, the SUite of Nonlinear and DIfferential/ALgebraic equation Solvers (SUNDIALS) has been redesigned to better enable the use of application-specific and third-party algebraic solvers and data structures. Throughout this work, we have adhered to specific guiding principles that minimized the impact to current users while providing maximum flexibility for later evolution of solvers and data structures. The redesign was done through creation of new classes for linear and nonlinear solvers, enhancements to the vector class, and the creation of modern Fortran interfaces that leverage interoperability features of the Fortran 2003 standard. The vast majority of this work has been performed "behind-the-scenes," with minimal changes to the user interface and no reduction in solver capabilities or performance. However, these changes now allow advanced users to create highly customized solvers that exploit their problem structure, enabling SUNDIALS use on extreme-scale, heterogeneous computational architectures.

Experimental, numerical, and analytical velocity spectra in turbulent quantum fluid

Abstract: Turbulence in superfluid helium is unusual and presents a challenge to fluid dynamicists because it consists of two coupled, interpenetrating turbulent fluids: the first is inviscid with quantized vorticity, and the second is viscous with continuous vorticity. Despite this double nature, the observed spectra of the superfluid turbulent velocity at sufficiently large length scales are similar to those of ordinary turbulence. We present experimental, numerical, and theoretical results that explain these similarities, and illustrate the limits of our present understanding of superfluid turbulence at smaller scales.