>[g8\dZddlmZddlZddlZddlmZddlm Z m Z m Z m Z m Z mZmZmZmZmZmZmZddlmZmZmZddlmZmZmZmZej<j?djAddd Z!ej<j?djAddd Z"d e!jF_$d e"jF_$e d edzZ%e d e de dd e ddzd e d ze%d e dd g zd eddzd eddzd eddzd eddzd edd gdzded dgdzdeddgdzeddedde d e d dze d g ed g gZ&eD] Z'e'd vre&jQed!e'""ejRjUd#e&d$Z+ejRjUd#e&Dcgc]}e,|ee fs|c}d%Z-ejRjUd#e&Dcgc] }|e%k7r| c}d&Z.ejRjUd#e&d'Z/d(Z0d)Z1ejRjUd#e&Dcgc]}|jes|c}d*Z3ejRjUd#e&d+Z4d,Z5d-Z6ejRjUd#e&d.Z7ejRjUd#e&d/Z8d0Z9ejRjUd#e&d1Z:ejRjUd#e&d2Z;ejRjUd#e&d3Z.eval_kernel_for_thetaTs$..u5 % 0绽|=)rr8shaper9rr)r2r; K_gradientr=K_gradient_approxs` r<test_kernel_gradientrDJs6]F1D1MAz   A !''!* ,, ,   A !''!* ,, ,   A &,,"4"4Q"7 77 7 'v||5JER $5q9r>ct|}|j}|td\}}t|jj j 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..qs!AY(r>c$|jdS)NrH)endswithrJs r<rLz#test_kernel_theta..qs 1::i;Pr>c34K|]}|jywN)name).0hyperparameters r< z$test_kernel_theta..ssP9O~~""9Osr%rHrrr.*+)rr9r8r __class____init__ parametersvaluesrQmapfiltersethyperparameters enumeratenploggetattr get_paramsrArrIrsetattr)r2r9_rB init_signpargs theta_varsirSparams kernel_class new_kernelK_gradient_news r<test_kernel_thetarp^s6]F 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:9K:cnt|}|t}|tt}t||dyNr)rr8r)r2K_autoK_crosss r<test_auto_vs_crossrus-6]F AYFQlG+r>ct|}tj|t}|jt}t ||dyrr)rradiagr8r)r2 K_call_diagK_diags r<test_kernel_diagrzs6 6]F''&)$K [[^F VQ/r>cttddztdtdzttdtdzttddzty)Nrr/r)rrr8rWr>r< test_kernel_operator_commutativer|sVSC+ncCHna-@As3x+nc#hna-@Ar>c8dtddgz}|t}tj}|dddfxxdzcc<dtd|z}t||tj}|dddfxxdzcc<dtd|z}t|||jt j dz|_t|jt j gdt|jjd d gy) Nrr"rrr@rr)g@r/r)r/r)) rr8copyrr9rarbrk2r)r2r;X1K1X2K2s r<test_kernel_anisotropicrs 3Sz? "Fq A Bq!tHMH xs3x| B2 Bq!tHMH xs3x| B2<<"&&)+FLv||RVVO%<=vyy--Sz:r>ct|}|ttdz}t|dtj|y)Nr)rr)rr8rrarw)r2r;s r<test_kernel_stationaryrs36]Fq!a%A$,r>c t|}t|tr%|j|jjk(sJt|t r>|j|j jxs|jjk(sJyyrP)r isinstancer requires_vector_inputr2r k1rr2s r<test_kernel_input_typersv 6]F&.)++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|}t|}||k(sJt|t|k7sJ|j|jk(sJt||yrP)ridrdr)r2 kernel_cloneds r<test_kernel_clonersi 6]F&MM ] "" " f:M* ** *    -":":"< << <  6r>ct|}d}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_paramsr s 6]FF&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_kernelr(s((#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>ct|}|tk7r)|t}tt|}t |||tt }ttt |}t ||y)N)r1)rkernel_rbf_plus_whiter8rrY)r2rrs r<test_kernel_versus_pairwiserEs\ 6]F&&AY a /!"b) 1B !Qv .Bb"%r>c t|}d}|j}|jD]}tdt |j r|j dk(r2|j }|dkDrBttj|j|||z||j||z }ttj|j|||j|dz }d}d}|jD]}tdt |j r|j dk(r2|j }|dkDr[|jdi|j|g|zittj|j|||z|g|z||z }|jdi|j|ittj|j|||dz }y)Nrstringr%r rW) rrdr_rtyper n_elementsrrarr9rQr)r2indexrlrSsizevalues r<test_set_get_paramsrXs 6]F E    F 00 h^%:%: ; <$$/(( !8 v||EEDL9:F>CVCVc0t|}t|yrP)rreprrs r<test_repr_kernelsrs 6]F 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_kernelrs= 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 RandomStatenormalr8rflags writeablerkernelsr1appendmark 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/9d F * * ,,7+0,0B;* GFG&v/C/C/EvGF--7+ , (.7+7,7"7+H,H:1:7+&,&$7+','T7+,W vnGsN N$ N) 5N)