Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Added dm nondiag test #114

Open
wants to merge 11 commits into
base: master
Choose a base branch
from
Open

Added dm nondiag test #114

Changes from 8 commits
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
478eaf6
Added dm nondiag test
Jyvol Jun 15, 2021
53547f9
Update test_models.py
Jyvol Jun 21, 2021
73fb4f1
Update test_models.py
Jyvol Jun 22, 2021
6508dbd
Update test_models.py
Jyvol Jun 22, 2021
46e90e8
Update test_models.py
Jyvol Jun 22, 2021
874ced1
Nondiag fix
Jyvol Jun 23, 2021
e160403
Update test_models.py
Jyvol Jul 14, 2021
a6b24f3
Update test_models.py
Jyvol Jul 14, 2021
8c37444
Update test_models.py
Jyvol Jul 15, 2021
0b93865
Update test_models.py
Jyvol Jul 20, 2021
734329a
Update test_models.py
Jyvol Jul 22, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
83 changes: 83 additions & 0 deletions tests/test_models.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -44,6 +44,12 @@ def nodmx_psrs(caplog):
psrs.append(pickle.load(fin))

return psrs
def psr_name(nodmx_psrs):
psn = []
for p in psr_names:
with open(datadir+'/{0}_ng9yr_nodmx_DE436_epsr.pkl'.format(p),'rb'):
psn.append(p)
return psn
paulthebaker marked this conversation as resolved.
Show resolved Hide resolved

def test_model_singlepsr_noise(nodmx_psrs,caplog):
# caplog.set_level(logging.CRITICAL)
Expand Down Expand Up @@ -83,6 +89,83 @@ def test_model_singlepsr_noise_dip_cusp(nodmx_psrs,caplog):
x0 = {pname:p.sample() for pname,p in zip(m.param_names, m.params)}
m.get_lnlikelihood(x0)

def test_model_singlepsr_noise_dm_nondiag(nodmx_psrs,caplog):
paulthebaker marked this conversation as resolved.
Show resolved Hide resolved
# caplog.set_level(logging.CRITICAL)
ii = 1
mn=models.model_singlepsr_noise(nodmx_psrs[ii], dm_var=True,
dm_type='gp', dmgp_kernel ='nondiag',
dm_nondiag_kernel ='dmx_like')
assert hasattr(mn,'get_lnlikelihood')
x0 = {pname:p.sample() for pname,p in zip(mn.param_names, mn.params)}
assert psr_name[ii]+'_dm_gp_log10_sigma' in mn.param_names
paulthebaker marked this conversation as resolved.
Show resolved Hide resolved
mn.get_lnlikelihood(x0)

mn=models.model_singlepsr_noise(nodmx_psrs[1], dm_var=True,
dm_type='gp', dmgp_kernel ='nondiag',
dm_nondiag_kernel ='sq_exp')
assert hasattr(mn,'get_lnlikelihood')
x0 = {pname:p.sample() for pname,p in zip(mn.param_names, mn.params)}
mn.get_lnlikelihood(x0)

mn=models.model_singlepsr_noise(nodmx_psrs[1], dm_var=True,
dm_type='gp', dmgp_kernel ='nondiag',
dm_nondiag_kernel ='sq_exp_rfband')
assert hasattr(mn,'get_lnlikelihood')
x0 = {pname:p.sample() for pname,p in zip(mn.param_names, mn.params)}
mn.get_lnlikelihood(x0)

mn=models.model_singlepsr_noise(nodmx_psrs[1], dm_var=True,
dm_type='gp', dmgp_kernel ='nondiag',
dm_nondiag_kernel ='periodic')
assert hasattr(mn,'get_lnlikelihood')
x0 = {pname:p.sample() for pname,p in zip(mn.param_names, mn.params)}
mn.get_lnlikelihood(x0)

mn=models.model_singlepsr_noise(nodmx_psrs[1], dm_var=True,
dm_type='gp', dmgp_kernel ='nondiag',
dm_nondiag_kernel ='periodic_rfband')
assert hasattr(mn,'get_lnlikelihood')
x0 = {pname:p.sample() for pname,p in zip(mn.param_names, mn.params)}
mn.get_lnlikelihood(x0)

mn=models.model_singlepsr_noise(nodmx_psrs[1], dm_var=True,
dm_type='gp', dmgp_kernel ='nondiag',
dm_nondiag_kernel ='None')
paulthebaker marked this conversation as resolved.
Show resolved Hide resolved
assert hasattr(mn,'get_lnlikelihood')
x0 = {pname:p.sample() for pname,p in zip(mn.param_names, mn.params)}
mn.get_lnlikelihood(x0)

def test_model_singlepsr_noise_dm_diag(nodmx_psrs,caplog):
# caplog.set_level(logging.CRITICAL)
mn=models.model_singlepsr_noise(nodmx_psrs[1], dm_var=True,
dm_type='gp', dmgp_kernel ='diag',
dm_psd ='powerlaw')
assert hasattr(mn,'get_lnlikelihood')
x0 = {pname:p.sample() for pname,p in zip(mn.param_names, mn.params)}
mn.get_lnlikelihood(x0)

mn=models.model_singlepsr_noise(nodmx_psrs[1], dm_var=True,
dm_type='gp', dmgp_kernel ='diag',
dm_psd ='turnover')
assert hasattr(mn,'get_lnlikelihood')
x0 = {pname:p.sample() for pname,p in zip(mn.param_names, mn.params)}
mn.get_lnlikelihood(x0)

mn=models.model_singlepsr_noise(nodmx_psrs[1], dm_var=True,
dm_type='gp', dmgp_kernel ='diag',
dm_psd ='tprocess')
assert hasattr(mn,'get_lnlikelihood')
x0 = {pname:p.sample() for pname,p in zip(mn.param_names, mn.params)}
mn.get_lnlikelihood(x0)

mn=models.model_singlepsr_noise(nodmx_psrs[1], dm_var=True,
dm_type='gp', dmgp_kernel ='diag',
dm_psd ='tprocess_adapt')
assert hasattr(mn,'get_lnlikelihood')
x0 = {pname:p.sample() for pname,p in zip(mn.param_names, mn.params)}
mn.get_lnlikelihood(x0)


def test_model_singlepsr_noise_chrom_nondiag(nodmx_psrs,caplog):
# caplog.set_level(logging.CRITICAL)
m=models.model_singlepsr_noise(nodmx_psrs[1], dm_var=True,
Expand Down