Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Efficient and robust gradient enhanced Kriging emulators.

Technical Report ·
DOI:https://doi.org/10.2172/1096451· OSTI ID:1096451
%E2%80%9CNaive%E2%80%9D or straight-forward Kriging implementations can often perform poorly in practice. The relevant features of the robustly accurate and efficient Kriging and Gradient Enhanced Kriging (GEK) implementations in the DAKOTA software package are detailed herein. The principal contribution is a novel, effective, and efficient approach to handle ill-conditioning of GEK's %E2%80%9Ccorrelation%E2%80%9D matrix, RN%CC%83, based on a pivoted Cholesky factorization of Kriging's (not GEK's) correlation matrix, R, which is a small sub-matrix within GEK's RN%CC%83 matrix. The approach discards sample points/equations that contribute the least %E2%80%9Cnew%E2%80%9D information to RN%CC%83. Since these points contain the least new information, they are the ones which when discarded are both the easiest to predict and provide maximum improvement of RN%CC%83's conditioning. Prior to this work, handling ill-conditioned correlation matrices was a major, perhaps the principal, unsolved challenge necessary for robust and efficient GEK emulators. Numerical results demonstrate that GEK predictions can be significantly more accurate when GEK is allowed to discard points by the presented method. Numerical results also indicate that GEK can be used to break the curse of dimensionality by exploiting inexpensive derivatives (such as those provided by automatic differentiation or adjoint techniques), smoothness in the response being modeled, and adaptive sampling. Development of a suitable adaptive sampling algorithm was beyond the scope of this work; instead adaptive sampling was approximated by omitting the cost of samples discarded by the presented pivoted Cholesky approach.
Research Organization:
Sandia National Laboratories (SNL-NM), Albuquerque, NM (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA)
DOE Contract Number:
AC04-94AL85000
OSTI ID:
1096451
Report Number(s):
SAND2013-7022; 471013
Country of Publication:
United States
Language:
English

Similar Records

Modified Cholesky factorizations in interior-point algorithms for linear programming.
Journal Article · Thu Dec 31 23:00:00 EST 1998 · SIAM J. Optimization · OSTI ID:937894

A novel method for modeling the neutron time of flight (nTOF) detector response in current mode to inertial confinement fusion experiments.
Technical Report · Sun Sep 01 00:00:00 EDT 2013 · OSTI ID:1096470

Kriging and its relation to least squares
Technical Report · Wed Oct 31 23:00:00 EST 1984 · OSTI ID:7204187

Related Subjects