Search NASA⌕ Search

SEARCH · Search NASA

Results for “arraying algorithms”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Software for Sparse Tensor Decomposition on Emerging Computing Architectures

In this paper, we develop software for decomposing sparse tensors that is portable to and performant on a variety of multicore, manycore, and GPU computing architectures. The result is a single code whose performance matches optimized architecture-specific implementations. The key to a portable approach is to determine multiple levels of parallelism that can be mapped in different ways to different architectures, and we explain how to do this for the matricized tensor times Khatri--Rao product (MTTKRP), which is the key kernel in canonical polyadic tensor decomposition. Our implementation leverages the Kokkos framework, which enables a single code to achieve high performance across multiple architectures that differ in how they approach fine-grained parallelism. We also introduce a new construct for portable thread-local arrays, which we call compile-time polymorphic arrays. Not only are the specifics of our approaches and implementation interesting for tuning tensor computations, but they also provide a roadmap for developing other portable high-performance codes. As a last step in optimizing performance, we modify the MTTKRP algorithm itself to do a permuted traversal of tensor nonzeros to reduce atomic-write contention. Lastly, we test the performance of our implementation on 16- and 68-core Intel CPUs and the K80 and P100 NVIDIA GPUs, showing that we are competitive with state-of-the-art architecture-specific codes while having the advantage of being able to run on a variety of architectures.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

An Algorithm for Synthesizing Mass and Stiffness Matrices from Experimental Vibration Modes

An algorithm is described for synthesizing the mass and stiffness matrices from experimentally derived modal data in a way that preserves the physical significance of the individual mass and stiffness elements. The mass and stiffness matrices are derived for a rollup solar array example, and are then used to define the modal response of a modified array.

Ross, R. G., Jr.↗
NASA NTRS1985

PISCES: An environment for parallel scientific computation

The parallel implementation of scientific computing environment (PISCES) is a project to provide high-level programming environments for parallel MIMD computers. Pisces 1, the first of these environments, is a FORTRAN 77 based environment which runs under the UNIX operating system. The Pisces 1 user programs in Pisces FORTRAN, an extension of FORTRAN 77 for parallel processing. The major emphasis in the Pisces 1 design is in providing a carefully specified virtual machine that defines the run-time environment within which Pisces FORTRAN programs are executed. Each implementation then provides the same virtual machine, regardless of differences in the underlying architecture. The design is intended to be portable to a variety of architectures. Currently Pisces 1 is implemented on a network of Apollo workstations and on a DEC VAX uniprocessor via simulation of the task level parallelism. An implementation for the Flexible Computing Corp. FLEX/32 is under construction. An introduction to the Pisces 1 virtual computer and the FORTRAN 77 extensions is presented. An example of an algorithm for the iterative solution of a system of equations is given. The most notable features of the design are the provision for several granularities of parallelism in programs and the provision of a window mechanism for distributed access to large arrays of data.

Pratt, T. W.↗
NASA NTRS2010

Reconfigurable Hardware for Compressing Hyperspectral Image Data

High-speed, low-power, reconfigurable electronic hardware has been developed to implement ICER-3D, an algorithm for compressing hyperspectral-image data. The algorithm and parts thereof have been the topics of several NASA Tech Briefs articles, including Context Modeler for Wavelet Compression of Hyperspectral Images (NPO-43239) and ICER-3D Hyperspectral Image Compression Software (NPO-43238), which appear elsewhere in this issue of NASA Tech Briefs. As described in more detail in those articles, the algorithm includes three main subalgorithms: one for computing wavelet transforms, one for context modeling, and one for entropy encoding. For the purpose of designing the hardware, these subalgorithms are treated as modules to be implemented efficiently in field-programmable gate arrays (FPGAs). The design takes advantage of industry- standard, commercially available FPGAs. The implementation targets the Xilinx Virtex II pro architecture, which has embedded PowerPC processor cores with flexible on-chip bus architecture. It incorporates an efficient parallel and pipelined architecture to compress the three-dimensional image data. The design provides for internal buffering to minimize intensive input/output operations while making efficient use of offchip memory. The design is scalable in that the subalgorithms are implemented as independent hardware modules that can be combined in parallel to increase throughput. The on-chip processor manages the overall operation of the compression system, including execution of the top-level control functions as well as scheduling, initiating, and monitoring processes. The design prototype has been demonstrated to be capable of compressing hyperspectral data at a rate of 4.5 megasamples per second at a conservative clock frequency of 50 MHz, with a potential for substantially greater throughput at a higher clock frequency. The power consumption of the prototype is less than 6.5 W. The reconfigurability (by means of reprogramming) of the FPGAs makes it possible to effectively alter the design to some extent to satisfy different requirements without adding hardware. The implementation could be easily propagated to future FPGA generations and/or to custom application-specific integrated circuits.

Aranki, Nazeeh↗
DOE OSTI2025

Evaluation and Development of Phase Array Ultrasonic Testing (PAUT) System for Additively Manufactured Parts

This research focuses on the application of advanced ultrasonic testing techniques developed by The Phased Array Company (TPAC) for inspecting defects in additive manufacturing (AM) parts. Traditionally, X-ray computed tomography is the standard for inspecting AM components. Although, the long inspection and analysis time, along with relatively high cost make implementation difficult. Thus, an alternative nondestructive evaluation (NDE) approach is necessary to support quality assurance efforts within the field of AM. TPAC is recognized as a leader in ultrasonic testing innovation, deploying sophisticated algorithms such as Total Focusing Method (TFM) and Phased Wave Imaging (PWI) for ultrasonic data processing and interpretation. This work will explore how the TFM and PWI algorithms can assist defect detection within polymer AM parts. The AM field is seeking novel NDE methods to provide support within quality control and assurance efforts. Advanced ultrasonics inspection have the potential to fulfill this need.

36 MATERIALS SCIENCE↗
NASA NTRS2015

An Integrated 0-1 Hour First-Flash Lightning Nowcasting, Lightning Amount and Lightning Jump Warning Capability

Lightning one of the most dangerous weather-related phenomena, especially as many jobs and activities occur outdoors, presenting risk from a lightning strike. Cloud-to-ground (CG) lightning represents a considerable safety threat to people at airfields, marinas, and outdoor facilities-from airfield personnel, to people attending outdoor stadium events, on beaches and golf courses, to mariners, as well as emergency personnel. Holle et al. (2005) show that 90% of lightning deaths occurred outdoors, while 10% occurred indoors despite the perception of safety when inside buildings. Curran et al. (2000) found that nearly half of fatalities due to weather were related to convective weather in the 1992-1994 timeframe, with lightning causing a large component of the fatalities, in addition to tornadoes and flash flooding. Related to the aviation industry, CG lightning represents a considerable hazard to baggage-handlers, aircraft refuelers, food caterers, and emergency personnel, who all become exposed to the risk of being struck within short time periods while convective storm clouds develop. Airport safety protocols require that ramp operations be modified or discontinued when lightning is in the vicinity (typically 16 km), which becomes very costly and disruptive to flight operations. Therefore, much focus has been paid to nowcasting the first-time initiation and extent of lightning, both of CG and of any lightning (e.g, in-cloud, cloud-to-cloud). For this project three lightning nowcasting methodologies will be combined: (1) a GOESbased 0-1 hour lightning initiation (LI) product (Harris et al. 2010; Iskenderian et al. 2012), (2) a High Resolution Rapid Refresh (HRRR) lightning probability and forecasted lightning flash density product, such that a quantitative amount of lightning (QL) can be assigned to a location of expected LI, and (3) an algorithm that relates Pseudo-GLM data (Stano et al. 2012, 2014) to the so-called "lightning jump" (LJ) methodology (Shultz et al. 2011) to monitor lightning trends and to anticipate/forecast severe weather (hail > or =2.5 cm, winds > or =25 m/s, tornadoes). The result will be a time-continuous algorithm that uses GOES satellite, radar fields, and HRRR model fields to nowcast first-flash LI and QL, and subsequently monitors lightning trends on a perstorm basis within the LJ algorithm for possible severe weather occurrence out to > or =3 hours. The LI-QL-LJ product will also help prepare the operational forecast community for Geostationary Lightning Mapper (GLM) data expected in late 2015, as these data are monitored for ongoing convective storms. The LI-QL-LJ product will first predict where new lightning is highly probable using GOES imagery of developing cumulus clouds, followed by n analysis of NWS (dual-polarization) radar indicators (reflectivity at the -10 C altitude) of lightning occurrence, to increase confidence that LI is immanent. Once lightning is observed, time-continuous lightning mapping array and Pseudo-GLM observations will be analyzed to assess trends and the severe weather threat as identified by trends in lightning (i.e. LJs). Additionally, 5- and 15-min GOES imagery will then be evaluated on a per-storm basis for overshooting and other cloud-top features known to be associated with severe storms. For the processing framework, the GOES-R 0-1 hour convective initiation algorithm's output will be developed within the Warning Decision Support System - Integrated Information (WDSS-II) tracking tool, and merged with radar and lightning (LMA/Psuedo-GLM) datasets for active storms. The initial focus of system development will be over North Alabama for select lightning-active days in summer 2014, yet will be formed in an expandable manner. The lightning alert tool will also be developed in concert with National Weather Service (NWS) forecasters to meet their needs for real-time, accurate first-flash LI and timing, as well as anticipated lightning trends, amounts, continuation and cessation, so to provide key situational awareness and decision support information. The NASA Short-term Prediction Research and Transition (SPoRT) Center will provide important logistical and collaborative support and training, involving interactions with the NWS and broader user community.

Mecikalski, John↗
DOE OSTI2021

Unsupervised Azimuth Estimation of Solar Arrays in Low-Resolution Satellite Imagery through Semantic Segmentation and Hough Transform

This paper explains the use of a convolutional neural network (CNN) to segment solar panels in a satellite image containing solar arrays, and extract associated metadata from the arrays. A novel unsupervised technique is introduced to estimate the azimuth of each individual solar panel from the predicted mask of the convolutional neural network. This pipeline was developed with the aim of extracting necessary metadata for a solar installation, using only a set of latitude–longitude coordinates. Azimuth prediction results for 669 individual solar installations associated with 387 sites located across the United States are provided. A mean average error and median average error of 21.65 degrees and 1.0 degrees were obtained, respectively, when predicting the azimuth of the solar fleet data set, with about 80% of the results within an error of zero degrees of the ground truth azimuth value and about 85% within an error of 25 degrees. The predicted azimuth was then used to estimate the energy conversion of the solar arrays. Results show a 90.9 and 90.6 R-squared value for estimating alternating current (AC) and direct current (DC) energy, respectively, and a mean absolute percentage error (MAPE) of 1.70% in estimating the alternating current (AC) energy using the fully automated algorithm.

14 SOLAR ENERGY↗
DOE OSTI2021

The Illumination of Thunderclouds by Lightning: 4. Volumetric Thunderstorm Imagery

Optical instruments such as the Geostationary Lightning Mapper (GLM) detect lightning based on transient changes in cloud illumination. The horizontal location of lightning is determined from the coordinates of the pixels on the imaging array illuminated during the flash. However, the vertical position of the lightning pulses (approximated by GLM “groups”) below the cloud top cannot be routinely measured from a single space-based instrument. In our prior work, we have developed a machine learning algorithm that can infer optical source altitude for a given pulse based on how the optical energy is distributed across the group footprint and the local Advanced Baseline Imager Cloud-Top Height (CTH). In this fourth part of our thundercloud illumination study, we leverage these source altitudes to generate volumetric GLM imagery of a Colombia thunderstorm. We find that 3D versions of the current GLM meteorological imagery products (that describe thunderstorm kinematics) and thundercloud imagery products (that depict how the flashes appear from space) provide additional insights into lightning activity in the thunderstorm that are lost in the vertical integration used to generate the current 2D GLM gridded products. This new volumetric imaging capability provides a more comprehensive picture of where lightning occurs in the storm, how its physical characteristics vary across three-dimensional space, and how its optical emissions interact with surrounding the cloud medium.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Accelerating Advanced Light Source Science Through Multi-Facility HPC Workflows

Synchrotron light sources support a wide array of techniques to investigate materials, often producing complex, high-volume data that challenge traditional workflows. At the Advanced Light Source (ALS), we developed infrastructure to move microtomography data over ESnet to ALCF and NERSC, where CPU- and GPU-based algorithms generate 3D reconstructed volumes of experimental samples. We employ two data movement and reconstruction models: real-time processing as data streams directly to NERSC compute nodes, and automated file transfer to NERSC and ALCF file systems. The streaming pipeline provides users with feedback in under ten seconds, while the file-based workflow produces high-quality reconstructions suitable for deeper analysis in 20-30 minutes. This infrastructure enables users to utilize HPC resources without direct access to backend systems. We plan to extend this architecture to more endstations, supporting our beamline scientists and users.

Abramov, David↗
DOE OSTI2019

Neutron Scattering Instrumentation Research and Development for High Spatial and Temporal Resolution Imaging at Oak Ridge National Laboratory

Neutron scattering facilities like those available at ORNL are enabling the discovery and development of advanced materials include those to be used for next generation energy storage systems. In spite of all the success that neutron imaging has enjoyed in recent years, even higher-impact research is limited by the current temporal and spatial resolution of neutron detection devices. In particular, the spatial resolution of all neutron imaging methods in use and under development is fundamentally limited by the variance introduced by the charged particles emitted from neutron absorption. Thus, we are investigating, designing and demonstrating an array of ~1 micron-sized glass scintillating fibers arranged via photonics-crystal-fiber patterning techniques-- with particle tracking (a must) so that position resolution can be finer than than that allowed by the charged particle ranges. The processes, knowledge and algorithms developed will allow for spatially (1 micron) and time resolved (<100 ns), high throughput neutron transmission imaging in large areas at reasonable costs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1988

Digital and optical shape representation and pattern recognition; Proceedings of the Meeting, Orlando, FL, Apr. 4-6, 1988

The present conference discusses topics in pattern-recognition correlator architectures, digital stereo systems, geometric image transformations and their applications, topics in pattern recognition, filter algorithms, object detection and classification, shape representation techniques, and model-based object recognition methods. Attention is given to edge-enhancement preprocessing using liquid crystal TVs, massively-parallel optical data base management, three-dimensional sensing with polar exponential sensor arrays, the optical processing of imaging spectrometer data, hybrid associative memories and metric data models, the representation of shape primitives in neural networks, and the Monte Carlo estimation of moment invariants for pattern recognition.

Juday, Richard D.↗
NASA NTRS1988

Recursive least squares estimation and Kalman filtering by systolic arrays

One of the most promising new directions for high-throughput-rate problems is that based on systolic arrays. In this paper, using the matrix-decomposition approach, a systolic Kalman filter is formulated as a modified square-root information filter consisting of a whitening filter followed by a simple least-squares operation based on the systolic QR algorithm. By proper skewing of the input data, a fully pipelined time and measurement update systolic Kalman filter can be achieved with O(n squared) processing cells, resulting in a system throughput rate of O (n).

Chen, M. J.↗
NASA NTRS1989

A temperature and pressure controlled calibration system for pressure sensors

A data acquisition and experiment control system capable of simulating temperatures from -184 to +220 C and pressures either absolute or differential from 0 to 344.74 kPa is developed to characterize silicon pressure sensor response to temperature and pressure. System software is described that includes sensor data acquisition, algorithms for numerically derived thermal offset and sensitivity correction, and operation of the environmental chamber and pressure standard. This system is shown to be capable of computer interfaced cryogenic testing to within 1 C and 34.47 Pa of single channel or multiplexed arrays of silicon pressure sensors.

Chapman, John J.↗
NASA NTRS1992

Real time thermal imaging for analysis and control of crystal growth by the Czochralski technique

A real time thermal imaging system with temperature resolution better than +/- 0.5 C and spatial resolution of better than 0.5 mm has been developed. It has been applied to the analysis of melt surface thermal field distributions in both Czochralski and liquid encapsulated Czochralski growth configurations. The sensor can provide single/multiple point thermal information; a multi-pixel averaging algorithm has been developed which permits localized, low noise sensing and display of optical intensity variations at any location in the hot zone as a function of time. Temperature distributions are measured by extraction of data along a user selectable linear pixel array and are simultaneously displayed, as a graphic overlay, on the thermal image.

Wargo, M. J.↗
NASA NTRS1991

Acquisition and tracking performance measurements for a high speed area array detector system

A proof-of-concept (POC) demonstration system has been developed which demonstrates acquisition, tracking and point-ahead angle sensing for a space optical communications terminal utilizing a single high speed area array detector. The detector is the 128 x 128 pixel Kodak HS-40 photodiode array. It has 64 parallel readout channels and can operate at frames rates up to 40,000 frames/sec with rms readout noise of 20 photoelectrons. A windowing scheme and special purpose digital signal processing electronics are employed to implement acquisition and tracking algorithms. The system operates at greater than 1 kHz sample (frame) rates. Acquisition can be performed in as little as 30 milliseconds with less than 1 picowatt of 0.85 micron beacon power on the detector. At the same power level, the rms tracking accuracy is approximately 1/16 pixel. Results of system analysis and measurements using the POC system are presented.

Short, R. C.↗
NASA NTRS2007

Digital Signal Processing Techniques for the GIFTS SM EDU

The Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) Sensor Module (SM) Engineering Demonstration Unit (EDU) is a high resolution spectral imager designed to measure infrared (IR) radiance using a Fourier transform spectrometer (FTS). The GIFTS instrument employs three Focal Plane Arrays (FPAs), which gather measurements across the long-wave IR (LWIR), short/mid-wave IR (SMWIR), and visible spectral bands. The raw interferogram measurements are radiometrically and spectrally calibrated to produce radiance spectra, which are further processed to obtain atmospheric profiles via retrieval algorithms. This paper describes several digital signal processing (DSP) techniques involved in the development of the calibration model. In the first stage, the measured raw interferograms must undergo a series of processing steps that include filtering, decimation, and detector nonlinearity correction. The digital filtering is achieved by employing a linear-phase even-length FIR complex filter that is designed based on the optimum equiripple criteria. Next, the detector nonlinearity effect is compensated for using a set of pre-determined detector response characteristics. In the next stage, a phase correction algorithm is applied to the decimated interferograms. This is accomplished by first estimating the phase function from the spectral phase response of the windowed interferogram, and then correcting the entire interferogram based on the estimated phase function. In the calibration stage, we first compute the spectral responsivity based on the previous results and the ideal Planck blackbody spectra at the given temperatures, from which, the calibrated ambient blackbody (ABB), hot blackbody (HBB), and scene spectra can be obtained. In the post-calibration stage, we estimate the Noise Equivalent Spectral Radiance (NESR) from the calibrated ABB and HBB spectra. The NESR is generally considered as a measure of the instrument noise performance, and can be estimated as the standard deviation of calibrated radiance spectra from multiple scans. To obtain an estimate of the FPA performance, we developed an efficient method of generating pixel performance assessments. In addition, a random pixel selection scheme is developed based on the pixel performance evaluation. This would allow us to perform the calibration procedures on a random pixel population that is a good statistical representation of the entire FPA. The design and implementation of each individual component will be discussed in details.

Tian, Jialin↗
NASA NTRS2009

A Survey of New Trends in Symbolic Execution for Software Testing and Analysis

Symbolic execution is a well-known program analysis technique which represents values of program inputs with symbolic values instead of concrete (initialized) data and executes the program by manipulating program expressions involving the symbolic values. Symbolic execution has been proposed over three decades ago but recently it has found renewed interest in the research community, due in part to the progress in decision procedures, availability of powerful computers and new algorithmic developments. We provide a survey of some of the new research trends in symbolic execution, with particular emphasis on applications to test generation and program analysis. We first describe an approach that handles complex programming constructs such as input data structures, arrays, as well as multi-threading. We follow with a discussion of abstraction techniques that can be used to limit the (possibly infinite) number of symbolic configurations that need to be analyzed for the symbolic execution of looping programs. Furthermore, we describe recent hybrid techniques that combine concrete and symbolic execution to overcome some of the inherent limitations of symbolic execution, such as handling native code or availability of decision procedures for the application domain. Finally, we give a short survey of interesting new applications, such as predictive testing, invariant inference, program repair, analysis of parallel numerical programs and differential symbolic execution.

Pasareanu, Corina S.↗
NASA NTRS2019

Optical Geolocation for Small Unmanned Aerial Systems

This paper presents an airborne optical geolocation system using four optical targets to provide position and attitude estimation for a sUAS supporting the NASA Acoustic Research Mission (ARM), where the goal is to reduce nuisance airframe noise during approach and landing. A large precision positioned microphone array captures the airframe noise for multiple passes of a Gulfstream III aircraft. For health monitoring of the microphone array, the Acoustic Calibration Vehicle (ACV) sUAS completes daily flights with an onboard speaker emitting tones at frequencies optimized for determining microphone functionality. An accurate position estimate of the ACV relative to the array is needed for microphone health monitoring. To this end, an optical geolocation system using a downward facing camera mounted to the ACV was developed. The 3D positioning of the ACV is computed using the pinhole camera model. A novel optical geolocation algorithm first detects the targets, then a recursive algorithm tightens the localization of the targets. Finally, the position of the sUAS is computed using the image coordinates of the targets, the 3D world coordinates of the targets, and the camera matrix. A Real-Time Kinematic GPS system is used to compare the optical geolocation system.

Dolph, Chester V.↗