phoebe2
phoebe2 copied to clipboard
Error when using ELLC in backend with GPs
PHOEBE Can't find gp_exclude_phases_enabled qualifier when using ELLC in backend with GPs
ValueError: 0 results found for twig: 'None', {'qualifier': 'gp_exclude_phases_enabled', 'dataset': 'lc01', 'check_visible': False, 'check_default': False, 'check_advanced': False, 'check_single': False}
""" import phoebe import numpy as np
def test_gps(verbose=False, plot=False): # Make fake data b = phoebe.default_binary()
b.add_dataset('lc', compute_times=phoebe.linspace(0,5.,501))
#Make a set of fast compute options
b.add_compute('ellc', compute='fast_compute')
b.set_value_all('ld_mode', value='manual')
b.set_value('irrad_method', compute='fast_compute', value='none')
b.set_value_all('distortion_method', compute='fast_compute', value='sphere')
b.run_compute(compute='fast_compute')
times = b.get_value(qualifier='times', context='model')
fluxes = b.get_value(qualifier='fluxes', context='model') + np.random.normal(size=times.shape) * 0.07 + 0.2*np.sin(times)
sigmas = np.ones_like(fluxes) * 0.05
#Upload the fake data to PHOEBE
b.add_dataset('lc', dataset='lc01', times=times, fluxes=fluxes, sigmas=sigmas, overwrite=True)
b.set_value_all('ld_mode', value='manual')
b.run_compute(model='withoutGPs', compute='fast_compute')
#Make a model with GPs
b.add_gaussian_process('celerite2', dataset='lc01', kernel='sho')
print (b['gp_exclude_phases_enabled'])
b['gp_exclude_phases_enabled'] = False
# Compute model in phase space
b.flip_constraint('compute_phases', solve_for='compute_times')
b.set_value('compute_phases', phoebe.linspace(0,1,101))
b.run_compute(model='withGPs', compute='fast_compute')
print ("with GPs")
if name == 'main': logger = phoebe.logger(clevel='INFO') test_gps(verbose=True, plot=True) """