Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel clustering”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

Data Structures for Extreme Scale Computing

As computing problems of national importance grow, the government meets the increased demand by funding the development of ever larger systems. The overarching goal of the work supported in part by this grant is to increase efficiency of programming and performing computations on these large computing systems. In past work, we have demonstrated that some of these computations once thought to require expensive hardware designs and/or complex, special-purpose programming may be executed efficiently on low-cost commodity cluster computing systems using a general-purpose “latency-tolerant” programming framework. One important developed application of the ideas underlying this framework is graph database technology supporting social network pattern matching used by US intelligence agencies to more quickly identify potential terrorist threats. This database application has been spun out by the Pacific Northwest National Laboratory, a Department of Energy Laboratory, into a commercial start-up, Trovares Inc. We explore an alternative application of the same underlying ideas to a well-studied challenge arising in engineering: solving unstructured sparse linear equations. Solving these equations is key to predicting the behavior of large electronic circuits before they are fabricated. Predicting that behavior ahead of fabrication means that designs can optimized and errors corrected ahead of the expense of manufacture.

42 ENGINEERING↗

Galaxy Alignments with Surrounding Structure in the Sloan Digital Sky Survey

Using data from the Sloan Digital Sky Survey Legacy Survey, we study the alignment of luminous galaxies with spectroscopic data with the surrounding larger-scale structure as defined by galaxies with only photometric data. We find that galaxies from the red sequence have a statistically significant tendency for their apparent long axes to align parallel to the projected surrounding structure. Red galaxies more luminous than the median of our sample (M r < –21.78) have a mean alignment angle $\langle$Φ$\rangle$ < 45°, indicating preferred parallel alignment, at a significance level >4.5σ on projected scales 0.1 Mpc < r p ≤ 7.5 Mpc. Fainter red galaxies have $\langle$Φ$\rangle$ < 45° at a significance level >4.3σ at scales 1 Mpc < r p < 3 Mpc. At a projected scale r p = 3.0 Mpc, the mean alignment angle decreases steadily with increasing luminosity for red galaxies with M r ≲ –22.5, reaching $\langle$Φ$\rangle$ = 40$^{°}_{.}$49 ± 0$^{°}_{.}$56 for the most luminous 1% (M r ~ –23.57). Galaxies from the blue sequence show no statistically significant tendency for their axes to align with larger-scale structure, regardless of galaxy luminosity. Galaxies in higher-density regions do not show a statistically significant difference in the mean alignment angle from galaxies in lower-density regions; this holds true for the faint blue, luminous blue, faint red, and luminous red subsets.

79 ASTRONOMY AND ASTROPHYSICS↗

Considerations on the Use of Custom Accelerators for Big Data Analytics

Accelerators, including Graphic Processing Units (GPUs) for gen- eral purpose computation, many-core designs with wide vector units (e.g., Intel Phi), have become a common component of many high performance clusters. The appearance of more stable and reliable tools tools that can automatically convert code written in high-level specifications with annotations (such as C or C++) to hardware de- scription languages (High-Level Synthesis - HLS), is also setting the stage for a broader use of reconfigurable devices (e.g., Field Pro- grammable Gate Arrays - FPGAs) in high performance system for the implementation of custom accelerators, helped by the fact that new processors include advanced cache-coherent interconnects for these components. In this chapter, we briefly survey the status of the use of accelerators in high performance systems targeted at big data analytics applications. We argue that, although the progress in the use of accelerators for this class of applications has been sig- nificant, differently from scientific simulations there still are gaps to close. This is particularly true for the ”irregular” behaviors exhibited by no-SQL, graph databases. We focus our attention on the limits of HLS tools for data analytics and graph methods, and discuss a new architectural template that better fits the requirement of this class of applications. We validate the new architectural templates by mod- ifying the Graph Engine for Multithreaded System (GEMS) frame- work to support accelerators generated with such a methodology, and testing with queries coming from the Lehigh University Benchmark (LUBM). The architectural template enables better supporting the task and memory level parallelism present in graph methods by sup- porting a new control model and a enhanced memory interface. We show that out solution allows generating parallel accelerators, pro- viding speed ups with respect to conventional HLS flows. We finally draw conclusions and present a perspective on the use of reconfig- urable devices and Design Automation tools for data analytics.

Data Analytics↗

Carbon Nanotubes in Water: MD Simulations of Internal and External Flow, Self Organization

We have developed computational tools, based on particle codes, for molecular dynamics (MD) simulation of carbon nanotubes (CNT) in aqueous environments. The interaction of CNTs with water is envisioned as a prototype for the design of engineering nano-devices, such as artificial sterocillia and molecular biosensors. Large scale simulations involving thousands of water molecules are possible due to our efficient parallel MD code that takes long range electrostatic interactions into account. Since CNTs can be considered as rolled up sheets of graphite, we expect the CNT-water interaction to be similar to the interaction of graphite with water. However, there are fundamental differences between considering graphite and CNTs, since the curvature of CNTs affects their chemical activity and also since capillary effects play an important role for both dynamic and static behaviour of materials inside CNTs. In recent studies Gordillo and Marti described the hydrogen bond structure as well as time dependent properties of water confined in CNTs. We are presenting results from the development of force fields describing the interaction of CNTs and water based on ab-initio quantum mechanical calculations. Furthermore, our results include both water flows external to CNTs and the behaviour of water nanodroplets inside heated CNTs. In the first case (external flows) the hydrophobic behaviour of CNTs is quantified and we analyze structural properties of water in the vicinity of CNTs with diagnostics such as hydrogen bond distribution, water dipole orientation and radial distribution functions. The presence of water leads to attractive forces between CNTs as a result of their hydrophobicity. Through extensive simulations we quantify these attractive forces in terms of the number and separation of the CNT. Results of our simulations involving arrays of CNTs indicate that these exhibit a hydrophobic behaviour that leads to self-organising structures capable of trapping water clusters. In the second case (internal flows) we study the behaviour of water droplets confined inside CNTs. Constant temperature simulations allow us to capture structural properties such as the contact angles and density profiles of the equilibrated drops. By heating and subsequently cooling of the CNT, we are able to measure the evaporation and the condensation rate of the entrapped water.

Jaffe, Richard L.↗

Life cycle of streaks in the buffer layer of wall-bounded turbulence

In this work, streaks in the buffer layer of wall-bounded turbulence are tracked in time to study their life cycle. Spatially and temporally resolved direct numerical simulation data are used to analyze the strong wall-parallel movements conditioned to low-speed streamwise flow. The analysis of the streaks shows that there is a clear distinction between wall-attached and detached streaks, and that the wall-attached streaks can be further categorized into streaks that are contained in the buffer layer and the ones that reach the outer region. The results reveal that streaks are born in the buffer layer, coalescing with each other to create larger streaks that are still attached to the wall. Once the streak becomes large enough, it starts to meander due to the large streamwise-to-wall-normal aspect ratio, and consequently the elongation in the streamwise direction, which makes it more difficult for the streak to be oriented strictly in the streamwise direction. While the continuous interaction of the streaks allows the superstructure to span extremely long temporal and length scales, individual streak components are relatively small and short-lived. Tall-attached streaks eventually split into wall-attached and wall-detached components. These wall-detached streaks have a strong wall-normal velocity away from the wall, similar to ejections or bursts observed in the literature. Conditionally averaging the flow fields to these split events show that the detached streak has not only a larger wall-normal velocity compared to the wall-attached counterpart, it also has a larger (less negative) streamwise velocity, similar to the velocity field at the tip of a vortex cluster.

42 ENGINEERING↗

Accelerated Simulation of Air Pollution Using NVIDIA RAPIDS

Atmospheric chemistry models are a central tool to study and forecast the impact of air pollution on the environment, vegetation, and human health. However, the numerical simulation of chemical kinetics is computationally expensive due to the stiffness of the system of ordinary differential equations that describes atmospheric chemistry. Here we present an alternative approach to the computation of atmospheric chemistry based on machine learning. Our training data set is produced using the NASA Goddard Earth Observing System (GEOS) model with GEOS-Chem chemistry, run on the NASA Center for Climate Simulation (NCCS) Discover supercomputing cluster on 384 Intel Xeon Haswell cores. This model spends more than 50% of total run time on solving atmospheric chemistry. The data set contains as input features the air pollution concentrations before solving the differential equations, together with some key physical parameters such as temperature and sun intensity. As target variables we define the air pollution concentrations after solving the differential equations. Using Dask-cuDF and Dask-XGBoost on the NVIDIA RAPIDS platform on 8 Tesla V100 GPUs, we generate from this training set gradient boosted decision tree models that can reproduce the simulation of chemical kinetics. We do this on the NCCS Advanced Data Analytics Platform (ADAPT) science cloud environment. Our application takes full advantage of recent advances in Dask-XGBoost, such as multi-node and multi-GPU scaling for distributed training with large data sets. The increase in training data size enabled by this is critical to capture the full range of chemical environments encountered across the globe and all annual seasons.The boosted tree models offer good predictability and show many of the features of the full chemistry reference simulation. Further improvements can be achieved through mass balance considerations and by accounting for error correlations. We incorporate the boosted tree models into the GEOS reference model using XGBoost's C API. This enables a seamless integration of the GPU trained models into GEOS-Chem, which is written in Fortran and optimized for use in a massively parallel CPU environment. We show the benefits of this approach and discuss the potential speedup of this machine learning accelerated atmospheric chemistry model.

Keller, Christoph A.↗

Machine learning–driven multiscale modeling reveals lipid-dependent dynamics of RAS signaling proteins

Significance Here we present an unprecedented multiscale simulation platform that enables modeling, hypothesis generation, and discovery across biologically relevant length and time scales to predict mechanisms that can be tested experimentally. We demonstrate that our predictive simulation-experimental validation loop generates accurate insights into RAS-membrane biology. Evaluating over 100,000 correlated simulations, we show that RAS–lipid interactions are dynamic and evolving, resulting in: 1) a reordering and selection of lipid domains in realistic eight-lipid bilayers, 2) clustering of RAS into multimers correlating with specific lipid fingerprints, 3) changes in the orientation of the RAS G-domain impacting its ability to interact with effectors, and 4) demonstration that RAS–RAS G-domain interfaces are nonspecific in these putative signaling domains.

59 BASIC BIOLOGICAL SCIENCES↗

Equatorial Electronic Structure in the Uranyl Ion: Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4

Electric field gradient (EFG) tensors in the equatorial plane of the linear UO 2 2+ ion have been measured by nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) experiments and computed by relativistic Kohn–Sham methods with and without environment embedding for Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4 . This approach expands the possibilities for probing the electronic structure in uranyl complexes beyond the strongly covalent U–O bonds. The combined analyses find that one of the two largest principal EFG tensor components at the halogen sites points along the U–X bond (X = Cl, Br), and the second is parallel to the UO 2 2+ ion; in Cs 2 UO 2 Cl 4 , the components are nearly equal in magnitude, whereas in Cs 2 UO 2 Br 4 , due to short-range bromide–cesium interactions, the equatorial component is dominant for one pair of Br sites and the axial component is larger for the second pair. Additionally, the directions and relative magnitudes of the field gradient principal axes are found to be sensitive to the σ and π electron donation by the ligands and the model of the environment. Chlorine-35 NQR spectra of 235 U-depleted and 235 U-enriched Cs 2 UO 2 Cl 4 exhibited no uranium-isotope-dependent shift, but the resonance of the depleted sample displayed a 58% broader line width.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identification of 2-Hydroxyacyl-CoA Synthases with High Acyloin Condensation Activity for Orthogonal One-Carbon Bioconversion

One-carbon (C1) compounds are emerging as cost-effective and potentially carbon-negative feedstocks for biomanufacturing, which require efficient, versatile metabolic platforms for the synthesis of value-added products. Synthetic formyl-CoA elongation (FORCE) pathways allow diverse product synthesis from C1 compounds via iterative C1 elongation, operating independently from the host metabolism with reduced engineering complexity and improved theoretical yields. However, a major bottleneck was identified as the suboptimal kinetics of the core C1–C1 condensation enzyme, 2-hydroxyacyl-CoA synthase (HACS), catalyzing the acyloin condensation reaction between formaldehyde and formyl-CoA. Furthermore, we used a combinatorial approach of bioprospecting and rational protein engineering to identify multiple HACS variants with significantly improved activities toward C1 substrates. Sequence and structure alignment of the active variants elucidated the key regions for the catalytic function, which were targeted for mutagenesis, leading to improved catalytic efficiency. In parallel, a consecutive round of bioprospecting for homologs with high similarity with active variants revealed a highly active HACS variant exhibiting up to 7-fold improvement in catalytic efficiency (k cat /K M ) and 14-fold improvement in the FORCE pathway flux in vivo compared to the previous reports. Upon further optimization of the downstream pathway, the orthogonal C1-to-product bioconversion system showed a metabolic flux of up to 700 μM glycolate OD –1 h –1 (2.1 mmol gDCW –1 h –1 ) and an industrially relevant glycolate titer, rate, and yield of 5.2 g L –1 (67.8 mM), 0.22 g L –1 h –1 , and 94% carbon yield, respectively.

2-hydroxyacyl-CoA synthase↗

Washboard Terrain on Pluto

Washboard texture or patterning consists of fields of parallel to sub-parallel ridges typically spaced ~1-2 km crest to crest and a few 100 m in amplitude (Fig. 4a in Moore et al., 2016, Science, 351, 1284-1293). For the most part, underlying topography can be easily discerned. We will refer to discrete, well-bounded patches of these landforms as Washboard Terrain (WT). WT is observed to occur along the rim, and just beyond the rim, of Sputnik basin from the West to NNW. Where it is seen in high-resolution data, it has clearly defined limits, beyond which it would be able to be seen if it were there. WT doesn't occur at very low latitudes or very high latitudes (ranging from 22degN to 62degN). WT seems to occur most conspicuously on relatively level, gently sloping terrain. It is restricted to elevations between approximately 2 km to less than +1.5 km (i.e. not at high elevations). The most noticeable regional aspect of the area in which WT occurs is the sinuous valley network, which is suspected to have been formed, or at least substantially modified, by glaciation. WT also appears to occur mainly on an intermediate-albedo reddish material, where seen in enhanced color data. Where it occurs in level terrain, WT tends to trend ENE - there doesn't seem to be a strong local control of its orientation in response to valley drainage directions. WT can display a greater range of orientations where it occurs in higher-relief (not higher elevation) settings such as spurs. WT appears superposed on very ancient landscapes, but is itself cratered locally by clusters of small (approximately 1-3 km) craters, which may be secondaries. This implies that WT may be intermediate in age. Of several working hypotheses, we currently provisionally favor that WT may be akin to terrestrial recessional moraines (or de Geer moraines) associated with the retreat of a higher stand of N2 glaciation that once overfilled Sputnik basin. These putative moraine features may owe their spacing to superseasonal retreat on Milankovitch timescales of approximately 1 Ma. If this hypothesis has validity, then perhaps the intermediate-albedo reddish material may be akin to ground moraine deposits.

Moore, Jeffrey M.↗

Quasi-heterogeneous efficient 3-D discrete ordinates CANDU calculations using Attila

In this paper, 3-D quasi-heterogeneous large scale parallel Attila calculations of a generic CANDU test problem consisting of 42 complete fuel channels and a perpendicular to fuel reactivity device are presented. The solution method is that of discrete ordinates SN and the computational model is quasi-heterogeneous, i.e. fuel bundle is partially homogenized into five homogeneous rings consistently with the DRAGON code model used by the industry for the incremental cross-section generation. In calculations, the HELIOS-generated 45 macroscopic cross-sections library was used. This approach to CANDU calculations has the following advantages: 1) it allows detailed bundle (and eventually channel) power calculations for each fuel ring in a bundle, 2) it allows the exact reactivity device representation for its precise reactivity worth calculation, and 3) it eliminates the need for incremental cross-sections. Our results are compared to the reference Monte Carlo MCNP solution. In addition, the Attila SN method performance in CANDU calculations characterized by significant up scattering is discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Magnetic Reconnection: Theoretical and Observational Perspectives: Preface

Magnetic reconnection is a fundamental plasma-physical process by which energy stored in a magnetic field is converted, often explosively, into heat and the kinetic energy of the charged particles that constitute the plasma. It occurs in a variety of astrophysical settings, ranging from the solar corona to pulsar magnetospheres and winds, as well as in laboratory fusion experiments, where it is responsible for sawtooth crashes. First proposed by R.G. Giovanelli in the late I 940s as the mechanism responsible for solar flares, magnetic reconnection was invoked at the beginning of the space age to explain not just solar flares but also the transfer of energy, mass, and momentum from the solar wind to Earth's magnetosphere and the subsequent storage and release of the transferred energy in the magnetotai\. During the half century or so that has followed the seminal theoretical works by J.W. Dungey, P.A. Sweet, E.N. Parker, and H.E. Petschek, in-situ measurements by Earth-orbiting satellites and remote-sensing observations of the solar corona have provided a growing body of evidence for the occurrence of reconnection at the Sun, in the solar wind, and in the near-Earth space environment. The last thirty years have also seen the development of laboratory reconnection experiments at a number of institutions. In parallel with the efforts of experimentalists in both space and laboratory plasma physics, theorists have investigated, analytically and with the help of increasingly powerful MHD, hybrid, and kinetic numerical simulations, the structure of the diffusion region, the factors controlling the rate, onset, and cessation of reconnection, and the detailed physics that enables the demagnetization of the ions and electrons and the topological reconfiguration of the magnetic field. Moreover, the scope of theoretical reconnection studies has been extended well beyond solar system and laboratory plasmas to include more exotic astrophysical plasma systems whose strong (10(exp 14)-10(exp 15) G) magnetic fields require that models of reconnection in these systems incorporate quantum electrodynamical, special relativistic, and radiative effects. The papers collected in this topical issue of Space Science Reviews cover different aspects of recent theoretical and observational work on magnetic reconnection in solar and space physics, astrophysics, and laboratory plasma physics. They derive from presentations given at a workshop on magnetic reconnection held in the Yosemite National Park, February 8-12,2010. The intent of the workshop was to stimulate, through a combination of tutorial talks, shorter focused talks, and extensive informal discussions, an interdisciplinary dialogue among members of the different research communities working on the problem of magnetic reconnection. One of the motivating considerations for holding the workshop was its relevance to NASA's Magnetospheric Multiscale (MMS) mission, scheduled for launch in 2014. The four identically instrumented MMS spacecraft are designed to study reconnect ion in Earth's magnetosphere and, specifically, to probe the electron diffusion region in order to determine the microphysical processes that enable the change in the topology of the magnetic field. Building on the achievements of the multi spacecraft Cluster and THEMIS missions, MMS will use the magnetosphere as an astrophysical plasma laboratory in which to test, through in-situ measurement of the plasma, energetic particles, and electric and magnetic fields, various models and theories that have emerged during the past twenty years, a period of extraordinarily productive theoretical and observational work.

Lewis, W. S.↗

HELICAL MAGNETIC FIELDS IN THE NGC 1333 IRAS 4A PROTOSTELLAR OUTFLOWS

We present Submillimeter Array polarization observations of the CO J = 3–2 line toward NGC 1333 IRAS 4A. The CO Stokes I maps at an angular resolution of ∼1″ reveal two bipolar outflows from the binary sources of NGC 1333 IRAS 4A. The kinematic features of the CO emission can be modeled by wind-driven outflows at ∼20° inclined from the plane of the sky. Close to the protostars the CO polarization, at an angular resolution of ∼2.″3, has a position angle approximately parallel to the magnetic field direction inferred from the dust polarizations. The CO polarization direction appears to vary smoothly from an hourglass field around the core to an arc-like morphology wrapping around the outflow, suggesting a helical structure of magnetic fields that inherits the poloidal fields at the launching point and consists of toroidal fields at a farther distance of outflow. The helical magnetic field is consistent with the theoretical expectations for launching and collimating outflows from a magnetized rotating disk. Considering that the CO polarized emission is mainly contributed from the low-velocity and low-resolution data, the helical magnetic field is likely a product of the wind–envelope interaction in the wind-driven outflows. The CO data reveal a PA of ∼30° deflection in the outflows. The variation in the CO polarization angle seems to correlate with the deflections. We speculate that the helical magnetic field contributes to ∼10° deflection of the outflows by means of Lorentz force.

79 ASTRONOMY AND ASTROPHYSICS↗

A 2000 Solar Mass Rotating Molecular Disk Around NGC 6334A

We present millimeter and centimeter wave spectroscopic observations of the H II region NGC 6334A. We have mapped the source in several transitions of CO, CS, and NH3. The molecular emission shows a distinct flattened structure in the east-west direction. This structure is probably a thick molecular disk or torus (2.2 x 0.9 pc) responsible for the bipolarity of the near-infrared (NIR) and radio continuum emission which extends in two "lobes" to the north and south of the shell-like H II region. The molecular disk is rotating from west to east (omega approximately equals 2.4 km/s.pc) about an axis approximately parallel to the radio and NIR emission lobes. By assuming virial equilibrium, we find that the molecular disk contains approximately 2000 solar mass. Single-component gas excitation model calculations show that the molecular gas in the disk is warmer and denser (T(sub k) approximately equals 60 K, n approximately equals 3000/cc) than the gas to the north and south (T(sub k) approximately equals 50 K, n approximately equals 400/cc). High resolution (approximately 5 sec) NH3 (3, 3) images of NGC 6334A reveal several small (approximately 0.1 pc) clumps, one of which lies southwest of the radio continuum shell, and is spatially coincident with a near-infrared source, IRS 20. A second NH3 clump is coincident with an H2O maser and the center of a molecular outflow. The dense gas tracers, CS J = 5 approaches 4 and 7 approaches 6, peak near IRS 20 and the H2O maser, not at NGC 6334A. IRS 20 has a substantial far-infrared (FIR) luminosity L(sub FIR) approximately 10(exp 5) solar luminosity, which indicates the presence of an O 7.5 star but has no detected radio continuum (F(sub 6 cm) < 0.02 Jy). The combination of dense gas, a large FIR luminosity and a lack of radio continuum can best be explained if IRS 20 is a protostar. A third clump of NH3 emission lies to the west of IRS 20 but is not associated with any other molecular or continuum features. The star formation activity in the region has moved west of NGC 6334A to IRS 20 and the H2O maser position. We suggest that NGC 6334A, IRS 20, and the H2O maser spot are part of a "protocluster" of stars which is condensing from the massive molecular disk. The similarity between the structure around NGC 6334A and other large (r approximately 1 pc), massive (M approximately 10(exp 3) solar mass), rotating disks (K3-50A and G10.6-0.4) suggests that this may be a common mechanism by which open clusters form.

Kraemer, Kathleen E.↗

HUBBLE TARANTULA TREASURY PROJECT: UNRAVELING TARANTULA'S WEB. I. OBSERVATIONAL OVERVIEW AND FIRST RESULTS

The Hubble Tarantula Treasury Project (HTTP) is an ongoing panchromatic imaging survey of stellar populations in the Tarantula Nebula in the Large Magellanic Cloud that reaches into the sub-solar mass regime (<0.5 M{sub Sun }). HTTP utilizes the capability of the Hubble Space Telescope to operate the Advanced Camera for Surveys and the Wide Field Camera 3 in parallel to study this remarkable region in the near-ultraviolet, optical, and near-infrared spectral regions, including narrow-band H{alpha} images. The combination of all these bands provides a unique multi-band view. The resulting maps of the stellar content of the Tarantula Nebula within its main body provide the basis for investigations of star formation in an environment resembling the extreme conditions found in starburst galaxies and in the early universe. Access to detailed properties of individual stars allows us to begin to reconstruct the temporal and spatial evolution of the stellar skeleton of the Tarantula Nebula over space and time on a sub-parsec scale. In this first paper we describe the observing strategy, the photometric techniques, and the upcoming data products from this survey and present preliminary results obtained from the analysis of the initial set of near-infrared observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Avoiding excess computation in asynchronous evolutionary algorithms

Abstract Asynchronous evolutionary algorithms are becoming increasingly popular as a means of making full use of many processors while solving computationally expensive search and optimization problems. These algorithms excel at keeping large clusters fully utilized, but may sometimes inefficiently sample an excess of fast‐evaluating solutions at the expense of higher‐quality, slow‐evaluating ones. We have previously introduced a steady‐state parent selection strategy, SWEET (“Selection whilE EvaluaTing”), that sometimes selects individuals that are still being evaluated and allows them to reproduce early. We perform a takeover‐time analysis that confirms that this strategy gives slow‐evaluating individuals that have higher fitnesses an increased ability to multiply in the population. We also find that SWEET appears effective at improving optimization performance on problems in which solution quality is positively correlated with evaluation time. We evaluate our approach on six simulated real‐valued optimization problems and three real‐world applications: an autonomous vehicle controller problem that involves tuning a spiking neural network and two adversarial EA problems. We further evaluate SWEET versus a basic asynchronous process in a simulated setting. We present evidence that SWEET outperforms basic asynchronous processes in a use‐case in which performance is positively correlated with evaluation time, and performs comparably (and often better) than basic asynchronous processes in several use‐cases where performance is negatively correlated with evaluation time. That said, in the cases where performance and evaluation time are negatively correlated the variance of outcomes for SWEET is notably high.

97 MATHEMATICS AND COMPUTING↗

Award ER25844: Minimizing System Noise Effects for Extreme-Scale Scientific Simulation Through Function Delegation

In software running on distributed computing clusters, time spent on communication between nodes in the cluster can be a significant portion of the overall computation time; background operating system tasks and other computational noise on the nodes of the system can have a significant impact on the amount of time this communication takes, especially on large systems. The research completed in this period has improved understanding of when such noise will have a significant impact. Specifically, it was demonstrated that not just noise on the nodes, but also noise on the network between nodes can have a significant impact on computation time. It was also demonstrated that noise patterns matter more than noise intensity: very regular noise can cause less disruption than lighter (on average) but less regular noise. It was also demonstrated that the effect of noise is more prominent as the speed of the network between nodes is increased. Furthermore, a tracing tool, Netgauge, was improved via our work, and a system simulator, LogGOPSim, was developed; they can be used by application developers to improve performance of their program and by system designers to mitigate the effects of noise by adjusting the noise characteristics of the operating system. Both have been made freely available as open source programs. In the course of developing these tools, we demonstrated weaknesses in existing methodologies for modeling communication, and we introduced a more detailed model, LogGOPS, for simulating systems. Not only were the deleterious effects of noise explored but we have also offered solutions. Our studies of simulations of system noise have led to specific recommendations on tuning systems to mitigate noise. We have also improved existing approaches to mitigating noise. Non-blocking collective communication avoids the effects of noise by letting communication continue simultaneously with computation (thus being non-blocking), so that the delays in communication introduced by noise have a smaller impact on overall computation time. Potentially, noise can be reduced much further by offloading communication tasks to a separate processing element than the operating system is using. We have improved our library LibNBC, which provides an implementation of non-blocking collectives, via this work. During this research, our proposal to include non-blocking collectives (which used LibNBC as a reference implementation) in the upcoming MPI-3 standard was accepted. As MPI is a ubiquitous and important standard for communication in parallel computing, this demonstrates a certain acceptance of the practicality and desirability of non-blocking collectives. Now that non-blocking collectives are a part of the standard we can expect to see optimized platform-specific implementations of non-blocking collectives. Also as part of this work we have also developed a language GOAL (Global Operation Assembly Language) that can be used as a starting point for defining languages to express optimized communication algorithms.

97 MATHEMATICS AND COMPUTING↗

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗