optic.py 文件源码

def apexes(self):
        '''Returns the positions of the apexes of HR and AR as a tuple.'''
        if self.HRK == 0.:
            apex1 = self.HRCenter
        else:
            theta1 = np.arcsin(self.Dia * self.ARK/2.)\
                if np.abs(self.Dia * self.ARK/2.) < 1. else np.pi/2.
            apex1 = self.HRCenter - (1-np.cos(theta1))*self.HRNorm/self.HRK

        if self.ARK == 0.:
            apex2 = self.ARCenter
        else:
            theta2 = np.arcsin(self.Dia * self.ARK/(2.*np.cos(self.Wedge)))\
                if np.abs(self.Dia * self.ARK/(2.*np.cos(self.Wedge))) < 1.\
                else np.pi/2.
            apex2 = self.ARCenter - (1-np.cos(theta2))*self.ARNorm/self.ARK

        return apex1, apex2