Picture of boy being examining by doctor at a tuberculosis sanatorium

Understanding our future through Open Access research about our past...

Strathprints makes available scholarly Open Access content by researchers in the Centre for the Social History of Health & Healthcare (CSHHH), based within the School of Humanities, and considered Scotland's leading centre for the history of health and medicine.

Research at CSHHH explores the modern world since 1800 in locations as diverse as the UK, Asia, Africa, North America, and Europe. Areas of specialism include contraception and sexuality; family health and medical services; occupational health and medicine; disability; the history of psychiatry; conflict and warfare; and, drugs, pharmaceuticals and intoxicants.

Explore the Open Access research of the Centre for the Social History of Health and Healthcare. Or explore all of Strathclyde's Open Access research...

Image: Heart of England NHS Foundation Trust. Wellcome Collection - CC-BY.

Efficient Gaussian process based on BFGS updating and logdet approximation

Leithead, W.E. and Zhang, Y. and Leith, D.J. (2005) Efficient Gaussian process based on BFGS updating and logdet approximation. In: Proceedings of the 16th IFAC World Congress, 2005. UNSPECIFIED, p. 217. ISBN 978-3-902661-75-3

Full text not available in this repository. Request a copy from the Strathclyde author

Abstract

Gaussian process (GP) is a Bayesian nonparametric regression model, showing good performance in various applications. However, its hyperparameterestimation procedure suffers from numerous covariance-matrix inversions of prohibitively O(N3) operations. In this paper, we propose using the quasi-Newton BFGS O(N2)-operation formula to update recursively the inverse of covariance matrix at every iteration. As for the involved log det computation, a power-series expansion based approximation and compensation scheme is proposed with only 50N2 operations. A number of numerical tests are performed based on the 2D- sinusoidal regression example and the Wiener-Hammerstein identification example. It is shown that by using the proposed implementation, more than 80% O(N3) operations are eliminated, and the speedup of 5 - 9 can be achieved.