0PhL UddlZddlZddlmZmZddlZddlZddlm Z ddl m Z m Z ddlmZddlmZmZmZmZmZddlmZmZmZmZmZddlmZegZeeed<egZ eeed <ee zZ!egZ"eeed <egZ#eeed <e"e#zZ$ d=d Z%dZ&d\Z'Z(e)e'e(zdz Z*ej+,dddgddgfddgddgfddgddgfgdZ-dZ.dZ/dZ0d Z1d!Z2ej+,d"e!d#Z3ej+,d"ed$Z4d%Z5d&Z6d'Z7ej+,d(ed)Z8ej+,d(ed*Z9ej+,d(ed+Z:ej+,d(ed,Z;ej+,d(ed-Zd0Z?ej+,d(eej+,d1e$d2Z@ej+,d(eej+,d3d4ej+,d5d4ej+,d1e$ej+,d6d7d8gd9ZAej+,d1e$ej+,d:ejBejBfejCejCfejDejCfejEejCffd;ZFej+,d1e$d<ZGdS)>N)AnyList)DataDimensionalityWarningNotFittedError)euclidean_distances)GaussianRandomProjectionSparseRandomProjection_gaussian_random_matrix_sparse_random_matrixjohnson_lindenstrauss_min_dim)assert_allcloseassert_allclose_dense_sparseassert_almost_equalassert_array_almost_equalassert_array_equal)COO_CONTAINERSall_sparse_random_matrixall_dense_random_matrixall_SparseRandomProjectionall_DenseRandomProjectioncsrc6tj|}|||||||||ff||f}|||S|S)zMake some random data with uniformly located non zero entries with Gaussian distributed values; `sparse_format` can be `"csr"` (default) or `None` (in which case a dense array is returned). sizeshape)nprandom RandomStaterandnrandintasformattoarray) coo_container n_samples n_features n_nonzeros random_state sparse_formatrngdata_coos d/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/sklearn/tests/test_random_projection.pymake_sparse_random_datar-$s )   - -C} IIj ! ! IJ 77 JZ 88  *%   H   ///!!!cVtj|s|S|SN)spissparser#)matrixs r,densifyr4As) ;v    ~~r.) gY@zn_samples, epsdn?皙?Z皙?2i皙?ctjt5t||ddddS#1swxYwYdS)Neps)pytestraises ValueErrorr )r%rBs r,test_invalid_jl_domainrFQs z " "::%iS9999::::::::::::::::::s 9==c"tjt5tddgzddgzdddn #1swxYwYttjddd t jdd dS) Nr7r9rAr5)r5r5r?)rCrDrEr rrr!fullr.r,test_input_size_jl_min_dimrN^s z " "@@%a3%iQ#Y????@@@@@@@@@@@@@@@" !Rh//RWXs5K5KsAAAcgd}|D]B\}}tjt5|||dddn #1swxYwYCdS)N))rr)rJ)rJrP)rJr)rPrrCrDrE random_matrixinputs n_componentsr&s r,check_input_size_random_matrixrVjs 8 8 8F$*44 j ]: & & 4 4 M, 3 3 3 4 4 4 4 4 4 4 4 4 4 4 4 4 4 444s >A A cPgd}|D]\}}|||j||fksJdS)N))rJ)rXrJ)rXrX)rJrJrrRs r,check_size_generatedrYqsb - - -F$*   j}\:66<  A        r.ct|ddd}tdtj|dtdtj|ddS)Ni'rJrr(rH?)r4rrmeanlinalgnorm)rSAs r,check_zero_mean_and_unit_normrazsb  eQQ77788AaQ///c29>>!#4#4a88888r.cd\}}dD]A}tjt5||||dddn #1swxYwYBdS)N)rXr5)gr=r:densityrQ)rSrUr&rds r,%check_input_with_sparse_random_matrixres$L*#EE ]: & & E E M, G D D D D E E E E E E E E E E E E E E EEEs>A A rSc`t|t|t|dSr0)rVrYra)rSs r,$test_basic_property_of_random_matrixrgs4#=111'''!-00000r.cnt|tj|d}t|dS)Nr\rc)re functoolspartialra)rSrandom_matrix_denses r,+test_basic_property_of_sparse_random_matrixrls;)-888#+M3GGG!"566666r.cd}d}t||d}tdtj|dttj|dd|z ddS) Nr7r6rr[r=rIrJddof)r rrr]var)rUr&r`s r,test_gaussian_random_matrixrqsh LJ jqIIIAc271::q111bfQQ///\1A1EEEEEr.c vd}d}dD]}d|z }t|||d}t|}tj|}tj|tj|z |vsJtj| tj|z |vsJ|dkrtj|dksJn d |vsJtj|d ksJt tj|d kdd|z z d t tj|tj|tj|z kdd|zz d t tj|tj| tj|z kdd|zz d t tj|d kd dd|z z dz|z d t tj|tj|tj|z kd ddd|zz z dzd|zz d t tj|tj| tj|z kd ddd|zz z dzd|zz d dS) Nr7i)g333333?r\rJr)rdr(r\rIr=rH)decimalrn) r r4runiquesqrtrrr]rp)rUr&rdsr`valuess r,test_sparse_random_matrixrxsLJ, ,  K ! *gA    AJJ1wqzzBGL111V;;;; {RW\222f<<<< c>>76??a'''''&====76??a'''' BGAH--q1q5y!DDDD GAbgl&;&;;; < * * B B  & 1 1 1 ; ;D A A A B B B B B B B B B B B B B B BBBs A33A7 :A7 ct|ddd|d}tD]Y}|dd}d}tjt|5||dddn #1swxYwYZdS) Nr6r7r%r&r'r(r)rzr<)rUrBz~eps=0.100000 and n_samples=1000 lead to a target dimension of 5920 which is larger than the original space with n_features=100)match)r-r|rCrDrEr})r$rrrrp expected_msgs r,.test_too_many_samples_to_find_a_safe_embeddingrs "'    D1  6s ; ; ;   ]:\ : : :   FF4LLL               s A++A/ 2A/ ct|ddddd}d}t|d}|}|d k}||}tD]}|d |d }||}t|d}|}||}||z } | d |zksJd |z | ksJdS) Nii:rrr?T)squaredr=rz)rUrBr(rJ)r-rravelr| fit_transformmaxmin) r$rrBoriginal_distances non_identicalrr projectedprojected_distancesdistances_ratios r,(test_random_projection_embedding_qualityrs3 "    D C,T4@@@+1133&#-M,M:0//  6s K K K$$T** 1)TJJJ177992-@-0BB ""$$q3w....3w,,.......!//r.ct|tttdd}t|tttdd}tD]}|ddd}||t ||tj sJt ||tj sJ|ddd}||}t ||tj sJtj ||sJdS)Nrrrr5T)rU dense_outputr(F) r-r%r&r'rr} isinstancerrndarrayr1r2)r$ dense_data sparse_dataSparseRandomProjrs r,+test_SparseRandomProj_output_representationr7sU( J* K766 2Dq Q Q Q z",,z22BJ?????",,{33RZ@@@@@ 2EPQ R R R VVJ  ",,z22BJ?????{2<< 4455555566r.ct|ttt|d}tD]}|ddd|}|jdksJ|jdksJ|tvr#|j dksJt|j dd|j j dtfksJ||}|j tdfksJ||}t|||dd }||}t||t#jt&5||ddd d fdddn #1swxYwY|tvrt|d d d}||} | j td fksJ|j j d tfksJ|j jdksJd|j jksJdS)NrrzrrK)rUr(rBr8gQ?rI)r(rBrJrXr7gMbP?)rUrdr(sU)r-r%r&r'r|r}rU n_components_rrdrdensity_ components_rrrrrCrDrEnnz) r$rrrr projected_1 projected_2rp2 projected_3rs r,2test_correct_RandomProjection_dimensions_embeddingr[s #'    D1&+&+  6s K K K O OPT U U&((((3&&&& 9 9 9:''''  T1 5 5 5~#Z'88888ll4((  Y$44444ll4(( ; 444A3777''-- ; 444]: & & ' ' LLaaa1f & & & ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 9 9 9!!sEPQRRRB((..I?y#&66666>'C+<<<<<>%++++*****M&+&+s"E66E: =E: cd}d}t|dz }t|||||d}tD]U}tjt 5||dz|dddn #1swxYwYVdS)NrXrrJr{)intr-r|rCwarnsrr})r$rr&r%r'rrs r,1test_warning_n_components_greater_than_n_featuresrsJIZ!^$$J "'    D1DD \3 4 4 D D  *q. 9 9 9 = =d C C C D D D D D D D D D D D D D D DDDs#A==B B cd}d}t|dz }t|||||d}t|||||d}tD]v}|dd|}|dd|} t t |jt | jwdS) NrrXrrrrHrJ)rUr()rr-r|r}rr4r) r$rr&r%r'rrrrp_dense rp_sparses r,test_works_with_sparse_datarsJIZ!^$$J(' J*' K1  ##CCCGG SS$$!!DDDHHUU ! H( ) )793H+I+I      r.c4tdddksJdS)zyTest Johnson-Lindenstrauss for small eps. Regression test for #17111: before #19374, 32-bit systems would fail. r7h㈵>rAlJWN)r rMr.r,"test_johnson_lindenstrauss_min_dimrs( )$ 7 7 7< G G G G G Gr.random_projection_clsc~t|ttt|d}|d}|||}|jtj fdt|j Dt}t||dS)NrrIr{cg|]}| SrMrM).0iclass_name_lowers r, z.s%RRRa  !a ! !RRRr.)dtype)r-r%r&r'r}get_feature_names_out__name__lowerrarrayrangerobjectr)r$rrrrandom_projection names_outexpected_names_outrs @r,(test_random_projection_feature_names_outrs #'    D.-1===$!7799I,5;;==RRRR51B1P+Q+QRRR y"455555r.r%)rI r5 r6r&compute_inverse_componentsTFcd}||||}t|||||zdzdz|d}t|||||zdzdz|d} || fD]} tj5tjddt || } dddn #1swxYwY|r(t |d sJ|j} | j||fksJ| | } | j| jksJ| | }t | d r| } t| |d d dS)Nr5)rUrr(r7rJ)r'r(r)rignorez>The number of components is higher than the number of features)messagecategoryinverse_components_r#gHz>g|=)rtolatol) r-warningscatch_warningsfilterwarningsrrhasattrrrinverse_transformrr#r )r$r%r&rrrrUrX_denseX_csrXrinv_componentsprojected_backprojected_agains r,test_inverse_transformrsL--!#=' &z)S014' G $z)S014'    Eu KK  $ & & ; ;  #T2     *77::I ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; & F,.CDD D D D.BN!'J +EEEEE*< > 1 1 1 1   . . . . . .r.cd}tjd}|dd}|d}|d}||tj}||tj}t|||t|j |j dS)Nrrrrrr[)r) rrrrrrfloat32float64r rr)rrr*rrp_32rp_64 projection_32 projection_64s r,,test_random_projection_numerical_consistencyr:s D )   # #C TA ! !q 1 1 1E ! !q 1 1 1E''(<(<==M''(<(<==MM=t<<<< !2E4EFFFFFr.)Nr)HrirtypingrrnumpyrrC scipy.sparsesparser1sklearn.exceptionsrrsklearn.metricsrsklearn.random_projectionrr r r r sklearn.utils._testingr rrrrsklearn.utils.fixesrr__annotations__rall_random_matrixrrr|r-r4r%r&rr'mark parametrizerFrNrVrYrarergrlrqrxrrrrrrrrrrrrrint32int64rrrMr.r,rs[ HHHHHHHH///////.....'<&=$s)===&=%>c>>>,/FF)?(@DI@@@(@'A49AAA14MM """":   # : SZ'%/ 0 0  sc3Z cS#J cS#J::: 444   999EEE*;<<11=<1*BCC77DC7 F F F1 1 1 rFFF.99 B B:9 B.99:9*.99"/"/:9"/J.99 6 6:9 6F.991+1+:91+h.99DD:9D(.99  :9 :HHH.9902FGG66HG:96..99&:;;';<<02FGG5e}EE7K7KFEHG=<<;:9 7Kt02FGG! RZ RZ 2: 2:  / /HG /02FGG G GHG G G Gr.