The use of the L-curve in the regularization of discrete ill-posed problems
SIAM Journal on Scientific Computing
Using the L-curve for determining optimal regularization parameters
Numerische Mathematik
Comparison of approximate methods for handling hyperparameters
Neural Computation
Evaluating derivatives: principles and techniques of algorithmic differentiation
Evaluating derivatives: principles and techniques of algorithmic differentiation
Inverse Problem Theory and Methods for Model Parameter Estimation
Inverse Problem Theory and Methods for Model Parameter Estimation
Gaussian Processes for Machine Learning (Adaptive Computation and Machine Learning)
Gaussian Processes for Machine Learning (Adaptive Computation and Machine Learning)
Coarse-gradient Langevin algorithms for dynamic data integration and uncertainty quantification
Journal of Computational Physics - Special issue: Uncertainty quantification in simulation science
Reservoir Simulation: Mathematical Techniques in Oil Recovery (CBMS-NSF Regional Conference Series in Applied Mathematics)
Optimal Choice of AR and MA Parts in Autoregressive Moving Average Models
IEEE Transactions on Pattern Analysis and Machine Intelligence
Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images
IEEE Transactions on Pattern Analysis and Machine Intelligence
An iterative stochastic ensemble method for parameter estimation of subsurface flow models
Journal of Computational Physics
Hi-index | 31.45 |
A Hybrid Nested Sampling (HNS) algorithm is proposed for efficient Bayesian model calibration and prior model selection. The proposed algorithm combines, Nested Sampling (NS) algorithm, Hybrid Monte Carlo (HMC) sampling and gradient estimation using Stochastic Ensemble Method (SEM). NS is an efficient sampling algorithm that can be used for Bayesian calibration and estimating the Bayesian evidence for prior model selection. Nested sampling has the advantage of computational feasibility. Within the nested sampling algorithm, a constrained sampling step is performed. For this step, we utilize HMC to reduce the correlation between successive sampled states. HMC relies on the gradient of the logarithm of the posterior distribution, which we estimate using a stochastic ensemble method based on an ensemble of directional derivatives. SEM only requires forward model runs and the simulator is then used as a black box and no adjoint code is needed. The developed HNS algorithm is successfully applied for Bayesian calibration and prior model selection of several nonlinear subsurface flow problems.