Commit e8079ae3 authored by Martin Reinecke's avatar Martin Reinecke
Browse files

adjust notebook

parent 2fedb50b
......@@ -260,11 +260,11 @@
mask = np.full(s_space.shape, 1.)
mask[l:h] = 0
mask = ift.Field.from_global_data(s_space, mask)
R = ift.DiagonalOperator(mask).chain(HT)
n = n.to_global_data().copy()
n = n.to_global_data_rw()
n[l:h] = 0
n = ift.Field.from_global_data(s_space, n)
d = R(sh) + n
```
......@@ -299,11 +299,11 @@
# Get signal data and reconstruction data
s_data = s.to_global_data()
m_data = HT(m).to_global_data()
m_var_data = m_var.to_global_data()
uncertainty = np.sqrt(m_var_data)
d_data = d.to_global_data().copy()
d_data = d.to_global_data_rw()
# Set lost data to NaN for proper plotting
d_data[d_data == 0] = np.nan
```
......@@ -360,11 +360,11 @@
mask = np.full(s_space.shape, 1.)
mask[l:h,l:h] = 0.
mask = ift.Field.from_global_data(s_space, mask)
R = ift.DiagonalOperator(mask).chain(HT)
n = n.to_global_data().copy()
n = n.to_global_data_rw()
n[l:h, l:h] = 0
n = ift.Field.from_global_data(s_space, n)
curv = Curvature(R=R, N=N, Sh=Sh)
D = curv.inverse
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment