Ë
    âQ(hù  ã                   óP   — d dl Z d dlZd dlmZ d dlmZ d„ Zd„ Zd„ Z	d„ Z
d„ Zd	„ Zy)
é    N)Úspecial)Úprimes_from_2_toc                 ó>   — t        t        j                  | «      «      S ©N)r   ÚmathÚceil)Úns    úY/var/www/skyplay_api_hub/venv/lib/python3.12/site-packages/scipy/stats/tests/data/_mvt.pyÚ_primesr      s   € ô œDŸI™I a›LÓ)Ð)ó    c                 ó.   — t        j                  || «      S r   )r   Úgammaincinv)ÚaÚbs     r
   Ú_gaminvr      s   € ô ×Ñ˜q !Ó$Ð$r   c                 óx  — t        dt        j                  |«      «      }t        |||z  ||z  «      \  }}}	t	        |«      }
d}t        j
                  | |z  «      }t        j                  |«      }d}d}t        j                  t        d|
z  t        j                  |
dz   «      z  dz  «      «      }|dd…t        j                  f   }d}d}t        |«      D �]Ç  }|j                  «       }t        j                  |
|f«      }t        |
«      D �]W  }t        j                  dt        j                  ||   t        j                   d|dz   «      z  |j#                  «       z   d«      z  dz
  «      }|dk(  r-|}|dkD  rTt        j                  dt%        ||dz  «      z  «      }n.t'        |||z  z   «      }||dxxx ||d…|dz
  |…f   |z  z  ccc ||dd…f   }|j                  «       }||   z  |z
  }|j                  «       }|	|   |z  |z
  }d||dk  <   t        |«      d	k  }t)        ||   «      ||<   d||dk  <   t        |«      d	k  }t)        ||   «      ||<   ||z
  }||z  }�ŒZ t        j*                  |«      |z
  |dz   z  }||z   }|dz
  |z  |dz   z  |dz  z   }�ŒÊ t        j                  |«      }||fS )
a  Estimates the multivariate t CDF using randomized QMC

    Parameters
    ----------
    m : int
        The number of points
    nu : float
        Degrees of freedom
    sigma : ndarray
        A 2D positive semidefinite covariance matrix
    a : ndarray
        Lower integration limits
    b : ndarray
        Upper integration limits.
    rng : Generator
        Pseudorandom number generator

    Returns
    -------
    p : float
        The estimated CDF.
    e : float
        An absolute error estimate.

    é   é
   r   é   é   Né   i÷ÿÿÿé	   )Úmaxr   ÚsqrtÚ_chlrpsÚlenr   ÚnpÚonesr   ÚlogÚnewaxisÚrangeÚcopyÚzerosÚabsÚmodÚarangeÚrandomr   Ú_PhinvÚ_PhiÚmean) ÚmÚnuÚsigmar   r   ÚrngÚsnÚchÚazÚbzr	   ÚNÚPÚonÚpÚeÚpsÚqÚcÚdcÚSÚvpÚsÚiÚxÚrÚyÚsiÚaiÚdÚbiÚtls                                    r
   Ú_qsimvtvrH      s§  € ôP 
ˆQ”—	‘	˜"“Ó	€B¬W°U¸A¸b¹DÀ!ÀBÁ$Ó-G¡
  B¨ÜˆE‹
€A˜�A¤§	¡	¨!¨A©#£˜A´R·W±W¸Q³Z°ÀQÀÈAÈÜ	�‰”˜˜1™œTŸX™X a¨¡c›]Ñ*¨1Ñ,Ó-Ó	.€B°B²q¼"¿*¹*°}Ñ4E°ð 	€A�4ˆbÜ�1‹Xó OˆØ�W‰W‹YˆœBŸH™H a¨ VÓ,˜Ü�q“ó 	%ˆAÜ—‘�qœŸ™  !¡¤R§Y¡Y¨q°!°A±#Ó%6Ñ 6¸¿¹»Ñ EÀqÓIÑIÈ!ÑKÓLˆAØ�AŠvØ�Ø˜’6ÜŸ™ ¤'¨!¨R°©TÓ"2Ñ 2Ó3‘Aä˜1˜q ™t™8Ó$�Ø�!�"“˜˜A™B  !¡ A ˜I™¨Ñ*Ñ*“Ø�1’a�4‘ˆB˜bŸg™g›i˜!¨b°©e°A©g¸©l¨ÀÇÁÃ	¸AÐPRÐSTÑPUÐVWÑPWÐZ\ÑP\È2ØˆAˆb�B‰h‰K¤# b£'¨A¡+˜R´t¸B¸r¹F³|¨q°©uØˆAˆb�B‰h‰K¤# b£'¨A¡+˜R´t¸B¸r¹F³|¨q°©uØ�Q‘ˆB˜R "™Wšð	%ô �W‰W�R‹[˜1‰_˜q 1™uÑ%ˆ¨1¨q©5 q°q¸1±u¸a±iÀÀQÁÑ6GÈ!ÈQÉ$Ñ6N²!ðOô  	�	‰	�!‹€AØˆaˆ4€Kr   c                 ó,   — t        j                  | «      S r   )r   Úndtr)Úzs    r
   r)   r)   v   s   € Ü�<‰<˜‹?Ðr   c                 ó,   — t        j                  | «      S r   )r   Úndtri)r6   s    r
   r(   r(   z   s   € Ü�=‰=˜ÓÐr   c           	      ó   — d}t        j                  | j                  «      j                  }t	        | «      }| j                  «       }|j                  «       }|j                  «       }t        j                  t        j                  t        j                  |«      d«      «      }	t        |«      D ]S  }
|	|
   dkD  sŒ|dd…|
fxx   |	|
   z  cc<   ||
dd…fxx   |	|
   z  cc<   ||
xx   |	|
   z  cc<   ||
xx   |	|
   z  cc<   ŒU t        j                  |df«      }t        j                  dt        j                  z  «      }t        |«      D �]»  }|}d}d}d}t        ||«      D ]‚  }
||
|
f   |kD  sŒt        j                  t        ||
|
f   d«      «      }|
dkD  r||
d|…f   |d| z  }||
   |z
  |z  }||
   |z
  |z  }t        |«      t        |«      z
  }||k  sŒy|}|}|}|}|
}Œ„ ||kD  rà|||g   |||g<   |||g   |||g<   |||f   |||f<   ||d|…f   j                  «       }||d|…f   ||d|…f<   |||d|…f<   ||dz   d…|f   j                  «       }||dz   d…|f   ||dz   d…|f<   |||dz   d…|f<   ||dz   |…|f   j                  «       }|||dz   |…f   j                  ||dz   |…|f<   |j                  |||dz   |…f<   |||dz   z  kD  �r||||f<   d|||dz   d…f<   t        |dz   |«      D ]O  }
||
|f   |z  ||
|f<   ||
|dz   |
dz   …f   ||
|f   ||dz   |
dz   …|f   j                  z  z
  ||
|dz   |
dz   …f<   ŒQ t!        |«      |kD  rCt        j"                  dz   dz  «      t        j"                  dz   dz  «      z
  ||z  z  ||<   n z   dz  ||<   |dk  r|||<   n
|dkD  r|||<   ||d|dz   …fxx   |z  cc<   ||xx   |z  cc<   ||xx   |z  cc<   nd||d…|f<   ||   ||   z   dz  ||<   �Œ¾ |||fS )z�
    Computes permuted and scaled lower Cholesky factor c for R which may be
    singular, also permuting and scaling integration limit vectors a and b.
    g»½×Ùß|Û=r   Nr   r   iöÿÿÿr   )r   ÚfinfoÚdtypeÚepsr   r"   r   ÚmaximumÚdiagr!   r#   r   Úpir   r)   ÚTr$   Úexp)ÚRr   r   ÚeprQ   r	   r:   ÚapÚbprE   r?   rB   ÚsqtpÚkÚimÚckkÚdemr>   ÚciirD   rF   ÚdeÚamÚbmÚts                            r
   r   r   ~   s¶  € ð
 
€BÜ
�(‰(�1—7‘7Ó
×
Ñ
€CäˆA‹€A�A—F‘F“H� 1§6¡6£8˜b°!·&±&³(¨RÄÇÁÌÏ
É
ÔSU×SZÑSZÐ[\ÓS]Ð_`ÓHaÓ@b¸AÜ�1‹Xò )ˆØˆQ‰4�!‹8ØŠa�ˆd‹G�q˜‘t‰O‹G˜Q˜q¢!˜t›W¨¨!©™_›WØˆq‹E�Q�q‘T‰M‹E˜2˜a›5 A a¡D™=œ5ð)ô 	�‰�!�Q�Ó€A¤§¡¨1¬T¯W©W©9Ó!5˜$ä�1‹Xó ˆØˆ�a�˜q˜ a !Ü�q˜!“ò 	BˆAØ��A�‰w˜‹}Ü—i‘i¤ A a¨ d¡G¨Q£Ó0�Ø�q’5˜a  2 A 2 ™h¨¨2¨A¨Ñ.˜!Ø˜‘e˜A‘g˜s‘]�¨"¨Q©%°©'°3© B¼TÀ"»XÄdÈ2ÃhÑ=N¸Ø˜“9Ø�C R˜s¨b¨°r°"À¹2ð	Bð �Š6Ø˜a ˜W™+ˆB��Aˆw‰K°R¸¸B¸±[ r¨2¨q¨'¡{ÈaÐPQÐSTÐPTÉgÀ!ÀBÈÀFÁ)Ø�"�b�q�b�&‘	—‘Ó ˆA¨a°°2°A°2°©h ! B¨¨¨ F¡)À1¸¸!¸R¸a¸R¸%¹Ø�"�Q‘$‘%˜�)‘×!Ñ!Ó#ˆA°A°b¸±d±e¸Q°h±K Q r¨!¡t¡u¨b y¡\ÈqÀÀ2ÀaÁ4Á5È!À8ÁØ�!�A‘#�b�&˜!�)‘×!Ñ!Ó#ˆA°A°b¸!¸A¹#¸b¸&°j±M·O±O Q q¨¡s¨2 v¨q y¡\ÐUV×UXÑUXÀQÀrÈ1ÈQÉ3ÈrÈ6ÀzÁ]Ø��Q�q‘S‘‹>ØˆAˆa�ˆd‰G¨˜1˜Q  !¡¡˜W™:Ü˜1˜Q™3 “]ò _�Ø˜A˜q˜D™' #™+��!�Q�$‘°q¸¸A¸a¹CÀÀ!Á¸G¸±}ÀqÈÈAÈÁwÈqÐQRÐSTÑQTÐUVÐWXÑUXÐQXÐZ[ÐQ[É}ÏÉÑG^Ñ7^ q¨¨A¨a©C°°!±¨G¨¢}ð_ä�3‹x˜"Š}ÜŸ™  A¡˜v a™xÓ(¬2¯6©6°2°q±5°&¸±(Ó+;Ñ;ÀÀSÁÑI��!’à˜R™ 1‘}��!‘Ø˜’8Ø�A�a’DØ˜"’WØ�A�a‘DØˆa��!�A‘#�ˆg‹J˜#Ñ‹J˜r !›u¨™|›u¨R°«U°c©\¬UàˆAˆa‰b�!ˆe‰H " Q¡%¨"¨Q©%¡-°Ñ!2˜!˜A™$Ùð;ð< ˆb�"ˆ9Ðr   )r   Únumpyr   Úscipyr   Úscipy.stats._qmcr   r   r   rH   r)   r(   r   © r   r
   ú<module>ri      s1   ðÛ Û Ý Ý -ò*ò%ò_òFòó-r   