Skip to content

Commit 0000000

Browse files
committed
stokes I fix
1 parent dfcc22b commit 0000000

File tree

2 files changed

+5
-6
lines changed

2 files changed

+5
-6
lines changed

src/polarizedtypes.jl

Lines changed: 3 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -19,11 +19,10 @@ function uvtable_values_to(::Type{CoherencyMatrix}, uvtbl)
1919
end
2020

2121
function uvtable_values_to(::Type{IPol}, uvtbl)
22-
grs = @p uvtbl group_vg((;_.datetime, _.freq_spec, _.spec))
22+
grs = @p uvtbl filter(is_parallel_hands(_.stokes)) group_vg((;_.datetime, _.freq_spec, _.spec))
2323
return map(grs) do gr
24-
par_hands = @p gr filter(x -> is_parallel_hands(x.stokes))
25-
length(par_hands) (1, 2) || error("expected 1 or 2 parallel-hand Stokes per group, got $(length(par_hands))")
26-
val = mean(x -> x.value, par_hands)
24+
length(gr) (1, 2) || error("expected 1 or 2 parallel-hand Stokes per group, got $(length(gr))")
25+
val = mean(x -> x.value, gr)
2726
@p let
2827
first(gr)
2928
@set __.value = val

test/runtests.jl

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -297,12 +297,12 @@ end
297297
@test length(ipol_single_ll) == 1
298298
@test ipol_single_ll[1].value == 7.0
299299

300-
# IPol: error when no parallel hands
300+
# IPol: when no parallel hands
301301
uvtbl_no_par = [
302302
(datetime=1, freq_spec=100.0, spec=1, stokes=:LR, value=12.0),
303303
(datetime=1, freq_spec=100.0, spec=1, stokes=:RL, value=13.0),
304304
]
305-
@test_throws "expected 1 or 2 parallel-hand" VLBI.uvtable_values_to(IPol, uvtbl_no_par)
305+
@test isempty(VLBI.uvtable_values_to(IPol, uvtbl_no_par))
306306
end
307307

308308
@testitem "_" begin

0 commit comments

Comments
 (0)