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

[WIP] Full VCS modelling #160

Open
wants to merge 16 commits into
base: master
Choose a base branch
from
+608 −0
Open

[WIP] Full VCS modelling #160

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
16 commits
Select commit Hold shift + click to select a range
2159f42
Add initial code for full VCS modelling
e-koch Aug 17, 2017
74133f0
Add integral of F_eps needed for normalization
e-koch Aug 31, 2017
68b9e93
Actually already wrote this function in F_eps_norm
e-koch Sep 1, 2017
fd52391
Fixes for C_eps
e-koch Sep 2, 2017
16836f5
Several critical VCS fixes
e-koch Sep 8, 2017
7d0ce9e
Various fixes for Dz, Ceps, and window functions
e-koch Sep 14, 2017
33b8431
Use analytic solution for integral in C_eps for shallow spectra
e-koch Sep 14, 2017
df36e2f
Add analytic solution for steep emissivity spectrum
e-koch Sep 14, 2017
dbb5ed5
Add analytical version of Int1 and Int2; use flint for a fast general…
e-koch Sep 15, 2017
b58d023
Various updates and changes to the VCS main model functions
e-koch Sep 15, 2017
e43955a
Simplify Dz a bit
e-koch Sep 15, 2017
98d46ca
Use hyp1f1 from flint; another Dz function to test
e-koch Sep 16, 2017
afa3239
Missing a factor of pi in gaussian autocorrelation function
e-koch Sep 16, 2017
94dffc5
Use Dz_simp for now (not completely sure it's right, though); normali…
e-koch Sep 16, 2017
c45060f
Add in a mu parameter for the thermal line width
e-koch Sep 18, 2017
ef4f363
Try imposing bounds on theta in the Gaussian fit
e-koch Sep 20, 2017
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
3 changes: 3 additions & 0 deletions turbustat/statistics/pca/width_estimate.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -400,6 +400,9 @@ def fit_2D_gaussian(xmat, ymat, z):
'y_mean': True}) + \
astropy_models.Const2D(amplitude=[np.percentile(z, 10)])

# Impose restrictions on theta
g.theta_0.bounds = (-np.pi, np.pi)

fit_g = fitting.LevMarLSQFitter()
output = fit_g(g, xmat, ymat, z)
cov = fit_g.fit_info['param_cov']
Expand Down
Loading