def _set_single_mass(self, new_mass):
"""
Initializes total_mass and epsilon if only one mass componenet
is defined (i.e. a point lens).
"""
if isinstance(new_mass, u.Quantity):
if new_mass.unit.physical_type == 'dimensionless':
new_mass *= u.solMass
elif new_mass.unit.physical_type != 'mass':
msg = 'wrong physical_type of new total_mass: {:}'
raise ValueError(msg.format(new_mass.unit.physical_type))
self._total_mass = new_mass
else:
self._total_mass = new_mass * u.solMass
self._epsilon = np.array([1.0])
评论列表
文章目录