Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
ift
NIFTy
Commits
2cc60994
Commit
2cc60994
authored
Jul 12, 2017
by
Theo Steininger
Browse files
Merge branch 'fix_demos' into 'master'
fix demos See merge request
!162
parents
92422402
a0388ae2
Pipeline
#14806
passed with stages
in 14 minutes and 36 seconds
Changes
3
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
demos/__init__.py
deleted
100644 → 0
View file @
92422402
demos/critical_filtering.py
View file @
2cc60994
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
=
1
0
e-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
=
1
0
e-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
)
...
...
demos/wiener_filter_advanced.py
View file @
2cc60994
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
=
1
0
e-5
,
convergence_tolerance
=
1e-5
,
preconditioner
=
None
)
# Setting starting position
m0
=
Field
(
h_space
,
val
=
.
0
)
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment