def linear_interp_basis_dm(toas, freqs, dt=30 * const.day):
# get linear interpolation basis in time
U, avetoas = utils.linear_interp_basis(toas, dt=dt)
# scale with radio frequency
Dm = (1400 / freqs)**2
return U * Dm[:, None], avetoas
def linear_interp_basis_chromatic(toas, freqs, dt=30 * 86400, idx=4):
"""Linear interpolation basis in time with nu^-4 scaling"""
# get linear interpolation basis in time
U, avetoas = utils.linear_interp_basis(toas, dt=dt)
# scale with radio frequency
Dm = (1400 / freqs)**idx
return U * Dm[:, None], avetoas
def linear_interp_basis_sw_dm(toas, freqs, planetssb, pos_t, dt=7*86400):
# get linear interpolation basis in time
U, avetoas = utils.linear_interp_basis(toas, dt=dt)
# scale with radio frequency
theta, R_earth = theta_impact(planetssb, pos_t)
dm_sol_wind = dm_solar(1.0, theta, R_earth)
dt_DM = dm_sol_wind * 4.148808e3 / (freqs**2)
return U * dt_DM[:, None], avetoas
def linear_interp_basis_freq(freqs, df=64):
return utils.linear_interp_basis(freqs, dt=df)
def linear_interp_basis_freq(freqs, df=64):
"""Linear interpolation in radio frequency"""
return utils.linear_interp_basis(freqs, dt=df)