Science/psdlag-agn
1
0
Fork 0
mirror of https://asciireactor.com/otho/psdlag-agn.git synced 2024-11-21 21:45:06 +00:00
Code Issues Packages Projects Releases Wiki Activity

changing argument handling

Browse Source
This commit is contained in:
caes 2017-01-31 21:38:20 -05:00
parent abd1587aec
commit 3075130e20
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
scripts/analyze_lightcurve.py
View File

@ -74,14 +74,14 @@ lc1_time, lc1_strength, lc1_strength_err = np.loadtxt(args[0],
# for pylab: errorbar(t1,l1,yerr=l1e,fmt='o') # for pylab: errorbar(t1,l1,yerr=l1e,fmt='o')
# Used throughout # Used throughout
initial_args = np.ones(nfq)
## initialize the psd class for multiple light curves ## ## initialize the psd class for multiple light curves ##
P1 = clag.clag('psd10r', P1 = clag.clag('psd10r',
[lc1_time], [lc1_strength], [lc1_strength_err], [lc1_time], [lc1_strength], [lc1_strength_err],
dt, fqL) dt, fqL)
ref_psd, ref_psd_err = clag.optimize(P1, initial_args) ref_psd = np.ones(nfq)
ref_psd, ref_psd_err = clag.optimize(P1, ref_psd)
ref_psd, ref_psd_err = clag.errors(P1, ref_psd, ref_psd_err) ref_psd, ref_psd_err = clag.errors(P1, ref_psd, ref_psd_err)
## plot ## ## plot ##
@ -91,7 +91,8 @@ ref_psd, ref_psd_err = clag.errors(P1, ref_psd, ref_psd_err)
# Load second light curve # Load second light curve
lc2_time, lc2_strength, lc2_strength_err = np.loadtxt(args[1],skiprows=1).T lc2_time, lc2_strength, lc2_strength_err = np.loadtxt(args[1],skiprows=1).T
P2 = clag.clag('psd10r', [lc2_time], [lc2_strength], [lc2_strength_err], dt, fqL) P2 = clag.clag('psd10r', [lc2_time], [lc2_strength], [lc2_strength_err], dt, fqL)
echo_psd, echo_psd_err = clag.optimize(P2, initial_args) echo_psd = np.ones(nfq)
echo_psd, echo_psd_err = clag.optimize(P2, echo_psd)
echo_psd, echo_psd_err = clag.errors(P2, echo_psd, echo_psd_err) echo_psd, echo_psd_err = clag.errors(P2, echo_psd, echo_psd_err)