def __init__(self,*args,**kwargs):
"""
NAME:
__init__
PURPOSE:
initialize an actionAngleAxi object
INPUT:
Either:
a) R,vR,vT
b) Orbit instance: initial condition used if that's it, orbit(t)
if there is a time given as well
pot= potential or list of potentials (planarPotentials)
verticalPot= the vertical Potential
gamma= (default=1.) replace Lz by Lz+gamma Jz in effective potential (if there is no vertical potential, this is set to zero)
OUTPUT:
HISTORY:
2010-12-01 - Written - Bovy (NYU)
"""
actionAngle.__init__(self,*args,**kwargs)
if not kwargs.has_key('pot'): #pragma: no cover
raise IOError("Must specify pot= for actionAngleAxi")
self._pot= kwargs['pot']
if kwargs.has_key('verticalPot'):
kwargs.pop('pot')
actionAngleVertical.__init__(self,*args,pot=kwargs['verticalPot'],
**kwargs)
if kwargs.has_key('gamma'):
self._gamma= kwargs['gamma']
else:
self._gamma= 1.
else:
self._gamma= 0.
return None
def __init__(self,*args,**kwargs):
"""
NAME:
__init__
PURPOSE:
initialize an actionAngleAxi object
INPUT:
Either:
a) R,vR,vT
b) Orbit instance: initial condition used if that's it, orbit(t)
if there is a time given as well
pot= potential or list of potentials (planarPotentials)
verticalPot= the vertical Potential
gamma= (default=1.) replace Lz by Lz+gamma Jz in effective potential (if there is no vertical potential, this is set to zero)
OUTPUT:
HISTORY:
2010-12-01 - Written - Bovy (NYU)
"""
actionAngle.__init__(self,*args,**kwargs)
if not kwargs.has_key('pot'):
raise IOError("Must specify pot= for actionAngleAxi")
self._pot= kwargs['pot']
if kwargs.has_key('verticalPot'):
kwargs.pop('pot')
actionAngleVertical.__init__(self,*args,pot=kwargs['verticalPot'],
**kwargs)
if kwargs.has_key('gamma'):
self._gamma= kwargs['gamma']
else:
self._gamma= 1.
else:
self._gamma= 0.
return None
def __init__(self,*args,**kwargs):
"""
NAME:
__init__
PURPOSE:
initialize an actionAngle object
INPUT:
Either:
a) R,vR,vT
b) Orbit instance: initial condition used if that's it, orbit(t)
if there is a time given as well
OUTPUT:
HISTORY:
2010-07-11 - Written - Bovy (NYU)
"""
actionAngle.__init__(self,*args,**kwargs)
if kwargs.has_key('verticalPot'):
actionAngleVertical.__init__(self,*args,pot=kwargs['verticalPot'],
**kwargs)
return None
def __init__(self, *args, **kwargs):
"""
NAME:
__init__
PURPOSE:
initialize an actionAnglePower object, power-law is vc=vo(r)^beta
INPUT:
Either:
a) R,vR,vT
b) Orbit instance: initial condition used if that's it, orbit(t)
if there is a time given as well
beta= power-law index
OUTPUT:
HISTORY:
2010-07-11 - Written - Bovy (NYU)
"""
actionAngle.__init__(self, *args, **kwargs)
if kwargs.has_key('verticalPot'):
actionAngleVertical.__init__(self,
*args,
pot=kwargs['verticalPot'],
**kwargs)
if not kwargs.has_key('beta'):
raise IOError("Must specify beta= for actionAnglePower")
self._beta = kwargs['beta']
if self._beta < 0.:
self._signbeta = -1
else:
self._signbeta = 1
self._absbeta = nu.fabs(self._beta)
(rperi, rap) = self.calcRapRperi()
self._e = (rap - rperi) / (rap + rperi)
if self._e == 0.:
self._X = 0.
self._Y = 0.
else:
self._X= (1.-self._e)**2./4./self._e\
*(((1.+self._e)/(1.-self._e))**(2.*self._beta+2.)-1.)
self._Y= (1.+self._e)**2./4./self._e\
*(1.-((1.-self._e)/(1.+self._e))**(2.*self._beta+2.))
return None
def __init__(self, *args, **kwargs):
"""
NAME:
__init__
PURPOSE:
initialize an actionAngle object
INPUT:
Either:
a) R,vR,vT
b) Orbit instance: initial condition used if that's it, orbit(t)
if there is a time given as well
OUTPUT:
HISTORY:
2010-07-11 - Written - Bovy (NYU)
"""
actionAngle.__init__(self, *args, **kwargs)
if kwargs.has_key('verticalPot'):
actionAngleVertical.__init__(self,
*args,
pot=kwargs['verticalPot'],
**kwargs)
return None
def __init__(self,*args,**kwargs):
"""
NAME:
__init__
PURPOSE:
initialize an actionAnglePower object, power-law is vc=vo(r)^beta
INPUT:
Either:
a) R,vR,vT
b) Orbit instance: initial condition used if that's it, orbit(t)
if there is a time given as well
beta= power-law index
OUTPUT:
HISTORY:
2010-07-11 - Written - Bovy (NYU)
"""
actionAngle.__init__(self,*args,**kwargs)
if kwargs.has_key('verticalPot'):
actionAngleVertical.__init__(self,*args,pot=kwargs['verticalPot'],
**kwargs)
if not kwargs.has_key('beta'):
raise IOError("Must specify beta= for actionAnglePower")
self._beta= kwargs['beta']
if self._beta < 0.:
self._signbeta= -1
else:
self._signbeta= 1
self._absbeta= nu.fabs(self._beta)
(rperi,rap)= self.calcRapRperi()
self._e= (rap-rperi)/(rap+rperi)
if self._e == 0.:
self._X= 0.
self._Y= 0.
else:
self._X= (1.-self._e)**2./4./self._e\
*(((1.+self._e)/(1.-self._e))**(2.*self._beta+2.)-1.)
self._Y= (1.+self._e)**2./4./self._e\
*(1.-((1.-self._e)/(1.+self._e))**(2.*self._beta+2.))
return None