format=logging_format,
handlers=[RichHandler(rich_tracebacks=True)],
)
else:
logging.basicConfig(
level=logging.INFO,
format=logging_format,
handlers=[RichHandler(rich_tracebacks=True)],
)
logging.info("Input Arguments:-")
for arg, value in sorted(vars(values).items()):
logging.info("%s: %r", arg, value)
your_object = Your(file=values.files)
max_delay = your_object.dispersion_delay(dms=np.max(values.dm))
dispersion_delay_samples = np.ceil(max_delay /
your_object.your_header.tsamp)
logging.info(f"Max Dispersion delay = {max_delay} s")
logging.info(f"Max Dispersion delay = {dispersion_delay_samples} samples")
if your_object.your_header.nspectra < 2**18:
nsamps_gulp = your_object.your_header.nspectra
else:
nsamps_gulp = int(
np.max([(2**np.ceil(np.log2(dispersion_delay_samples))), 2**18]))
if values.channel_end > 0:
c_max = values.channel_end
else:
c_max = None
if args.verbose:
logging.basicConfig(filename=log_filename,
level=logging.DEBUG,
format=logging_format)
logging.debug("Using debug mode")
else:
logging.basicConfig(filename=log_filename,
level=logging.INFO,
format=logging_format)
logging.info("Input Arguments:-")
for arg, value in sorted(vars(args).items()):
logging.info("%s: %r", arg, value)
your_object = Your(file=args.files)
max_delay = your_object.dispersion_delay(dms=np.max(args.dm))
dispersion_delay_samples = np.ceil(max_delay /
your_object.your_header.tsamp)
logging.info(f"Max Dispersion delay = {max_delay} s")
logging.info(f"Max Dispersion delay = {dispersion_delay_samples} samples")
nsamps_gulp = int(
np.max([(2**np.ceil(np.log2(dispersion_delay_samples))), 2**18]))
your_dada = dada.YourDada(your_object)
your_dada.setup()
if args.apply_savgol:
bandpass = your_object.bandpass(nspectra=8192)
mask = savgol_filter(bandpass,
your_object.your_header.foff,
