Skip to content

Commit

Permalink
Removing unneeded comments.
Browse files Browse the repository at this point in the history
  • Loading branch information
kmacdonald-stsci committed Oct 17, 2023
1 parent f644b9c commit 080393c
Showing 1 changed file with 7 additions and 37 deletions.
44 changes: 7 additions & 37 deletions src/stcal/ramp_fitting/ols_fit.py
Original file line number Diff line number Diff line change
Expand Up @@ -677,10 +677,10 @@ def ols_ramp_fit_single(

if ramp_data.drop_frames1 is None:
ramp_data.drop_frames1 = 0
dbg_print(" ---------------- Entering C Code ----------------")
print(" ---------------- Entering C Code ----------------")
image_info, integ_info, opt_info = ols_slope_fitter(
ramp_data, gain_2d, readnoise_2d, weighting, save_opt)
dbg_print(" ---------------- Return C Code ----------------")
print(" ---------------- Return C Code ----------------")

c_end = time.time()
print("=" * 80)
Expand All @@ -691,10 +691,10 @@ def ols_ramp_fit_single(
p_start = time.time()

print("\n");
dbg_print(" ---------------- Entering Python Code ----------------")
print(" ---------------- Entering Python Code ----------------")
image_info, integ_info, opt_info = ols_ramp_fit_single_python(
ramp_data, buffsize, save_opt, readnoise_2d, gain_2d, weighting)
dbg_print(" ---------------- Return Python Code ----------------")
print(" ---------------- Return Python Code ----------------")

# XXX end python time
p_end = time.time()
Expand Down Expand Up @@ -810,9 +810,9 @@ def c_python_time_comparison(c_start, c_end, p_start, p_end):
c_diff = c_end - c_start
p_diff = p_end - p_start
c_div_p = c_diff / p_diff * 100.
dbg_print(f"{c_diff = }")
dbg_print(f"{p_diff = }")
dbg_print(f"{c_div_p = :.4f}%")
print(f"{c_diff = }")
print(f"{p_diff = }")
print(f"{c_div_p = :.4f}%")


def discard_miri_groups(ramp_data):
Expand Down Expand Up @@ -1250,13 +1250,6 @@ def ramp_fit_compute_variances(ramp_data, gain_2d, readnoise_2d, fit_slopes_ans)
var_p4[num_int, :, med_rates <= 0.] = 0.
var_both4[num_int, :, :, :] = var_r4[num_int, :, :, :] + var_p4[num_int, :, :, :]

'''
if num_int == 0:
print("=" * 80)
dbg_print(f"var_both4 = integration 0\n{var_both4[0, :, 0, 0]}")
print("=" * 80)
'''

inv_var_both4[num_int, :, :, :] = 1. / var_both4[num_int, :, :, :]

# Want to retain values in the 4D arrays only for the segments that each
Expand Down Expand Up @@ -1315,12 +1308,6 @@ def ramp_fit_compute_variances(ramp_data, gain_2d, readnoise_2d, fit_slopes_ans)
ramp_data.groupdq = groupdq
ramp_data.pixeldq = inpixeldq

'''
print("=" * 80)
dbg_print(f"var_both4 = integration 0\n{var_both4[0, :, 0, 0]}")
print("=" * 80)
'''

return var_p3, var_r3, var_p4, var_r4, var_both4, var_both3, inv_var_both4, \
s_inv_var_p3, s_inv_var_r3, s_inv_var_both3

Expand Down Expand Up @@ -2695,14 +2682,6 @@ def fit_next_segment_long_end_of_ramp(
if len(g_pix) > 0:
inv_var[g_pix] += 1.0 / variance[g_pix]

'''
print("=" * 80)
dbg_print(f"{ramp_data.current_integ = }")
dbg_print(f"{variance = }")
dbg_print(f"{inv_var = }")
print("=" * 80)
'''

# Append results to arrays
opt_res.append_arr(num_seg, g_pix, intercept, slope, sig_intercept,
sig_slope, inv_var, save_opt)
Expand Down Expand Up @@ -3166,15 +3145,6 @@ def fit_double_read(mask_2d, wh_pix_2r, data_masked, slope_s, intercept_s,
sig_slope_s[pixel_ff] = np.sqrt(2) * rn
sig_intercept_s[pixel_ff] = np.sqrt(2) * rn

'''
if len(wh_pix_2r[0]) > 0:
print("=" * 80)
dbg_print(f"{mask_2d[:, 0] = }")
dbg_print(f"{rn = }")
dbg_print(f"{variance_s = }")
print("=" * 80)
'''

return slope_s, intercept_s, variance_s, sig_slope_s, sig_intercept_s


Expand Down

0 comments on commit 080393c

Please sign in to comment.