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.

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

Update example online and make quickstart.py downloadable #112

Open
mtazzari opened this issue Dec 13, 2017 · 1 comment
Open

Update example online and make quickstart.py downloadable #112

mtazzari opened this issue Dec 13, 2017 · 1 comment

Comments

@mtazzari
Copy link
Owner

mtazzari commented Dec 13, 2017

Remove from GaussianProfile() the following lines altogether:

# convert to radians
sigma *= arcsec     # [rad]	
Rmin *= arcsec      # [rad]
dR *= arcsec        # [rad]

and move them to lnpostfn():

 # convert to rad
inc *= deg
PA *= deg
dRA *= arcsec
dDec *= arcsec
sigma *= arcsec     # [rad]	
Rmin *= arcsec      # [rad]
dR *= arcsec        # [rad]

 # compute the model brightness profile
f = GaussianProfile(f0, sigma, Rmin, dR, nR)
chi2 = chi2Profile(f, Rmin, dR, nxy, dxy, u, v, Re, Im, w, inc=inc, PA=PA, dRA=dRA, dDec=dDec)

and re-define:

def GaussianProfile(f0, sigma, Rmin, dR, nR):
	""" Gaussian brightness profile. """
	  
	# radial grid
	R = np.linspace(Rmin, Rmin + dR*nR, nR, endpoint=False)

	return f0 * np.exp(-0.5*(R/sigma)**2)

Check the quickstart.py exactly reflects the online example.

@mtazzari mtazzari added this to the 1.0.3 milestone Dec 13, 2017
@mtazzari mtazzari self-assigned this Dec 13, 2017
@mtazzari mtazzari reopened this Apr 25, 2018
@mtazzari
Copy link
Owner Author

Still to do

@mtazzari mtazzari modified the milestones: 1.0.3, 1.1 Apr 25, 2018
@fredRos fredRos removed this from the 1.1 milestone May 15, 2018
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants