§
    fŠtjù  ã                   óR   — 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dS )
é    N)Úspecial)Úprimes_from_2_toc                 óD   — t          t          j        | ¦  «        ¦  «        S ©N)r   ÚmathÚceil)Úns    úY/var/www/html/CA-Chatbot/venv/lib/python3.11/site-packages/scipy/stats/tests/data/_mvt.pyÚ_primesr      s   € õ �DœI a™LœLÑ)Ô)Ð)ó    c                 ó,   — t          j        || ¦  «        S r   )r   Úgammaincinv)ÚaÚbs     r
   Ú_gaminvr      s   € õ Ô˜q !Ñ$Ô$Ð$r   c                 ó  — 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 �]
}|                     ¦   «         }t          j        |
|f¦  «        }t          |
¦  «        D �]–}t          j        dt          j        ||         t          j        d|dz   ¦  «        z  |                     ¦   «         z   d¦  «        z  dz
  ¦  «        }|dk    r1|}|dk    r(t          j        dt%          ||dz  ¦  «        z  ¦  «        }n9t'          |||z  z   ¦  «        }||d…xx         ||d…|dz
  |…f         |z  z  cc<   ||dd…f         }|                     ¦   «         }||         |z  |z
  }|                     ¦   «         }|	|         |z  |z
  }d||dk    <   t          |¦  «        d	k     }t)          ||         ¦  «        ||<   d||dk    <   t          |¦  «        d	k     }t)          ||         ¦  «        ||<   ||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Ô-G¡
  B¨ÝˆE‰
Œ
€A˜�A¥¤	¨!¨A©#¡¤˜AµR´W¸Q±Z´Z°ÀQÀÈAÈÝ	Œ•˜˜1™�TœX a¨¡c™]œ]Ñ*¨1Ñ,Ñ-Ô-Ñ	.Ô	.€B°B°q°q°q½"¼*°}Ô4E°ð 	€A�4ˆbÝ�1‰XŒXð Oñ OˆØ�WŠW‰YŒYˆ�BœH a¨ VÑ,Ô,˜Ý�q‘”ð 	%ñ 	%ˆAÝ”�q�œ  !¤¥R¤Y¨q°!°A±#Ñ%6Ô%6Ñ 6¸¿º¹¼Ñ EÀqÑIÔIÑIÈ!ÑKÑLÔLˆAØ�AŠvˆvØ�Ø˜’6�6Ýœ ¥'¨!¨R°©TÑ"2Ô"2Ñ 2Ñ3Ô3�Aøå˜1˜q ™t™8Ñ$Ô$�Ø�!�"�"��”˜˜A˜B˜B  !¡ A ˜Iœ¨Ñ*Ñ*��‘Ø�1�a�a�a�4”ˆB˜bŸgšg™iœ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�R‰[Œ[˜1‰_˜q 1™uÑ%ˆ¨1¨q©5 q°q¸1±u¸a±iÀÀQÁÑ6GÈ!ÈQÉ$Ñ6N°!±!ÝŒ	�!‰Œ€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	          | ¦  «        }|                      ¦   «         }|                     ¦   «         }|                     ¦   «         }t          j        t          j        t          j        |¦  «        d¦  «        ¦  «        }	t          |¦  «        D ]n}
|	|
         dk    r`|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<   Œot          j
        |df¦  «        }t          j        dt          j        z  ¦  «        }t          |¦  «        D �]}|}d}d}d}t          ||¦  «        D ]¢}
||
|
f         |k    r’t          j        t          ||
|
f         d¦  «        ¦  «        }|
dk    r||
d|…f         |d|…         z  }||
         |z
  |z  }||
         |z
  |z  }t          |¦  «        t          |¦  «        z
  }||k    r
|}|}|}|}|
}Œ£||k    rý|||g         |||g<   |||g         |||g<   |||f         |||f<   ||d|…f                              ¦   «         }||d|…f         ||d|…f<   |||d|…f<   ||dz   d…|f                              ¦   «         }||dz   d…|f         ||dz   d…|f<   |||dz   d…|f<   ||dz   |…|f                              ¦   «         }|||dz   |…f         j        ||dz   |…|f<   |j        |||dz   |…f<   |||dz   z  k    �r+||||f<   d|||dz   d…f<   t          |dz   |¦  «        D ]V}
||
|f         |z  ||
|f<   ||
|dz   |
dz   …f         ||
|f         ||dz   |
dz   …|f         j        z  z
  ||
|dz   |
dz   …f<   ŒWt!          |¦  «        |k    rAt          j        |dz   dz  ¦  «        t          j        |dz   dz  ¦  «        z
  ||z  z  ||<   n"||z   dz  ||<   |dk     r|||<   n|dk    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   ~   s5  € ð
 
€BÝ
Œ(�1”7Ñ
Ô
Ô
€CåˆA‰Œ€A�A—F’F‘H”H� 1§6¢6¡8¤8˜b°!·&²&±(´(¨RÅÄÍÌ
ÕSUÔSZÐ[\ÑS]ÔS]Ð_`ÑHaÔHaÑ@bÔ@b¸AÝ�1‰XŒXð )ð )ˆØˆQŒ4�!Š8ˆ8Øˆaˆaˆa�ˆdˆGˆGŒG�q˜”t‰OˆGˆG‰G˜Q˜q ! ! !˜t˜W˜WœW¨¨!¬™_˜W˜W™WØˆqˆEˆEŒE�Q�q”T‰MˆEˆE‰E˜2˜a˜5˜5œ5 A a¤D™=˜5˜5™5øÝ
Œ�!�Q�ÑÔ€A¥¤¨1­T¬W©9Ñ!5Ô!5˜$å�1‰XŒXð ñ ˆØˆ�a�˜q˜ a !Ý�q˜!‘”ð 	Bð 	BˆAØ��A�Œw˜Š}ˆ}Ý”i¥ A a¨ d¤G¨Q¡¤Ñ0Ô0�Ø�q’5�5˜a  2 A 2 œh¨¨2¨A¨2¬Ñ.˜!Ø˜”e˜A‘g˜s‘]�¨"¨Q¬%°©'°3© B½TÀ"¹X¼XÅdÈ2ÁhÄhÑ=N¸Ø˜’9�9Ø�C R˜s¨b¨°r°"À¸2øØ�Š6ˆ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°e¸Q°h´K Q r¨!¡t u u¨b y¡\ÈqÀÀ2ÀaÁ4À5À5È!À8ÁØ�!�A‘#�b�&˜!�)”×!Ò!Ñ#Ô#ˆA°A°b¸!¸A¹#¸b¸&°j´M´O Q q¨¡s¨2 v¨q y¡\ÐUVÔ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Œx˜"Š}ˆ}Ýœ  A¡˜v a™xÑ(Ô(­2¬6°2°q±5°&¸±(Ñ+;Ô+;Ñ;ÀÀSÁÑI��!‘�à˜R™ 1‘}��!‘Ø˜’8�8Ø�A�a‘D�DØ˜"’W�WØ�A�a‘DØˆa��!�A‘#�ˆgˆJˆJŒJ˜#ÑˆJˆJ‰J˜r !˜u˜uœu¨™|˜u˜u™u¨R°¨U¨U¬U°c©\¨U¨U©U¨UàˆAˆaˆbˆ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      s¦   ðØ €€€Ø Ð Ð Ð Ø Ð Ð Ð Ð Ð Ø -Ð -Ð -Ð -Ð -Ð -ð*ð *ð *ð%ð %ð %ð_ð _ð _ðFð ð ðð ð ð-ð -ð -ð -ð -r   