new feature request: iterative Gaussian fitting

[kmzzhang]

Fastest inference for ultra precise problem is to fit multivariate Gaussian to reweighted surrogate posterior