
    tKg+                     <    d dl Zd dlmZ d dlmZ d dlmZ d Z	d Z
y)    N)assert_allclose)_sph_harm_allc                     d } d }d }d }| |||g}g d}g d}t        j                  ddt         j                  z        }t        j                  dt         j                        }t        j                  ||      \  }}t	        |||      D ]:  \  }	}}t        t        j                  ||||       |	||      d	d	d
| d| d       < y )Nc                 T    dt        j                  dt         j                  z        z  S )N      ?   )npsqrtpithetaphis     e/home/alanp/www/video.onchill/myenv/lib/python3.12/site-packages/scipy/special/tests/test_sph_harm.pyY00z!test_first_harmonics.<locals>.Y00   s    2771RUU7###    c                     dt        j                  ddt         j                  z  z        z  t        j                  d| z        z  t        j                  |      z  S )Nr         y             r	   r
   r   expsinr   s     r   Yn11z"test_first_harmonics.<locals>.Yn11   sA    2771ag;''s5y(99"&&+EEr   c                     dt        j                  dt         j                  z        z  t        j                  |      z  S )Nr   r   )r	   r
   r   cosr   s     r   Y01z!test_first_harmonics.<locals>.Y01   s+    2771RUU7##BFF3K//r   c                     dt        j                  ddt         j                  z  z        z  t        j                  d| z        z  t        j                  |      z  S )Ng      r   r   y              ?r   r   s     r   Y11z!test_first_harmonics.<locals>.Y11   sA    BGGAqwK((5)99"&&+EEr   )r   r   r   )r   r   r   r   r   r   gV瞯<zY^_z
 incorrect)rtolatolerr_msg)r	   linspacer   meshgridzipr   scsph_harm)
r   r   r   r   harmsmnr   r   harms
             r   test_first_harmonicsr,      s    $F0F $S!EAAKK1RUU7#E
++a
CUC(JE3%A&
aAq%5UC(""$QCq: 6	8 'r   c                     d} t        j                  ddt         j                  z        }t        j                  dt         j                        }t        d| z  | ||      }dD ]  }dD ]  }||k  rt	        j
                  ||||      }nd}t         j                  j                  |||f   |d       ||k  rt	        j
                  | |||      }nd}t         j                  j                  || |f   |d         y )N2   r   r   )r   r   r      
      r.   gh㈵>)r    )r	   r#   r   r   r&   r'   testingr   )n_maxr   r   y_actualr*   r)   	y_desireds          r   test_all_harmonicsr6   )   s    EKK1ruu9%E
++a
CQYuc:H%)AQKK1eS9		JJ&&x1~y&OQKKAuc:		JJ&&xA	%&P * &r   )numpyr	   numpy.testingr   scipy.specialspecialr&   scipy.special._basicr   r,   r6    r   r   <module>r=      s     )  .8DQr   