Search NASA⌕ Search

SEARCH · Search NASA

Results for “Randomized methods”

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 343 records · Page 19

A Novel Pruning Method for Convolutional Neural Networks Based off Identifying Critical Filters

Convolutional Neural Networks (CNNs) are one of the most extensively used tools in machine learning, but they are still not well understood and in many cases they are over-parameterized, leading to slow inference and impeding their deployment on low-power devices. In the last few years, many methods for decreasing the number of parameters in a network by pruning its output channels have been suggested, but a very recent work has argued that random pruning of channels performs on-par with state-of-the-art pruning methods. While random and other pruning methods might be effectively used for lowering the number of parameters in a CNN, none of these methods can be used to gain any further understanding of the model that the CNN has built. In this work, we propose a novel method for pruning a network, that at the same time can lead to a better understanding of what the individual filters of the network learn about the data. The method proposed aims to keep only the filters that are "important" for a class. We define a filter as important for a class if its removal has the highest negative impact on the accuracy for that class. We demonstrate that our method is better than random pruning on two networks used on the EMNIST and CIFAR10 datasets. By analyzing the important filters, we find that the important filters in the pruned networks learn features which are more general across classes. We demonstrate the importance and applicability of that observation in two transfer-learning tasks.

Dimovska, Mihaela↗

Theoretically informed Monte Carlo simulation of liquid crystals by sampling of alignment-tensor fields

A theoretically informed coarse-grained Monte Carlo method is proposed for studying liquid crystals. The free energy functional of the system is described in the framework of the Landau-de Gennes formalism. The alignment field and its gradients are approximated by finite differences, and the free energy is minimized through a stochastic sampling technique. The validity of the proposed method is established by comparing the results of the proposed approach to those of traditional free energy minimization techniques. Its usefulness is illustrated in the context of three systems, namely, a nematic liquid crystal confined in a slit channel, a nematic liquid crystal droplet, and a chiral liquid crystal in the bulk. It is found that for systems that exhibit multiple metastable morphologies, the proposed Monte Carlo method is generally able to identify lower free energy states that are often missed by traditional approaches. Importantly, the Monte Carlo method identifies such states from random initial configurations, thereby obviating the need for educated initial guesses that can be difficult to formulate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The isotropy problem of sub-ankle ultra high energy cosmic rays

We study the time dependent propagation of sub-ankle ultra high energy cosmic rays (UHECRs) originating from point-like Galactic sources. We show that drift in the Galactic magnetic field (GMF) may play an important role in the propagation of UHECRs and their measured anisotropy, particularly when the transport is anisotropic. To fully account for the discreteness of UHECR sources in space and time, a Monte Carlo method is used to randomly place sources in the Galaxy. The low anisotropy measured by Auger is not generally characteristic of the theoretical models, given that the sources are distributed in proportion to the star formation rate, but it can possibly be understood as (1) intermittency effects due to the discrete nature of the sources or, with extreme parameters, (2) a cancellation of drift current along a current sheet with outward radial diffusive flux. We conclude that it is possible to interpret the Galactic sub-ankle CR flux as being due entirely to intermittent discrete Galactic sources distributed in proportion to star formation, but only with a probability of roughly 35%, of which the spectrum is in accord with observations about 30% of the time. An alternative explanation for the low anisotropy may be that they are mostly extragalactic and/or heavy.

79 ASTRONOMY AND ASTROPHYSICS↗

Preserving correlations between trajectories for efficient path sampling

Importance sampling of trajectories has proved a uniquely successful strategy for exploring rare dynamical behaviors of complex systems in an unbiased way. Carrying out this sampling, however, requires an ability to propose changes to dynamical pathways that are substantial, yet sufficiently modest to obtain reasonable acceptance rates. Satisfying this requirement becomes very challenging in the case of long trajectories, due to the characteristic divergences of chaotic dynamics. Here, we examine schemes for addressing this problem, which engineer correlation between a trial trajectory and its reference path, for instance using artificial forces. Our analysis is facilitated by a modern perspective on Markov chain Monte Carlo sampling, inspired by non-equilibrium statistical mechanics, which clarifies the types of sampling strategies that can scale to long trajectories. Viewed in this light, the most promising such strategy guides a trial trajectory by manipulating the sequence of random numbers that advance its stochastic time evolution, as done in a handful of existing methods. In cases where this “noise guidance” synchronizes trajectories effectively, as the Glauber dynamics of a two-dimensional Ising model, we show that efficient path sampling can be achieved for even very long trajectories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparison of Machine Learning Algorithms for Natural Gas Identification with Mixed Potential Electrochemical Sensor Arrays

Mixed-potential electrochemical sensor arrays consisting of indium tin oxide (ITO), La 0.87 Sr 0.13 CrO 3 , Au, and Pt electrodes can detect the leaks from natural gas infrastructure. Algorithms are needed to correctly identify natural gas sources from background natural and anthropogenic sources such as wetlands or agriculture. We report for the first time a comparison of several machine learning methods for mixture identification in the context of natural gas emissions monitoring by mixed potential sensor arrays. Random Forest, Artificial Neural Network, and Nearest Neighbor methods successfully classified air mixtures containing only CH 4 , two types of natural gas simulants, and CH 4 +NH 3 with >98% identification accuracy. The model complexity of these methods were optimized and the degree of robustness against overfitting was determined. Finally, these methods are benchmarked on both desktop PC and single-board computer hardware to simulate their application in a portable internet-of-things sensor package. The combined results show that the random forest method is the preferred method for mixture identification with its high accuracy (>98%), robustness against overfitting with increasing model complexity, and had less than 10 ms training time and less than 0.1 ms inference time on single-board computer hardware.

03 NATURAL GAS↗

WE-AB-207A-04: Random Undersampled Cone Beam CT: Theoretical Analysis and a Novel Reconstruction Method

Purpose: Reducing x-ray exposure and speeding up data acquisition motived studies on projection data undersampling. It is an important question that for a given undersampling ratio, what the optimal undersampling approach is. In this study, we propose a new undersampling scheme: random-ray undersampling. We will mathematically analyze its projection matrix properties and demonstrate its advantages. We will also propose a new reconstruction method that simultaneously performs CT image reconstruction and projection domain data restoration. Methods: By representing projection operator under the basis of singular vectors of full projection operator, matrix representations for an undersampling case can be generated and numerical singular value decomposition can be performed. We compared properties of matrices among three undersampling approaches: regular-view undersampling, regular-ray undersampling, and the proposed random-ray undersampling. To accomplish CT reconstruction for random undersampling, we developed a novel method that iteratively performs CT reconstruction and missing projection data restoration via regularization approaches. Results: For a given undersampling ratio, random-ray undersampling preserved mathematical properties of full projection operator better than the other two approaches. This translates to advantages of reconstructing CT images at lower errors. Different types of image artifacts were observed depending on undersampling strategies, which were ascribed to the unique singular vectors of the sampling operators in the image domain. We tested the proposed reconstruction algorithm on a Forbid phantom with only 30% of the projection data randomly acquired. Reconstructed image error was reduced from 9.4% in a TV method to 7.6% in the proposed method. Conclusion: The proposed random-ray undersampling is mathematically advantageous over other typical undersampling approaches. It may permit better image reconstruction at the same undersampling ratio. The novel algorithm suitable for this random-ray undersampling was able to reconstruct high-quality images.

60 APPLIED LIFE SCIENCES↗

Incorporating Experimental Information in the Total Monte Carlo Methodology Using File Weights

Some criticism has been directed towards the Total Monte Carlo method because experimental information has not been taken into account in a statistically well-founded manner. In this work, a Bayesian calibration method is implemented by assigning weights to the random nuclear data files and the method is illustratively applied to a few applications. In some considered cases, the estimated nuclear data uncertainties are significantly reduced and the central values are significantly shifted. The study suggests that the method can be applied both to estimate uncertainties in a more justified way and in the search for better central values. Some improvements are however necessary; for example, the treatment of outliers and cross-experimental correlations should be more rigorous and random files that are intended to be prior files should be generated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

UQ Toolkit v. 3.0

The Uncertainty Quantification (UQ) Toolkit is a software library for the characterization and propagation of uncertainties in computational models. This library provides Bayesian inference tools to infer uncertain model parameters, Bayesian compressive sensing methods for discovering sparse representations of high-dimensional input-output response surfaces, methods for constructing Karhunen-Loève representations of stochastic processes, and global sensitivity analysis tools used to compute Sobol indices in order to characterize the importance of uncertain inputs or the interactions between them. The basis for many of these methods relies on representing random variables with Polynomial Chaos expansions (PCEs). This library implements several spectral basis function types (e.g. Hermite basis functions in terms of Gaussian random variables or Legendre basis functions in terms of uniform random variables) that can be used to represent random variables with PCEs. For the propagation of uncertainty, tools are provided to propagate PCEs that describe the input uncertainty through the computational model using either intrusive methods (Galerkin projection of equations onto basis functions) or non-intrusive methods (perform deterministic operation at sampled values of the random values and project the obtained results onto basis functions).

Sargsyan, Khachik↗

A method for solving stochastic equations by reduced order models and local approximations

A method is proposed for solving equations with random entries, referred to as stochastic equations (SEs). The method is based on two recent developments. The first approximates the response surface giving the solution of a stochastic equation as a function of its random parameters by a finite set of hyperplanes tangent to it at expansion points selected by geometrical arguments. The second approximates the vector of random parameters in the definition of a stochastic equation by a simple random vector, referred to as stochastic reduced order model (SROM), and uses it to construct a SROM for the solution of this equation. The proposed method is a direct extension of these two methods. It uses SROMs to select expansion points, rather than selecting these points by geometrical considerations, and represents the solution by linear and/or higher order local approximations. The implementation and the performance of the method are illustrated by numerical examples involving random eigenvalue problems and stochastic algebraic/differential equations. The method is conceptually simple, non-intrusive, efficient relative to classical Monte Carlo simulation, accurate, and guaranteed to converge to the exact solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multilevel Hierarchical Decomposition of Finite Element White Noise with Application to Multilevel Markov Chain Monte Carlo

In this work we develop a new hierarchical multilevel approach to generate Gaussian random field realizations in an algorithmically scalable manner that is well suited to incorporating into multilevel Markov chain Monte Carlo (MCMC) algorithms. This approach builds off of other partial differential equation (PDE) approaches for generating Gaussian random field realizations; in particular, a single field realization may be formed by solving a reaction-diffusion PDE with a spatial white noise source function as the right-hand side. While these approaches have been explored to accelerate forward uncertainty quantification tasks, e.g., multilevel Monte Carlo, the previous constructions are not directly applicable to multilevel MCMC frameworks which build fine-scale random fields in a hierarchical fashion from coarse-scale random fields. Our new hierarchical multilevel method relies on a hierarchical decomposition of the white noise source function in $L^2$ which allows us to form Gaussian random field realizations across multiple levels of discretization in a way that fits into multilevel MCMC algorithmic frameworks. After presenting our main theoretical results and numerical scaling results to showcase the utility of this new hierarchical PDE method for generating Gaussian random field realizations, this method is tested on a four-level MCMC algorithm to explore its feasibility.

algebraic multigrid↗

Acute Toxicity and Quality of Life in Patients With Prostate Cancer Treated With Protons or Carbon Ions in a Prospective Randomized Phase II Study—The IPI Trial

Purpose: The purpose of this study was to compare safety and feasibility of proton therapy with that of carbon ion therapy in hypofractionated raster-scanned irradiation of the prostate, in a prospective randomized phase 2 trial. Methods and Materials: In this trial, 92 patients with localized prostate cancer were enrolled. Patients were randomized to receive either proton therapy (arm A) or carbon ion therapy (arm B) and treated with a total dose of 66 Gy(relative biological effectiveness [RBE]) administered in 20 fractions (single dose of 3.3 Gy[RBE]). Patients were stratified by the use of antihormone therapy. Primary endpoint was the combined assessment of safety and feasibility. Secondary endpoints were specific toxicities, prostate-specific antigen progression-free survival (PFS), overall survival (OS), and quality of life (QoL). Results: Ninety-one patients completed therapy and have had a median follow-up of 22.3 months. Among acute genitourinary toxicities, grade 1 cystitis rates were 34.1% (39.1% in A; 28.9% in B) and 17.6% grade 2 (21.7% in A; 13.3% in B). Seven patients (8%) required urinary catheterization during treatment due to urinary retention, 5 of whom were in arm A. Regarding acute gastrointestinal toxicities, 2 patients treated with protons developed grade 3 rectal fistulas. Grade 1 radiation proctitis occurred in 12.1% (13.0% in A; 11.1% in B) and grade 2 in 5.5% (8.7% in A; 2.2% in B). No statistically significant differences in toxicity profiles between arms were found. Reduced QoL was evident mainly in fatigue, pain, and urinary symptoms during therapy and 6 weeks thereafter. All European Organization for Research and Treatment of Cancer QLQ-C30 and -PR25 scores improved during follow-up. Conclusions: Hypofractionated irradiation using either carbon ions or protons results in comparable acute toxicities and QoL parameters. We found that hypofractionated particle irradiation is feasible and may be safe. Due to the occurrence of gel in the rectal wall and the consecutive occurrence of 2 rectal fistulas, we stopped using the insertion of spacer gel. Longer follow-up is necessary for evaluation of PFS and OS. (Ion Prostate Irradiation (IPI); (NCT01641185); (ClinicalTrials.gov).)

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Stochastic Analysis of Levee Stability Subject to Variable Seepage Conditions

© ASCE. Levee stability is highly influenced by seepage. Specifically, hydraulic conductivity distribution within a levee influences pore pressure distribution and controls the drained strength of the soil. In this study the influence of hydraulic conductivity and blanket layer thickness on failure probability is evaluated within the context of reliability analyses that also include soil strength and unit weight as random variables. First-Order Reliability Method (FORM) is used to evaluate reliability, rank random variables by importance and to obtain sensitivity of the solution to each random variable and its distribution parameters. Stability is computed using Spencer's method of slices coupled to a finite element seepage code to directly evaluate pore pressure. In addition, response surface method solutions are compared to the direct reliability solution to assess response surface accuracy. The results show that blanket layer thickness is more important than most strength parameters and that the uncertainty in hydraulic conductivity is less important to the variance in safety factor for the case of a low-permeability blanket layer. Numerical challenges caused by implementing a finite element limit-state function are discussed. Response surface methods are found to give a reasonable approximation to the direct reliability solution when the design point is between response surface fitting points.

Lanzafame, Robert↗

Technical Note: Introduction of variance component analysis to setup error analysis in radiotherapy

Purpose: The purpose of this technical note is to introduce variance component analysis to the estimation of systematic and random components in setup error of radiotherapy. Methods: Balanced data according to the one-factor random effect model were assumed. Results: Analysis-of-variance (ANOVA)-based computation was applied to estimate the values and their confidence intervals (CIs) for systematic and random errors and the population mean of setup errors. The conventional method overestimates systematic error, especially in hypofractionated settings. The CI for systematic error becomes much wider than that for random error. The ANOVA-based estimation can be extended to a multifactor model considering multiple causes of setup errors (e.g., interpatient, interfraction, and intrafraction). Conclusions: Variance component analysis may lead to novel applications to setup error analysis in radiotherapy.

60 APPLIED LIFE SCIENCES↗