Improved maximum entropy analysis with an extended search space
Improved maximum entropy analysis with an extended search space
复制标题
DOI:
10.1016/j.jcp.2012.12.023
复制
发表时间:
2011-10
期刊:
影响因子:
--
通讯作者:
A. Rothkopf
中科院分区:
文献类型:
--
作者:
A. Rothkopf
The standard implementation of the Maximum Entropy Method (MEM) follows Bryan [1] and deploys a Singular Value Decomposition (SVD) to limit the dimensionality of the underlying solution space apriori. Here we present arguments based on the shape of the SVD basis functions and numerical evidence from a mock data analysis, which show that the correct Bayesian solution is not in general recovered with this approach. As a remedy we propose to extend the search basis systematically, which will eventually recover the full solution space and the correct solution. In order to adequately approach problems where an exponentially damped kernel is used, we provide an open-source implementation, using the C/C++ language that utilizes high precision arithmetic adjustable at run-time [2]. The LBFGS algorithm is included in the code in order to attack problems without the need to resort to a particular search space restriction.