A nonlinear preconditioner for optimum experimental design problems

Abstract

We show how to efficiently compute A-optimal experimental designs, which are formulated in terms of the minimization of the trace of the covariance matrix of the underlying regression process, using quasi-Newton sequential quadratic programming methods. In particular, we introduce a modification of the problem that leads to significantly faster convergence. To derive this modification, we model each iteration in terms of an initial experimental design that is to be improved, and show that the absolute condition number of the model problem grows without bounds as the quality of the initial design improves. As a remedy, we devise a preconditioner that ensures that the absolute condition number will instead stay uniformly bounded. Using numerical experiments, we study the effect of this reformulation on relevant cases of the general problem class, and find that it leads to significant improvements in both stability and convergence behavior.

Keywords

• 90C30
• 90C55
• 62K99