Skip to content

Commit

Permalink
test refined gains without Winv in symbolic expression
Browse files Browse the repository at this point in the history
  • Loading branch information
landmanbester committed Dec 17, 2023
1 parent 6150144 commit 5ade14d
Show file tree
Hide file tree
Showing 2 changed files with 6 additions and 6 deletions.
4 changes: 2 additions & 2 deletions pfb/operators/gridder.py
Original file line number Diff line number Diff line change
Expand Up @@ -649,8 +649,8 @@ def image_data_products(uvw,
uvw=uvw,
freq=freq,
vis=vis,
wgt=wgt,
# wgt=imwgt, # data already naturally weighted
# wgt=wgt,
wgt=imwgt, # data already naturally weighted
mask=mask,
npix_x=nx, npix_y=ny,
pixsize_x=cellx, pixsize_y=celly,
Expand Down
8 changes: 4 additions & 4 deletions pfb/utils/stokes.py
Original file line number Diff line number Diff line change
Expand Up @@ -344,15 +344,15 @@ def stokes_funcs(data, jones, product, pol):
[0, 1.0, 1.0j, 0],
[0, 1.0, -1.0j, 0],
[1.0, 0, 0, -1.0]])
Tinv = T.inv()
# Tinv = T.inv()

# Full Stokes weights
W = T.H * Mpq.H * Sinv * Mpq * T
Winv = Tinv * Mpqinv * S * Mpqinv.H * Tinv.H
# Winv = Tinv * Mpqinv * S * Mpqinv.H * Tinv.H

# Full Stokes coherencies
C = Winv * (T.H * (Mpq.H * (Sinv * Vpq)))
# C = T.H * (Mpq.H * (Sinv * Vpq))
# C = Winv * (T.H * (Mpq.H * (Sinv * Vpq)))
C = T.H * (Mpq.H * (Sinv * Vpq))

if product == literal('I'):
i = 0
Expand Down

0 comments on commit 5ade14d

Please sign in to comment.