Skip to content
Snippets Groups Projects
Commit 9f471cba authored by Martin Reinecke's avatar Martin Reinecke
Browse files

adjustments

parent d490d56c
Branches
Tags
No related merge requests found
Pipeline #
...@@ -41,7 +41,7 @@ if __name__ == "__main__": ...@@ -41,7 +41,7 @@ if __name__ == "__main__":
# Setting up the noise covariance and drawing a random noise realization # Setting up the noise covariance and drawing a random noise realization
noiseless_data = R(mock_signal) noiseless_data = R(mock_signal)
noise_amplitude = noiseless_data.std()/signal_to_noise noise_amplitude = noiseless_data.val.std()/signal_to_noise
N = ift.DiagonalOperator( N = ift.DiagonalOperator(
ift.Field.full(data_domain, noise_amplitude**2)) ift.Field.full(data_domain, noise_amplitude**2))
noise = ift.Field.from_random( noise = ift.Field.from_random(
...@@ -56,9 +56,9 @@ if __name__ == "__main__": ...@@ -56,9 +56,9 @@ if __name__ == "__main__":
inverter = ift.ConjugateGradient(controller=ctrl) inverter = ift.ConjugateGradient(controller=ctrl)
energy = ift.library.LogNormalWienerFilterEnergy(m0, data, R, energy = ift.library.LogNormalWienerFilterEnergy(m0, data, R,
N, S, inverter=inverter) N, S, inverter=inverter)
# minimizer = ift.VL_BFGS(controller=ctrl2, max_history_length=20) #minimizer = ift.VL_BFGS(controller=ctrl2, max_history_length=20)
minimizer = ift.RelaxedNewton(controller=ctrl2) minimizer = ift.RelaxedNewton(controller=ctrl2)
# minimizer = ift.SteepestDescent(controller=ctrl2) #minimizer = ift.SteepestDescent(controller=ctrl2)
me = minimizer(energy) me = minimizer(energy)
m = ht(me[0].position) m = ht(me[0].position)
......
...@@ -81,7 +81,7 @@ if __name__ == "__main__": ...@@ -81,7 +81,7 @@ if __name__ == "__main__":
data_domain = R.target data_domain = R.target
noiseless_data = R(mock_signal) noiseless_data = R(mock_signal)
noise_amplitude = noiseless_data.std()/signal_to_noise noise_amplitude = noiseless_data.val.std()/signal_to_noise
# Setting up the noise covariance and drawing a random noise realization # Setting up the noise covariance and drawing a random noise realization
ndiag = ift.Field.full(data_domain, noise_amplitude**2) ndiag = ift.Field.full(data_domain, noise_amplitude**2)
N = ift.DiagonalOperator(ndiag) N = ift.DiagonalOperator(ndiag)
......
...@@ -42,7 +42,7 @@ if __name__ == "__main__": ...@@ -42,7 +42,7 @@ if __name__ == "__main__":
data_domain = R.target[0] data_domain = R.target[0]
noiseless_data = R(mock_signal) noiseless_data = R(mock_signal)
noise_amplitude = noiseless_data.std()/signal_to_noise noise_amplitude = noiseless_data.val.std()/signal_to_noise
# Setting up the noise covariance and drawing a random noise realization # Setting up the noise covariance and drawing a random noise realization
ndiag = ift.Field.full(data_domain, noise_amplitude**2) ndiag = ift.Field.full(data_domain, noise_amplitude**2)
N = ift.DiagonalOperator(ndiag) N = ift.DiagonalOperator(ndiag)
......
...@@ -42,7 +42,7 @@ if __name__ == "__main__": ...@@ -42,7 +42,7 @@ if __name__ == "__main__":
noiseless_data = R(sh) noiseless_data = R(sh)
signal_to_noise = 1. signal_to_noise = 1.
noise_amplitude = noiseless_data.std()/signal_to_noise noise_amplitude = noiseless_data.val.std()/signal_to_noise
N = ift.DiagonalOperator(ift.Field.full(s_space, noise_amplitude**2)) N = ift.DiagonalOperator(ift.Field.full(s_space, noise_amplitude**2))
n = ift.Field.from_random(domain=s_space, n = ift.Field.from_random(domain=s_space,
random_type='normal', random_type='normal',
......
...@@ -36,7 +36,7 @@ if __name__ == "__main__": ...@@ -36,7 +36,7 @@ if __name__ == "__main__":
R = Instrument*ht R = Instrument*ht
noiseless_data = R(sh) noiseless_data = R(sh)
signal_to_noise = 1. signal_to_noise = 1.
noise_amplitude = noiseless_data.std()/signal_to_noise noise_amplitude = noiseless_data.val.std()/signal_to_noise
N = ift.DiagonalOperator(ift.Field.full(s_space, noise_amplitude**2)) N = ift.DiagonalOperator(ift.Field.full(s_space, noise_amplitude**2))
n = ift.Field.from_random(domain=s_space, n = ift.Field.from_random(domain=s_space,
random_type='normal', random_type='normal',
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment