An O(log N) Algorithm for Massively Parallel Molecular Dynamics Simulations
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
In this paper new factorization techniques for computation of the Operational Space Mass Matrix and its inverse are developed.
The National Airspace System is a highly complex system of systems within which a number of participants with widely varying business and operating models exist. From the airspace user's perspective, a means by which to operate flights in a more flexible and efficient manner is highly desired to meet their business objectives. From the air navigation service provider's viewpoint, there is a need for increasing the capacity of the airspace, while maintaining or increasing the levels of efficiency and safety that currently exist in order to meet the charter under which they operate. Enhancing the communication between airspace operators and users is essential in order to meet these demands. In the spring of 2015, a prototype system that implemented an airborne tool to optimize en-route flight paths for fuel and time savings was designed and tested. The system utilized in-flight Internet as a high-bandwidth data link to facilitate collaborative decision making between the flight deck and an airline dispatcher. The system was tested and demonstrated in a laboratory environment, as well as in-situ. Initial results from these tests indicate that this system is not only feasible, but could also serve as a growth path and testbed for future air traffic management concepts that rely on shared situational awareness through data exchange and electronic negotiation between multiple entities operating within the National Airspace System.
There is an overwhelming desire to create and enhance communication mechanisms between entities that operate within the National Airspace System. Furthermore, airlines are always extremely interested in increasing the efficiency of their flights. An innovative system prototype was developed and tested that improves collaborative decision making without modifying existing infrastructure or operational procedures within the current Air Traffic Management System. This system enables collaboration between flight crew and airline dispatchers to share and assess optimized flight routes through an Internet connection. Using a sophisticated medium-fidelity flight simulation environment, a rapid-prototyping development, and a unified modeling language, the software was designed to ensure reliability and scalability for future growth and applications. Ensuring safety and security were primary design goals, therefore the software does not interact or interfere with major flight control or safety systems. The system prototype demonstrated an unprecedented use of in-flight Internet to facilitate effective communication with Airline Operations Centers, which may contribute to increased flight efficiency for airlines.
The EMU (Extravehicular Mobility Unit) Garments and associated hardware were evaluated to determine if they are qualified for use in the first Lunar Landing Mission, Apollo XI (s/c 107/LM-5). It has been determined that the subject equipment is qualified for this application as summarized in this report.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
This highlight details the key takeaways from a project that utilized NLR's Fleet Research, Energy Data, and Insights (FleetREDI) data analysis pipeline. National Laboratory of the Rockies researchers developed and demonstrated low-cost, open-source Arduino data loggers with 3D-printed cases that are compatible with global navigational systems and built with components available ubiquitously worldwide, enabling cost-effective collection and analysis of fleet operational data. Validated on an overseas transit bus fleet, NLR analysis showed that, with sufficient charging opportunities, 90% of observed duty cycles could be accomplished by electric buses with no modifications to operations.
BACKGROUND. The purposes of the present study were (1) to establish normal values for the regression of log(power) on log(frequency) for, RR-interval fluctuations in healthy middle-aged persons, (2) to determine the effects of myocardial infarction on the regression of log(power) on log(frequency), (3) to determine the effect of cardiac denervation on the regression of log(power) on log(frequency), and (4) to assess the ability of power law regression parameters to predict death after myocardial infarction. METHODS AND RESULTS. We studied three groups: (1) 715 patients with recent myocardial infarction; (2) 274 healthy persons age and sex matched to the infarct sample; and (3) 19 patients with heart transplants. Twenty-four-hour RR-interval power spectra were computed using fast Fourier transforms and log(power) was regressed on log(frequency) between 10(-4) and 10(-2) Hz. There was a power law relation between log(power) and log(frequency). That is, the function described a descending straight line that had a slope of approximately -1 in healthy subjects. For the myocardial infarction group, the regression line for log(power) on log(frequency) was shifted downward and had a steeper negative slope (-1.15). The transplant (denervated) group showed a larger downward shift in the regression line and a much steeper negative slope (-2.08). The correlation between traditional power spectral bands and slope was weak, and that with log(power) at 10(-4) Hz was only moderate. Slope and log(power) at 10(-4) Hz were used to predict mortality and were compared with the predictive value of traditional power spectral bands. Slope and log(power) at 10(-4) Hz were excellent predictors of all-cause mortality or arrhythmic death. To optimize the prediction of death, we calculated a log(power) intercept that was uncorrelated with the slope of the power law regression line. We found that the combination of slope and zero-correlation log(power) was an outstanding predictor, with a relative risk of > 10, and was better than any combination of the traditional power spectral bands. The combination of slope and log(power) at 10(-4) Hz also was an excellent predictor of death after myocardial infarction. CONCLUSIONS. Myocardial infarction or denervation of the heart causes a steeper slope and decreased height of the power law regression relation between log(power) and log(frequency) of RR-interval fluctuations. Individually and, especially, combined, the power law regression parameters are excellent predictors of death of any cause or arrhythmic death and predict these outcomes better than the traditional power spectral bands.
Binary heaps and binary search trees merge efficiently. We introduce a new amortized analysis that allows us to prove the cost of merging either binary heaps or balanced binary trees is O(l), in the amortized sense. The standard set of other operations (create, insert, delete, extract minimum, in the case of binary heaps, and balanced binary trees, as well as a search operation for balanced binary trees) remain with a cost of O(log n). For binary heaps implemented as arrays, we show a new merge algorithm that has a single operation cost for merging two heaps, a and b, of O(absolute value of a + min(log absolute value of b log log absolute value of b. log absolute value of a log absolute value of b). This is an improvement over O(absolute value of a + log absolute value of a log absolute value of b). The cost of the new merge is so low that it can be used in a new structure which we call shadow heaps. to implement the insert operation to a tunable efficiency. Shadow heaps support the insert operation for simple priority queues in an amortized time of O(f(n)) and other operations in time O((log n log log n)/f (n)), where 1 less than or equal to f (n) less than or equal to log log n. More generally, the results here show that any data structure with operations that change its size by at most one, with the exception of a merge (aka meld) operation, can efficiently amortize the cost of the merge under conditions that are true for most implementations of binary heaps and search trees.
The spectral flux of gravitational radiation incident on earth from the Galactic W UMa binaries, unevolved binaries, cataclysmic binaries, neutron star binaries, black hole-neutron star binaries, and close white dwarf binary (CWDB) systems is calculated. The peak values for the strain amplitude hv produced by the binaries are: log hv = -17.0/sq rt Hz at log v = -4.40 Hz for W UMa's; log hv = -16.5/sq rt Hz at log v = -5.50 Hz for unevolved binaries; log hv = -18.2/sq rt Hz at log v = -5.15 Hz for neutron star binaries; log hv = -17.0/sq rt Hz at log v = -4.7 Hz for black hole-neutron star binaries; log hv = -18.0 /sq rt Hz at log v = -4.10 Hz for cataclysmic binaries, and log hv = -27.0/sq rt Hz at log v = -4.0 Hz for CWDBs. The gravitational flux at ultralow frequencies is emitted mainly by unevolved binaries. The integrated flux incident on earth is about 2.4 x 10 to the -9th ergs/sq cm/s.
The known biases affecting the KONUS gamma burst cumulative number distribution are discussed and categorized as intrinsic (unavoidable) and artificial (avoidable). Artificial biases are removed to obtain a revised, more faithful representation of the true KONUS log N-log Pc, where Pc is in counts (1/0.25 s). A low Pc turnover is still apparent in the revised log N-log Pc. Conflicting claims that the KONUS log N-log Pc is compatible with a nearby disk or an isotropic, homogeneous distribution are tested using the revised log N-log Pc and a rigorous treatment of the KONUS intrinsic gamma background bias. The claims are shown to be improbable. Other possibilities and problms concerning log N-log Pc are discussed, but the need for a complete explanation of the low Pc turnover remains.
Historically, tremendous effort has been put into processing and displaying mission health and safety telemetry data; and relatively little attention has been paid to extracting information from missions time-tagged event log messages. Todays missions may log tens of thousands of messages per day and the numbers are expected to dramatically increase as satellite fleets and constellations are launched, as security monitoring continues to evolve, and as the overall complexity of ground system operations increases. The logs may contain information about orbital events, scheduled and actual observations, device status and anomalies, when operators were logged on, when commands were resent, when there were data drop outs or system failures, and much much more. When dealing with distributed space missions or operational fleets, it becomes even more important to systematically analyze this data. Several advanced information systems technologies make it appropriate to now develop analytic capabilities which can increase mission situational awareness, reduce mission risk, enable better event-driven automation and cross-mission collaborations, and lead to improved operations strategies: Industry Standard for Log Messages. The Object Management Group (OMG) Space Domain Task Force (SDTF) standards organization is in the process of creating a formal standard for industry for event log messages. The format is based on work at NASA GSFC. Open System Architectures. The DoD, NASA, and others are moving towards common open system architectures for mission ground data systems based on work at NASA GSFC with the full support of the commercial product industry and major integration contractors. Text Analytics. A specific area of data analytics which applies statistical, linguistic, and structural techniques to extract and classify information from textual sources. This presentation describes work now underway at NASA to increase situational awareness through the collection of non-telemetry mission operations information into a common log format and then providing display and analytics tools to provide in-depth assessment of the log contents. The work includes: Common interface formats for acquiring time-tagged text messages Conversion of common files for schedules, orbital events, and stored commands to the common log format Innovative displays to depict thousands of messages on a single display Structured English text queries against the log message data store, extensible to a more mature natural language query capability Goal of speech-to-text and text-to-speech additions to create a personal mission operations assistant to aid on-console operations. A wide variety of planned uses identified by the mission operations teams will be discussed.
We develop a framework to study the relation between the stellar mass of a galaxy and the total mass of its host dark matter halo using galaxy clustering and galaxy-galaxy lensing measurements. We model a wide range of scales, roughly from $\sim 100 \; {\rm kpc}$ to $\sim 100 \; {\rm Mpc}$, using a theoretical framework based on the Halo Occupation Distribution and data from Year 3 of the Dark Energy Survey (DES) dataset. The new advances of this work include: 1) the generation and validation of a new stellar mass-selected galaxy sample in the range of $\log M_\star/M_\odot \sim 9.6$ to $\sim 11.5$; 2) the joint-modeling framework of galaxy clustering and galaxy-galaxy lensing that is able to describe our stellar mass-selected sample deep into the 1-halo regime; and 3) stellar-to-halo mass relation (SHMR) constraints from this dataset. In general, our SHMR constraints agree well with existing literature with various weak lensing measurements. We constrain the free parameters in the SHMR functional form $\log M_\star (M_h) = \log(εM_1) + f\left[ \log\left( M_h / M_1 \right) \right] - f(0)$, with $f(x) \equiv -\log(10^{αx}+1) + δ[\log(1+\exp(x))]^γ/ [1+\exp(10^{-x})]$, to be $\log M_1 = 11.506^{+0.325}_{-0.404}$, $\log ε= -1.632^{+0.306}_{-0.181}$, $α= -1.638^{+0.108}_{-0.099}$, $γ= 0.596^{+0.251}_{-0.210}$ and $δ= 3.810^{+2.045}_{-1.811}$. The inferred average satellite fraction is within $\sim 5-35\%$ for our fiducial results and we do not see any clear trends with redshift or stellar mass. Furthermore, we find that the inferred average galaxy bias values follow the generally expected trends with stellar mass and redshift. Our study is the first SHMR in DES in this mass range, and we expect the stellar mass sample to be of general interest for other science cases.
The double-ring planetary nebula IC 1297 resembles NGC 7662 in appearance, although it is of much lower surface brightness. What is remarkable is the great strength of the dielectronic recombination O V line. Although this line is seen as a P Cygni feature in a number of planetary nebulae, it is in those instances accompanied by a strong continuum and other easily recognized features of stellar origin. No star is visible on CCD images of IC 1297. Optical region measurements are supplemented by IUE observations. The following logarithmic abundance values are found: log N(He) = 11.065; log N(forbidden C) = 8.6; log N(N) = 8.1; log N(O) = 8.74; log N(Ne) = 8.16; log N(S) = 7.0; log N(Cl) = 5.4; log N(Ar) = 6.2. The nebula shows no dramatic pattern of nucleogenesis events.
In accordance with NASA Planetary Protection (PP) policy requirements, flight project hardware may be required to undergo microbial reduction processes to prevent the forward contamination of target planetary bodies with Earth organisms. Heat microbial reduction (HMR) is the most commonly employed modality used at JPL for reducing the microbial bioburden on flight hardware. In 2013, longstanding HMR specifications were abandoned, and revised specifications were developed which integrated the latest findings on bacterial spore heat resistivity. Revised decimal reduction values (D-values) for time-temperature lethality curves (110 °C to 200 °C) were developed to account for “hardy” bacterial spores that exhibit greater heat resistance than previously understood. Presented here is a comparative analysis of the revised NASA HMR specifications against empirical data compiled from recent JPL studies, and peer-reviewed, published literature. Bacillus sp. strain ATCC 29669 displayed high heat resistance, and this strain’s 4-log heat lethality curve was comparable to the revised 4-log specification. Spores of Bacillus atrophaeus ATCC 9372 displayed less heat resistance, and exhibited D-values which were less than the revised 3-log microbial reduction specifications. Extrapolations indicate that the current 6-log reduction credit applied to 350 °C for 1 hr. and 500 °C for 0.5 sec. is highly conservative. Projections indicate that a 10- to 18-log reduction of both hardy and non-hardy spores may be achievable at bakeouts of 350 °C for 1 hr. The findings reported here indicate the revised NASA HMR specifications from 110 °C to 200 °C are appropriate for achieving 4-log and 6-log reductions with hardy spore populations; however, for non-hardy spore populations, or for temperatures above 200 °C, the specifications are exceedingly conservative.
In recent years, there has been widespread adoption of machine learning-based approaches to automate the solving of partial differential equations (PDEs). Among these approaches, Gaussian processes (GPs) and kernel methods have garnered considerable interest due to their flexibility, robust theoretical guarantees, and close ties to traditional methods. They can transform the solving of general nonlinear PDEs into solving quadratic optimization problems with nonlinear, PDE-induced constraints. However, the complexity bottleneck lies in computing with dense kernel matrices obtained from pointwise evaluations of the covariance kernel, and its partial derivatives, a result of the PDE constraint and for which fast algorithms are scarce. The primary goal of this paper is to provide a near-linear complexity algorithm for working with such kernel matrices. We present a sparse Cholesky factorization algorithm for these matrices based on the near-sparsity of the Cholesky factor under a novel ordering of pointwise and derivative measurements. The near-sparsity is rigorously justified by directly connecting the factor to GP regression and exponential decay of basis functions in numerical homogenization. We then employ the Vecchia approximation of GPs, which is optimal in the Kullback-Leibler divergence, to compute the approximate factor. This enables us to compute ϵ-approximate inverse Cholesky factors of the kernel matrices with complexity O(N log d (N/ϵ)) in space and O(N log 2d (N/ϵ)) in time. We integrate sparse Cholesky factorizations into optimization algorithms to obtain fast solvers of the nonlinear PDE. We numerically illustrate our algorithm’s near-linear space/time complexity for a broad class of nonlinear PDEs such as the nonlinear elliptic, Burgers, and Monge-Ampère equations. In summary, we provide a fast, scalable, and accurate method for solving general PDEs with GPs and kernel methods.