Skip to content

Commit c76684f

Browse files
committed
doc: fix blowing up in deblending tutorial
1 parent e297055 commit c76684f

File tree

1 file changed

+2
-3
lines changed

1 file changed

+2
-3
lines changed

tutorials/deblending.py

Lines changed: 2 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -128,9 +128,8 @@
128128
Op = Bop * Sop
129129

130130
# Compute max eigenvalue (we do this explicitly to be able to run this fast)
131-
Op1 = pylops.LinearOperator(Op.H * Op, explicit=False)
132-
X = np.random.rand(Op1.shape[0], 1).astype(Op1.dtype)
133-
maxeig = sp_lobpcg(Op1, X=X, maxiter=5, tol=1e-10)[0][0]
131+
Op1 = Op.H * Op
132+
maxeig = np.abs(Op1.eigs(1, niter=5, ncv=5, tol=5e-2))[0]
134133
alpha = 1.0 / maxeig
135134

136135
# Deblend

0 commit comments

Comments
 (0)