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

fix demos

parent 92422402
No related branches found
No related tags found
1 merge request!162fix demos
Pipeline #
from nifty import *
from nifty.library.wiener_filter import WienerFilterEnergy
from nifty.library.critical_filter import CriticalPowerEnergy
......@@ -113,7 +112,7 @@ if __name__ == "__main__":
print (x, iteration)
minimizer1 = RelaxedNewton(convergence_tolerance=10e-2,
minimizer1 = RelaxedNewton(convergence_tolerance=1e-2,
convergence_level=2,
iteration_limit=3,
callback=convergence_measure)
......@@ -124,7 +123,7 @@ if __name__ == "__main__":
max_history_length=3)
# Setting starting position
flat_power = Field(p_space,val=10e-8)
flat_power = Field(p_space,val=1e-8)
m0 = flat_power.power_synthesize(real_signal=True)
t0 = Field(p_space, val=log(1./(1+p_space.kindex)**2))
......@@ -136,7 +135,7 @@ if __name__ == "__main__":
distribution_strategy=distribution_strategy)
# Initializing the nonlinear Wiener Filter energy
map_energy = WienerFilterEnergy(position=m0, d=d, R=R, N=N, S=S0, inverter=inverter)
map_energy = WienerFilterEnergy(position=m0, d=d, R=R, N=N, S=S0)
# Solving the Wiener Filter analytically
D0 = map_energy.curvature
m0 = D0.inverse_times(j)
......
from nifty import *
import plotly.offline as pl
import plotly.graph_objs as go
from nifty.library.wiener_filter import *
from mpi4py import MPI
comm = MPI.COMM_WORLD
......@@ -103,7 +103,7 @@ if __name__ == "__main__":
#
inverter = ConjugateGradient(convergence_level=3,
convergence_tolerance=10e-5,
convergence_tolerance=1e-5,
preconditioner=None)
# Setting starting position
m0 = Field(h_space, val=.0)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment