← Back
Editing: _pade.cpython-311.pyc
� d�c# � �@ � d dl mZmZmZmZmZmZ d dlmZ dgZ dd�Z dS )� )�zeros�asarray�eye�poly1d�hstack�r_)�linalg�padeNc �J � t | � � } |�*t | � � dz |z }|dk rt d� � �|dk rt d� � �||z }|t | � � dz k rt d� � �| d|dz � } t |dz |dz | j �� � }t |dz |f| j �� � }t d|dz � � D ] }| d|� ddd� ||d|�f< � t |dz |dz � � D ] }| ||z |� ddd� ||dd�f< �!t ||f� � }t j || � � }|d|dz � } t d ||dz d� f } t | ddd� � � t | ddd� � � fS ) a� Return Pade approximation to a polynomial as the ratio of two polynomials. Parameters ---------- an : (N,) array_like Taylor series coefficients. m : int The order of the returned approximating polynomial `q`. n : int, optional The order of the returned approximating polynomial `p`. By default, the order is ``len(an)-1-m``. Returns ------- p, q : Polynomial class The Pade approximation of the polynomial defined by `an` is ``p(x)/q(x)``. Examples -------- >>> import numpy as np >>> from scipy.interpolate import pade >>> e_exp = [1.0, 1.0, 1.0/2.0, 1.0/6.0, 1.0/24.0, 1.0/120.0] >>> p, q = pade(e_exp, 2) >>> e_exp.reverse() >>> e_poly = np.poly1d(e_exp) Compare ``e_poly(x)`` and the Pade approximation ``p(x)/q(x)`` >>> e_poly(1) 2.7166666666666668 >>> p(1)/q(1) 2.7179487179487181 N� r z.Order of q <m> must be smaller than len(an)-1.z&Order of p <n> must be greater than 0.z0Order of q+p <m+n> must be smaller than len(an).)�dtype���g �?) r �len� ValueErrorr r r �ranger r �solver r ) �an�m�n�N�Akj�Bkj�row�C�pq�p�qs �9/usr/lib/python3/dist-packages/scipy/interpolate/_pade.pyr r s� � �N ����B��y� ��G�G�a�K�!�O�� �q�5�5� � M�N�N� N��1�u�u��A�B�B�B� �A��A��3�r�7�7�1�9�}�}��K�L�L�L� �D�Q�q�S�D��B� �a��c�1�Q�3�b�h� '� '� '�C� ��1��a���� )� )� )�C��Q��!��}�}� *� *�� �T�c�T�(�D�D�b�D�)�)��C����H� � ��Q�q�S�!�A�#��� ,� ,���#�a�%��)�}�d�d��d�+�+��C����E� � ���S�z���A� � �a�� � �B� �4�A�a�C�4��A� �3��1�Q�3�4�4��=��A� �!�D�D�b�D�'�?�?�F�1�T�T�r�T�7�O�O� +�+� )N) �numpyr r r r r r �scipyr �__all__r � r r �<module>r$ sl �� 9� 9� 9� 9� 9� 9� 9� 9� 9� 9� 9� 9� 9� 9� 9� 9� � � � � � � �(��<,� <,� <,� <,� <,� <,r
Save File
Cancel