← Back
Editing: _rotation_groups.cpython-311.pyc
� d�cF � �D � d dl Zd dlmZ d� Zd� Zd� Zd d�Zd d�Z d d �Z dS ) � N)�goldenc �J � t | � � � � � }d}dt z }t dz }t j g | | | dg�| | | dg�| | d| g�| | d| g�| | | dg�| | | dg�| | d| g�| | d| g�| d| | g�| d| | g�| d| | g�| d| | g�| | d| g�| | d| g�| | | dg�| | | dg�| | d| g�| | d| g�| | | dg�| | | dg�| d| | g�| d| | g�| d| | g�| d| | g�| | | dg�| | | dg�| | d| g�| | d| g�| | | dg�| | | dg�| | d| g�| | d| g�| d| | g�| d| | g�| d| | g�| d| | g�d| | | g�d| | | g�d| | | g�d| | | g�d| | | g�d| | | g�d| | | g�d| | | g�d| | | g�d| | | g�d| | | g�d| | | g�� � }| � t j ||f� � � � S )N� �?� r )� tetrahedral�as_quat�phi�np�array� from_quat� concatenate)�cls�g1�a�b�c�g2s �J/usr/lib/python3/dist-packages/scipy/spatial/transform/_rotation_groups.py� icosahedralr s � � �S� � � !� !� #� #�B� �A� �c� �A� �a��A� �� /$�Q�B���Q�B��?� /$��B���Q�B��?�/$��B���A��r�?�/$� �B���A��r�?�/$� �B���Q�B��?� /$� �B���Q�B��?� /$� �B���A��r�?� /$� �B���A��r�?�/$� �B��A�2��r�?�/$� �B��A�2��r�?�/$� �B��A�2��r�?�/$� �B��A�2��r�?�/$� �B���A��r�?�/$� �B���A��r�?�/$� �B���Q�B��?� /$� �B���Q�B��?� /$� �B���A��r�?�!/$�" �B���A��r�?�#/$�$ �B���Q�B��?�%/$�& �B���Q�B��?�'/$�( �B��A�2��r�?�)/$�* �B��A�2��r�?�+/$�, �B��A�2��r�?�-/$�. �B��A�2��r�?�//$�0 �B���Q�B��?�1/$�2 �B���Q�B��?�3/$�4 �B���A��r�?�5/$�6 �B���A��r�?�7/$�8 �B���Q�B��?�9/$�: �B���Q�B��?�;/$�< �B���A��r�?�=/$�> �B���A��r�?�?/$�@ �B��A�2��r�?�A/$�B �B��A�2��r�?�C/$�D �B��A�2��r�?�E/$�F �B��A�2��r�?�G/$�H ��r�A�2��r�?�I/$�J ��r�A�2��r�?�K/$�L ��r�A�2��r�?�M/$�N ��r�A�2��r�?�O/$�P ��r�A�2��r�?�Q/$�R ��r�A�2��r�?�S/$�T ��r�A�2��r�?�U/$�V ��r�A�2��r�?�W/$�X ��r�A�2��r�?�Y/$�Z ��r�A�2��r�?�[/$�\ ��r�A�2��r�?�]/$�^ ��r�A�2��r�?�_/$� / %� / %�B�` �=�=�� ��R��1�1� 2� 2�2� c � � t | � � � � � }t j d� � dz }t j | dd| gd| d| gdd| | gdd| | gd| d| g| dd| gd| | dgd| | dg| d| dg| d| dg| | ddg| | ddgg � � }| � t j ||f� � � � S )Nr r )r r r �sqrtr r r �r r r r s r � octahedralr = s � � �S� � � !� !� #� #�B� ��� � �Q��A� ��Q�B��1�q�b�>���r�1�q�b�>��q�1�"�q�b�>��q�1�"�q�b�>���r�1�q�b�>��B��1�q�b�>���r�A�2�q�>���r�A�2�q�>��B��A�2�q�>��B��A�2�q�>��B���A�q�>��B���A�q�>� #� $� $�B� �=�=�� ��R��1�1� 2� 2�2r c �( � t j d� � }d}t j || | | g|| | | g|| | | g|| | | g|| | | g|| | | g|| | | g|| | | gg� � }| � t j ||f� � � � S )N� r )r �eyer r r r s r r r O s� � � �����B� �A� ��A��r�A�2��r�?���r�A�2��r�?���r�A�2��r�?���r�A�2��r�?���r�A�2��r�?���r�A�2��r�?���r�A�2��r�?���r�A�2��r�?�$� %� %�B� �=�=�� ��R��1�1� 2� 2�2r r c �� � t | ||� � � � � }t j dt j |d�� � }t j t j t j |� � t j |� � t j |� � g� � j z }t j ||d�� � }| � t j ||f� � � � S )Nr F��endpoint� ��axis) �cyclic� as_rotvecr �linspace�pi�vstack�zeros�cos�sin�T�roll� from_rotvecr )r �nr# r �thetas�rvr s r �dicyclicr2 ] s� � � ��Q�� � � '� '� )� )�B� �[��B�E�1�u� 5� 5� 5�F� ����B�H�Q�K�K���� � ���� � �H�I�I�K� K�B� ���T�� "� "� "�B� �?�?�2�>�2�r�(�3�3� 4� 4�4r c � � t j ddt j z |d�� � }t j |t j |� � t j |� � g� � j }| � t j ||d�� � � � S )Nr r Fr r! r" )r r& r' r( r) r, r. r- )r r/ r# r0 r1 s r r$ r$ f sl � � �[��A���I�q�5� 9� 9� 9�F� ��F�B�H�Q�K�K���!���5� 6� 6� 8�B� �?�?�2�7�2�t�!�4�4�4� 5� 5�5r �Zc �� � t |t � � st d� � �g d�}||vr%t dd� |� � z � � �|dv r|}d}nY|d d� dv r>|dd � � � � r"|d d� }t |dd � � � }nt d� � �|dk rt d � � �d � |� � � � � }|d k rt | � � S |d k rt | � � S |d k rt | � � S |dk rt | ||�� � S |dk rt | ||�� � S J �)Nz!`group` argument must be a string)�x�y�z�X�Yr4 z`axis` must be one of z, )�I�Or, r! )�C�Dz0`group` must be one of 'I', 'O', 'T', 'Dn', 'Cn'z Group order must be positive�xyzr; r< r, r> r" r= ) � isinstance�str� ValueError�join�isdigit�int�index�lowerr r r r2 r$ )r �groupr# �permitted_axes�symbol�orders r � create_grouprL l s� � � �e�S� !� !� >��<�=�=�=�3�3�3�N� �>�!�!��1�D�I�I�n�4M�4M�M�N�N�N� �� � ������ �r��r��j� � �U�1�2�2�Y�%6�%6�%8�%8� ��r��r����E�!�"�"�I������K�L�L�L� �q�y�y��7�8�8�8� �;�;�t�z�z�|�|� $� $�D� ��}�}��3� � � � �3� � ��#��� � �3� � ��3� � � � �3� � ���U��.�.�.�.� �3� � ��c�5�t�,�,�,�,�� r )r )r4 ) �numpyr �scipy.constantsr r r r r r2 r$ rL � r r �<module>rP s� �� � � � � )� )� )� )� )� )�53� 53� 53�p3� 3� 3�$ 3� 3� 3� 5� 5� 5� 5�6� 6� 6� 6� � � � � � r
Save File
Cancel