U _{j$'@s`dZddlZddlZddlZddlmZddlmZddl m Z ddl m Z m Z mZddl mZddlmZdd lmZmZd d Zeedd d jZedddddgjZejeedkedZeeZ ej!ej"edZ#de#e dk<de#e dke dk@<de#e dk<e dddZ$e dde$e dddedde dddgZ%dd e%DZ&ej'(d!e%d"d#Z)d$d%Z*ej'(d!e&d&d'Z+ej'(d!e%d(d)Z,ej'(d!e%d*d+Z-ej'(d!e&d,d-Z.ej'(d!e%d.d/Z/d0d1Z0ej'(d!e&d2d3Z1ej'(d!e%d4d5Z2ej'(d!e%d6d7Z3d8d9Z4ej'(d:d!e die5d;fgdz,Testing for Gaussian process classification N) approx_fprime)ConvergenceWarning)GaussianProcessClassifier)RBFCompoundKernel WhiteKernel)ConstantKernel) MiniSeqKernel)assert_almost_equalassert_array_equalcCs t|S)N)npsin)xrf/var/www/html/bwave/python/venv/lib/python3.8/site-packages/sklearn/gaussian_process/tests/test_gpc.pyfsr @g@g@g@g@Zdtypegffffffֿgffffff??fixed length_scalelength_scale_boundsg?)r)MbP?@@g{Gz?Y@cCsg|]}|tkr|qSr) fixed_kernel).0kernelrrr /sr$r#cCs<t|dtt}t|t|tdddfdkdS)Nr#r?)rfitXyr predict predict_probar#gpcrrrtest_predict_consistent2sr.cCs`dddg}tdddg}tdd}t|d||}t||||ddd fd kdS) NAZABBTFr)Zbaseline_similarity_boundsr%rr&)r arrayr rr'r r*r+)r(r)r#r-rrr"test_predict_consistent_structured9s   r2cCs4t|dtt}||jj||jks0tdS)Nr%)rr'r(r)log_marginal_likelihoodkernel_thetaAssertionErrorr,rrrtest_lml_improvingBsr7cCs0t|dtt}t||jj|ddS)Nr%)rr'r(r)r r3r4r5r,rrrtest_lml_precomputedKs  r9cCsJt|dtt}tj|jjjtj d}|j |ddt |jj|ddS)Nr%rF)Z clone_kernelr8) rr'r(r)r Zonesr4r5shapeZfloat64r3r )r#r-Z input_thetarrrtest_lml_without_cloning_kernelTsr;cCszt|dtt}||jjd\}}tt |dk|jj|jj dddfkB|jj|jj dddfkBsvt dS)Nr%T-C6?rr) rr'r(r)r3r4r5r allabsboundsr6)r#r-lml lml_gradientrrrtest_converged_to_local_maximum^s rBcsJt|dtt|jd\}}t|jfddd}t||ddS)Nr%Tcs |dS)NF)r3)r5r-rrsz#test_lml_gradient..g|=)rr'r(r)r3r5rr )r#r@rAZlml_gradient_approxrrCrtest_lml_gradientls rGc Csd\}}tj|}|||dd}t|jddtd|jdddk}tddtd g|d g|d }tj }t d D]F}t |||d  ||} | | j j} | |ttjjkst| }qdS)N)rrr)ZaxisrFrrrr)r<r r)r#n_restarts_optimizerZ random_state)r random RandomStateZrandnr sumCrinfrangerr'r3r4r5ZfinfoZfloat32epsr6) global_random_seedZ n_samplesZ n_featuresrngr(r)r#Zlast_lmlrJgpr@rrrtest_random_startsys* ,  rUcsFfdd}t||d}|tt||jj||jksBtdS)Nc stj}|||dd}}tdD]\}t|td|dddftd|dddf}||dd}||kr&||}}q&||fS)NF)Z eval_gradientrrr)r rKrLrPZ atleast_1duniformmaximumminimum) Zobj_funcZ initial_thetar?rSZ theta_optZfunc_min_r5rrRrr optimizers   2  z(test_custom_optimizer..optimizer)r#r\)rr'r(y_mcr3r4r5r6)r#rRr\r-rr[rtest_custom_optimizers    r^cCsPt|d}|tt|t}t|dd|t}t t |d|dS)Nr%r) rr'r(r]r+X2r rMr*r r Zargmax)r#r-y_probZy_predrrrtest_multi_classs     racCsPt|d}|ttt|dd}|tt|t}|t}t||dS)Nr%r)r#Zn_jobs)rr'r(r]r+r_r )r#r-Zgpc_2r`Zy_prob_2rrrtest_multi_class_n_jobss      rbc Cstddgd}t|d}d}tjt|d|ttW5QRXtddgdtdd gd}t|d}t j d d }t d |ttt |d kst t|djtst |djjddkst t|djtst |djjddkst W5QRXttd }tddgddgd}t|d}t j d d }t d ||tt |d ks^t t|djtstt |djjddkst t|djtst |djjddkst W5QRXdS)Ngh㈵>r)rr%zThe optimal value found for dimension 0 of parameter length_scale is close to the specified upper bound 0.001. Increasing the bound and calling fit again may find a better value.match)Znoise_level_boundsrgj@T)recordalwaysrrzThe optimal value found for dimension 0 of parameter k1__noise_level is close to the specified upper bound 0.001. Increasing the bound and calling fit again may find a better value.rzThe optimal value found for dimension 0 of parameter k2__length_scale is close to the specified lower bound 1000.0. Decreasing the bound and calling fit again may find a better value.rrg$@r rzThe optimal value found for dimension 0 of parameter length_scale is close to the specified upper bound 100.0. Increasing the bound and calling fit again may find a better value.zThe optimal value found for dimension 1 of parameter length_scale is close to the specified upper bound 100.0. Increasing the bound and calling fit again may find a better value.)rrpytestZwarnsrr'r(r)rwarningscatch_warnings simplefilterlenr6 issubclasscategorymessageargsr Ztile) r#r-Zwarning_messageZ kernel_sumZgpc_sumreZX_tileZ kernel_dimsZgpc_dimsrrrtest_warning_boundssR          rpzparams, error_type, err_msgz!kernel cannot be a CompoundKernelc Cs4tf|}tj||d|ttW5QRXdS)z0Check that expected error are raised during fit.rcN)rrgZraisesr'r(r))paramsZ error_typeerr_msgr-rrrtest_gpc_fit_errors rs)7__doc__rhnumpyr rgZscipy.optimizerZsklearn.exceptionsrZsklearn.gaussian_processrZ sklearn.gaussian_process.kernelsrrrrrNZ4sklearn.gaussian_process.tests._mini_sequence_kernelr Zsklearn.utils._testingr r rZ atleast_2dZlinspaceTr(r_r1Zravelintr)ZfXemptyr:r]r!ZkernelsZnon_fixed_kernelsmarkZ parametrizer.r2r7r9r;rBrGrUr^rarbrp ValueErrorrsrrrrsp                         M