Skip to content
Snippets Groups Projects
Commit 83d9e58f authored by Philipp Arras's avatar Philipp Arras
Browse files

Use traditional settings for getting started. No preconditioning here

Although it is working
parent fc889da1
Branches
Tags
1 merge request!333Operator spectra
Pipeline #61919 passed
...@@ -103,8 +103,7 @@ if __name__ == '__main__': ...@@ -103,8 +103,7 @@ if __name__ == '__main__':
data = signal_response(mock_position) + N.draw_sample() data = signal_response(mock_position) + N.draw_sample()
# Minimization parameters # Minimization parameters
ic_sampling = ift.AbsDeltaEnergyController(0.5, convergence_level=5, ic_sampling = ift.GradientNormController(iteration_limit=100)
iteration_limit=100)
ic_newton = ift.GradInfNormController( ic_newton = ift.GradInfNormController(
name='Newton', tol=1e-7, iteration_limit=35) name='Newton', tol=1e-7, iteration_limit=35)
minimizer = ift.NewtonCG(ic_newton) minimizer = ift.NewtonCG(ic_newton)
...@@ -129,7 +128,7 @@ if __name__ == '__main__': ...@@ -129,7 +128,7 @@ if __name__ == '__main__':
# Draw new samples to approximate the KL five times # Draw new samples to approximate the KL five times
for i in range(5): for i in range(5):
# Draw new samples and minimize KL # Draw new samples and minimize KL
KL = ift.MetricGaussianKL(mean, H, N_samples, napprox=20) KL = ift.MetricGaussianKL(mean, H, N_samples)
KL, convergence = minimizer(KL) KL, convergence = minimizer(KL)
mean = KL.position mean = KL.position
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment