Skip to content

Commit ebad100

Browse files
committed
fix syntax
1 parent e61e7da commit ebad100

3 files changed

Lines changed: 30 additions & 11 deletions

File tree

litebird_sim/mapmaking/pair_differencing.py

Lines changed: 11 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -141,10 +141,14 @@ def _accumulate_pair_differenced_samples_and_build_nobs_matrix(
141141
info_pix[1, 0] += pair_cos * pair_sin * inv_sigma2
142142
info_pix[1, 1] += pair_sin * pair_sin * inv_sigma2
143143

144-
145-
for cur_sample_t, cur_sample_b, cur_pix_idx, cur_psi_t, cur_psi_b, cur_t_mask in zip(
146-
tod_t, tod_b, pix_t, psi_t, psi_b, t_mask
147-
):
144+
for (
145+
cur_sample_t,
146+
cur_sample_b,
147+
cur_pix_idx,
148+
cur_psi_t,
149+
cur_psi_b,
150+
cur_t_mask,
151+
) in zip(tod_t, tod_b, pix_t, psi_t, psi_b, t_mask):
148152
if cur_t_mask:
149153
pair_sample = cur_sample_t - cur_sample_b
150154
pair_cos = np.cos(2 * cur_psi_t) - np.cos(2 * cur_psi_b)
@@ -198,7 +202,9 @@ def _build_nobs_matrix(
198202
)
199203

200204
for t_idx, b_idx in detector_pairs:
201-
if not np.array_equal(pixidx_all[t_idx, cur_t_mask], pixidx_all[b_idx, cur_t_mask]):
205+
if not np.array_equal(
206+
pixidx_all[t_idx, cur_t_mask], pixidx_all[b_idx, cur_t_mask]
207+
):
202208
raise ValueError(
203209
f"Observation {obs_idx}: detectors {t_idx} and {b_idx} do not "
204210
"observe the same map pixels for the selected time samples."

litebird_sim/simulations.py

Lines changed: 10 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -2216,7 +2216,9 @@ def make_pair_differenced_map(
22162216
self.check_valid_splits(detector_split, time_split)
22172217

22182218
if append_to_report and MPI_COMM_WORLD.rank == 0:
2219-
template_file_path = get_template_file_path("report_pair_differenced_map.md")
2219+
template_file_path = get_template_file_path(
2220+
"report_pair_differenced_map.md"
2221+
)
22202222
with template_file_path.open("rt") as inpf:
22212223
markdown_template = "".join(inpf.readlines())
22222224
self.append_to_report(
@@ -2264,7 +2266,9 @@ def make_pair_differenced_map_splits(
22642266
self.check_valid_splits(detector_splits, time_splits)
22652267

22662268
if append_to_report and MPI_COMM_WORLD.rank == 0:
2267-
template_file_path = get_template_file_path("report_pair_differenced_map_splits.md")
2269+
template_file_path = get_template_file_path(
2270+
"report_pair_differenced_map_splits.md"
2271+
)
22682272
with template_file_path.open("rt") as inpf:
22692273
markdown_template = "".join(inpf.readlines())
22702274
self.append_to_report(
@@ -2300,7 +2304,10 @@ def make_pair_differenced_map_splits(
23002304

23012305
filenames.append(
23022306
self.write_healpix_map(
2303-
Path(file), mapp, column_names=names, coord=result.coordinate_system.name
2307+
Path(file),
2308+
mapp,
2309+
column_names=names,
2310+
coord=result.coordinate_system.name,
23042311
)
23052312
)
23062313

test/test_pair_differencing.py

Lines changed: 9 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -7,7 +7,9 @@
77
from litebird_sim.mapmaking.pair_differencing import make_pair_differenced_map
88

99

10-
def _build_observation(detectors: list[dict]) -> tuple[lbs.Observation, np.ndarray, int]:
10+
def _build_observation(
11+
detectors: list[dict],
12+
) -> tuple[lbs.Observation, np.ndarray, int]:
1113
n_samples = 4
1214
theta, phi = hp.pix2ang(1, 0)
1315
psi = np.array([0.0, np.pi / 8.0, np.pi / 4.0, 3.0 * np.pi / 8.0])
@@ -71,7 +73,9 @@ def test_make_pair_differenced_map_recovers_qu():
7173
unseen_mask[pix_idx] = False
7274
assert np.all(result.binned_map[:, unseen_mask] == hp.UNSEEN)
7375

74-
pair_weight = 0.5 * (detectors[0]["net_ukrts"] ** 2 + detectors[1]["net_ukrts"] ** 2)
76+
pair_weight = 0.5 * (
77+
detectors[0]["net_ukrts"] ** 2 + detectors[1]["net_ukrts"] ** 2
78+
)
7579
design_matrix = np.column_stack(
7680
(
7781
np.cos(2 * psi_t) - np.cos(2 * psi_b),
@@ -219,4 +223,6 @@ def test_simulation_make_pair_differenced_map_splits(tmp_path):
219223
)
220224

221225
assert "full_full" in results
222-
assert np.allclose(results["full_full"].binned_map[:, pix_idx], np.array([q_true, u_true]))
226+
assert np.allclose(
227+
results["full_full"].binned_map[:, pix_idx], np.array([q_true, u_true])
228+
)

0 commit comments

Comments
 (0)