Skip to content

Commit

Permalink
remove debug statements in clark clean
Browse files Browse the repository at this point in the history
  • Loading branch information
landmanbester committed Nov 30, 2024
1 parent 90f853c commit 10d1190
Showing 1 changed file with 2 additions and 5 deletions.
7 changes: 2 additions & 5 deletions pfb/deconv/clark.py
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@
import pyscilog
log = pyscilog.get_logger('CLARK')

# @njit(nogil=True, cache=True) # parallel=True,
@njit(nogil=True, cache=True) # parallel=True,
def subtract(A, psf, Ip, Iq, xhat, nxo2, nyo2):
'''
Subtract psf centered at location of xhat
Expand All @@ -25,7 +25,7 @@ def subtract(A, psf, Ip, Iq, xhat, nxo2, nyo2):
return A


# @njit(nogil=True, cache=True) # parallel=True,
@njit(nogil=True, cache=True) # parallel=True,
def subminor(A, psf, Ip, Iq, model, wsums, gamma=0.05, th=0.0, maxit=10000):
"""
Run subminor loop in active set
Expand Down Expand Up @@ -83,9 +83,6 @@ def subminor(A, psf, Ip, Iq, model, wsums, gamma=0.05, th=0.0, maxit=10000):
# mask = (left_x & right_x) & (left_y & right_y)

mask = (np.abs(Idelp) <= nxo2) & (np.abs(Idelq) <= nyo2)
assert mask.size == Ip.size
assert mask.all()
# import ipdb; ipdb.set_trace()
# Ipp, Iqq = psf[:, nxo2 - Ip[mask], nyo2 - Iq[mask]]
A = subtract(A[:, mask], psf,
Idelp[mask], Idelq[mask],
Expand Down

0 comments on commit 10d1190

Please sign in to comment.