U
    Y»{j>] ã                   @   sÆ
  d Z ddlZddlZddlZddlmZ ej ej	j
¡sFejddd� dd„ Zd	d
„ Zejded�Zdd„ Zdd„ Zdd„ Zdd„ Zdd„ Zdd„ Zdd„ Zdd„ Zdd„ Zdd „ Zd!d"„ Zd#d$„ Zd%d&„ Zd'd(„ Zd)d*„ Zd+d,„ Z d-d.„ Z!d/d0„ Z"d1d2„ Z#d3d4„ Z$d5d6„ Z%d7d8„ Z&d9d:„ Z'd;d<„ Z(d=d>„ Z)d?d@„ Z*dAdB„ Z+dCdD„ Z,dEdF„ Z-dGdH„ Z.dIdJ„ Z/dKdL„ Z0dMdN„ Z1dOdP„ Z2ej3j4ejdQk dRdS�dTdU„ ƒZ5dVdW„ Z6dXdY„ Z7dZd[„ Z8d\d]„ Z9d^d_„ Z:d`da„ Z;dbdc„ Z<ddde„ Z=dfdg„ Z>dhdi„ Z?djdk„ Z@dldm„ ZAdndo„ ZBdpdq„ ZCdrds„ ZDdtdu„ ZEdvdw„ ZFdxdy„ ZGdzd{„ ZHd|d}„ ZId~d„ ZJd€d�„ ZKd‚dƒ„ ZLd„d…„ ZMd†d‡„ ZNdˆd‰„ ZOdŠd‹„ ZPdŒd�„ ZQdŽd�„ ZRd�d‘„ ZSd’d“„ ZTd”d•„ ZUd–d—„ ZVd˜d™„ ZWdšd›„ ZXdœd�„ ZYdždŸ„ ZZd d¡„ Z[d¢d£„ Z\d¤d¥„ Z]d¦d§„ Z^d¨d©„ Z_dªd«„ Z`d¬d­„ Zad®d¯„ Zbd°d±„ Zcd²d³„ Zdd´dµ„ Zed¶d·„ Zfd¸d¹„ Zgdºd»„ Zhd¼d½„ Zid¾d¿„ ZjdÀdÁ„ ZkdÂdÃ„ ZldÄdÅ„ ZmdÆdÇ„ ZndÈdÉ„ ZodÊdË„ ZpdÌdÍ„ ZqdÎdÏ„ ZrdÐdÑ„ ZsdÒdÓ„ Ztej3j4ejdQk dRdS�dÔdÕ„ ƒZudÖd×„ ZvdØdÙ„ ZwdÚdÛ„ ZxdÜdÝ„ ZydÞdß„ Zzdàdá„ Z{dâdã„ Z|dädå„ Z}dædç„ Z~dèdé„ Zdêdë„ Z€dìdí„ Z�dîdï„ Z‚dðdñ„ Zƒdòdó„ Z„dôdõ„ Z…död÷„ Z†dødù„ Z‡dúdû„ Zˆdüdý„ Z‰dþdÿ„ ZŠej3j4ejdQk dRdS��d �d„ ƒZ‹�d�d„ ZŒ�d�d„ Z��d�d„ ZŽ�d�d	„ Z��d
�d„ Z��d�d„ Z‘�d�d„ Z’�d�d„ Z“�d�d„ Z”�d�d„ Z•�d�d„ Z–�d�d„ Z—�d�d„ Z˜�d�d„ Z™�d�d„ Zš�d �d!„ Z›�d"�d#„ Zœ�d$�d%„ Z��d&�d'„ Zž�d(�d)„ ZŸ�d*�d+„ Z �d,�d-„ Z¡�d.�d/„ Z¢�d0�d1„ Z£�d2�d3„ Z¤�d4�d5„ Z¥�d6�d7„ Z¦�d8�d9„ Z§�d:�d;„ Z¨�d<�d=„ Z©�d>�d?„ Zª�d@�dA„ Z«�dB�dC„ Z¬�dD�dE„ Z­�dF�dG„ Z®�dH�dI„ Z¯�dJ�dK„ Z°�dL�dM„ Z±�dN�dO„ Z²�dP�dQ„ Z³�dR�dS„ Z´�dT�dU„ Zµ�dV�dW„ Z¶�dX�dY„ Z·�dZ�d[„ Z¸�d\�d]„ Z¹�d^�d_„ Zº�d`�da„ Z»�db�dc„ Z¼�dd�de„ Z½�df�dg„ Z¾�dh�di„ Z¿�dj�dk„ ZÀ�dl�dm„ ZÁ�dn�do„ ZÂ�dp�dq„ ZÃ�dr�ds„ ZÄ�dt�du„ ZÅ�dv�dw„ ZÆ�dx�dy„ ZÇ�dz�d{„ ZÈ�d|�d}„ ZÉ�d~�d„ ZÊ�d€�d�„ ZË�d‚�dƒ„ ZÌ�d„�d…„ ZÍ�d†�d‡„ ZÎ�dˆ�d‰„ ZÏ�dŠ�d‹„ ZÐ�dŒ�d�„ ZÑ�dŽ�d�„ ZÒ�d��d‘„ ZÓ�d’�d“„ ZÔ�d”�d•„ ZÕ�d–�d—„ ZÖ�d˜�d™„ Z×�dš�d›„ ZØ�dœ�d�„ ZÙ�dž�dŸ„ ZÚ�d �d¡„ ZÛ�d¢�d£„ ZÜ�d¤�d¥„ ZÝ�d¦�d§„ ZÞ�d¨�d©„ Zß�dª�d«„ Zà�d¬�d­„ Zá�d®�d¯„ Zâ�d°�d±„ Zã�d²�d³„ Zä�d´�dµ„ Zå�d¶�d·„ Zæ�d¸�d¹„ Zç�dº�d»„ Zè�d¼�d½„ Zé�d¾�d¿„ Zê�dÀ�dÁ„ Zë�dÂ�dÃ„ Zì�dÄ�dÅ„ Zí�dÆ�dÇ„ Zî�dÈ�dÉ„ Zï�dÊ�dË„ Zð�dÌ�dÍ„ Zñ�dÎ�dÏ„ Zò�dÐ�dÑ„ Zó�dÒ�dÓ„ Zô�dÔ�dÕ„ Zõ�dÖ�d×„ Zö�dØ�dÙ„ Z÷�dÚ�dÛ„ Zø�dÜ�dÝ„ Zù�dÞ�dß„ Zú�dà�dá„ Zû�dâ�dã„ Zü�dä�då„ Zý�dæ�dç„ Zþ�dè�dé„ Zÿ�dê�dë„ �Z �dì�dí„ �Z�dî�dï„ �Z�dð�dñ„ �Z�dò�dó„ �Z�dô�dõ„ �Z�dö�d÷„ �Zej3j4ejdQk dRdS��dø�dù„ ƒ�Zej3j4ej�dúk �dûdS��dü�dý„ ƒ�Zej3j4ejdQk dRdS��dþ�dÿ„ ƒ�Z	�d �d„ �Z
dS (  aê  
Basic tests of the ERFA library routines.

These are just the tests bundled with ERFA itself, in ``t_erfa_c.c``,
but translated to python to make sure that the code returns the same numbers.

The viv and vid functions emulate the corresponding functions in ``t_erfa_c.c``.

The tests are skipped if a system library is used that does not match in
version number, since in that case the precise numbers may have changed
(e.g., due to the earth orientation changes between 1.7.2 and 1.7.3).
é    N)Úufuncz3skipping for system library with mismatched versionT)Zallow_module_levelc                 C   s*   | |ks&t |› d|› d|› d| › �ƒ‚dS )zValidate an integer result.ú	 failed: ú want ú got N)ÚAssertionError)ZivalZivalokÚfuncÚtestÚ_© r
   úT/var/www/html/bwave/python/venv/lib/python3.8/site-packages/erfa/tests/test_ufunc.pyÚviv!   s    r   c                 C   sX   | | }|dksTt |ƒt |ƒksTt|› d|› d|d›d| d›dt || ƒd›d�
ƒ‚d	S )
zValidate a double result.ç        r   r   z.20gr   z (1/z.3gú)N)Úabsr   )ÚvalZvalokÚdvalr   r   r	   Úar
   r
   r   Úvvd&   s    .ÿr   r
   )Zdtypec                  C   sx   t  dd¡\} }t| d dddtƒ t|d dddtƒ t|d	 d
ddtƒ t|d dddtƒ t|d dddtƒ d S )Né   gÃõ(\�Â@r   ó   +Úa2afÚsé†   Ú0é   é   Ú1é   é   Ú2é   iê%  Ú3)Ú
erfa_ufuncr   r   Ústatus)r   Zidmsfr
   r
   r   Ú	test_a2af2   s    r$   c                  C   sx   t  dd¡\} }t| d dddtƒ t|d dddtƒ t|d	 d
ddtƒ t|d dddtƒ t|d dddtƒ d S )Nr   gvÃ¶EÀr   ó   -Úa2tfr   é   r   r   r   r   r   é   r   r    iT  r!   )r"   r&   r   r#   ©r   Úihmsfr
   r
   r   Ú	test_a2tf;   s    r+   c                  C   s²   t  dt¡} t  dt¡}t  dt¡}d| d< d| d< d| d< d|d< d	|d< d
|d< d}d}t | |||¡}t|d ddddtƒ t|d ddddtƒ t|d ddddtƒ d S )N©r    g›Ý—°Klè¿r   gVpºølzã¿r   g`/¾ñì¾Ë¿r   g=Ö÷õö>g²M?¥ÿ[¿g™w'UTA¿gÓì…äoþï?gÔVYýÿÿï?gd�
Õkè¿çê-�™—q=Úabr   gK{†iìzã¿r   g}•Z]Ü¿Ë¿r!   )ÚnpÚemptyÚfloatr"   r.   r   r#   )ZpnatÚvZpprr   Úbm1r
   r
   r   Útest_abD   s    r4   c                  C   sF   d} d}d}t  | ||¡\}}t|ddddtƒ t|dddd	tƒ d S )
Ng      @çš™™™™™ñ?gffffffæ?g<U�üâ?ç›+¡†›„=Úae2hdÚhgRåV@¼Ãî?Úd)r"   r7   r   r#   )r   ÚeÚpr8   r9   r
   r
   r   Ú
test_ae2hdV   s    r<   c                  C   s:   t  dddd¡\} }t| ddddtƒ t|d	dd
tƒ d S )Nú-é-   é   g333333;@gÍ`€[
Bé¿r-   Úaf2ar   r   Új)r"   r@   r   r#   r   ©r   rA   r
   r
   r   Ú	test_af2a_   s    rC   c                   C   s   t t d¡ddddtƒ d S )Ngš™™™™™¹¿g²ÆÝí”»@r-   ÚanpÚ )r   r"   rD   r#   r
   r
   r
   r   Útest_anpe   s    rF   c                   C   s   t t d¡ddddtƒ d S )Ng      Àg1Zˆ¨öC@r-   ÚanpmrE   )r   r"   rG   r#   r
   r
   r
   r   Ú	test_anpmi   s    rH   c                  C   s²  t  dtj¡} t  dt¡}d}d}d| d d< d| d d	< d
| d d< d| d d< d| d d	< d| d d< d|d< d|d	< d|d< t ||| |¡}| t j¡}t|j	ddddt
ƒ t|jd ddddt
ƒ t|jd	 ddddt
ƒ t|jd d
dddt
ƒ t|jd ddddt
ƒ t|jd	 ddddt
ƒ t|jd dddd t
ƒ t|jd!ddd"t
ƒ t|jd d#d$dd%t
ƒ t|jd	 d&d$dd't
ƒ t|jd d(d$dd)t
ƒ t|jd*ddd+t
ƒ t|jd d d,d-dd.t
ƒ t|jd	 d d-d-dd/t
ƒ t|jd d d-d-dd0t
ƒ t|jd d	 d-d-dd1t
ƒ t|jd	 d	 d,d-dd2t
ƒ t|jd d	 d-d-dd3t
ƒ t|jd d d-d-dd4t
ƒ t|jd	 d d-d-dd5t
ƒ t|jd d d,d-dd6t
ƒ d S )7Nr
   r,   ç   À2½BAç\ùò#ú¬Ù?ç¿îtç‰×ì?r;   r   çq+É©¹¶Ú¿r   çžÀxm*Ç¿r   ç‚��+l~?r2   çÃÎÿ0£ÆŒ?ç°SXÙNòx?ç„Ya-Pèì?çuçÕ•Ú¿ç¹ÉT3þÇ¿çe= *|M)@ç•dyáý¥=ÚapcgÚpmtr-   úeb(1)úeb(2)úeb(3)çð�8«›ì?úeh(1)ç�¾€ YOÚ¿úeh(2)çëœ�~àÏÆ¿úeh(3)ç1t¨Ý*ð?Úemç«)(µt}?ç¼‰Ø—²Òœ<úv(1)çU û;èE?úv(2)çKç½.q?úv(3)çEúgýÿÿï?r3   ç      ð?r   úbpn(1,1)úbpn(2,1)úbpn(3,1)úbpn(1,2)úbpn(2,2)úbpn(3,2)úbpn(1,3)úbpn(2,3)úbpn(3,3))r/   r0   r"   Údt_pvr1   rV   ÚviewÚrecarrayr   rW   r#   ÚebÚehrb   r2   r3   Úbpn)ÚebpvÚehpÚdate1Údate2Úastromr
   r
   r   Ú	test_apcgm   sH    r€   c                  C   s4  d} d}t  | |¡}| tj¡}t|jddddtƒ t|jd dd	dd
tƒ t|jd dd	ddtƒ t|jd dd	ddtƒ t|j	d dd	ddtƒ t|j	d dd	ddtƒ t|j	d dd	ddtƒ t|j
dd	ddtƒ t|jd ddddtƒ t|jd ddddtƒ t|jd ddddtƒ t|jd d	dd!tƒ t|jd d d"d#dd$tƒ t|jd d d#d#dd%tƒ t|jd d d#d#dd&tƒ t|jd d d#d#dd'tƒ t|jd d d"d#dd(tƒ t|jd d d#d#dd)tƒ t|jd d d#d#dd*tƒ t|jd d d#d#dd+tƒ t|jd d d"d#dd,tƒ d S )-NrI   rJ   rT   rU   Úapcg13rW   r   çùaPç‰×ì?r-   rX   r   çBUÔ©¹¶Ú¿rY   r   çÒˆ™wm*Ç¿rZ   ç�þŒ8«›ì?r\   çI.w YOÚ¿r^   çÐ©Î~àÏÆ¿r`   ç@”~¨Ý*ð?rb   çbµt}?rd   re   çy#;èE?rg   ç?IÔ½.q?ri   rj   r3   rk   r   rl   rm   rn   ro   rp   rq   rr   rs   rt   )r"   r�   rv   r/   rw   r   rW   r#   rx   ry   rb   r2   r3   rz   )r}   r~   r   r
   r
   r   Útest_apcg13”   s2    rŒ   c               	   C   sÄ  t  dtj¡} t  dt¡}d}d}d| d d< d| d d	< d
| d d< d| d d< d| d d	< d| d d< d|d< d|d	< d|d< d}d}d}t ||| ||||¡}| t j¡}t|j	ddddt
ƒ t|jd ddddt
ƒ t|jd	 ddddt
ƒ t|jd d
dddt
ƒ t|jd ddddt
ƒ t|jd	 d ddd!t
ƒ t|jd d"ddd#t
ƒ t|jd$ddd%t
ƒ t|jd d&d'dd(t
ƒ t|jd	 d)d'dd*t
ƒ t|jd d+d'dd,t
ƒ t|jd-ddd.t
ƒ t|jd d d/ddd0t
ƒ t|jd	 d d1ddd2t
ƒ t|jd d dddd3t
ƒ t|jd d	 d4ddd5t
ƒ t|jd	 d	 d6ddd7t
ƒ t|jd d	 dddd8t
ƒ t|jd d d9ddd:t
ƒ t|jd	 d d;ddd<t
ƒ t|jd d d=ddd>t
ƒ d S )?Nr
   r,   rI   rJ   rK   r;   r   rL   r   rM   r   rN   r2   rO   rP   rQ   rR   rS   ç¸Å+:áU?çç]‰´þ¾çPÜN0j`>rT   rU   ÚapcirW   r-   rX   rY   rZ   r[   r\   r]   r^   r_   r`   ra   rb   rc   rd   re   rf   rg   rh   ri   rj   r3   ç!Å1þÿï?rl   çÇ:Äšºj>rm   rn   ç}€pªgH¾ro   çÅÿÿÿï?rp   rq   çKÅj:áU¿rr   ç
3MÇ´þ>rs   ç+3Š1þÿï?rt   )r/   r0   r"   ru   r1   r�   rv   rw   r   rW   r#   rx   ry   rb   r2   r3   rz   )r{   r|   r}   r~   ÚxÚyr   r   r
   r
   r   Ú	test_apci°   sN    rš   c                  C   sJ  d} d}t  | |¡\}}| tj¡}t|jddddtƒ t|jd dd	dd
tƒ t|jd dd	ddtƒ t|jd dd	ddtƒ t|j	d dd	ddtƒ t|j	d dd	ddtƒ t|j	d dd	ddtƒ t|j
dd	ddtƒ t|jd ddddtƒ t|jd ddddtƒ t|jd ddddtƒ t|jd d	dd!tƒ t|jd d d"d	dd#tƒ t|jd d d$d	dd%tƒ t|jd d d&d	dd'tƒ t|jd d d(d	dd)tƒ t|jd d d*d	dd+tƒ t|jd d d,d	dd-tƒ t|jd d d.d	dd/tƒ t|jd d d0d	dd1tƒ t|jd d d2d	dd3tƒ t|d4d	dd5tƒ d S )6NrI   rJ   rT   rU   Úapci13rW   r   r‚   r-   rX   r   rƒ   rY   r   r„   rZ   r…   r\   r†   r^   r‡   r`   rˆ   rb   r‰   rd   re   rŠ   rg   r‹   ri   rj   r3   g-¦¾Uþÿï?rl   g3œM¨Q$f>rm   g•MÀ\¥T?rn   g‚FD‰y‰K¾ro   g	Ýÿÿÿï?rp   gCéµ®þ¥÷¾rq   gT4BÀ\¥T¿rr   g.EpWü¥÷>rs   g<²›Uþÿï?rt   çJ÷IÐ	Ãg¿Úeo)r"   r›   rv   r/   rw   r   rW   r#   rx   ry   rb   r2   r3   rz   )r}   r~   r   r�   r
   r
   r   Útest_apci13Ú   s4    rž   c                  C   s®  t  dtj¡} t  dt¡}d}d}d| d d< d| d d	< d
| d d< d| d d< d| d d	< d| d d< d|d< d|d	< d|d< d}d}d}d}d}d}	d}
d}d}d}d}d}t ||| |||||||	|
|||||¡}| t j¡}t|j	dd d!d"t
ƒ t|jd d#d$d!d%t
ƒ t|jd	 d&d$d!d't
ƒ t|jd d(d$d!d)t
ƒ t|jd d*d$d!d+t
ƒ t|jd	 d,d$d!d-t
ƒ t|jd d.d$d!d/t
ƒ t|jd0d$d!d1t
ƒ t|jd d2d3d!d4t
ƒ t|jd	 d5d3d!d6t
ƒ t|jd d7d3d!d8t
ƒ t|jd9d$d!d:t
ƒ t|jd d d;d$d!d<t
ƒ t|jd	 d d=d$d!d>t
ƒ t|jd d dd$d!d?t
ƒ t|jd d	 d@d$d!dAt
ƒ t|jd	 d	 dBd$d!dCt
ƒ t|jd d	 dd$d!dDt
ƒ t|jd d dEd$d!dFt
ƒ t|jd	 d dGd$d!dHt
ƒ t|jd d dId$d!dJt
ƒ t|jdKd$d!dLt
ƒ t|jdMdNd!dOt
ƒ t|jdPdNd!dQt
ƒ t|jdRd$d!dSt
ƒ t|jdTd$d!dUt
ƒ t|jddd!dVt
ƒ t|jdWd$d!dXt
ƒ t|jddYd!dZt
ƒ t|jdd[d!d\t
ƒ d S )]Nr
   r,   ç   @ ½BAçkbÍ
ï?çr~{g,ï¿r;   r   ç‘QË¿r   çUÉñÛx}·¿r   çà²PNÙm?r2   ç¤oW™�¿ç.'ªâ'f{¿çÓD†ò‘&ï¿çßœ±2‘ÇÊ¿ç4V¡”8·¿r�   rŽ   r�   çJÄç�Ì)	@çoQ„¾ãà¿çôÌ�ÔÜÀó¿ç     d¥@çaÆÁGd—�>ç¡‚#Y¤¾>çvÇ’5é™À½ç@: #~f*?çÊÐ?·¼±�¾çr)<¬E�*@rU   ÚapcorW   g×‹8�,ï¿r-   rX   g”±«ÍÛË¿rY   g£ÃÚ�€·¿rZ   g¡o(ï¿r\   g'ÌrŒÈÊ¿r^   g#Ð�YÅ;·¿r`   g6õâ{�þï?rb   g3!«1ûËõ>rd   re   ghØ¿&×y¿rg   gšm³+IA¿ri   ç¹žTýÿÿï?r3   r‘   rl   r’   rm   rn   r“   ro   r”   rp   rq   r•   rr   r–   rs   r—   rt   çàb„¾ãà¿Úalongçc ³>ç—ÔFFõg<Úxplç©Â¬”b¸>Úyplç˜®e¿W5î¿ÚsphiçV‚²HÕ?ÚcphiÚdiurabç’+âïÜð@ÚeralçVçž¯Ò<Úrefaç¬CÒÑ]r2<Úrefb)r/   r0   r"   ru   r1   r´   rv   rw   r   rW   r#   rx   ry   rb   r2   r3   rz   r·   rº   r¼   r¾   rÀ   rÁ   rÃ   rÅ   rÇ   )r{   r|   r}   r~   r˜   r™   r   ÚthetaÚelongÚphiÚhmÚxpÚypÚsprÅ   rÇ   r   r
   r
   r   Ú	test_apco÷   sr    (rÏ   c                  C   sL  d} d}d}d}d}d}d}d}d	}d
}	d}
d}t  | |||||||||	|
|¡\}}}| tj¡}t|jddddtƒ t|jd ddddtƒ t|jd ddddtƒ t|jd ddddtƒ t|j	d ddddtƒ t|j	d ddddtƒ t|j	d dddd tƒ t|j
d!ddd"tƒ t|jd d#d$dd%tƒ t|jd d&d$dd'tƒ t|jd d(d$dd)tƒ t|jd*ddd+tƒ t|jd d d,ddd-tƒ t|jd d d.ddd/tƒ t|jd d d0ddd1tƒ t|jd d d2ddd3tƒ t|jd d d4ddd5tƒ t|jd d d6ddd7tƒ t|jd d d8ddd9tƒ t|jd d d:ddd;tƒ t|jd d d<ddd=tƒ t|jd>ddd?tƒ t|jd@dAddBtƒ t|jdCdAddDtƒ t|jdEdddFtƒ t|jdGdddHtƒ t|jddddItƒ t|jdJdddKtƒ t|jdLdMddNtƒ t|jdOdPddQtƒ t|dRdSddTtƒ t|dddUtƒ d S )VNrŸ   ç¿ïQ!ï?çÓ–x@ÙÃ?r«   r¬   r­   r®   r¯   ç     Ø†@çš™™™™™)@çáz®Gáâ?çš™™™™™á?gþ&<¬E�*@rU   Úapco13rW   r   gÂˆ]8�,ï¿r-   rX   r   gë,¿ÍÛË¿rY   r   goÅö�€·¿rZ   gÊo­o(ï¿r\   g_‰ŒÈÊ¿r^   gß¡YÅ;·¿r`   gTØ{�þï?rb   gM[ó1ûËõ>rd   re   gÓŽå'×y¿rg   gy©—+IA¿ri   g¸žTýÿÿï?r3   gÅ1þÿï?rl   gàzÄšºj>rm   g>›g:áU?rn   gÖ	ªgH¾ro   r”   rp   gKK‰´þ¾rq   gÑš¦:áU¿rr   g6 yÇ´þ>rs   g!3Š1þÿï?rt   r¶   r·   r¸   r¹   rº   r»   r¼   r½   r¾   r¿   rÀ   rÁ   ç³€•ðÜð@rÃ   ç4±;"~f*?rÄ   rÅ   çßÌ\¸¼±�¾rÆ   rÇ   çl­ËŒŒ¾h¿r6   r�   rA   )r"   rÖ   rv   r/   rw   r   rW   r#   rx   ry   rb   r2   r3   rz   r·   rº   r¼   r¾   rÀ   rÁ   rÃ   rÅ   rÇ   r   )Úutc1Úutc2Údut1rÉ   rÊ   rË   rÌ   rÍ   ÚphpaÚtcÚrhÚwlr   r�   rA   r
   r
   r   Útest_apco133  s\    &râ   c                  C   s
  t  dtj¡} t  dtj¡}t  dt¡}d}d}d| d d< d| d d	< d
| d d< d| d d< d| d d	< d| d d< d|d d< d|d d	< d|d d< d|d d< d|d d	< d|d d< d|d< d|d	< d|d< t ||| ||¡}| t j¡}t|j	ddddt
ƒ t|jd ddddt
ƒ t|jd	 d ddd!t
ƒ t|jd d"ddd#t
ƒ t|jd d$ddd%t
ƒ t|jd	 d&ddd't
ƒ t|jd d(ddd)t
ƒ t|jd*ddd+t
ƒ t|jd d,d-dd.t
ƒ t|jd	 d/d-dd0t
ƒ t|jd d1d-dd2t
ƒ t|jd3ddd4t
ƒ t|jd d d	ddd5t
ƒ t|jd	 d dddd6t
ƒ t|jd d dddd7t
ƒ t|jd d	 dddd8t
ƒ t|jd	 d	 d	ddd9t
ƒ t|jd d	 dddd:t
ƒ t|jd d dddd;t
ƒ t|jd	 d dddd<t
ƒ t|jd d d	ddd=t
ƒ d S )>Nr
   r,   rŸ   r    gq=
ø<Ár;   r   g…ëQ¸_0Ar   g
×£Ð2âVÁr   g$Ø¸PBSÀr2   g~8Hˆò©`ÀgIœK±'á¸?r¡   r¢   r£   r¤   r¥   r¦   r§   r¨   r©   r³   rU   ÚapcsrW   gJ‹8�,ï¿r-   rX   g¨À«ÍÛË¿rY   g_ÑÚ�€·¿rZ   g¡o(ï¿r\   g´ÛrŒÈÊ¿r^   gIÞ�YÅ;·¿r`   g¥ôâ{�þï?rb   geœ»1ûËõ>rd   re   gæ¾&×y¿rg   g>p³+IA¿ri   rµ   r3   rl   rm   rn   ro   rp   rq   rr   rs   rt   )r/   r0   r"   ru   r1   rã   rv   rw   r   rW   r#   rx   ry   rb   r2   r3   rz   )Úpvr{   r|   r}   r~   r   r
   r
   r   Ú	test_apcsd  sV    rå   c                  C   sŒ  t  dtj¡} d}d}d| d d< d| d d< d	| d d
< d| d d< d| d d< d| d d
< t ||| ¡}| t j¡}t|jddddt	ƒ t|j
d ddddt	ƒ t|j
d ddddt	ƒ t|j
d
 ddddt	ƒ t|jd ddddt	ƒ t|jd ddddt	ƒ t|jd
 ddddt	ƒ t|jd ddd!t	ƒ t|jd d"d#dd$t	ƒ t|jd d%d#dd&t	ƒ t|jd
 d'd#dd(t	ƒ t|jd)ddd*t	ƒ t|jd d dddd+t	ƒ t|jd d dddd,t	ƒ t|jd
 d dddd-t	ƒ t|jd d dddd.t	ƒ t|jd d dddd/t	ƒ t|jd
 d dddd0t	ƒ t|jd d
 dddd1t	ƒ t|jd d
 dddd2t	ƒ t|jd
 d
 dddd3t	ƒ d S )4Nr
   rI   rJ   g¤p=J6ÏWÁr;   r   g%�•Ar   g¤p=
@3Ár   g¥ž¡¼C=Àr2   g pvkYp|ÀgQÉ›mþ@›?rT   rU   Úapcs13rW   gÝ„h2×ì?r-   rX   gÇ['i®¶Ú¿rY   gL¦µŸ³*Ç¿rZ   g²îà™›ì?r\   g#bmúŒOÚ¿r^   g¸£Š\ÐÆ¿r`   g#Yô¶*ð?rb   gXý¤Zp?rd   re   g¦Â†Îà?rg   gä§»Ê1q?ri   gmRyýÿÿï?r3   rl   rm   rn   ro   rp   rq   rr   rs   rt   )r/   r0   r"   ru   ræ   rv   rw   r   rW   r#   rx   ry   rb   r2   r3   rz   )rä   r}   r~   r   r
   r
   r   Útest_apcs13’  s@    rç   c                  C   sP   t  dtj¡ t j¡} d| _d}t || ¡} |  t j¡} t| j	ddddt
ƒ d S )Nr
   çX9´Èv¾ó?gƒÀÊ¡E¶@gÙÎ÷Sã¥@r-   ÚaperrW   )r/   Úzerosr"   Údt_eraASTROMrv   rw   r·   ré   r   rÃ   r#   )r   rÈ   r
   r
   r   Ú	test_aperµ  s    rì   c                  C   sV   t  dtj¡ t j¡} d| _d}d}t ||| ¡} |  t j¡} t| j	ddddt
ƒ d S )	Nr
   rè   rI   rJ   g°vA§‡
@r-   Úaper13rW   )r/   rê   r"   rë   rv   rw   r·   rí   r   rÃ   r#   )r   Zut11Zut12r
   r
   r   Útest_aper13¾  s    rî   c            
      C   s  d} d}d}d}d}d}d}d}d	}t  | ||||||||¡	}	|	 tj¡}	t|	jd
dddtƒ t|	jddddtƒ t|	j	ddddtƒ t|	j
ddddtƒ t|	jddddtƒ t|	jddddtƒ t|	jddddtƒ t|	jddddtƒ t|	jd	dddtƒ d S )Nr°   rª   r«   r¬   r­   r®   r¯   r±   r²   r¶   r-   Úapior·   r¸   r¹   rº   r»   r¼   r½   r¾   r¿   rÀ   çTh÷§;¡>rÁ   rÂ   rÃ   rÄ   rÅ   rÆ   rÇ   )r"   rï   rv   r/   rw   r   r·   r#   rº   r¼   r¾   rÀ   rÁ   rÃ   rÅ   rÇ   )
rÎ   rÈ   rÉ   rÊ   rË   rÌ   rÍ   rÅ   rÇ   r   r
   r
   r   Ú	test_apioÈ  s(    rñ   c                  C   s(  d} d}d}d}d}d}d}d}d	}d
}	d}
d}t  | |||||||||	|
|¡\}}| tj¡}t|jddddtƒ t|jddddtƒ t|j	ddddtƒ t|j
ddddtƒ t|jddddtƒ t|jddddtƒ t|jddddtƒ t|jdddd tƒ t|jd!d"dd#tƒ t|d$dd%tƒ d S )&NrŸ   rÐ   rÑ   r«   r¬   r­   r®   r¯   rÒ   rÓ   rÔ   rÕ   r¶   r-   Úapio13r·   r¸   r¹   rº   r»   r¼   r½   r¾   r¿   rÀ   rð   rÁ   r×   rÃ   rØ   rÄ   rÅ   rÙ   rÆ   rÇ   r   rA   )r"   rò   rv   r/   rw   r   r·   r#   rº   r¼   r¾   rÀ   rÁ   rÃ   rÅ   rÇ   r   )rÛ   rÜ   rÝ   rÉ   rÊ   rË   rÌ   rÍ   rÞ   rß   rà   rá   r   rA   r
   r
   r   Útest_apio13ß  s0    $ró   c            
   
   C   sd   d} d}d}d}d}d}d}d}t  | |||||||¡\}}	t|d	d
ddtƒ t|	dd
ddtƒ d S )Nç®Gáz®@çƒÀÊ¡EÆ?çñhãˆµøä>çñhãˆµøÔ>çš™™™™™¹?ç     €K@rI   rJ   ç]áÎV®@r-   Úatcc13Úraça,g´GÆ?Úda)r"   rû   r   r#   )
ÚrcÚdcÚprÚpdÚpxÚrvr}   r~   rü   rþ   r
   r
   r   Útest_atcc13ú  s    r  c               	   C   sr   d} d}t  | |¡\}}d}d}d}d}d}d}	t  ||||||	|¡\}
}t|
d	d
ddtƒ t|dd
ddtƒ d S )NrI   rJ   rô   rõ   rö   r÷   rø   rù   rú   r-   Úatccqrü   rý   rþ   )r"   r›   r  r   r#   )r}   r~   r   r�   rÿ   r   r  r  r  r  rü   rþ   r
   r
   r   Ú
test_atccq  s    r  c               
   C   sx   d} d}d}d}d}d}d}d}t  | |||||||¡\}}	}
t|d	d
ddtƒ t|	dd
ddtƒ t|
ddddtƒ d S )Nrô   rõ   rö   r÷   rø   rù   rI   rJ   çS­ 8T®@r-   Úatci13ÚriçÁ	>ÙÍ"Æ?Údirœ   r6   r�   )r"   r	  r   r#   )rÿ   r   r  r  r  r  r}   r~   r
  r  r�   r
   r
   r   Útest_atci13  s    r  c               	   C   sr   d} d}t  | |¡\}}d}d}d}d}d}d}	t  ||||||	|¡\}
}t|
d	d
ddtƒ t|dd
ddtƒ d S )NrI   rJ   rô   rõ   rö   r÷   rø   rù   r  r-   Úatciqr
  r  r  )r"   r›   r  r   r#   )r}   r~   r   r�   rÿ   r   r  r  r  r  r
  r  r
   r
   r   Ú
test_atciq&  s    r  c               
   C   s
  t  dtj¡ t j¡} d}d}t ||¡\}}d}d}d}d}d}	d	}
d
| d _d| d _d| d j	d d< d| d j	d d< d| d j	d d< d| d j	d d< d| d j	d d< d| d j	d d< d| d _d| d _d| d j	d d< d| d j	d d< d| d j	d d< d| d j	d d< d| d j	d d< d| d j	d d< d| d _d | d _d!| d j	d d< d"| d j	d d< d#| d j	d d< d$| d j	d d< d%| d j	d d< d&| d j	d d< t 
|||||	|
|| ¡\}}t|d'd(d)d*tƒ t|d+d(d)d,tƒ d S )-Nr    rI   rJ   rô   rõ   rö   r÷   rø   rù   ç‡M…ôë¹2?r   çLÎaã§�ô=çfÄ¢u–=Àr;   çKµ&H2pÀr   çyj`qY±ÿ¿r   ç1bûÞ!+i?r2   çÆ2"*¨«p¿ç›dÇÆµ]¿çCèû«EO?çßA:ÜÅ)>çç.Á]�žç?çqØ	qÝ‹@ç¢9úN‚ÿ?ç'ká™Põ~¿çˆ¨¶%ËT?çŒíMóåÚG?rk   çTäqs*Ù>ç BÔ&VG¿ç¯‹Ÿ÷záb¿ç´À�'WXQ¿çËÎO•\dÚ>ç@>„ž4–¾ç¢=7õ™å“¾g”Ø‡rT®@r-   Úatciqnr
  g8§2OÎ"Æ?r  )r/   r0   r"   Údt_eraLDBODYrv   rw   r›   ÚbmÚdlrä   r'  r   r#   )Úbr}   r~   r   r�   rÿ   r   r  r  r  r  r
  r  r
   r
   r   Útest_atciqn5  sJ    





r,  c                  C   sZ   d} d}t  | |¡\}}d}d}t  |||¡\}}t|ddddtƒ t|d	ddd
tƒ d S )NrI   rJ   rô   rõ   gÀDÎ®@r-   Úatciqzr
  gÑ£þÒ¼ Æ?r  )r"   r›   r-  r   r#   )r}   r~   r   r�   rÿ   r   r
  r  r
   r
   r   Útest_atciqz]  s    r.  c                  C   s  d} d}d}d}d}d}d}d}d	}d
}	d}
d}d}d}d}d}d}d}t  | |||||||||	|
|||||||¡\}}}}}}}t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|dd dd!tƒ t|d"dd#tƒ d S )$Nrô   rõ   rö   r÷   rø   rù   rŸ   rÐ   rÑ   r«   r¬   r­   r®   r¯   rÒ   rÓ   rÔ   rÕ   gÅ†s0>¯·?r-   Úatco13Úaobgg°6÷Ç…ö?Úzobg›Âø¸·¿Úhobgdq¿º
ùÅ?ÚdobgG‘[�®@ÚrobrÚ   r6   r�   r   rA   )r"   r/  r   r#   r   )rÿ   r   r  r  r  r  rÛ   rÜ   rÝ   rÉ   rÊ   rË   rÌ   rÍ   rÞ   rß   rà   rá   r0  r1  r2  r3  r4  r�   rA   r
   r
   r   Útest_atco13h  s4    :r5  c                  C   s`   d} d}d}d}t  | |||¡\}}}t|ddddtƒ t|d	ddd
tƒ t|ddddtƒ d S )NçU° 8T®@çV÷=ÙÍ"Æ?rI   rJ   çdäÎV®@r-   Úatic13rÿ   çûN,g´GÆ?r   rœ   r6   r�   )r"   r9  r   r#   )r
  r  r}   r~   rÿ   r   r�   r
   r
   r   Útest_atic13…  s    r;  c                  C   sZ   d} d}t  | |¡\}}d}d}t  |||¡\}}t|ddddtƒ t|d	ddd
tƒ d S )NrI   rJ   r6  r7  r8  r-   Úaticqrÿ   r:  r   )r"   r›   r<  r   r#   )r}   r~   r   r�   r
  r  rÿ   r   r
   r
   r   Ú
test_aticq�  s    r=  c            	      C   sò  t  dtj¡ t j¡} d}d}t ||¡\}}d}d}d| d _d| d _d	| d j	d
 d< d| d j	d
 d< d| d j	d
 d< d| d j	d d< d| d j	d d< d| d j	d d< d| d _d| d _d| d j	d
 d< d| d j	d
 d< d| d j	d
 d< d| d j	d d< d| d j	d d< d| d j	d d< d| d _d| d _d| d j	d
 d< d| d j	d
 d< d| d j	d
 d< d | d j	d d< d!| d j	d d< d"| d j	d d< t 
|||| ¡\}}t|d#d$d%d&tƒ t|d'd$d%d(tƒ d S ))Nr    rI   rJ   gÃDÎ®@gg‘þÒ¼ Æ?r  r   r  r  r;   r  r   r  r   r  r2   r  r  r  r  r  r  r  r  r  r  rk   r   r!  r"  r#  r$  r%  r&  g¨†×A®@r-   Z	eraAtciqnrÿ   gdOò€¡EÆ?r   )r/   r0   r"   r(  rv   rw   r›   r)  r*  rä   Zaticqnr   r#   )	r+  r}   r~   r   r�   r
  r  rÿ   r   r
   r
   r   Útest_aticqn›  sB    





r>  c                  C   sÖ   d} d}d}d}d}d}d}d}d	}d
}	d}
d}d}d}t  | |||||||||	|
|||¡\}}}}}}t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|dddtƒ d S )Nr6  r7  rŸ   rÐ   rÑ   r«   r¬   r­   r®   r¯   rÒ   rÓ   rÔ   rÕ   ç¦ 0�£·?r-   Úatio13r0  çl¸Èý-†ö?r1  çÓ)…¬·¿r2  ç9+Ÿ"hüÅ?r3  çJáA®@r4  r   rA   )r"   r@  r   r#   r   )r
  r  rÛ   rÜ   rÝ   rÉ   rÊ   rË   rÌ   rÍ   rÞ   rß   rà   rá   r0  r1  r2  r3  r4  rA   r
   r
   r   Útest_atio13¿  s*    0rE  c                  C   sÒ   d} d}d}d}d}d}d}d}d	}d
}	d}
d}t  | |||||||||	|
|¡\}}d}d}t  |||¡\}}}}}t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ d S )NrŸ   rÐ   rÑ   r«   r¬   r­   r®   r¯   rÒ   rÓ   rÔ   rÕ   r6  r7  r?  r-   Úatioqr0  rA  r1  rB  r2  rC  r3  rD  r4  )r"   rò   rF  r   r#   )rÛ   rÜ   rÝ   rÉ   rÊ   rË   rÌ   rÍ   rÞ   rß   rà   rá   r   rA   r
  r  r0  r1  r2  r3  r4  r
   r
   r   Ú
test_atioq×  s*    $rG  c                  C   sl  d} d}d}d}d}d}d}d}d	}d
}	d}
d}d}d}t  d||| |||||||||	|
|¡\}}}t|ddddtƒ t|ddddtƒ t|dddtƒ d}d}t  d||| |||||||||	|
|¡\}}}t|ddddtƒ t|ddddtƒ t|dddtƒ d}d}t  d ||| |||||||||	|
|¡\}}}t|dddd!tƒ t|dddd"tƒ t|ddd#tƒ d S )$NrŸ   rÐ   rÑ   r«   r¬   r­   r®   r¯   rÒ   rÓ   rÔ   rÕ   çDqáA®@çÌGŸ"hüÅ?ÚRg¼?Íý­@r-   Úatoc13zR/rcgj»‹ë0KÆ?zR/dcr   zR/jçeº~)…¬·¿ÚHgþ?Íý­@zH/rczH/dczH/jçÍæ.�£·?çu»Èý-†ö?ÚAzA/rczA/dczA/j)r"   rK  r   r#   r   )rÛ   rÜ   rÝ   rÉ   rÊ   rË   rÌ   rÍ   rÞ   rß   rà   rá   Úob1Úob2rÿ   r   rA   r
   r
   r   Útest_atoc13ï  s<    ,,,rS  c                  C   sl  d} d}d}d}d}d}d}d}d	}d
}	d}
d}d}d}t  d||| |||||||||	|
|¡\}}}t|ddddtƒ t|ddddtƒ t|dddtƒ d}d}t  d||| |||||||||	|
|¡\}}}t|ddddtƒ t|ddddtƒ t|dddtƒ d}d}t  d ||| |||||||||	|
|¡\}}}t|dddd!tƒ t|dddd"tƒ t|ddd#tƒ d S )$NrŸ   rÐ   rÑ   r«   r¬   r­   r®   r¯   rÒ   rÓ   rÔ   rÕ   rH  rI  rJ  ç_}S8T®@r-   Úatoi13úR/riç—Ú”>Î"Æ?úR/dir   zR/JrL  rM  ç¢‚S8T®@úH/riúH/dizH/JrN  rO  rP  úA/riúA/dizA/J)r"   rU  r   r#   r   )rÛ   rÜ   rÝ   rÉ   rÊ   rË   rÌ   rÍ   rÞ   rß   rà   rá   rQ  rR  r
  r  rA   r
   r
   r   Útest_atoi13  s<    ,,,r^  c                  C   s  d} d}d}d}d}d}d}d}d	}d
}	d}
d}t  | |||||||||	|
|¡\}}d}d}t  d|||¡\}}t|ddddtƒ t|ddddtƒ d}d}t  d|||¡\}}t|ddddtƒ t|ddddtƒ d}d}t  d|||¡\}}t|ddddtƒ t|ddddtƒ d S ) NrŸ   rÐ   rÑ   r«   r¬   r­   r®   r¯   rÒ   rÓ   rÔ   rÕ   rH  rI  rJ  rT  r-   ÚatoiqrV  rW  rX  rL  rM  rY  rZ  r[  rN  rO  rP  r\  r]  )r"   rò   r_  r   r#   )rÛ   rÜ   rÝ   rÉ   rÊ   rË   rÌ   rÍ   rÞ   rß   rà   rá   r   rA   rQ  rR  r
  r  r
   r
   r   Ú
test_atoiq1  s8    $r`  z1.24z+numpy < 1.24 do not support no-input ufuncs)Úreasonc                  C   sH   t  ¡ \} }}t| ddddtƒ t|ddddtƒ t|ddddtƒ d S )	NgGMÅè.‹¾r-   Úbi00Údpsibig]ã9íÉ¿a¾Údepsbigx/$( s¾Údra)r"   rb  r   r#   )rc  rd  re  r
   r
   r   Ú	test_bi00P  s    rf  c                  C   sø  t  dt¡} t  dt¡}t  dt¡}t dd¡\} }}t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ t|d d ddddtƒ t|d d	 dddd tƒ t|d d d!ddd"tƒ t|d	 d d#ddd$tƒ t|d	 d	 d%ddd&tƒ t|d	 d d'ddd(tƒ t|d d d)ddd*tƒ t|d d	 d+ddd,tƒ t|d d d-ddd.tƒ t|d d d/ddd0tƒ t|d d	 d1ddd2tƒ t|d d d3ddd4tƒ t|d	 d d5ddd6tƒ t|d	 d	 d7ddd8tƒ t|d	 d d9ddd:tƒ t|d d d;ddd<tƒ t|d d	 d=ddd>tƒ t|d d d?ddd@tƒ d S )AN©r    r    ç   @€OBAçéH.ÿyè@r   çÌÿÿÿÿÿï?r-   Úbp00Úrb11r   çb/$( s¾rd   Úrb12r   çâOiµ0 u>Úrb13ç÷c0( s>Úrb21çåÿÿÿÿÿï?Úrb22çLã9íÉ¿a>Úrb23çøEó¿0 u¾Úrb31çÝ¹É¿a¾Úrb32çÞÿÿÿÿÿï?Úrb33giR«ÿÿï?Úrp11géú÷9Ô~L?r6   Úrp12g»ÿ‘ýÃ8?Úrp13gÊïú8Ô~L¿Úrp21gx¬ 5ÿÿï?Úrp22gÏñ4!±†¾Úrp23g¢œ¢•ýÃ8¿Úrp31gÐnÿ+j…¾Úrp32gñ¥ªÙÿÿï?Úrp33çÈ˜¯ÿÿï?Úrbp11çaÝ1<~L?Úrbp12çÞw³³WÅ8?Úrbp13ç£¼7<~L¿Úrbp21çG"	5ÿÿï?Úrbp22çI²â—C‚¾Úrbp23çù*°˜WÅ8¿Úrbp31çaÌK?2Ù‰¾Úrbp32çyv¦Ùÿÿï?Úrbp33)r/   r0   r1   r"   rk  r   r#   ©ÚrbÚrpÚrbpr
   r
   r   Ú	test_bp00X  s>    rœ  c                  C   sø  t  dt¡} t  dt¡}t  dt¡}t dd¡\} }}t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ t|d d ddddtƒ t|d d	 ddddtƒ t|d d d ddd!tƒ t|d	 d d"ddd#tƒ t|d	 d	 d$ddd%tƒ t|d	 d d&ddd'tƒ t|d d d(ddd)tƒ t|d d	 d*ddd+tƒ t|d d d,ddd-tƒ t|d d d.ddd/tƒ t|d d	 d0ddd1tƒ t|d d d2ddd3tƒ t|d	 d d4ddd5tƒ t|d	 d	 d6ddd7tƒ t|d	 d d8ddd9tƒ t|d d d:ddd;tƒ t|d d	 d<ddd=tƒ t|d d d>ddd?tƒ d S )@Nrg  rh  ri  r   rj  r-   Úbp06rl  r   çu&Ö7 s¾r6   rn  r   ç‰­&0 u>rp  ç3Ö'Ê7 s>rr  rs  rt  ç¢@.Œ§¿a>rv  çéõ610 u¾rx  çZ7ÑX§¿a¾rz  r{  r|  gžU«ÿÿï?r}  gTØôÓ~L?r~  gÇrÛÅýÃ8?r  gt$ÖóÓ~L¿r€  gQ° 5ÿÿï?r�  gêñÍ´s´†¾r‚  g§ºÊýÃ8¿rƒ  gQ2‘Øf…¾r„  gM¥ªÙÿÿï?r…  çœ¯ÿÿï?r‡  çä?ì;~L?r‰  çŸÁoèWÅ8?r‹  çL¬ò;~L¿r�  ç'&	5ÿÿï?r�  çºËë)òF‚¾r‘  ç¾\„ÍWÅ8¿r“  çF+w×Õ‰¾r•  çÖu¦Ùÿÿï?r—  )r/   r0   r1   r"   r�  r   r#   r˜  r
   r
   r   Ú	test_bp06z  s>    r­  c                  C   s®   t  dt¡} d| d d< d| d d< d| d d< d| d d< d	| d d< d
| d d< d| d d< d| d d< d| d d< t | ¡\}}t|ddddtƒ t|ddddtƒ d S )Nrg  ç¸¥%øÿï?r   çÇ•çO�d¿r   ç .±�ÃêQ¿r   çN¬¢rg�d?çs=­[ùÿï?çqìóù ?ç¯sÃ¦SêQ?çüGCEPr¿çMjŒ¾þÿï?r-   Úbpn2xyr˜   r™   )r/   r0   r1   r"   r·  r   r#   )Úrbpnr˜   r™   r
   r
   r   Útest_bpn2xyœ  s    r¹  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  ç     =ê@r   g‘ÿ÷¥ÿÿï?r-   Úc2i00aÚ11r   gikÐÿ�->Ú12r   çrl‡úB¿Ú13g4ÍÐ3Ò™Y¾Ú21gæé�ÿÿÿï?Ú22g]JèY¿Ú23çFcBúB?Ú31ç+Áp&Z?Ú32çxéˆ¥ÿÿï?Ú33)r/   r0   r1   r"   r»  r   r#   ©Úrc2ir
   r
   r   Útest_c2i00a¬  s    rË  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  rº  r   gåø¥ÿÿï?r-   Úc2i00br¼  r   g{Òÿ�->r½  r   gÔãÝžúB¿r¿  gò7O>¿™Y¾rÀ  g~ê�ÿÿÿï?rÁ  gt¬{K¿rÂ  ç	Û˜žúB?rÄ  çðIÒ¹K?rÆ  çe÷ˆ¥ÿÿï?rÈ  )r/   r0   r1   r"   rÌ  r   r#   rÉ  r
   r
   r   Útest_c2i00bº  s    rÐ  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  rº  r   g‹ÿ÷¥ÿÿï?r-   Úc2i06ar¼  r   g¼úÌ4->r½  r   g¬ŒúB¿r¿  g<‹kÉ™Y¾rÀ  gê�ÿÿÿï?rÁ  gÞJn�T¿rÂ  ç±ƒØúB?rÄ  ça’”ÍT?rÆ  ç«éˆ¥ÿÿï?rÈ  )r/   r0   r1   r"   rÑ  r   r#   rÉ  r
   r
   r   Útest_c2i06aÈ  s    rÕ  c                  C   s€  t  dt¡} t  dt¡}d| d d< d| d d< d| d d< d| d d< d	| d d< d
| d d< d| d d< d| d d< d| d d< t dd| ¡}t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d dddd tƒ d S )!Nrg  r®  r   r¯  r   r°  r   r±  r²  r³  r´  rµ  r¶  rh  ri  g}`
¿þÿï?r-   Úc2ibpnr¼  gçàÖ9ž0¾r½  g¡§SêQ¿r¿  gÍuò ù5k>rÀ  gÊ	‚ÿÿÿï?rÁ  geÝè>Or?rÂ  rÄ  rÆ  rÈ  )r/   r0   r1   r"   rÖ  r   r#   )r¸  rÊ  r
   r
   r   Útest_c2ibpnÖ  s*    r×  c                  C   s  t  dt¡} d}d}t dd||¡} t| d d ddd	d
tƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ d S )Nrg  çmýúB?çoB¶´T?rh  éèÑ  r   ç‰ÿ÷¥ÿÿï?r-   Úc2ixyr¼  r   géÑÿ�->r½  r   gvTúB¿r¿  g{´Ò¿Ë™Y¾rÀ  çê�ÿÿÿï?rÁ  g¬í�vT¿rÂ  rÄ  rÆ  çªéˆ¥ÿÿï?rÈ  )r/   r0   r1   r"   rÜ  r   r#   )rÊ  r˜   r™   r
   r
   r   Ú
test_c2ixyî  s    rß  c                  C   s  t  dt¡} d}d}d}t |||¡} t| d d dddd	tƒ t| d d
 ddddtƒ t| d d ddddtƒ t| d
 d ddddtƒ t| d
 d
 ddddtƒ t| d
 d ddddtƒ t| d d ddddtƒ t| d d
 ddddtƒ t| d d ddddtƒ d S )Nrg  rØ  rÙ  ç�2ØØ>3J¾r   rÛ  r-   Úc2ixysr¼  r   guå³?÷->r½  r   g‡TúB¿r¿  g†ÂÑùÎ™Y¾rÀ  rÝ  rÁ  g]Þ�vT¿rÂ  rÄ  rÆ  rÞ  rÈ  )r/   r0   r1   r"   rá  r   r#   )rÊ  r˜   r™   r   r
   r
   r   Útest_c2ixysþ  s    râ  c                  C   sZ   t  dt¡} d| d< d| d< d| d< t | ¡\}}t|dd	d
dtƒ t|dd	d
dtƒ d S )Nr,   ç      Y@r   ç      IÀr   ç      9@r   çO»ag¬Ý¿r6   Úc2srÈ   ç¢�6è�(Ì?rÊ   )r/   r0   r1   r"   rç  r   r#   )r;   rÈ   rÊ   r
   r
   r   Útest_c2s  s    ré  c                  C   s&  t  dt¡} d}d}d}d}d}d}t ||||||¡} t| d d ddd	d
tƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ d S )Nrg  rh  rº  çË4Êæ‘>ç¹Ÿ.6¿>r   gÔ‚L¨,Ç¿r-   Úc2t00ar¼  r   gQ]Ž·¤xï?r½  r   gíqrX/?r¿  gÚŠ¿]¤xï¿rÀ  gøz¸,Ç¿rÁ  gÜÕ:cH×B?rÂ  g6¨Ž$ëB?rÄ  g g�>{Å?rÆ  çÙÞ¦ÿÿï?rÈ  )r/   r0   r1   r"   rì  r   r#   ©Úrc2tÚttaÚutaÚttbÚutbrÌ   rÍ   r
   r
   r   Útest_c2t00a  s"    rô  c                  C   s&  t  dt¡} d}d}d}d}d}d}t ||||||¡} t| d d ddd	d
tƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ d S )Nrg  rh  rº  rê  rë  r   gŒËS¨,Ç¿r-   Úc2t00br¼  r   gÐŽ·¤xï?r½  r   g*hd�]/?r¿  g�B¿]¤xï¿rÀ  g ä�¸,Ç¿rÁ  g9¹º×F×B?rÂ  g< å²"ëB?rÄ  g¶áñÑlÅ?rÆ  g¹ì¦ÿÿï?rÈ  )r/   r0   r1   r"   rõ  r   r#   rî  r
   r
   r   Útest_c2t00b-  s"    rö  c                  C   s&  t  dt¡} d}d}d}d}d}d}t ||||||¡} t| d d ddd	d
tƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ d S )Nrg  rh  rº  rê  rë  r   g½pL¨,Ç¿r-   Úc2t06ar¼  r   gÊ]Ž·¤xï?r½  r   gŽki[/?r¿  g©‹¿]¤xï¿rÀ  g|zz¸,Ç¿rÁ  g‰QTH×B?rÂ  g¢È$$ëB?rÄ  g@8´åuÅ?rÆ  gß¦ÿÿï?rÈ  )r/   r0   r1   r"   r÷  r   r#   rî  r
   r
   r   Útest_c2t06aA  s"    rø  c                  C   sü  t  dt¡} t  dt¡}t  dt¡}d| d d< d| d d< d| d d< d| d d< d	| d d< d
| d d< d| d d< d| d d< d| d d< d}d|d d< d|d d< d|d d< d|d d< d|d d< d|d d< d|d d< d|d d< d|d d< t | ||¡}t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d d ddd!tƒ t|d d d"ddd#tƒ t|d d d$ddd%tƒ t|d d d&ddd'tƒ t|d d d(ddd)tƒ t|d d d*ddd+tƒ d S ),Nrg  g�ÿ÷¥ÿÿï?r   g  \û�->r   r¾  r   g ØÀ3Ò™Y¾gåé�ÿÿÿï?g}kJèY¿rÃ  çŒÏp&Z?rÇ  gðâ-âšü?çÛþÿÿÿÿï?ç+Nå7‚®½ç—4Êæ‘>çE�g¡¯=çÃÿÿÿÿï?çó¸Ÿ.6¿¾çL ,Êæ‘¾ç‹4¹Ÿ.6¿>çøÁÿÿÿÿï?gÐ�L¨,Ç¿r-   Úc2tcior¼  g\]Ž·¤xï?r½  g¤þprX/?r¿  gåŠ¿]¤xï¿rÀ  gô~z¸,Ç¿rÁ  g
Ö:cH×B?rÂ  g7¨Ž$ëB?rÄ  gcu�>{Å?rÆ  rí  rÈ  )r/   r0   r1   r"   r  r   r#   )rÊ  Úrpomrï  Zerar
   r
   r   Útest_c2tcioU  s@    r  c                  C   sü  t  dt¡} t  dt¡}t  dt¡}d| d d< d| d d< d| d d< d| d d< d	| d d< d
| d d< d| d d< d| d d< d| d d< d}d|d d< d|d d< d|d d< d|d d< d|d d< d|d d< d|d d< d|d d< d|d d< t | ||¡}t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d d ddd!tƒ t|d d d"ddd#tƒ t|d d d$ddd%tƒ t|d d d&ddd'tƒ t|d d d(ddd)tƒ t|d d d*ddd+tƒ d S ),Nrg  ç7Éýÿï?r   çÉÅ»ä‚ÖU¿r   ç±_iÖ¥ùB¿r   çý�2§iÖU?ç^£«"þÿï?çmkó{¿rÃ  rù  rÇ  g@+*ƒü?rú  rû  rü  rý  rþ  rÿ  r   r  r  g5°X¨,Ç¿r-   Úc2teqxr¼  gÜÍ�·¤xï?r½  gy™rX/?r¿  gdû¾]¤xï¿rÀ  gY­†¸,Ç¿rÁ  g2¡:cH×B?rÂ  g±¨Ž$ëB?rÄ  g9^�>{Å?rÆ  rí  rÈ  )r/   r0   r1   r"   r  r   r#   )r¸  r  rï  Úgstr
   r
   r   Útest_c2teqxx  s@    r  c            	   
   C   s2  t  dt¡} d}d}d}d}d}d}d}d}t ||||||||¡} t| d d d	d
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ d S )Nrg  rh  rº  çð5–Y.Ú?ç÷”Ž2ä¾rê  rë  r   gf%Ý_7Ç¿r-   Úc2tper¼  r   g1×.`	àì?r½  r   g ÉL'üÙ¿r¿  gô^M="xï¿rÀ  g7ËwP?Å¿rÁ  gøóÕà†¶²?rÂ  gX×3‘¤\Q?rÄ  g@b¸ótÙ?rÆ  g+‘Â\í?rÈ  )r/   r0   r1   r"   r  r   r#   )	rï  rð  rñ  rò  ró  ÚdepsÚdpsirÌ   rÍ   r
   r
   r   Ú
test_c2tpe›  s&    r  c            	   
   C   s2  t  dt¡} d}d}d}d}d}d}d}d}t ||||||||¡} t| d d d	d
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ d S )Nrg  rh  rº  rØ  rÙ  rê  rë  r   gvL¨,Ç¿r-   Úc2txyr¼  r   g‰]Ž·¤xï?r½  r   gØæò [/?r¿  gh‹¿]¤xï¿rÀ  g€z¸,Ç¿rÁ  g?TH×B?rÂ  g]B[$ëB?rÄ  gDèÕÌuÅ?rÆ  g
ß¦ÿÿï?rÈ  )r/   r0   r1   r"   r  r   r#   )	rï  rð  rñ  rò  ró  r˜   r™   rÌ   rÍ   r
   r
   r   Ú
test_c2txy±  s&    r  c                  C   sL   t  ddd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
ddtƒ d S )NéÓ  é   r   rh  r   Úcal2jdÚdjm0g    àÆé@Údjmr   rA   )r"   r  r   r#   r   )r  r  rA   r
   r
   r   Útest_cal2jdÇ  s    r  c                  C   s€   t  dt¡} t  dt¡}d| d< d| d< d| d< t | ¡}t|d ddd	d
tƒ t|d ddd	dtƒ t|d ddd	dtƒ d S )Nr,   ç333333Ó?r   ç333333ó?r   ç      Àr   r   Úcpr   r   r!   )r/   r0   r1   r"   r   r   r#   )r;   Úcr
   r
   r   Útest_cpÎ  s    
r"  c                  C   s  t  dtj¡} t  dtj¡}d| d d< d| d d< d| d d< d	| d
 d< d| d
 d< d| d
 d< t | ¡}t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d d	dddtƒ t|d d ddddtƒ t|d d ddddtƒ d S )Nr
   r  r;   r   r  r   r  r   ç      à¿r2   çÍÌÌÌÌÌ@çÍÌÌÌÌÌì?r   ÚcpvÚp1Úp2Úp3Úv1Úv2Úv3)r/   r0   r"   ru   r&  r   r#   )rä   r!  r
   r
   r   Útest_cpvÚ  s    
r-  c                  C   s|  t  dt¡} t  dt¡}d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< t | ¡}t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ d S )Nrg  ç       @r   ç      @r   r   ç      @ç      @r   Úcrr¼  r½  r¿  rÀ  rÁ  rÂ  rÄ  rÆ  rÈ  )r/   r0   r1   r"   r2  r   r#   )Úrr!  r
   r
   r   Útest_crì  s*    
r4  c                  C   s®   t  dddd¡\} }}}}t| dddtƒ t|ddd	tƒ t|d
ddtƒ t|d dddtƒ t|d dddtƒ t|d dddtƒ t|d dddtƒ t|dddtƒ d S )NÚUTCé   rh  gJëÿ¿/è@éÊ  Úd2dtfr™   r  Úmor   r9   r   é   r8   r   é;   Úmr   é<   r   r    i5  ÚfrA   )r"   r8  r   r#   )ÚiyÚimÚidr*   rA   r
   r
   r   Ú
test_d2dtf  s    rB  c                  C   sx   t  dd¡\} }t| d dddtƒ t|d dddtƒ t|d	 d
ddtƒ t|d dddtƒ t|d dddtƒ d S )Nr   g¸V<Ýšï¿r   r%   Úd2tfr   r:  r   r   é*   r   r   r?   r   r    i  r!   )r"   rC  r   r#   r)   r
   r
   r   Ú	test_d2tf  s    rE  c                  C   s¦   t  dddd¡\} }t| ddddtƒ t|ddd	tƒ t  d
ddd¡\} }t| ddddtƒ t|dddtƒ t  dddd¡\} }t| ddddtƒ t|dddtƒ d S )Nr  r  r   r   g      @@ÚdatZd1r   Új1iØ  é   g     €@@Zd2Új2iá  é	   g     €B@Zd3Új3)r"   rF  r   r#   r   )ZdeltatrA   r
   r
   r   Útest_dat  s    rL  c                  C   s*   t  dddddd¡} t| ddd	d
tƒ d S )Ng   À¯BAç°rh‘í|¿?g®,g~è?g”ö_˜@g¢E¶ó=•µ@g     ì¨@gå—­@úT¿rÄ   ÚdtdbrE   )r"   rN  r   r#   )rN  r
   r
   r   Ú	test_dtdb%  s    rO  c               	   C   sF   t  ddddddd¡\} }}t| | dd	d
dtƒ t|dd
dtƒ d S )Nr5  r7  r  r   r:  r;  g¡JÍhN@g¬ÿ??°BAç�íµ ÷Æ°>Údtf2dÚur   rA   )r"   rQ  r   r#   r   ©Úu1Úu2rA   r
   r
   r   Ú
test_dtf2d*  s    rV  c                  C   sL   d} d}d}d}t  | |||¡\}}t|ddddtƒ t|d	ddd
tƒ d S )NrI   rJ   çffffff@gÍÌÌÌÌÌì¿gãi¶DC"@r6   Úeceq06ÚdrgT† ×ñó¿Údd)r"   rX  r   r#   )r}   r~   r*  ÚdbrY  rZ  r
   r
   r   Útest_eceq060  s    r\  c                  C   s  t  dt¡} d}d}t ||¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rI   rJ   r   g/hûöÿï?r6   Úecm06Úrm11r   g“6ƒ ù,g¿Úrm12r   guíq„|#T¿Úrm13gG™²~ÃDi?Úrm21gìv¨H\í?Úrm22gaT8·tÙ?Úrm23gÒ¤‰v
àÄ>Úrm31gÓÆœ/¿tÙ¿Úrm32gt•On\í?Úrm33)r/   r0   r1   r"   r]  r   r#   )Úrmr}   r~   r
   r
   r   Ú
test_ecm06:  s    rh  c                  C   s.   d} d}t  dd| |¡}t|ddddtƒ d S )	Nr  r  rh  rº  g!ØåEÒ†â¾rÆ   Úee00rE   )r"   ri  r   r#   )Úepsar  Úeer
   r
   r   Ú	test_ee00J  s    rl  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  rº  g¼ò–*Ò†â¾rÆ   Úee00arE   )r"   rm  r   r#   ©rk  r
   r
   r   Ú
test_ee00aQ  s    ro  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  rº  g““^¡‡â¾rÆ   Úee00brE   )r"   rp  r   r#   rn  r
   r
   r   Ú
test_ee00bV  s    rq  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  rº  g=&…†Ò†â¾rÄ   Úee06arE   )r"   rr  r   r#   rn  r
   r
   r   Ú
test_ee06a[  s    rs  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  rº  gaWLuc“!>ç#B’¡œÇ;Úeect00rE   )r"   ru  r   r#   )Zeectr
   r
   r   Útest_eect00`  s    rv  c                  C   s  t  d¡\} }}t|dddtƒ t  tj¡\} }}t|dddtƒ t| ddddtƒ t|d	d
ddtƒ t  tj¡\} }}t|dddtƒ t| ddddtƒ t|dd
ddtƒ t  tj¡\} }}t|dddtƒ t| ddddtƒ t|dd
ddtƒ t  d¡\} }}t|dddtƒ d S )Nr   éÿÿÿÿÚeformÚj0rG  g   @¦TXAç»½×Ùß|Û=Úa1g)áó„Zwk?rÆ   Úf1rI  Úa2g!6‡Zwk?Úf2g   À¥TXAÚa3gv8gÃIwk?Zf3r   rK  )	r"   rx  r   r#   ÚerfaÚWGS84r   ÚGRS80ÚWGS72)r   r>  rA   r
   r
   r   Ú
test_eforme  s     r„  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  rº  gÈšŒƒÖU¿rÄ   Úeo06arE   )r"   r…  r   r#   )r�   r
   r
   r   Ú
test_eo06ax  s    r†  c                  C   sž   t  dt¡} d| d d< d| d d< d| d d< d| d d< d	| d d< d
| d d< d| d d< d| d d< d| d d< d}t | |¡}t|ddddtƒ d S )Nrg  r  r   r  r   r  r   r	  r
  r  rÃ  rù  rÇ  rà  g–ÝïêƒÖU¿r6   ÚeorsrE   )r/   r0   r1   r"   r‡  r   r#   )Úrnpbr   r�   r
   r
   r   Ú	test_eors}  s    r‰  c                  C   s"   t  dd¡} t| ddddtƒ d S )NgœÄ èÕlBAgWÏIïËeÝ@gÕÿaw¬ùž@r-   ÚepbrE   )r"   rŠ  r   r#   )rŠ  r
   r
   r   Útest_epb�  s    r‹  c                  C   s:   d} t  | ¡\}}t|ddddtƒ t|ddddtƒ d S )Ng33333•ž@rh  ç•Ö&è.>Úepb2jdr  g–åÙ †�á@Úmjd)r"   r�  r   r#   )rŠ  r  r  r
   r
   r   Útest_epb2jd’  s    r�  c                  C   s"   t  dd¡} t| ddddtƒ d S )NiYh% g    €à¼ÀgÇF°
ïž@r-   ÚepjrE   )r"   r�  r   r#   )r�  r
   r
   r   Útest_epj™  s    r‘  c                  C   s:   d} t  | ¡\}}t|ddddtƒ t|ddddtƒ d S )Ng333333Ÿ@rh  rŒ  Úepj2jdr  gffffö˜è@rŽ  )r"   r’  r   r#   )r�  r  r  r
   r
   r   Útest_epj2jdž  s    r“  c                  C   sz  t  dtj¡} t  dtj¡}t dd¡\} }}t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t|d d ddddtƒ t|d d	 ddddtƒ t|d d ddddtƒ t|d	 d ddddtƒ t|d	 d	 ddddtƒ t|d	 d d ddd!tƒ t|ddd"tƒ d S )#Nr
   rh  gß%åÌpê@r   gÜ9jãºÒè¿r6   Úepv00zph(x)r   gÆt9¥ìéá?zph(y)r   g‰ŽàÉÏ?zph(z)g«á	f¨\†¿rÄ   zvh(x)g1ejÚŠ‰¿zvh(y)gÍ±Ò\?&v¿zvh(z)g{Ÿôd¯è¿zpb(x)g®6ð8êá?zpb(y)gFÛ¦-�Ï?zpb(z)gËÕÕ�\†¿zvb(x)gÁ‡W'b‡‰¿zvb(y)g’mßÀP#v¿zvb(z)rA   )r/   r0   r"   ru   r”  r   r#   r   )ZpvhZpvbrA   r
   r
   r   Ú
test_epv00¥  s     r•  c                  C   sL   d} d}d}d}t  | |||¡\}}t|ddddtƒ t|d	ddd
tƒ d S )Ng     J“@g   @ �BArè   ç/Ý$�•ï?gÑ3J®ëzõ?r6   Úeqec06r*  g\kÁöâ?r[  )r"   r—  r   r#   )r}   r~   rY  rZ  r*  r[  r
   r
   r   Útest_eqec06¸  s    r˜  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  g    @"ä@gKÿú·?r¹   Úeqeq94rE   )r"   r™  r   r#   )Zeqeqr
   r
   r   Útest_eqeq94Â  s    rš  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  ç    €Žê@g#ˆ“;¿Ù?r-   Úera00rE   )r"   rœ  r   r#   )rœ  r
   r
   r   Ú
test_era00Ç  s    r�  c                   C   s   t t d¡ddddtƒ d S )Nçš™™™™™é?gõ@�¸%ÿ?r-   Úfad03rE   )r   r"   rŸ  r#   r
   r
   r
   r   Ú
test_fad03Ì  s    r   c                   C   s   t t d¡ddddtƒ d S )Nrž  g”7ôXêû?r-   Úfae03rE   )r   r"   r¡  r#   r
   r
   r
   r   Ú
test_fae03Ð  s    r¢  c                   C   s   t t d¡ddddtƒ d S )Nrž  ge� Ð?r-   Úfaf03rE   )r   r"   r£  r#   r
   r
   r
   r   Ú
test_faf03Ô  s    r¤  c                   C   s   t t d¡ddddtƒ d S )Nrž  gr®¤(T@r-   Úfaju03rE   )r   r"   r¥  r#   r
   r
   r
   r   Útest_faju03Ø  s    r¦  c                   C   s   t t d¡ddddtƒ d S )Nrž  gA÷Ÿï‹‡@r-   Úfal03rE   )r   r"   r§  r#   r
   r
   r
   r   Ú
test_fal03Ü  s    r¨  c                   C   s   t t d¡ddddtƒ d S )Nrž  g[\º=è@r-   Úfalp03rE   )r   r"   r©  r#   r
   r
   r
   r   Útest_falp03à  s    rª  c                   C   s   t t d¡ddddtƒ d S )Nrž  gø‰&î<4
@r-   Úfama03rE   )r   r"   r«  r#   r
   r
   r
   r   Útest_fama03ä  s    r¬  c                   C   s   t t d¡ddddtƒ d S )Nrž  g,·s³Z«@r-   Úfame03rE   )r   r"   r­  r#   r
   r
   r
   r   Útest_fame03è  s    r®  c                   C   s   t t d¡ddddtƒ d S )Nrž  g‰zµ¢ @r-   Úfane03rE   )r   r"   r¯  r#   r
   r
   r
   r   Útest_fane03ì  s    r°  c                   C   s   t t d¡ddddtƒ d S )Nrž  gd¬Š;üäÀr-   Úfaom03rE   )r   r"   r±  r#   r
   r
   r
   r   Útest_faom03ð  s    r²  c                   C   s   t t d¡ddddtƒ d S )Nrž  g[š ú“?r-   Úfapa03rE   )r   r"   r³  r#   r
   r
   r
   r   Útest_fapa03ô  s    r´  c                   C   s   t t d¡ddddtƒ d S )Nrž  gÎ²;	~|@r-   Úfasa03rE   )r   r"   rµ  r#   r
   r
   r
   r   Útest_fasa03ø  s    r¶  c                   C   s   t t d¡ddddtƒ d S )Nrž  g.ö÷Âø¸@r-   Úfaur03rE   )r   r"   r·  r#   r
   r
   r
   r   Útest_faur03ü  s    r¸  c                   C   s   t t d¡ddddtƒ d S )Nrž  gôap62f@r-   Úfave03rE   )r   r"   r¹  r#   r
   r
   r
   r   Útest_fave03   s    rº  c                  C   s¨   d} d}d}d}d}d}t  | |||||¡\}}}}	}
}t|ddd	d
tƒ t|ddd	dtƒ t|ddd	dtƒ t|	ddd	dtƒ t|
ddd	dtƒ t|ddd	dtƒ d S )Ng¦¹ r]†³?g?UhûÏ2ò¿g“Gt">²ô>gËÀÙ–½×>gôýÔxé&Á?gffffff!@gšçnã¢k¶?r6   Úfk425Úr2000gÑZ´Ðò¿r-   Úd2000g
„19X|ô>r¹   Údr2000g ù¤ƒk¥×>rÆ   Údd2000ge½YR¦&Á?ç‚vIhÂ%<=Úp2000g¶œHXy!@Úv2000)r"   r»  r   r#   )Úr1950Úd1950Údr1950Údd1950Úp1950Úv1950r¼  r½  r¾  r¿  rÁ  rÂ  r
   r
   r   Ú
test_fk425  s     rÉ  c                  C   sF   d} d}d}t  | ||¡\}}t|ddddtƒ t|dd	dd
tƒ d S )NgÆ¹>Kh�?g©9ª««Î½¿g¯eáµŠž@gòW˜xØ›?rÄ   Úfk45zr¼  gã«¸H�¼¿rÀ  r½  )r"   rÊ  r   r#   )rÃ  rÄ  Úbepochr¼  r½  r
   r
   r   Ú
test_fk45z  s    rÌ  c                  C   s¨   d} d}d}d}d}d}t  | |||||¡\}}}}	}
}t|ddd	d
tƒ t|ddd	dtƒ t|ddd	dtƒ t|	ddd	dtƒ t|
ddd	dtƒ t|ddd	dtƒ d S )Ng�*LKêë?g#a_,ýè¿g3‚4îé,õ>g‚$Ö]sµÍ>g£#¹ü‡ôÃ?gHáz®·U@gÈ±äç¢ë?rÀ  Úfk524rÃ  g?pþ×0)è¿rÄ  gíƒ´"8õ>r¹   rÅ  gø�ÒwgÎ>rÆ   rÆ  gÙfíøÃ?rÇ  g{Âw´ò²U@rU   rÈ  )r"   rÍ  r   r#   )r¼  r½  r¾  r¿  rÁ  rÂ  rÃ  rÄ  rÅ  rÆ  rÇ  rÈ  r
   r
   r   Ú
test_fk524  s     rÎ  c                  C   s¨   d} d}d}d}d}d}t  | |||||¡\}}}}	}
}t|ddd	d
tƒ t|ddd	dtƒ t|ddd	dtƒ t|	ddd	dtƒ t|
ddd	dtƒ t|ddd	dtƒ d S )Nç¼wµâHü?ç¤˜’]¬Ò¿gä¶Fø¿‰¾g +ù3†Ø¾ç÷;úDØ?çffffffÀgA@™âHü?r6   Úfk52hrü   gªI¶î¬Ò¿Údecg(!»òTŠ¾g¬Ò¶OÉƒý;Údr5gó3£…Ø¾Údd5r  gâNµlffÀrU   r  )r"   rÓ  r   r#   )Úr5Úd5rÕ  rÖ  Úpx5Úrv5rà   ÚdhÚdrhÚddhÚpxhÚrvhr
   r
   r   Ú
test_fk52h-  s     rà  c                  C   sn   d} d}d}t  | ||¡\}}}}t|ddddtƒ t|dd	dd
tƒ t|ddddtƒ t|ddddtƒ d S )NgÌ·CáÃ×›?gÄœ«6›�¼¿g~7E�µŠž@gˆdv–g�?r6   Úfk54zrÃ  gÎþá*þÎ½¿rÀ  rÄ  gçØ¤sŒ?I¾rt  rÅ  gÃ„Ý¹¢V>rÆ  )r"   rá  r   r#   )r¼  r½  rË  rÃ  rÄ  rÅ  rÆ  r
   r
   r   Ú
test_fk54z=  s    râ  c                  C   sT  t  dt¡} t  dt¡}t ¡ \} }t| d d ddddtƒ t| d d d	d
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d ddddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d ddddtƒ t|d dd
ddtƒ t|d dd
ddtƒ t|d dd
dd tƒ d S )!Nrg  r,   r   gÀÿÿÿÿÿï?r6   Úfk5hipr¼  r   g›ÿ’eÍ}>r¹   r½  r   gª\¢…‰¯g>r¿  góÈÓeÍ}¾rÀ  gŸÿÿÿÿÿï?rÁ  gn¢›éåy¾rÂ  gˆ·æ‰¯g¾rÄ  gð…ù„éåy>rÆ  gÍÿÿÿÿÿï?rÈ  gçS	³ü¾Ús1gçS	³ü(>Ús2gŽØßÐ&->Zs3)r/   r0   r1   r"   rã  r   r#   )Zr5hZs5hr
   r
   r   Útest_fk5hipH  s    ræ  c                  C   sD   d} d}t  | |dd¡\}}t|ddddtƒ t|d	ddd
tƒ d S )NrÏ  rÐ  rh  ç    à™ê@g`æ�âHü?r-   Úfk5hzrü   gK½Qç¬Ò¿rÔ  )r"   rè  r   r#   )r×  rØ  rà   rÛ  r
   r
   r   Ú
test_fk5hz[  s
    ré  c                  C   s  t  dt¡} d}d}d}d}t ||||¡} t| d d ddd	d
tƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ t| d d ddd	dtƒ d S )Nrg  ç¥;9j£ÑÂ¾ç6Dä·.Ú?çÁDÑ"O¿ç½%nÅ·.Ú?r   r¤  r-   Úfw2mr¼  r   gä?ì;~L?r½  r   g“ÁoèWÅ8?r¿  g>¬ò;~L¿rÀ  r¨  rÁ  g
©à)òF‚¾rÂ  ç±\„ÍWÅ8¿rÄ  çËM‚×Õ‰¾rÆ  r¬  rÈ  )r/   r0   r1   r"   rî  r   r#   )r3  ÚgambÚphibÚpsiÚepsr
   r
   r   Ú	test_fw2mc  s    rõ  c                  C   sL   d} d}d}d}t  | |||¡\}}t|ddddtƒ t|d	ddd
tƒ d S )Nrê  rë  rì  rí  rï  r6   Úfw2xyr˜   rð  r™   )r"   rö  r   r#   )rñ  rò  ró  rô  r˜   r™   r
   r
   r   Ú
test_fw2xyu  s    r÷  c                  C   s@   d} d}t  | |¡\}}t|ddddtƒ t|ddddtƒ d S )	Nç‡DçJW@ç-DTû!é¿çFÖ8¼@r6   Úg2icrsrJ  çïª´?Ûò¿ÚD)r"   rû  r   r#   )r*  r[  rY  rZ  r
   r
   r   Útest_g2icrs  s
    rþ  c                  C   sp  t  dddg¡} t d| ¡\}}}}t|dddtƒ t tj| ¡\}}}}t|dddtƒ t|d	d
ddtƒ t|dd
ddtƒ t|ddddtƒ t tj	| ¡\}}}}t|dddtƒ t|d	d
ddtƒ t|dd
ddtƒ t|ddddtƒ t tj
| ¡\}}}}t|dddtƒ t|d	d
ddtƒ t|dd
ddtƒ t|ddddtƒ t d| ¡\}}}}t|dddtƒ d S )Nç    €„>Aç    `ãFAç    TAr   rw  Úgc2gdry  rG  ç›ö�Òsï?r6   Úe1g_ÌVÂ\ï?r'  gÿ�J
­¶t@ç:Œ0âŽyE>Zh1rI  Úe2g‰çXÂ\ï?r(  gÉi*U­¶t@Zh2rK  Ze3g;�0²\ï?r)  gl9±ãnÔt@Zh3r   Új4)r/   Úarrayr"   r  r   r#   r€  r�  r   r‚  rƒ  )Úxyzr:   r;   r8   rA   r
   r
   r   Ú
test_gc2gd‡  s(    r
  c                  C   sx   d} d}t  dddg¡}t | ||¡\}}}}t|dddtƒ t|d	d
ddtƒ t|dd
ddtƒ t|ddddtƒ d S )Nç    ¦TXAç™o3ËTwk?rÿ  r   r  r   Úgc2gderA   r  r6   r:   g·Ì¼¼\ï?r;   g ™HãåÅt@r  r8   )r/   r  r"   r  r   r#   r   )r   r>  r	  r:   r;   r8   rA   r
   r
   r   Útest_gc2gdež  s    r  c                  C   sœ  d} d}d}t  dt¡}t d| ||¡\}}t|dddtƒ t tj| ||¡\}}t|ddd	tƒ t	|d d
dddtƒ t	|d ddddtƒ t	|d ddddtƒ t tj
| ||¡\}}t|dddtƒ t	|d ddddtƒ t	|d ddddtƒ t	|d ddddtƒ t tj| ||¡\}}t|dddtƒ t	|d ddddtƒ t	|d ddddtƒ t	|d dddd tƒ t d!| ||¡\}}t|ddd"tƒ d S )#Nr$  r#  ç     ˆ£@r,   r   rw  Úgd2gcry  rG  gJp±#Æ[UÁçH¯¼šò×z>z1/1r   gx¢¼`�qAz2/1r   gÔ×?¼F3GÁz3/1rI  gÎÈ±#Æ[UÁz1/2gY½`�qAz2/2gï<¼F3GÁz3/2rK  g¯îÎ°Å[UÁz1/3gÈ±½ÇœqAz2/3gÊ§ÓWF3GÁz3/3r   r  )r/   r0   r1   r"   r  r   r#   r€  r�  r   r‚  rƒ  )r:   r;   r8   r	  rA   r
   r
   r   Ú
test_gd2gc©  s.    r  c                  C   sŒ   d} d}d}d}d}t  dt¡}t | ||||¡\}}t|ddd	tƒ t|d d
dddtƒ t|d ddddtƒ t|d ddddtƒ d S )Nr  r  r$  r#  r  r,   r   Úgd2gcerA   gh ¨êÅ[UÁr  r   r   g–Ç�qAr   r   gŒ"Ÿ†F3GÁr!   )r/   r0   r1   r"   r  r   r#   r   )r   r>  r:   r;   r8   r	  rA   r
   r
   r   Útest_gd2gceÃ  s    r  c                  C   s&   t  dddd¡} t| ddddtƒ d S )Nrh  rº  gžµ‘Æü?r-   Úgmst00rE   )r"   r  r   r#   ©rÈ   r
   r
   r   Útest_gmst00Ñ  s    r  c                  C   s&   t  dddd¡} t| ddddtƒ d S )Nrh  rº  gÙLzÆü?r-   Úgmst06rE   )r"   r  r   r#   r  r
   r
   r   Útest_gmst06Ö  s    r  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  rº  gÙ(Éü?r-   Úgmst82rE   )r"   r  r   r#   r  r
   r
   r   Útest_gmst82Û  s    r  c                  C   s&   t  dddd¡} t| ddddtƒ d S )Nrh  rº  gR (ƒü?r-   Úgst00arE   )r"   r  r   r#   r  r
   r
   r   Útest_gst00aà  s    r  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  rº  gTÁ‚ü?r-   Úgst00brE   )r"   r  r   r#   r  r
   r
   r   Útest_gst00bå  s    r  c                  C   s    t  dt¡} d| d d< d| d d< d| d d< d| d d< d	| d d< d
| d d< d| d d< d| d d< d| d d< t dddd| ¡}t|ddddtƒ d S )Nrg  r  r   r  r   r  r   r	  r
  r  rÃ  rù  rÇ  rh  rº  gRž(ƒü?r-   Úgst06rE   )r/   r0   r1   r"   r   r   r#   )rˆ  rÈ   r
   r
   r   Ú
test_gst06ê  s    r!  c                  C   s&   t  dddd¡} t| ddddtƒ d S )Nrh  rº  g–	ƒü?r-   Úgst06arE   )r"   r"  r   r#   r  r
   r
   r   Útest_gst06aù  s    r#  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  rº  g‰YŸ‚ü?r-   Úgst94rE   )r"   r$  r   r#   r  r
   r
   r   Ú
test_gst94þ  s    r%  c                  C   s¨   d} d}d}d}d}d}t  | |||||¡\}}}}	}
}t|ddd	d
tƒ t|ddd	dtƒ t|ddd	dtƒ t|	ddd	dtƒ t|
ddd	dtƒ t|ddd	dtƒ d S )NçÓKþºâHü?çä%ly¬Ò¿g„G®ì/Ç¾g‡÷Íº:ßØ¾rÑ  rÒ  g¬‚ÔÖâHü?rÀ  Úh2fk5rü   guHè¬Ò¿rÔ  g4Ý¦�&Ç¾rÆ   rÕ  gx3$G&àØ¾rÖ  r  gø]/offÀrU   r  )r"   r(  r   r#   )rà   rÛ  rÜ  rÝ  rÞ  rß  r×  rØ  rÕ  rÖ  rÙ  rÚ  r
   r
   r   Ú
test_h2fk5  s     r)  c                  C   sF   d} d}d}t  | ||¡\}}t|ddddtƒ t|dd	dd
tƒ d S )Nr5   r  r  gÍ�åª@rÀ  Úhd2aer   gÌð	ä:ŸÜ?r6   r:   )r"   r*  r   r#   )r8   r9   r;   r   r:   r
   r
   r   Ú
test_hd2ae  s    r+  c                  C   s0   d} d}d}t  | ||¡}t|ddddtƒ d S )Nr5   r  r  g¶àTèþ?rÀ  Úhd2paÚq)r"   r,  r   r#   )r8   r9   r;   r-  r
   r
   r   Ú
test_hd2pa  s
    r.  c                  C   sl   d} d}t  | |dd¡\}}}}t|ddddtƒ t|d	d
ddtƒ t|ddddtƒ t|ddddtƒ d S )Nr&  r'  rh  rç  g‹c.àâHü?rÀ  Úhfk5zrü   g:­ï¬Ò¿r6   rÔ  glÚHÐ\Ÿ2>çæ^ 9^;rÕ  g³Ì]ñq¾gQ²@³-(;rÖ  )r"   r/  r   r#   )rà   rÛ  r×  rØ  rÕ  rÖ  r
   r
   r   Ú
test_hfk5z$  s    r1  c                  C   s@   d} d}t  | |¡\}}t|ddddtƒ t|ddddtƒ d S )	Nrú  rü  rø  r6   Úicrs2gÚLrù  ÚB)r"   r2  r   r#   )rY  rZ  r*  r[  r
   r
   r   Útest_icrs2g.  s
    r5  c                  C   sn  t  dt¡} d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< t ¡ } t| d d d	d
ddtƒ t| d d d
d
ddtƒ t| d d d
d
ddtƒ t| d d d
d
ddtƒ t| d d d	d
ddtƒ t| d d d
d
ddtƒ t| d d d
d
ddtƒ t| d d d
d
ddtƒ t| d d d	d
ddtƒ d S )Nrg  r.  r   r/  r   r   r0  r1  rk   r   Úirr¼  r½  r¿  rÀ  rÁ  rÂ  rÄ  rÆ  rÈ  )r/   r0   r1   r"   r6  r   r#   ©r3  r
   r
   r   Útest_ir6  s(    r8  c                  C   st   d} d}t  | |¡\}}}}}t|dddtƒ t|dddtƒ t|ddd	tƒ t|d
dddtƒ t|dddtƒ d S )Nrh  ri  éÌ  Újd2calr™   r   r<  é
   r9   g§èH.ÿï?r  Úfdr   rA   )r"   r:  r   r#   r   )Údj1Údj2r?  r@  rA  r<  rA   r
   r
   r   Útest_jd2calN  s    r?  c                  C   s~   d} d}t  d| |¡\}}t|d dddtƒ t|d d	dd
tƒ t|d	 dddtƒ t|d dddtƒ t|dddtƒ d S )Nrh  ri  r   r   r9  Újdcalfr™   r   r   r<  r;  r9   r    i'  r>  rA   )r"   r@  r   r#   )r=  r>  ZiydmfrA   r
   r
   r   Útest_jdcalfY  s    rA  c                  C   sÞ   t  dt¡} t  dt¡}t  dt¡}t  dt¡}d}d| d< d| d< d| d< d|d< d|d< d|d< d	|d< d
|d< d|d< d}d}t || ||||¡}t|d ddddtƒ t|d ddddtƒ t|d ddddtƒ d S )Nr,   r  ç7:!SÂlè¿r   ç{3�…íyã¿r   ç1lˆ‚ý½Ë¿r   g£»iL‹è?g¿ÉœRQaã?gÅ3Á Ë?g?@1‘VÓ!@r  gÀSÂlè¿r-   Úldr   gÖ…›…íyã¿r   g=Ñ‚ý½Ë¿r!   )r/   r0   r1   r"   rE  r   r#   )r;   r-  r:   r'  r)  rb   Zdlimr
   r
   r   Útest_ldd  s(    rF  c                  C   s>  t  dtj¡ t j¡} t  dt¡}t  dt¡}t  dt¡}d| d _d| d _d| d j	d d< d| d j	d d	< d
| d j	d d< d| d j	d d< d| d j	d d	< d| d j	d d< d| d	 _d| d	 _d| d	 j	d d< d| d	 j	d d	< d| d	 j	d d< d| d	 j	d d< d| d	 j	d d	< d| d	 j	d d< d| d _d| d _d| d j	d d< d| d j	d d	< d| d j	d d< d| d j	d d< d| d j	d d	< d| d j	d d< d |d< d!|d	< d"|d< d#|d< d$|d	< d%|d< t 
| ||¡}t|d d&d'd(d)tƒ t|d	 d*d'd(d+tƒ t|d d,d'd(d-tƒ d S ).Nr    r,   r  r   r  r  r;   r  r   r  r   r  r2   r  r  r  r  r  r  r  r  r  r  rk   r   r!  r"  r#  r$  r%  r&  g¤ôzg,ï¿gÐ+ù3Ë¿gd™y}·¿rB  rC  rD  gYT¹TÂlè¿r-   Úldnr   gÖ/¾ƒíyã¿r   gö2†€ý½Ë¿r!   )r/   r0   r"   r(  rv   rw   r1   r)  r*  rä   rG  r   r#   )r+  ÚobÚscZsnr
   r
   r   Útest_ldn{  sL    





rJ  c                  C   s¬   t  dt¡} t  dt¡}t  dt¡}d| d< d| d< d| d< d|d< d	|d< d
|d< d}t | ||¡}t|d ddddtƒ t|d ddddtƒ t|d ddddtƒ d S )Nr,   rB  r   rC  r   rD  r   gMc—‚(ï¿g{ø®àÈÊ¿g ÕÆó99·¿g©Ò‹äoþï?gÄ˜ÇTÂlè¿r-   Úldsunr   gIÆ²ƒíyã¿r   g9i=€ý½Ë¿r!   )r/   r0   r1   r"   rK  r   r#   )r;   r:   r'  rb   r
   r
   r   Ú
test_ldsun¤  s    rL  c                  C   sF   d} d}d}t  | ||¡\}}t|ddddtƒ t|dddd	tƒ d S )
Nr  ç      ø?ç333333ã?gP3 
gô?r6   ÚlteceqrY  gÛLÈö�äï?rZ  )r"   rO  r   r#   )r�  r*  r[  rY  rZ  r
   r
   r   Útest_lteceqµ  s    rP  c                  C   s  t  dt¡} d}t |¡} t| d d ddddtƒ t| d d d	ddd
tƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ d S )Nrg  g     p§Àr   g«¢šAnÏÖ?r6   Últecmr^  r   g¾ìªe?Lë?r_  r   gÍáE=†dØ?r`  g3÷½�æí¿ra  g£ñÿâ>ÉÔ?rb  gFw éìÉÂ?rc  gÈ¾§·tW?rd  gBõÝ—#Ú¿re  g6ÉDŠ—5í?rf  )r/   r0   r1   r"   rQ  r   r#   )rg  r�  r
   r
   r   Ú
test_ltecm¾  s    
rR  c                  C   sF   d} d}d}t  | ||¡\}}t|ddddtƒ t|dddd	tƒ d S )
Nç     p—Àrè   r–  gÛmDRX à?r6   Últeqecr*  gZøì’·â?r[  )r"   rT  r   r#   )r�  rY  rZ  r*  r[  r
   r
   r   Útest_lteqecÍ  s    rU  c                  C   s  t  dt¡} d}t |¡} t| d d ddddtƒ t| d d d	ddd
tƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ d S )Nrg  ç‹lçû©
š@r   gàÍ§ åï?r6   Últpr}  r   gRšÚp
³?r~  r   gð{ž“ ?r  gÎeß:p
³¿r€  g–çÈ›Méï?r�  g·¿o\ÁS¿r‚  g]›Õ�“ ¿rƒ  gˆ‘ÄÅ‡ÃS¿r„  gï³ûï?r…  )r/   r0   r1   r"   rW  r   r#   )rš  r�  r
   r
   r   Útest_ltpÖ  s    
rX  c                  C   s  t  dt¡} d}t |¡} t| d d ddddtƒ t| d d d	ddd
tƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ d S )Nrg  rV  r   gÙ:© åï?r6   ÚltpbZrpb11r   g ±?ãn
³?Zrpb12r   gõA¡“ ?Zrpb13gúµIo
³¿Zrpb21gÇô¡žMéï?Zrpb22g‡›jx?ÁS¿Zrpb23gl]þ‰ “ ¿Zrpb31gT{¤Ê¨ÃS¿Zrpb32g0$˜
³ûï?Zrpb33)r/   r0   r1   r"   rY  r   r#   )Zrpbr�  r
   r
   r   Ú	test_ltpbå  s    
rZ  c                  C   s`   t  dt¡} d}t |¡} t| d ddddtƒ t| d d	ddd
tƒ t| d ddddtƒ d S )Nr,   rS  r   gœÒú!m@??r6   ÚltpeclZvec1r   g<=Ö8ïÙ¿Zvec2r   g#‘÷ÕBAí?Zvec3)r/   r0   r1   r"   r[  r   r#   )Zvecr�  r
   r
   r   Útest_ltpeclô  s    
r\  c                  C   s`   t  dt¡} d}t |¡} t| d ddddtƒ t| d d	ddd
tƒ t| d ddddtƒ d S )Nr,   g     ˆ£Àr   ghV5_ôÖ¿r6   ÚltpequZveq1r   g»Ð S³�É¿Zveq2r   g@”ë-í?Zveq3)r/   r0   r1   r"   r]  r   r#   )Zveqr�  r
   r
   r   Útest_ltpequý  s    
r^  c                  C   sº   t  dtj¡} t dd¡} t| d d dddd	tƒ t| d d
 ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d
 ddddtƒ t| d d ddddtƒ d S )Nr
   rh  çÍÌÌÌü{å@r;   r   gñËh"POe¿rU   Úmoon98úx 4r   g–×¾yeD?úy 4r   g¥	�„N1?úz 4r2   g2ø	@O ¿úxd 4g€’<=A¿úyd 4gµêæXb~&¿úzd 4)r/   r0   r"   ru   r`  r   r#   ©rä   r
   r
   r   Útest_moon98  s    rh  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  rº  r   ç? ùÿÿÿï?r-   Únum00ar¼  r   ç‡Í`ë‡â>r½  r   ç[ÔøSÐ>r¿  ç.ÂâûÕ‡â¾rÀ  çæ-‰ÿÿÿï?rÁ  çÍÐ;î–M¿rÂ  çÙ= ©¶Ð¾rÄ  ç�yÄ•M?rÆ  çC‰�ÿÿÿï?rÈ  )r/   r0   r1   r"   rj  r   r#   ©Úrmatnr
   r
   r   Útest_num00a  s    ru  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  rÚ  r   ç±Ÿùÿÿÿï?r-   Únum00br¼  r   ç¤¬É”ºˆâ>r½  r   çHèà¡Ð>r¿  çB¯þ.¥ˆâ¾rÀ  ç.‰ÿÿÿï?rÁ  çI ÿyˆM¿rÂ  ç</¹VjÐ¾rÄ  ç å"P‡M?rÆ  çÆ‰�ÿÿÿï?rÈ  )r/   r0   r1   r"   rw  r   r#   rs  r
   r
   r   Útest_num00b  s    r€  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  rÚ  r   ri  r-   Únum06ar¼  r   çíTÚë‡â>r½  r   çŠè‹ÌSÐ>r¿  çYFjuÖ‡â¾rÀ  çè-‰ÿÿÿï?rÁ  ç�¬²–M¿rÂ  çXY}¶Ð¾rÄ  ç]Ôéˆ•M?rÆ  çE‰�ÿÿÿï?rÈ  )r/   r0   r1   r"   r�  r   r#   rs  r
   r
   r   Útest_num06a-  s    rŠ  c                  C   s  t  dt¡} d}d}d}t |||¡} t| d d dddd	tƒ t| d d
 ddddtƒ t| d d ddddtƒ t| d
 d ddddtƒ t| d
 d
 ddddtƒ t| d
 d ddddtƒ t| d d ddddtƒ t| d d
 ddddtƒ t| d d ddddtƒ d S )Nrg  r  r  ç}ÁuY–M?r   ri  r-   Únumatr¼  r   gìñ|ë‡â>r½  r   g6÷ÙzSÐ>r¿  gNO2Ö‡â¾rÀ  gå-‰ÿÿÿï?rÁ  gÍ;î–M¿rÂ  g&G¦+¶Ð¾rÄ  g5yÄ•M?rÆ  rr  rÈ  )r/   r0   r1   r"   rŒ  r   r#   )rt  rj  r  r  r
   r
   r   Ú
test_numat;  s    r�  c                  C   s8   t  dd¡\} }t| ddddtƒ t|ddddtƒ d S )	Nrh  rº  çÑ”Ž2ä¾rÀ  Únut00ar  r‹  r  )r"   r�  r   r#   ©r  r  r
   r
   r   Útest_nut00aL  s    r‘  c                  C   s8   t  dd¡\} }t| ddddtƒ t|ddddtƒ d S )	Nrh  rº  çÎúço3ä¾rÀ  Únut00br  çE�,å‡M?r  )r"   r“  r   r#   r�  r
   r
   r   Útest_nut00bR  s    r•  c                  C   s8   t  dd¡\} }t| ddddtƒ t|ddddtƒ d S )	Nrh  rº  çÛqExŽ2ä¾rÀ  Únut06ar  ç�~æ–M?r  )r"   r—  r   r#   r�  r
   r
   r   Útest_nut06aX  s    r™  c                  C   s8   t  dd¡\} }t| ddddtƒ t|ddddtƒ d S )	Nrh  rº  g·(ÌOf9ä¾rÀ  Únut80r  g„qPçNI?r  )r"   rš  r   r#   r�  r
   r
   r   Ú
test_nut80^  s    r›  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  rº  r   gí›ùÿÿÿï?r-   Únutm80r¼  r   gÐh
	3Žâ>r½  r   g„WóÅÐ>r¿  gu«-¡Žâ¾rÀ  gÒW‰ÿÿÿï?rÁ  gLˆ{|OI¿rÂ  g¼<»±(Ð¾rÄ  g\ÞîQNI?rÆ  g#¶�ÿÿÿï?rÈ  )r/   r0   r1   r"   rœ  r   r#   rs  r
   r
   r   Útest_nutm80d  s    r�  c                   C   s   t t dd¡ddddtƒ d S )Nrh  r›  gÐ¬è•H.Ú?r6   Úobl06rE   )r   r"   rž  r#   r
   r
   r
   r   Ú
test_obl06r  s    rŸ  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  r›  g2ÞZyI.Ú?r6   Úobl80rE   )r"   r   r   r#   )Úeps0r
   r
   r   Ú
test_obl80v  s    r¢  c                  C   sP  t  dd¡\} }}}}}}}}}	}
}}}}}t| ddddtƒ t|ddddtƒ t|d	ddd
tƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|	ddddtƒ t|
ddddtƒ t|ddddtƒ t|ddddtƒ t|dddd tƒ t|d!ddd"tƒ t|d#ddd$tƒ d S )%Nrh  ç     §é@g±]'»’.Ú?r6   Úp06er¡  g)x{zÖE?Úpsiag~£±·’.Ú?ÚomagÈó­ý¨¢>Úbpag$?rIøÙ¾ÚbqagwmUÚÚ>Úpiag
i†j@ÚbpiagÈó4¿x.Ú?rj  g¡Lo $H·>Úchiagrù/ó%3?Úzag<wz»Õ4?ÚzetaagYþì†_1?Úthetaagà¸¸PÌËE?ÚpagWÕ§²t·>ÚgamgX-Ãx.Ú?rÊ   gã»oëŽÖE?ró  )r"   r¤  r   r#   )r¡  r¥  r¦  r§  r¨  r©  rª  rj  r«  r¬  r­  r®  r¯  r°  rÊ   ró  r
   r
   r   Ú	test_p06e{  s"    ,r±  c                  C   s$  t  dt¡} t  dtj¡}d| d< d| d< d| d< d	|d
 d< d|d
 d< d|d
 d< d|d d< d|d d< d|d d< t | ¡}t|d
 d ddddtƒ t|d
 d ddddtƒ t|d
 d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ d S )Nr,   r
   g      Ð?r   r  r   r/  r   r  r;   r  r#  r2   r$  r%  r   Úp2pvr'  r(  r)  r*  r+  r,  )r/   r0   r1   r"   ru   r²  r   r#   )r;   rä   r
   r
   r   Ú	test_p2pv�  s$    
r³  c                  C   sn   t  dt¡} d| d< d| d< d| d< t | ¡\}}}t|dd	d
dtƒ t|dd	d
dtƒ t|ddd
dtƒ d S )Nr,   rã  r   rä  r   rå  r   ræ  r-   Úp2srÈ   rè  rÊ   g°ÿ0¤\@rŒ  r3  )r/   r0   r1   r"   r´  r   r#   )r;   rÈ   rÊ   r3  r
   r
   r   Útest_p2s¤  s    rµ  c                  C   sj   t  dt¡} t  dt¡}d| d< d| d< d| d< d|d< d	|d< d|d< t | |¡}t|d
dddtƒ d S )Nr,   rk   r   rø   r   çš™™™™™É?r   ç      Àçü©ñÒMbP?gxë­²h×?r-   ÚpaprE   )r/   r0   r1   r"   r¹  r   r#   )r   r+  rÈ   r
   r
   r   Útest_pap¯  s    rº  c                  C   s6   d} d}d}d}t  | |||¡}t|ddddtƒ d S )	Nrk   rø   r¶  ç      ð¿gÇÐc5ÞËÀr-   ÚpasrE   )r"   r¼  r   r#   )ÚalÚapÚblÚbprÈ   r
   r
   r   Útest_pas¼  s    rÁ  c                  C   sL   t  dd¡\} }}t| ddddtƒ t|ddddtƒ t|d	ddd
tƒ d S )Nrh  ri  gËÙÓ°@¿r-   Úpb06Úbzetag·)HÑrœ7¿ÚbzgÃ¬MXÅ8¿Úbtheta)r"   rÂ  r   r#   )rÃ  rÄ  rÅ  r
   r
   r   Ú	test_pb06Å  s    rÆ  c                  C   sj   t  dt¡} t  dt¡}d| d< d| d< d| d< d|d< d|d< d|d< t | |¡}t|d	d
ddtƒ d S )Nr,   r.  r   r   r/  r   rk   r0  é   r-   ÚpdprE   )r/   r0   r1   r"   rÈ  r   r#   ©r   r+  Zadbr
   r
   r   Útest_pdpÌ  s    rÊ  c                  C   s`   t  dd¡\} }}}t| ddddtƒ t|dddd	tƒ t|d
dddtƒ t|ddddtƒ d S )Nrh  ri  g­;9j£ÑÂ¾rd   Úpfw06rñ  rë  r-   rò  gÐDÑ"O¿r6   Úpsibrí  rj  )r"   rË  r   r#   )rñ  rò  rÌ  rj  r
   r
   r   Ú
test_pfw06Ù  s
    rÍ  c                  C   sn  t  dtj¡} t ddd¡\} }t| d d ddddtƒ t| d d	 dddd
tƒ t| d d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ t|dddtƒ t ddd¡\} }t|dddtƒ t ddd¡\} }t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d	 d ddd!tƒ t| d d d"ddd#tƒ t|d	dd$tƒ t dd%d	¡\} }t| d d d&ddd'tƒ t| d d	 d(ddd)tƒ t| d d d*ddd+tƒ t| d d d,ddd-tƒ t| d d	 d.ddd/tƒ t| d d d0ddd1tƒ t|ddd2tƒ d S )3Nr
   rh  g    €„.Ar   r;   r   Úplan94zx 1r   zy 1r   zz 1r2   zxd 1zyd 1zzd 1rw  zj 1r;  zj 2i ûÿr    g¸¤Ì0%Éí?rU   zx 3g<HònÚÔ?zy 3gß©ß‘66Â?zz 3gâ�N¿ÆUz¿zxd 3g?„¼+ŠŽ?zyd 3gÞÜ9>z?zzd 3zj 3r_  g‡ŒÍÒ‘ÙÒ?ra  g±¯õacÏ¿rb  gN„&ùn­Ä¿rc  g±Ü{U½ôŒ?rd  g’°»sÂî“?re  gBÌÙÇÊI‚?rf  zj 4)r/   r0   r"   ru   rÎ  r   r#   r   )rä   rA   r
   r
   r   Útest_plan94á  s6    rÏ  c                  C   sD   t  dt¡} d| d< d| d< d| d< t | ¡}t|dd	d
dtƒ d S )Nr,   r  r   r  r   r  r   ç°^=jP@r-   ÚpmrE   )r/   r0   r1   r"   rÑ  r   r#   )r;   r3  r
   r
   r   Útest_pmÿ  s    
rÒ  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  ri  r   r†  r-   Úpmat00r¼  r   rˆ  r6   r½  r   rŠ  r¿  rŒ  rÀ  rŽ  rÁ  r�  rÂ  r’  rÄ  r”  rÆ  r–  rÈ  )r/   r0   r1   r"   rÓ  r   r#   ©r›  r
   r
   r   Útest_pmat00	  s    rÕ  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  ri  r   r¤  r-   Úpmat06r¼  r   r¥  r6   r½  r   r¦  r¿  r§  rÀ  r¨  rÁ  r©  rÂ  rª  rÄ  r«  rÆ  r¬  rÈ  )r/   r0   r1   r"   rÖ  r   r#   rÔ  r
   r
   r   Útest_pmat06	  s    r×  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  ri  r   g—õ£ÿÿï?r-   Úpmat76r¼  r   g^=×‹CL?r6   r½  r   gnã^Ä8?r¿  gSL×‹CL¿rÀ  gWzú4ÿÿï?rÁ  gAŸD“H†¾rÂ  g—Ãâ^Ä8¿rÄ  g²ðzæ[†¾rÆ  g@{©Ùÿÿï?rÈ  )r/   r0   r1   r"   rØ  r   r#   )Zrmatpr
   r
   r   Útest_pmat76$	  s    rÙ  c                  C   s¦   t  dt¡} t  dt¡}t  dt¡}d| d< d| d< d| d< d|d< d|d< d|d< t | |¡}t|d dd	d
dtƒ t|d dd	d
dtƒ t|d dd	d
dtƒ d S )Nr,   r.  r   r   r/  r   rk   r0  r-   Úpmpr   r»  r   r   )r/   r0   r1   r"   rÚ  r   r#   ©r   r+  Zambr
   r
   r   Útest_pmp2	  s    rÜ  c            	   
   C   sª   t  dt¡} t  dt¡}d}d}d}d}d}d}d}d	| d
< d| d< d| d< t |||||||| ¡}t|d
 ddddtƒ t|d ddddtƒ t|d ddddtƒ d S )Nr,   rè   çsh‘í|?é?rö   çñhãˆµøô¾ç{®Gáz„?ç      $@g     €!@r%  r   çš™™™™™Ù?r   rø   r   g~[Æc×ÌÍ?r-   Úpmpxr   gkÆÍ'”Hå?r   g=ãÃfo´æ?r!   )r/   r0   r1   r"   râ  r   r#   )	ZpobZpcorÿ   r   r  r  r  r  rW   r
   r
   r   Ú	test_pmpxB	  s     rã  c                  C   sÒ   d} d}d}d}d}d}d}d}d}d	}	t  | |||||||||	¡
\}
}}}}}}t|
d
dddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|dddtƒ d S )Nrè   rÝ  rö   rÞ  rß  rà  rh  g    ’›ç@g    +é@g=[¥„Ò¾ó?r-   ÚpmsafeÚra2gúY7ì>é?Údec2g›€®Îöä>Úpmr2gÜ ¼Ñøô¾rd   Úpmd2gÉ¿êàz„?Úpx2g|Ù�FõÄ$@rz  Úrv2r   rA   )r"   rä  r   r#   r   )Úra1Údec1Úpmr1Úpmd1Úpx1Úrv1Zep1aZep1bZep2aZep2brå  ræ  rç  rè  ré  rê  rA   r
   r
   r   Útest_pmsafeU	  s$    *rñ  c                  C   s–   t  dt¡} t  dt¡}d| d< d| d< d| d< t | ¡\}}t|dd	d
dtƒ t|d dd	d
dtƒ t|d dd	d
dtƒ t|d dd	d
dtƒ d S )Nr,   r  r   r  r   r  r   rÐ  r-   Úpnr3  g}\½ˆ»?rT  g}\½ˆÛ?rU  gÿ¬Œo®ì¿Úu3)r/   r0   r1   r"   rò  r   r#   )r;   rR  r3  r
   r
   r   Útest_pnj	  s    rô  c                  C   s  t  dt¡} t  dt¡}t  dt¡}t  dt¡}t  dt¡}d}d}t dd||¡\}} }}}}t|dddd	tƒ t| d
 d
 ddddtƒ t| d
 d ddddtƒ t| d
 d ddddtƒ t| d d
 ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d
 ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t|d
 d
 d ddd!tƒ t|d
 d d"d#dd$tƒ t|d
 d d%d#dd&tƒ t|d d
 d'd#dd(tƒ t|d d d)ddd*tƒ t|d d d+d#dd,tƒ t|d d
 d-d#dd.tƒ t|d d d/d#dd0tƒ t|d d d1ddd2tƒ t|d
 d
 d3ddd4tƒ t|d
 d d5d#dd6tƒ t|d
 d d7d#dd8tƒ t|d d
 d9d#dd:tƒ t|d d d;ddd<tƒ t|d d d=d#dd>tƒ t|d d
 d?d#dd@tƒ t|d d dAd#ddBtƒ t|d d dCdddDtƒ t|d
 d
 dEdddFtƒ t|d
 d dGdHddItƒ t|d
 d dJdHddKtƒ t|d d
 dLdHddMtƒ t|d d dNdddOtƒ t|d d dPdHddQtƒ t|d d
 dRdHddStƒ t|d d dTdHddUtƒ t|d d dVdddWtƒ t|d
 d
 dXdddYtƒ t|d
 d dZd#dd[tƒ t|d
 d d\d#dd]tƒ t|d d
 d^d#dd_tƒ t|d d d`dddatƒ t|d d dbd#ddctƒ t|d d
 ddd#ddetƒ t|d d dfd#ddgtƒ t|d d dhddditƒ d S )jNrg  çúço3ä¾ç+�,å‡M?rh  rº  çvƒÁoZ.Ú?r-   Úpn00rj  r   rj  rl  r   rm  rÆ   rn  r   ro  rp  rq  rr  rs  rt  ru  rv  rw  rx  ry  rz  r{  r|  ç±¢“Áýÿï?r}  çGÉŽ½FûU¿r6   r~  çÝVfeuC¿r  ç³Ñ»¾FûU?r€  ç°Ïþÿï?r�  ç åíé­À™¾r‚  ç¾Éü_uC?rƒ  çß¨Ë Ï¼š¾r„  ç òÃ¤ÿÿï?r…  ç¯	�Áýÿï?r‡  çH¿Ý¸’ûU¿r‰  çDEO|ÈC¿r‹  ç ãY¿’ûU?r�  ç£Âþÿï?r�  çO¾oÙ†—¾r‘  ç6¿t^ÈC?r“  çAmŽÎ}õœ¾r•  ç†fÊ¤ÿÿï?r—  rv  Úrn11g³®É”ºˆâ>rd   Úrn12gêà¡Ð>Úrn13gQ±þ.¥ˆâ¾Úrn21r{  Úrn22g/ ÿyˆM¿Úrn23g1¹VjÐ¾Úrn31g†å"P‡M?Úrn32r  Úrn33ç;VÉýÿï?Úrbpn11g«éSF�ÖU¿Úrbpn12g’]o¤ùB¿Úrbpn13g˜´ÝhÖU?Úrbpn21ç©ê«"þÿï?Úrbpn22g©ïä{¿Úrbpn23gÛ˜žúB?Úrbpn31gÖIÒ¹K?Úrbpn32rÏ  Úrbpn33)r/   r0   r1   r"   rø  r   r#   ©r™  rš  r›  Zrnr¸  r  r  rj  r
   r
   r   Ú	test_pn00w	  sl    r   c                  C   s$  t  dt¡} t  dt¡}t  dt¡}t  dt¡}t  dt¡}t dd¡\}}}} }}}}t|ddddtƒ t|dddd	tƒ t|d
dddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d d ddd!tƒ t|d d d"ddd#tƒ t|d d d$d%dd&tƒ t|d d d'd%dd(tƒ t|d d d)d%dd*tƒ t|d d d+ddd,tƒ t|d d d-d%dd.tƒ t|d d d/d%dd0tƒ t|d d d1d%dd2tƒ t|d d d3ddd4tƒ t|d d d5ddd6tƒ t|d d d7d%dd8tƒ t|d d d9d%dd:tƒ t|d d d;d%dd<tƒ t|d d d=ddd>tƒ t|d d d?d%dd@tƒ t|d d dAd%ddBtƒ t|d d dCd%ddDtƒ t|d d dEdddFtƒ t|d d dGdddHtƒ t|d d dId%ddJtƒ t|d d dKd%ddLtƒ t|d d dMd%ddNtƒ t|d d dOdddPtƒ t|d d dQd%ddRtƒ t|d d dSd%ddTtƒ t|d d dUd%ddVtƒ t|d d dWdddXtƒ t|d d dYdddZtƒ t|d d d[d%dd\tƒ t|d d d]d%dd^tƒ t|d d d_d%dd`tƒ t|d d dadddbtƒ t|d d dcd%dddtƒ t|d d ded%ddftƒ t|d d dgd%ddhtƒ t|d d didddjtƒ d S )kNrg  rh  rº  rŽ  r-   Úpn00ar  r‹  r  r÷  rj  r   rj  rl  r   rm  rd   rn  r   ro  rp  rq  rr  rs  rt  ru  rv  rw  rx  ry  rz  r{  r|  rù  r}  rú  r6   r~  rû  r  rü  r€  rý  r�  rþ  r‚  rÿ  rƒ  r   r„  r  r…  r  r‡  r  r‰  r  r‹  r  r�  r  r�  r  r‘  r  r“  r	  r•  r
  r—  ri  r  rk  r  rl  r  rm  r  rn  r  ro  r  rp  r  rq  r  rr  r  r  r  gÊÅ»ä‚ÖU¿r  g²_iÖ¥ùB¿r  gþ�2§iÖU?r  r
  r  g¿ïkó{¿r  rÃ  r  rÅ  r  rÇ  r  )r/   r0   r1   r"   r!  r   r#   r  r
   r
   r   Ú
test_pn00a°	  sl    r"  c                  C   s$  t  dt¡} t  dt¡}t  dt¡}t  dt¡}t  dt¡}t dd¡\}}}} }}}}t|ddddtƒ t|dddd	tƒ t|d
dddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d d ddd!tƒ t|d d d"ddd#tƒ t|d d d$d%dd&tƒ t|d d d'd%dd(tƒ t|d d d)d%dd*tƒ t|d d d+ddd,tƒ t|d d d-d%dd.tƒ t|d d d/d%dd0tƒ t|d d d1d%dd2tƒ t|d d d3ddd4tƒ t|d d d5ddd6tƒ t|d d d7d%dd8tƒ t|d d d9d%dd:tƒ t|d d d;d%dd<tƒ t|d d d=ddd>tƒ t|d d d?d%dd@tƒ t|d d dAd%ddBtƒ t|d d dCd%ddDtƒ t|d d dEdddFtƒ t|d d dGdddHtƒ t|d d dId%ddJtƒ t|d d dKd%ddLtƒ t|d d dMd%ddNtƒ t|d d dOdddPtƒ t|d d dQd%ddRtƒ t|d d dSd%ddTtƒ t|d d dUd%ddVtƒ t|d d dWdddXtƒ t|d d dYdddZtƒ t|d d d[d%dd\tƒ t|d d d]d%dd^tƒ t|d d d_d%dd`tƒ t|d d dadddbtƒ t|d d dcd%dddtƒ t|d d ded%ddftƒ t|d d dgd%ddhtƒ t|d d didddjtƒ d S )kNrg  rh  rº  r’  r-   Úpn00br  r”  r  r÷  rj  r   rj  rl  r   rm  rd   rn  r   ro  rp  rq  rr  rs  rt  ru  rv  rw  rx  ry  rz  r{  r|  rù  r}  rú  r6   r~  rû  r  rü  r€  rý  r�  rþ  r‚  rÿ  rƒ  r   r„  r  r…  r  r‡  r  r‰  r  r‹  r  r�  r  r�  r  r‘  r  r“  r	  r•  r
  r—  rv  r  rx  r  ry  r  rz  r  r{  r  r|  r  r}  r  r~  r  r  r  r  r  g¯éSF�ÖU¿r  g–]o¤ùB¿r  gœ´ÝhÖU?r  r  r  gÄïä{¿r  rÍ  r  rÎ  r  rÏ  r  )r/   r0   r1   r"   r#  r   r#   r  r
   r
   r   Ú
test_pn00bé	  sl    r$  c                  C   s  t  dt¡} t  dt¡}t  dt¡}t  dt¡}t  dt¡}d}d}t dd||¡\}} }}}}t|dddd	tƒ t| d
 d
 ddddtƒ t| d
 d ddddtƒ t| d
 d ddddtƒ t| d d
 ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d
 ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t|d
 d
 d ddd!tƒ t|d
 d d"ddd#tƒ t|d
 d d$ddd%tƒ t|d d
 d&ddd'tƒ t|d d d(ddd)tƒ t|d d d*ddd+tƒ t|d d
 d,ddd-tƒ t|d d d.ddd/tƒ t|d d d0ddd1tƒ t|d
 d
 d2ddd3tƒ t|d
 d d4ddd5tƒ t|d
 d d6ddd7tƒ t|d d
 d8ddd9tƒ t|d d d:ddd;tƒ t|d d d<ddd=tƒ t|d d
 d>ddd?tƒ t|d d d@dddAtƒ t|d d dBdddCtƒ t|d
 d
 dDdddEtƒ t|d
 d dFdddGtƒ t|d
 d dHdddItƒ t|d d
 dJdddKtƒ t|d d dLdddMtƒ t|d d dNdddOtƒ t|d d
 dPdddQtƒ t|d d dRdddStƒ t|d d dTdddUtƒ t|d
 d
 dVdddWtƒ t|d
 d dXdddYtƒ t|d
 d dZddd[tƒ t|d d
 d\ddd]tƒ t|d d d^ddd_tƒ t|d d d`dddatƒ t|d d
 dbdddctƒ t|d d dddddetƒ t|d d dfdddgtƒ d S )hNrg  rõ  rö  rh  rº  r  r-   Úpn06rj  r   rj  rl  r   rž  r6   rn  r   rŸ  rp  r   rr  rs  rt  r¡  rv  r¢  rx  r£  rz  r{  r|  çÁ²“Áýÿï?r}  çzð_FûU¿r~  çÿú¾euC¿r  çlC#aFûU?r€  ç$ÀÏþÿï?r�  ç‡B–¾™¾r‚  çö9`uC?rƒ  çÝªhb¿š¾r„  ç�òÃ¤ÿÿï?r…  ç´�Áýÿï?r‡  ç—Ü}[’ûU¿r‰  ç¤/¬|ÈC¿r‹  çü b’ûU?r�  ç"³Âþÿï?r�  çÝœíI„—¾r‘  çì‰µ^ÈC?r“  ç•îŸrøœ¾r•  çƒfÊ¤ÿÿï?r—  rv  r  gMÅ°ºˆâ>r  g-áà#Ð>r  g­oOJ¥ˆâ¾r  r{  r  g{œÿyˆM¿r  gš¹¹ØiÐ¾r  g:é"P‡M?r  r  r  g&fÉýÿï?r  g q½è€ÖU¿r  gŽpp¤ùB¿r  gLM«gÖU?r  g×ú«"þÿï?r  g;ÃÐß{¿r  g`³ÓŸúB?r  gÞ�œF?r  g�÷ˆ¥ÿÿï?r  )r/   r0   r1   r"   r%  r   r#   r  r
   r
   r   Ú	test_pn06"
  sl    r8  c                  C   s$  t  dt¡} t  dt¡}t  dt¡}t  dt¡}t  dt¡}t dd¡\}}}} }}}}t|ddddtƒ t|dddd	tƒ t|d
dddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d d ddd!tƒ t|d d d"ddd#tƒ t|d d d$ddd%tƒ t|d d d&ddd'tƒ t|d d d(ddd)tƒ t|d d d*ddd+tƒ t|d d d,ddd-tƒ t|d d d.ddd/tƒ t|d d d0ddd1tƒ t|d d d2ddd3tƒ t|d d d4ddd5tƒ t|d d d6ddd7tƒ t|d d d8ddd9tƒ t|d d d:ddd;tƒ t|d d d<ddd=tƒ t|d d d>ddd?tƒ t|d d d@dddAtƒ t|d d dBdddCtƒ t|d d dDdddEtƒ t|d d dFdddGtƒ t|d d dHdddItƒ t|d d dJdddKtƒ t|d d dLdddMtƒ t|d d dNdddOtƒ t|d d dPdddQtƒ t|d d dRdddStƒ t|d d dTdddUtƒ t|d d dVdddWtƒ t|d d dXdddYtƒ t|d d dZddd[tƒ t|d d d\ddd]tƒ t|d d d^ddd_tƒ t|d d d`dddatƒ t|d d dbdddctƒ t|d d dddddetƒ t|d d dfdddgtƒ t|d d dhddditƒ d S )jNrg  rh  rº  r–  r-   Úpn06ar  r˜  r  r  rj  r   rj  rl  r   rž  r6   rn  r   rŸ  rp  r   rr  rs  rt  r¡  rv  r¢  rx  r£  rz  r{  r|  r&  r}  r'  r~  r(  r  r)  r€  r*  r�  r+  r‚  r,  rƒ  r-  r„  r.  r…  r/  r‡  r0  r‰  r1  r‹  r2  r�  r3  r�  r4  r‘  r5  r“  r6  r•  r7  r—  ri  r  r‚  r  rƒ  r  r„  r  r…  r  r†  r  r‡  r  rˆ  r  r‰  r  gHÉýÿï?r  g?ßh†‚ÖU¿r  gÖ«×¥ùB¿r  gÿåHiÖU?r  g®³«"þÿï?r  g1ßuî{¿r  rÒ  r  rÓ  r  rÔ  r  )r/   r0   r1   r"   r9  r   r#   r  r
   r
   r   Ú
test_pn06a[
  sl    r:  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  ri  r   gNXF ÿÿï?r-   Úpnm00ar¼  r   gNûÓCoK?r6   r½  r   ga¢�`Ú7?r¿  ggÇ?doK¿rÀ  g@µ^Cÿÿï?rÁ  gSpß«?rÂ  g~×UWËÙ7¿rÄ  gT‹ÖÚ Ô¿rÆ  g¾$ûÛÿÿï?rÈ  )r/   r0   r1   r"   r;  r   r#   ©r¸  r
   r
   r   Útest_pnm00a”
  s    r=  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  ri  r   gÁF ÿÿï?r-   Úpnm00br¼  r   gUá!{GoK?r6   r½  r   g³‡¬cÚ7?r¿  gÄÇægoK¿rÀ  g[„^Cÿÿï?rÁ  g­¥÷Sþª?rÂ  g¨�ê‚ÎÙ7¿rÄ  gvº¦ZäÓ¿rÆ  g�ûÛÿÿï?rÈ  )r/   r0   r1   r"   r>  r   r#   r<  r
   r
   r   Útest_pnm00b¢
  s    r?  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  ri  r   g\F ÿÿï?r-   Úpnm06ar¼  r   g+óˆ‚CoK?r6   r½  r   g÷»´`Ú7?r¿  g¥‹PîcoK¿rÀ  gÀ¹^Cÿÿï?rÁ  g˜‰h««?rÂ  gý�†ŠËÙ7¿rÄ  gàè§ýÓ¿rÆ  gH$ûÛÿÿï?rÈ  )r/   r0   r1   r"   r@  r   r#   r<  r
   r
   r   Útest_pnm06a°
  s    rA  c                  C   s  t  dt¡} t dd¡} t| d d ddddtƒ t| d d	 d
dddtƒ t| d d ddddtƒ t| d	 d ddddtƒ t| d	 d	 ddddtƒ t| d	 d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ d S )Nrg  rh  ri  r   g¦‹: ÿÿï?r-   Úpnm80r¼  r   g€8ƒƒTpK?r6   r½  r   g0™AönÙ7?r¿  gK¹™ètpK¿rÀ  g-@PCÿÿï?rÁ  gó~ÒÆv§?rÂ  g¾ð$ßÙØ7¿rÄ  gmBñº\Ð¿rÆ  gQþÛÿÿï?rÈ  )r/   r0   r1   r"   rB  r   r#   )Zrmatpnr
   r
   r   Ú
test_pnm80¾
  s    rC  c                  C   s  t  dt¡} d}d}d}t |||¡} t| d d dddd	tƒ t| d d
 ddddtƒ t| d d ddddtƒ t| d
 d ddddtƒ t| d
 d
 ddddtƒ t| d
 d ddddtƒ t| d d ddddtƒ t| d d
 ddddtƒ t| d d ddddtƒ d S )Nrg  rê  rë  g>Oå7‚®½r   rú  r-   Úpom00r¼  r   rû  rd   r½  r   g–4Êæ‘>r¿  rý  rÀ  rþ  rÁ  rÿ  rÂ  r   rÄ  r  rÆ  r  rÈ  )r/   r0   r1   r"   rD  r   r#   )r  rÌ   rÍ   rÎ   r
   r
   r   Ú
test_pom00Ì
  s    rE  c                  C   s¦   t  dt¡} t  dt¡}t  dt¡}d| d< d| d< d| d< d|d< d|d< d|d< t | |¡}t|d dd	d
dtƒ t|d dd	d
dtƒ t|d dd	d
dtƒ d S )Nr,   r.  r   r   r/  r   rk   r0  r-   Úpppr   r1  r   ç      @r   )r/   r0   r1   r"   rF  r   r#   ©r   r+  Zapbr
   r
   r   Útest_pppÝ
  s    rI  c                  C   s¬   t  dt¡} t  dt¡}t  dt¡}d| d< d| d< d| d< d}d|d< d|d< d	|d< t | ||¡}t|d d
dddtƒ t|d ddddtƒ t|d ddddtƒ d S )Nr,   r.  r   r   r/  r   r1  rk   r0  rG  r-   Úppspr   g      1@r   g      7@r   )r/   r0   r1   r"   rJ  r   r#   )r   r+  Zapsbr   r
   r
   r   Ú	test_ppspí
  s    rK  c                  C   s8   t  dd¡\} }t| ddddtƒ t|ddddtƒ d S )	Nrh  rÚ  g}B	Åèew¾r0  Úpr00ÚdpsiprgzXTt¢ˆ?¾Údepspr)r"   rL  r   r#   )rM  rN  r
   r
   r   Ú	test_pr00þ
  s    rO  c                  C   s`   d} d}d}d}t  | |||¡\}}}t|ddddtƒ t|dddd	tƒ t|d
dddtƒ d S )Nrh  g    @@à@g     +é@gì+‰wsäv?r-   Úprec76Úzetag3=ú[uåv?Úzgx˜M>ùæs?rÈ   )r"   rP  r   r#   )Zep01Zep02Zep11Zep12rQ  rR  rÈ   r
   r
   r   Útest_prec76  s    rS  c                  C   s²   t  dtj¡} t  dt¡}d| d d< d| d d< d| d d	< d
| d d< d| d d< d| d d	< t | ¡}t|d ddddtƒ t|d ddddtƒ t|d	 ddddtƒ d S )Nr
   r,   r  r;   r   r  r   r  r   r#  r2   r$  r%  r   Úpv2pr   r   r!   )r/   r0   r"   ru   r1   rT  r   r#   ©rä   r;   r
   r
   r   Ú	test_pv2p  s    
rV  c                  C   sÜ   t  dtj¡} d| d d< d| d d< d| d d< d	| d
 d< d| d
 d< d| d
 d< t | ¡\}}}}}}t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ d S )Nr
   çZ:lpQåÜ¿r;   r   g:[1Â%­Ÿ?r   çßJÞ¥¬?r   g-~ˆÓë>r2   g\~à@Æ>g¼K"Xª‹Å>gƒ§-â•@r-   Úpv2srÈ   rM  rÊ   gÈv¾Ÿ/Ý?r3  çÐúJÖ—[à¾rd   Útdgc7…ø4åâ>r  goP± ƒËé¾Úrd)r/   r0   r"   ru   rY  r   r#   )rä   rÈ   rÊ   r3  r[  r  r\  r
   r
   r   Ú	test_pv2s  s    r]  c                  C   sö   t  dtj¡} t  dtj¡}t  dtj¡}d| d d< d| d d< d| d d< d| d	 d< d
| d	 d< d| d	 d< d|d d< d|d d< d|d d< d
|d	 d< d|d	 d< d|d	 d< t | |¡}t|d ddddtƒ t|d ddddtƒ d S )Nr
   r.  r;   r   r   r/  r   ç      @r2   r   r0  rk   ç       @ç      4@r-   Úpvdpvr   g      I@r   )r/   r0   r"   ru   Zdt_pvdpvra  r   r#   rÉ  r
   r
   r   Ú
test_pvdpv/  s$    rb  c                  C   sŒ   t  dtj¡} d| d d< d| d d< d| d d< d	| d
 d< d| d
 d< d| d
 d< t | ¡\}}t|ddddtƒ t|ddddtƒ d S )Nr
   r  r;   r   r  r   r  r   gÍÌÌÌÌÌÜ?r2   g      Ð¿r5   rÐ  r-   Úpvmr3  gTí� “nó?r   )r/   r0   r"   ru   rc  r   r#   )rä   r3  r   r
   r
   r   Útest_pvmD  s    rd  c                  C   sf  t  dtj¡} t  dtj¡}t  dtj¡}d| d d< d| d d< d| d d< d| d	 d< d
| d	 d< d| d	 d< d|d d< d|d d< d|d d< d|d	 d< d|d	 d< d|d	 d< t | |¡}t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ d S )Nr
   r.  r;   r   r   r/  r   r1  r2   r^  rk   r0  r-   Úpvmpvr¼  r»  rÀ  rÄ  r½  rÁ  rÆ  )r/   r0   r"   ru   re  r   r#   rÛ  r
   r
   r   Ú
test_pvmpvQ  s,    rf  c                  C   sf  t  dtj¡} t  dtj¡}t  dtj¡}d| d d< d| d d< d| d d< d| d	 d< d
| d	 d< d| d	 d< d|d d< d|d d< d|d d< d|d	 d< d|d	 d< d|d	 d< t | |¡}t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ d S )Nr
   r.  r;   r   r   r/  r   r1  r2   r^  rk   r0  r-   Úpvppvr'  r(  rG  r)  r_  r*  r+  r,  )r/   r0   r"   ru   rg  r   r#   rH  r
   r
   r   Ú
test_pvppvj  s,    rh  c                  C   sî   t  dtj¡} d| d d< d| d d< d| d d< d	| d
 d< d| d
 d< d| d
 d< t | ¡\}}}}}}}t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|ddddtƒ t|dddtƒ d S )Nr
   ç§ÜuÉìþ@r;   r   çMûÞ•± @r   ç2ï4ßðÁr   ç„‘%I­˜p¿r2   gItÍµ�y¿gêkbŠ£[ˆ?ç_Ö@ÊºE‘?r-   Úpvstarrü   ç©ªVzû€ñ¿rÔ  gÄëôÞÕ²ò¾rd   Úpmrg×ÍÍ’˜Ã>Úpmdçïá’ãNéç?r  g›´K›™™5ÀrU   r  rA   )r/   r0   r"   ru   rn  r   r#   r   ©rä   rü   rÔ  rp  rq  r  r  rA   r
   r
   r   Útest_pvstarƒ  s    rt  c               	   C   sà   t  dtj¡} d}d}d}d}d}d}d}t |||||||¡} t| d	 d
 ddddtƒ t| d	 d ddddtƒ t| d	 d ddddtƒ t| d d
 ddddtƒ t| d d ddddtƒ t| d d d
d
ddtƒ d S )Nr
   r.  ç      à?g     p§@rP  g�íµ ÷Æ ¾r  r1  r;   r   g~WPArö   Úpvtobzp(1)r   g-¡›KLAzp(2)r   goV³¾3GAzp(3)r2   gƒ–îWÝÇpÀrŒ  re   gÿû”�As@rg   ri   )r/   r0   r"   ru   rv  r   r#   )rä   rÉ   rÊ   rË   rÌ   rÍ   rÎ   rÈ   r
   r
   r   Ú
test_pvtob•  s    rw  c                  C   s  t  dtj¡} t  dtj¡}d| d d< d| d d< d| d d< d	| d
 d< d| d
 d< d| d
 d< t d| ¡}t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d
 d d	dddtƒ t|d
 d ddddtƒ t|d
 d ddddtƒ d S )Nr
   ç¨ÜuÉìþ@r;   r   çMûÞ•± @r   rk  r   rl  r2   çKtÍµ�y¿çìkbŠ£[ˆ?ç     Ð¦@çS„c(ìþ@rP  Úpvur'  çNÞÅ� @r  r(  ç¿Œó	ïÁr)  r-   r*  r+  r,  )r/   r0   r"   ru   r~  r   r#   )rä   Zupvr
   r
   r   Útest_pvu§  s    r�  c                  C   s´   t  dtj¡} t  dt¡}d| d d< d| d d< d| d d	< d
| d d< d| d d< d| d d	< t d| ¡}t|d ddddtƒ t|d ddddtƒ t|d	 ddddtƒ d S )Nr
   r,   rx  r;   r   ry  r   rk  r   rl  r2   rz  r{  r|  r}  rP  Úpvupr   r  r  r   r€  r!   )r/   r0   r"   ru   r1   r‚  r   r#   rU  r
   r
   r   Ú	test_pvup¹  s    rƒ  c                  C   sf  t  dtj¡} t  dtj¡}t  dtj¡}d| d d< d| d d< d| d d< d| d	 d< d
| d	 d< d| d	 d< d|d d< d|d d< d|d d< d
|d	 d< d|d	 d< d|d	 d< t | |¡}t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ d S )Nr
   r.  r;   r   r   r/  r   r^  r2   r   r0  rk   r_  r»  r-   Úpvxpvr'  ç      Àr(  r)  g       Àr*  g      BÀr+  g      6@r,  )r/   r0   r"   ru   r„  r   r#   ©r   r+  Zaxbr
   r
   r   Ú
test_pvxpvÈ  s,    r‡  c                  C   s¦   t  dt¡} t  dt¡}t  dt¡}d| d< d| d< d| d< d|d< d|d< d|d< t | |¡}t|d d	d
ddtƒ t|d dd
ddtƒ t|d dd
ddtƒ d S )Nr,   r.  r   r   r/  r   rk   r0  r»  r-   Úpxpr   r…  r   r!   )r/   r0   r1   r"   rˆ  r   r#   r†  r
   r
   r   Útest_pxpá  s    r‰  c                  C   sL   d} d}d}d}t  | |||¡\}}t|ddddtƒ t|d	d
ddtƒ d S )Ng      ‰@rà  r%  rá  gNÅ¤2é¯-?rÄ   ÚrefcorÅ   gÛÌzýup‘¾rÆ   rÇ   )r"   rŠ  r   r#   )rÞ   rß   rà   rá   rÅ   rÇ   r
   r
   r   Ú
test_refcoñ  s    r‹  c                  C   sÔ   t  dt¡} t  dt¡}d| d d< d| d d< d| d d< d	| d d< d
| d d< d| d d< d| d d< d| d d< d| d d< t | ¡}t|d ddddtƒ t|d ddddtƒ t|d ddddtƒ d S )Nrg  r,   r   r   gš™™™™™é¿r   g333333ã¿r   rž  g
×£p=
×¿g¸…ëQ¸Þ?rN  g{®Gázä¿r-   Úrm2vr   güŽpe•žö?r   gPië1Ç(þ¿r!   )r/   r0   r1   r"   rŒ  r   r#   )r3  Úwr
   r
   r   Ú	test_rm2vû  s    
rŽ  c                  C   s(  t  dt¡} t  dt¡}d| d< d| d< d| d< t | ¡}t|d d d	d
ddtƒ t|d d dd
ddtƒ t|d d dd
ddtƒ t|d d dd
ddtƒ t|d d dd
ddtƒ t|d d dd
ddtƒ t|d d dd
ddtƒ t|d d dd
ddtƒ t|d d dd
ddtƒ d S )Nr,   rg  r   r   gû¿Ud•žö?r   gíeW1Ç(þ¿r   gÿÍçdž æ¿r6   Úrv2mr¼  gOÓvSâ¿r½  gÚd|$'Û¿r¿  gOÓvSâ?rÀ  g/@{?±·¿rÁ  góÜÆò�8ê¿rÂ  gÚd|$'Û?rÄ  rÆ  gÇo «Ø?rÈ  )r/   r0   r1   r"   r�  r   r#   )r�  r3  r
   r
   r   Ú	test_rv2m  s    
r�  c                  C   sv  t  dt¡} d}d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d	| d d< t || ¡} t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d dd
ddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ d S )Nrg  çL´šdšÖ?r.  r   r/  r   r   r0  r1  r   Úrxr¼  r½  r¿  g-Ca\¶@r-   rÀ  g’8°qå	@rÁ  gˆÐü’@rÂ  gØ¡?&€åü?rÄ  gîßàˆ‰¯@rÆ  gèzë9t€@rÈ  )r/   r0   r1   r"   r’  r   r#   )r3  rÊ   r
   r
   r   Útest_rx  s*    r“  c                  C   sú   t  dt¡} t  dt¡}t  dt¡}d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d	| d d< d
|d< d|d< d|d< t | |¡}t|d ddddtƒ t|d ddddtƒ t|d ddddtƒ d S )Nrg  r,   r.  r   r/  r   r   r0  r1  r¶  rM  rø   rW  r-   Úrxpr   ç333333@r   çffffff@r!   )r/   r0   r1   r"   r”  r   r#   )r3  r;   rš  r
   r
   r   Útest_rxp7  s&    r—  c                  C   sˆ  t  dt¡} t  dtj¡}t  dtj¡}d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d	| d d< d
|d d< d|d d< d|d d< d|d d< d
|d d< d|d d< t | |¡}t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ d S )Nrg  r
   r.  r   r/  r   r   r0  r1  r¶  r;   rM  rø   r2   rW  r-   Úrxpvr¼  gffffff@r½  r•  rÀ  çÍÌÌÌÌÌ@rÁ  r–  rÄ  g333333@rÆ  )r/   r0   r1   r"   ru   r˜  r   r#   )r3  rä   Zrpvr
   r
   r   Ú	test_rxpvM  s2    rš  c                  C   sö  t  dt¡} t  dt¡}t  dt¡}d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d	|d d< d|d d< d|d d< d|d d< d	|d d< d	|d d< d|d d< d
|d d< d	|d d< t | |¡}t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ d S )Nrg  r.  r   r/  r   r   r0  r1  rk   r   r`  r-   Úrxrr¼  rG  r½  g      "@r¿  rÀ  r_  rÁ  g      &@rÂ  g      A@rÄ  rà  rÆ  g      .@rÈ  )r/   r0   r1   r"   r›  r   r#   )r   r+  Zatbr
   r
   r   Útest_rxri  s>    rœ  c                  C   sv  t  dt¡} d}d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d	| d d< t || ¡} t| d d d
dddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ d S )Nrg  r‘  r.  r   r/  r   r   r0  r1  ç¸ïçþ—¯ë?r-   Úryr¼  gÏ-‰a¡y÷?r½  gM~8°s È?r¿  r    rÀ  rÁ  rÂ  çóçþl @rÄ  g�T:›@rÆ  gÎÙ’‡@rÈ  )r/   r0   r1   r"   rž  r   r#   )r3  rÈ   r
   r
   r   Útest_ry‹  s*    r   c                  C   sv  t  dt¡} d}d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d	| d d< t || ¡} t| d d d
dddtƒ t| d d ddddtƒ t| d d d
dddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d d	dddtƒ d S )Nrg  r‘  r.  r   r/  r   r   r0  r1  g	XÝM‚/@r-   Úrzr¼  rŸ  r½  r¿  gð
{u¯(@rÀ  r�  rÁ  rÂ  rÄ  rÆ  rÈ  )r/   r0   r1   r"   r¡  r   r#   )r3  ró  r
   r
   r   Útest_rz£  s*    r¢  c                  C   s.   d} d}t  dd| |¡}t|ddddtƒ d S )	NrØ  rÙ  rh  rº  gBÚd83J¾rÆ   Ús00rE   )r"   r£  r   r#   ©r˜   r™   r   r
   r
   r   Útest_s00»  s    r¥  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  r£  gÄž}ÊL¾rÆ   Ús00arE   )r"   r¦  r   r#   ©r   r
   r
   r   Ú	test_s00aÂ  s    r¨  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  r£  gÿG’�ÊL¾rÆ   Ús00brE   )r"   r©  r   r#   r§  r
   r
   r   Ú	test_s00bÇ  s    rª  c                  C   s.   d} d}t  dd| |¡}t|ddddtƒ d S )	NrØ  rÙ  rh  rº  gÕƒ�±23J¾rÆ   Ús06rE   )r"   r«  r   r#   r¤  r
   r
   r   Útest_s06Ì  s    r¬  c                  C   s"   t  dd¡} t| ddddtƒ d S )Nrh  r£  g_ÁlùwÊL¾rÆ   Ús06arE   )r"   r­  r   r#   r§  r
   r
   r   Ú	test_s06aÓ  s    r®  c                  C   s^   t  dt¡} t dd¡} t| d ddddtƒ t| d	 d
dddtƒ t| d ddddtƒ d S )Nr,   g)í¾0@g+‡ÙÎ÷ï¿r   gGä§å,á¿r-   Ús2cr   r   g‡Ùÿ„Ü±?r   r   gL‚œ”æèê¿r!   )r/   r0   r1   r"   r¯  r   r#   )r!  r
   r
   r   Útest_s2cØ  s
    r°  c                  C   s`   t  dt¡} t ddd¡} t| d dddd	tƒ t| d
 ddddtƒ t| d ddddtƒ d S )Nr,   ç®Gáz®	ÀrM  çÉv¾Ÿ/Ý?r   rW  r-   Ús2pr˜   r   ç?[1Â%­Ÿ?r™   r   rX  rR  )r/   r0   r1   r"   r³  r   r#   ©r;   r
   r
   r   Útest_s2pà  s
    r¶  c                  C   sÂ   t  dtj¡} t dddddd¡} t| d d	 d
dddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ t| d d ddddtƒ d S )Nr
   r±  rM  r²  rZ  gb7…ø4åâ>gpP± ƒËé¾r;   r   rW  r-   Ús2pvr˜   r   r´  r™   r   rX  rR  r2   g-~ˆÓë>rd   Zvxg[~à@Æ>ZvygºK"Xª‹Å>Úvz)r/   r0   r"   ru   r·  r   r#   rg  r
   r
   r   Ú	test_s2pvè  s    r¹  c                  C   s  t  dtj¡} t  dtj¡}d}d}d| d d< d| d d< d	| d d
< d| d d< d| d d< d| d d
< t ||| ¡}t|d d ddddtƒ t|d d ddddtƒ t|d d
 ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d
 ddddtƒ d S )Nr
   r.  r/  r  r;   r   r  r   r  r   ru  r2   çffffff@gš™™™™™Ù¿rN  r-   Ús2xpvr'  ç333333@r(  r…  r)  rM  r*  gš™™™™™@r+  g333333ó¿r,  )r/   r0   r"   ru   r»  r   r#   )rä   Úspvrä  rå  r
   r
   r   Ú
test_s2xpvó  s"    r¾  c                  C   sj   t  dt¡} t  dt¡}d| d< d| d< d| d< d|d< d	|d< d|d< t | |¡}t|d
dddtƒ d S )Nr,   rk   r   rø   r   r¶  r   r·  r¸  gø �(â@r-   ÚsepprE   )r/   r0   r1   r"   r¿  r   r#   )r   r+  r   r
   r
   r   Ú	test_sepp  s    rÀ  c                  C   s6   d} d}d}d}t  | |||¡}t|ddddtƒ d S )	Nrk   rø   r¶  r·  g/q^1Æ@r6   ÚsepsrE   )r"   rÁ  r   r#   )r½  r¾  r¿  rÀ  r   r
   r
   r   Ú	test_seps  s    rÂ  c                   C   s   t t dd¡ddddtƒ d S )Nrh  r£  gBû¢VaW›½r-   Úsp00rE   )r   r"   rÃ  r#   r
   r
   r
   r   Ú	test_sp00  s    rÄ  c                  C   sÂ   d} d}d}d}d}d}t  | |||||dddd	¡
\}}}}	}
}}t|d
dddtƒ t|ddddtƒ t|ddddtƒ t|	ddddtƒ t|
ddddtƒ t|ddddtƒ t|dddtƒ d S )Nrm  ro  çQæôÞÕ²ò¾ç ÈÍ’˜Ã>rr  çš™™™™™5Àrh  g    `tè@rº  gæ…é ø‘?rÀ  Ústarpmrü   gÐ‰ü÷â€ñ¿rÔ  g?ÿr©û³ò¾r¹   rp  gùv‰,�Ã>rq  go]š¤Qêç?r  g>“Øs[™5ÀrU   r  r   rA   )r"   rÈ  r   r#   r   )rë  rì  rí  rî  rï  rð  rå  ræ  rç  rè  ré  rê  rA   r
   r
   r   Útest_starpm!  s    *rÉ  c                  C   sî   t  dtj¡} d}d}d}d}d}d}t ||||||¡\} }t| d d	 d
dddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d d	 ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t|d	ddtƒ d S )Nr
   rm  ro  rÅ  rÆ  rr  rÇ  r;   r   ri  rz  Ústarpvr¼  r   rj  r-   r½  r   rk  r¿  r2   g\\NG­˜p¿rÀ  rÀ  gaƒÍµ�y¿rÄ   rÁ  gŸAšˆ£[ˆ?rÂ  rA   )r/   r0   r"   ru   rÊ  r   r#   r   rs  r
   r
   r   Útest_starpv2  s    rË  c                  C   s†   t  dt¡} t  dt¡}d}d| d< d| d< d| d< t || ¡}t|d d	d
ddtƒ t|d dd
ddtƒ t|d dd
ddtƒ d S )Nr,   r.  r  r   r  r   r  r   rN  r   Úsxpr   r¼  r   r…  r!   )r/   r0   r1   r"   rÌ  r   r#   )r;   rÎ   r   r
   r
   r   Útest_sxpD  s    rÍ  c                  C   s  t  dtj¡} t  dtj¡}d}d| d d< d| d d< d| d d	< d
| d d< d| d d< d| d d	< t || ¡}t|d d ddddtƒ t|d d ddddtƒ t|d d	 ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d	 ddddtƒ d S )Nr
   r.  r  r;   r   r  r   r  r   ru  r2   gš™™™™™	@gffffffæ¿rN  r   Úsxpvr'  r¼  r(  r…  r)  rk   r*  gš™™™™™@r+  gffffffö¿r,  )r/   r0   r"   ru   rÎ  r   r#   )rä   r½  r   r
   r
   r   Ú	test_sxpvQ  s     rÏ  c                  C   sJ   t  dd¡\} }}t| ddddtƒ t|ddddtƒ t|d	dd
tƒ d S )Nç   @{¸BAçºXÀ7�ì?rP  ÚtaittÚt1çÑAÂðD’ì?r-   Út2r   rA   )r"   rÒ  r   r#   r   ©rÓ  rÕ  rA   r
   r
   r   Ú
test_taittd  s    r×  c                  C   sL   t  ddd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
ddtƒ d S )NrÐ  rÑ  ç46<U@ÀrP  Útaiut1rT  çgË�ÝŒì?r-   rU  r   rA   )r"   rÙ  r   r#   r   rS  r
   r
   r   Útest_taiut1k  s    rÛ  c                  C   sJ   t  dd¡\} }}t| ddddtƒ t|ddddtƒ t|d	dd
tƒ d S )NrÐ  rÑ  rP  ÚtaiutcrT  g!\{ÁŒì?r-   rU  r   rA   )r"   rÜ  r   r#   r   rS  r
   r
   r   Útest_taiutcr  s    rÝ  c                  C   sJ   t  dd¡\} }}t| ddddtƒ t|ddddtƒ t|d	dd
tƒ d S )NrÐ  gŒ‹ŒÖ�“ì?rP  ÚtcbtdbÚb1g‰¯KïD’ì?r-   Úb2r   rA   )r"   rÞ  r   r#   r   ©rß  rà  rA   r
   r
   r   Útest_tcbtdby  s    râ  c                  C   sJ   t  dd¡\} }}t| ddddtƒ t|ddddtƒ t|d	dd
tƒ d S )NrÐ  gL OqT’ì?rP  ÚtcgttrÓ  g`¾ŸðD’ì?r-   rÕ  r   rA   )r"   rã  r   r#   r   rÖ  r
   r
   r   Ú
test_tcgtt€  s    rä  c                  C   sJ   t  dd¡\} }}t| ddddtƒ t|ddddtƒ t|d	dd
tƒ d S )NrÐ  ça¯ïD’ì?rP  Útdbtcbrß  g=ðÖ�“ì?r-   rà  r   rA   )r"   ræ  r   r#   r   rá  r
   r
   r   Útest_tdbtcb‡  s    rç  c                  C   sL   t  ddd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
ddtƒ d S )NrÐ  rå  ç¯]ÚpX*¿rP  ÚtdbttrÓ  g›ïðD’ì?r-   rÕ  r   rA   )r"   ré  r   r#   r   rÖ  r
   r
   r   Ú
test_tdbttŽ  s    rê  c                  C   s:   t  dddd¡\} }t| ddddtƒ t|d	dd
tƒ d S )Nú+r   é:   g3333334@g#â»�ìÓô?r-   Útf2ar   r   rA   )r"   rí  r   r#   r   rB   r
   r
   r   Ú	test_tf2a•  s    rî  c                  C   s:   t  dddd¡\} }t| ddddtƒ t|d	dd
tƒ d S )Nú r:  é7   gÍÌÌÌÌÌ%@g5š!Ë–äï?r-   Útf2dr9   r   rA   )r"   rñ  r   r#   r   )r9   rA   r
   r
   r   Ú	test_tf2d›  s    rò  c            	      C   s†   d} d}d}d}t  | |||¡\}}}}}t|ddddtƒ t|d	ddd
tƒ t|ddddtƒ t|ddddtƒ t|dddtƒ d S )Nç¸…ëQ¸ž¿çìQ¸…ë±?çÍÌÌÌÌÌô?rM  gvµ!µ3Éû?rÀ  ÚtporsÚaz1gÕg!zßüö?Úbz1g•@#@Úaz2g`ð÷œ•
ù?Úbz2r   Ún)r"   rö  r   r#   r   )	ÚxiÚetarü   rÔ  r÷  rø  rù  rú  rû  r
   r
   r   Ú
test_tpors¡  s    rþ  c                  C   sì   t  dt¡} t  dt¡}t  dt¡}d}d}d}d}t ||¡} t ||| ¡\}}}t|d ddd	d
tƒ t|d ddd	dtƒ t|d ddd	dtƒ t|d ddd	dtƒ t|d ddd	dtƒ t|d ddd	dtƒ t|dd	dtƒ d S )Nr,   ró  rô  rõ  rM  r   g�=À4�—–¿rÄ   ÚtporvÚx1r   go¬›#¥ßÀ?r6   Úy1r   gãc`-¶ï?Zz1gÕ%âƒin¿rd   Zx2g„Å¢ãiÈq¿Úy2gW{ÉÝÿï?Zz2rû  )	r/   r0   r1   r"   r¯  rÿ  r   r#   r   )r2   Zvz1Zvz2rü  rý  rü   rÔ  rû  r
   r
   r   Ú
test_tporv®  s     r  c                  C   sL   d} d}d}d}t  | |||¡\}}t|ddddtƒ t|d	d
ddtƒ d S )Nró  rô  rº  rM  gª6 T¿Nè?r6   Útpstsrü   g1é¤¯a§ø?rÀ  rÔ  )r"   r  r   r#   )rü  rý  ÚrazÚdeczrü   rÔ  r
   r
   r   Ú
test_tpstsÁ  s    r  c                  C   sˆ   t  dt¡} t  dt¡}d}d}d}d}t ||¡} t ||| ¡}t|d ddd	d
tƒ t|d ddd	dtƒ t|d ddd	dtƒ d S )Nr,   ró  rô  rº  rM  r   g¹ògÉš8–?rÄ   Útpstvr˜   r   ghŒ§�•?r™   r   g¯q"¥Tüï?r6   rR  )r/   r0   r1   r"   r¯  r  r   r#   )r¸  r2   rü  rý  r  r  r
   r
   r   Ú
test_tpstvË  s    r	  c                  C   s^   d} d}d}d}t  | |||¡\}}}t|ddddtƒ t|d	ddd
tƒ t|dddtƒ d S )Nrõ  çÍÌÌÌÌÌø?rº  rM  ç"åsCéó‘¿rÄ   Útpxesrü  ç°A¦¾‡®?rý  r   rA   )r"   r  r   r#   r   )rü   rÔ  r  r  rü  rý  rA   r
   r
   r   Ú
test_tpxesÙ  s    r  c            	      C   sŠ   t  dt¡} t  dt¡}d}d}d}d}t ||¡} t ||¡}t | |¡\}}}t|dddd	tƒ t|d
dddtƒ t|dddtƒ d S )Nr,   rõ  r
  rº  rM  r  rÄ   Útpxevrü  r  rý  r   rA   )	r/   r0   r1   r"   r¯  r  r   r#   r   )	r2   r¸  rü   rÔ  r  r  rü  rý  rA   r
   r
   r   Ú
test_tpxevä  s    r  c                  C   s|  t  dt¡} t  dt¡}d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< t | ¡}t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ t|d d dd	d
dtƒ d S )Nrg  r.  r   r/  r   r   r0  r1  r   Útrr¼  r½  r¿  rÀ  rÁ  rÂ  rÄ  rÆ  rÈ  )r/   r0   r1   r"   r  r   r#   )r3  Úrtr
   r
   r   Útest_tró  s*    
r  c                  C   sú   t  dt¡} t  dt¡}t  dt¡}d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d	| d d< d
|d< d|d< d|d< t | |¡}t|d ddddtƒ t|d ddddtƒ t|d ddddtƒ d S )Nrg  r,   r.  r   r/  r   r   r0  r1  r¶  rM  rø   r™  r-   Útrxpr   r   çš™™™™™@r!   )r/   r0   r1   r"   r  r   r#   )r3  r;   Ztrpr
   r
   r   Ú	test_trxp  s&    r  c                  C   sˆ  t  dt¡} t  dtj¡}t  dtj¡}d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d	| d d< d
|d d< d|d d< d|d d< d|d d< d
|d d< d|d d< t | |¡}t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ t|d d ddddtƒ d S )Nrg  r
   r.  r   r/  r   r   r0  r1  r¶  r;   rM  rø   r2   r™  r-   Útrxpvr'  r  r•  r*  g333333@r+  gffffff@r,  )r/   r0   r1   r"   ru   r  r   r#   )r3  rä   Ztrpvr
   r
   r   Ú
test_trxpv!  s2    r  c                  C   sJ   t  dd¡\} }}t| ddddtƒ t|ddddtƒ t|d	dd
tƒ d S )NrÐ  rÑ  rP  Útttair{  g£ëí�*Œì?r-   r}  r   rA   )r"   r  r   r#   r   ©r{  r}  rA   r
   r
   r   Ú
test_tttai=  s    r  c                  C   sJ   t  dd¡\} }}t| ddddtƒ t|ddddtƒ t|d	dd
tƒ d S )NrÐ  rÑ  rP  ÚtttcgÚg1gƒOAG�ì?r-   Úg2r   rA   )r"   r  r   r#   r   )r  r  rA   r
   r
   r   Ú
test_tttcgD  s    r  c                  C   sL   t  ddd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
ddtƒ d S )NrÐ  rÔ  rè  rP  Útttdbrß  gI…‚ïD’ì?r-   rà  r   rA   )r"   r   r   r#   r   rá  r
   r
   r   Ú
test_tttdbK  s    r!  c                  C   sL   t  ddd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
ddtƒ d S )NrÐ  rÔ  çµ7øÂd6P@rP  Úttut1rT  rÚ  r-   rU  r   rA   )r"   r#  r   r#   r   rS  r
   r
   r   Ú
test_ttut1R  s    r$  c                  C   sL   t  ddd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
ddtƒ d S )NrÐ  ç pCÝŒì?rØ  rP  Úut1tair{  gr»À7�ì?r-   r}  r   rA   )r"   r&  r   r#   r   r  r
   r
   r   Útest_ut1taiY  s    r'  c                  C   sL   t  ddd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
ddtƒ d S )NrÐ  r%  r"  rP  Úut1ttrÓ  g‰æƒðD’ì?r-   rÕ  r   rA   )r"   r(  r   r#   r   rÖ  r
   r
   r   Ú
test_ut1tt`  s    r)  c                  C   sL   t  ddd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
ddtƒ d S )NrÐ  r%  ç‹ýe÷äaÕ?rP  Úut1utcrT  gÚ =ÁŒì?r-   rU  r   rA   )r"   r+  r   r#   r   rS  r
   r
   r   Útest_ut1utcg  s    r,  c                  C   sJ   t  dd¡\} }}t| ddddtƒ t|ddddtƒ t|d	dd
tƒ d S )NrÐ  çB/ÁŒì?rP  ÚutctairT  gÚ»À7�ì?r-   rU  r   rA   )r"   r.  r   r#   r   rS  r
   r
   r   Útest_utctain  s    r/  c                  C   sL   t  ddd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
ddtƒ d S )NrÐ  r-  r*  rP  Úutcut1rT  gˆp5ÝŒì?r-   rU  r   rA   )r"   r0  r   r#   r   rS  r
   r
   r   Útest_utcut1u  s    r1  c                  C   s8   t  dd¡\} }t| ddddtƒ t|dddd	tƒ d S )
Nrh  rº  rØ  rÄ   Úxy06r˜   gjB¶´T?rd   r™   )r"   r2  r   r#   )r˜   r™   r
   r
   r   Ú	test_xy06|  s    r3  c                  C   sL   t  dd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
dddtƒ d S )Nrh  rº  rÃ  r6   Úxys00ar˜   rÅ  rÄ   r™   gs!ØØ>3J¾rÆ   r   )r"   r4  r   r#   r¤  r
   r
   r   Útest_xys00a‚  s    r5  c                  C   sL   t  dd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
dddtƒ d S )Nrh  rº  rÍ  r6   Úxys00br˜   rÎ  rÄ   r™   g.|Vã+3J¾rÆ   r   )r"   r6  r   r#   r¤  r
   r
   r   Útest_xys00b‰  s    r7  c                  C   sL   t  dd¡\} }}t| ddddtƒ t|dddd	tƒ t|d
dddtƒ d S )Nrh  rº  rÒ  r6   Úxys06ar˜   rÓ  rÄ   r™   gÜ’ÔÎ23J¾rÆ   r   )r"   r8  r   r#   r¤  r
   r
   r   Útest_xys06a�  s    r9  c                  C   sr   t  dt¡} d| d< d| d< d| d< t ¡ } t| d ddd	d
tƒ t| d ddd	dtƒ t| d ddd	dtƒ d S )Nr,   r  r   r  r   r  r   r   Úzpr   r   r!   )r/   r0   r1   r"   r:  r   r#   rµ  r
   r
   r   Útest_zp—  s    r;  z1.21zneeds numpy >= 1.21c                  C   sþ   t  dtj¡} d| d d< d| d d< d| d d< d	| d
 d< d| d
 d< d| d
 d< t ¡ } t| d d ddddtƒ t| d d ddddtƒ t| d d ddddtƒ t| d
 d ddddtƒ t| d
 d ddddtƒ t| d
 d ddddtƒ d S )Nr
   r  r;   r   r  r   r  r   r#  r2   r$  r%  r   Úzpvr'  r(  r)  r*  r+  r,  )r/   r0   r"   ru   r<  r   r#   rg  r
   r
   r   Útest_zpv£  s    r=  c                  C   sn  t  dt¡} d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< d| d d< t ¡ } t| d d d	d	d
dtƒ t| d d d	d	d
dtƒ t| d d d	d	d
dtƒ t| d d d	d	d
dtƒ t| d d d	d	d
dtƒ t| d d d	d	d
dtƒ t| d d d	d	d
dtƒ t| d d d	d	d
dtƒ t| d d d	d	d
dtƒ d S )Nrg  r.  r   r/  r   r   r0  r1  r   ÚzrZ00Z01Z02Ú10r¼  r½  Z20rÀ  rÁ  )r/   r0   r1   r"   r>  r   r#   r7  r
   r
   r   Útest_zrµ  s(    r@  c                   C   s   t dkstdƒ‚d S )Nr   zSanity check failed!)r#   r   r
   r
   r
   r   Útest_statusÍ  s    rA  (  Ú__doc__Únumpyr/   Zpytestr€  r   r"   Ú__version__Ú
startswithÚversionZerfa_versionÚskipr   r   rê   Úintr#   r$   r+   r4   r<   rC   rF   rH   r€   rŒ   rš   rž   rÏ   râ   rå   rç   rì   rî   rñ   ró   r  r  r  r  r,  r.  r5  r;  r=  r>  rE  rG  rS  r^  r`  ÚmarkZxfailrf  rœ  r­  r¹  rË  rÐ  rÕ  r×  rß  râ  ré  rô  rö  rø  r  r  r  r  r  r"  r-  r4  rB  rE  rL  rO  rV  r\  rh  rl  ro  rq  rs  rv  r„  r†  r‰  r‹  r�  r‘  r“  r•  r˜  rš  r�  r   r¢  r¤  r¦  r¨  rª  r¬  r®  r°  r²  r´  r¶  r¸  rº  rÉ  rÌ  rÎ  rà  râ  ræ  ré  rõ  r÷  rþ  r
  r  r  r  r  r  r  r  r  r!  r#  r%  r)  r+  r.  r1  r5  r8  r?  rA  rF  rJ  rL  rP  rR  rU  rX  rZ  r\  r^  rh  ru  r€  rŠ  r�  r‘  r•  r™  r›  r�  rŸ  r¢  r±  r³  rµ  rº  rÁ  rÆ  rÊ  rÍ  rÏ  rÒ  rÕ  r×  rÙ  rÜ  rã  rñ  rô  r   r"  r$  r8  r:  r=  r?  rA  rC  rE  rI  rK  rO  rS  rV  r]  rb  rd  rf  rh  rt  rw  r�  rƒ  r‡  r‰  r‹  rŽ  r�  r“  r—  rš  rœ  r   r¢  r¥  r¨  rª  r¬  r®  r°  r¶  r¹  r¾  rÀ  rÂ  rÄ  rÉ  rË  rÍ  rÏ  r×  rÛ  rÝ  râ  rä  rç  rê  rî  rò  rþ  r  r  r	  r  r  r  r  r  r  r  r!  r$  r'  r)  r,  r/  r1  r3  r5  r7  r9  r;  r=  r@  rA  r
   r
   r
   r   Ú<module>   s  ÿ			'*<1.#	
($!!
""
##	

	

	
)						99999
"	
