{Kg5dZddlmZddlZddlZddlmZddlm Z m Z m Z m Z m Z mZmZmZmZmZmZmZddlmZmZmZddlmZmZmZmZej<j?djAddd Z!ej<j?djAddd Z"e d ed zZ#e d e de dd e ddzd e d ze#d e dd g zd eddzd eddzd eddzd eddzd edd gdzd ed dgdzdeddgdzeddedde d e d dze d g ed g gZ$eD] Z%e%dvre$jMed e%!"ejNjQd"e$d#Z)ejNjQd"e$Dcgc]}e*|ee fs|c}d$Z+ejNjQd"e$Dcgc] }|e#k7r| c}d%Z,ejNjQd"e$d&Z-d'Z.d(Z/ejNjQd"e$Dcgc]}|jas|c}d)Z1ejNjQd"e$d*Z2d+Z3d,Z4ejNjQd"e$d-Z5ejNjQd"e$d.Z6d/Z7ejNjQd"e$d0Z8ejNjQd"e$d1Z9ejNjQd"e$d2Z:d3Z;ycc}wcc}wcc}w)4z+Testing for kernels for Gaussian processes.) signatureN)clone) RBFCompoundKernelConstantKernel DotProductExponentiationExpSineSquaredKernelOperatorMaternPairwiseKernelRationalQuadratic WhiteKernel_approx_fprime)PAIRWISE_KERNEL_FUNCTIONSeuclidean_distancespairwise_kernels)assert_allcloseassert_almost_equalassert_array_almost_equalassert_array_equal))r@ length_scale@ noise_level)?r)length_scale_boundsg$@)constant_valuegQ?fixed)rr#r")rnu?@@)ralpha)r periodicity)sigma_0r) additive_chi2chi2?)gammametrickernelcdtd\}}|jdtjdk(sJ|jdtjdk(sJ|jdjjdk(sJfd}tj|d}t ||dy) NT eval_gradientrrrcFj|}|td}|S)NFr4)clone_with_thetaX)theta kernel_cloneKr2s o/home/alanp/www/video.onchill/myenv/lib/python3.12/site-packages/sklearn/gaussian_process/tests/test_kernels.pyeval_kernel_for_thetaz3test_kernel_gradient..eval_kernel_for_thetaOs$..u5 % 0绽|=)r8shaper9rr)r2r; K_gradientr=K_gradient_approxs` r<test_kernel_gradientrDFs1D1MAz   A !''!* ,, ,   A !''!* ,, ,   A &,,"4"4Q"7 77 7 'v||5JER $5q9r>c|j}|td\}}t|jjj j }|Dcgc]}|jdk7s|j }}tdtd|}td|jDt|k(sJt|jD]6\}} ||tjt|| jk(r6Jt|jD](\}} |j!} d| | jdz<|j} | di| } | td\}} |j"d | jj"d d zk(sJ|j"d | j"d d zk(sJ|d kDr1|d|| jd|k(sJt%|d d|f| d d|f|d zt'|jks||d zd| j|dk(sJt%|d |d zdf| d |df+t|jD]\}} tjd ||<||_t)t|| jd t+|| jdt)|j|tjdycc}w)NTr4selfc |dtd S)Nr_bounds)lenss r<z#test_kernel_theta..ks!AY(r>c$|jdS)NrH)endswithrJs r<rLz#test_kernel_theta..ks 1::i;Pr>c34K|]}|jywN)name).0hyperparameters r< z$test_kernel_theta..msP9O~~""9Osr%rHrrr.*+)r9r8r __class____init__ parametersvaluesrQmapfiltersethyperparameters enumeratenploggetattr get_paramsrArrIrsetattr)r2r9_rB init_signpargs theta_varsirSparams kernel_class new_kernelK_gradient_news r<test_kernel_thetarpYs LLE1D1MAz&**334??FFHI% :Iq6)9AFFID :(&1PRV*WJ P9O9OP PTWU   'v'='=>>Qx266'&.2E2E"FGGGG?'v'='=>>""$29~""Y./'' !+F+ 'q=>{{1~!1!1!7!7!:Q!>>>>"n&:&:1&=&AAAA q5!9 0 0! 44 44 z#rr'2N374K L q53v--. .Q>Z%5%5ab%99 99 z#q1uw,7QR9P Q!?&'v'='=>>66":a GFN,?,?@"E++R0FLLORVVBZ8 ?C ;s K/.K/cX|t}|tt}t||dyNr)r8r)r2K_autoK_crosss r<test_auto_vs_crossrus$AYFQlG+r>ctj|t}|jt}t||dyrr)radiagr8r)r2 K_call_diagK_diags r<test_kernel_diagrzs/''&)$K [[^F VQ/r>cttddztdtdzttdtdzttddzty)Nrr/r)rrr8rWr>r< test_kernel_operator_commutativer|sVSC+ncCHna-@As3x+nc#hna-@Ar>cLdtddgz}|t}tjt}|dddfxxdzcc<dtd|z}t ||tjt}|dddfxxdzcc<dtd|z}t |||j tj dz|_t|j tj gdt|jjd d gy) Nrr"rrr@rr)g@r/r)r/r)) rr8raarrayrr9rbrk2r)r2r;X1K1X2K2s r<test_kernel_anisotropicrs 3Sz? "Fq A !Bq!tHMH xs3x| B2 !Bq!tHMH xs3x| B2<<"&&)+FLv||RVVO%<=vyy--Sz:r>cp|ttdz}t|dtj|y)Nr)rr)r8rrarw)r2r;s r<test_kernel_stationaryrs* q!a%A$,r>c t|tr%|j|jjk(sJt|tr>|j|j jxs|j jk(sJyyrP) isinstancer requires_vector_inputr2r k1rr2s r<test_kernel_input_typerso&.)++v}}/R/RRRR&.)++ II + + Nvyy/N/N   *r>cttdg}|jrJttdtdg}|jsJy)Nrr rr)rrrrrs r<test_compound_kernel_input_typersK 9: ;F+++ + 93C;PQ RF  ' '' 'r>ctt|t|zD]3}|jdst||}t||}||k(r3Jy)Nhyperparameter_)r^dir startswithrc)kernel1kernel2attr attr_value1 attr_value2s r<check_hyperparameters_equalrsQCL3w</0 ??, -!'40K!'40K+- -- 1r>ct|}||k(sJt|t|k7sJ|j|jk(sJt||yrP)ridrdr)r2 kernel_cloneds r<test_kernel_clonersb&MM ] "" " f:M* ** *    -":":"< << <  6r>cd}t|}|j}ttf}d|vrt ||s|d}t j |r|d|d<||d<n|gdz|d<|dz|d<|jdi|t|}|j|jk(sJt|t|k7sJt||yyy)N)gh㈵>gj@rrr#rrW) rrdr rrraiterable set_paramsrr)r2boundsrrlisotropic_kernelsrkernel_cloned_clones r<"test_kernel_clone_after_set_paramsrsF&MM    F'):; 6;L(Mn- ;;| $%1!_F> ",2F( )&2^a%7F> ",2QJF( )   *6*#M2"--/=3K3K3MMMM%&"]*;;;;#M3FG)Nr>ctddt}ttj|tj tj dtjtttd }tddt}t||tdt}ttjdt}t||t||d }d D]?}t|dt}t||zdt}t||Ad }t|dt}tdt}t||d y)Nr'r/)r&rrF)squaredr"rr?)r"r'r(dr)decimal) r r8rrarwonesrAexprrinfr)r;K_absexpK_rbftinyr&rrlarges r<test_matern_kernelrs((#C(+Abggaj"''!''!**=>vv*1a??@H(#C(+Aa* !CS !! $E+"&&s+A.Aa'Au D ,Vr ,Q / 3VrDys 3A 6!"b) E +5s +A .B # q !Bb"a0r>c|tk7r)|t}tt|}t|||tt}ttt|}t||y)N)r1)kernel_rbf_plus_whiter8rrY)r2rrs r<test_kernel_versus_pairwiser9sS &&AY a /!"b) 1B !Qv .Bb"%r>c d}|j}|jD]}tdt|jr|jdk(r2|j }|dkDrBt tj|j|||z||j||z }t tj|j|||j|dz }d}d}|jD]}tdt|jr|jdk(r2|j }|dkDr[|jdi|j|g|zit tj|j|||z|g|z||z }|jdi|j|it tj|j|||dz }y)Nrstringr%r rW) rdr_rtyper n_elementsrrarr9rQr)r2indexrlrSsizevalues r<test_set_get_paramsrKs E    F 00 h^%:%: ; <$$/(( !8 v||EEDL9:F>CVCVct|yrP)reprrs r<test_repr_kernelsrus  Lr>ctddg}d}tjt|5|tdddy#1swYyxYw)Nr/rzeRationalQuadratic kernel only supports isotropic version, please use a single scalar for length_scale)match)rpytestraisesAttributeErrorr8)r2messages r<test_rational_quadratic_kernelr|s= S#J 7F " ~W 5q  6 5 5s AA )<__doc__inspectrnumpyrar sklearn.baser sklearn.gaussian_process.kernelsrrrrr r r r r rrrsklearn.metrics.pairwiserrrsklearn.utils._testingrrrrrandom RandomStatenormalr8rrkernelsr1appendmark parametrizerDrrprurzr|r is_stationaryrrrrrrrrrrrrs0r<rs1       II!##Aq&1II!##Aq&1- 0LLSJ'$'#4W ==#3 #C: &&&d @@&cc **&cc **&cc **&sCjS 11&sCjS 11&sCjS 113c25ssq cU ) ,(F ** NN>F;<( 7+:,:$ F6NN#CD .9.9b F * * ,,7+0,0B;* GFG&v/C/C/EvGF-- 7+ , (.7+ 7, 7 7+H,H81:7+&,&"7+&,&R7+, C tjGs%M= N =N N