o \i+@sdZddlZddlZddlZddlmZddlmZddl m Z ddl m Z m Z mZddl mZddlmZdd lmZmZd d Zeedd d jZegdjZejeedkedZeeZ ej!ej"edZ#de#e dk<de#e dke dk@<de#e dk<e dddZ$e dde$e dddedde dddgZ%dde%DZ&ej'(de%ddZ)d d!Z*ej'(de&d"d#Z+ej'(de%d$d%Z,ej'(de%d&d'Z-ej'(de&d(d)Z.ej'(de%d*d+Z/d,d-Z0ej'(de&d.d/Z1ej'(de%d0d1Z2ej'(de%d2d3Z3d4d5Z4ej'(d6de die5d7fgd8d9Z6ej'(de%d:d;Z7dd?Z9dS)@z+Testing for Gaussian process classificationN) approx_fprime)ConvergenceWarning)GaussianProcessClassifier)RBFCompoundKernel WhiteKernel)ConstantKernel) MiniSeqKernel)assert_almost_equalassert_array_equalcCs t|S)N)npsin)xr/var/www/www-root/data/www/176.119.141.140/sports-predictor/venv/lib/python3.10/site-packages/sklearn/gaussian_process/tests/test_gpc.pyfs r )@g@g@g@g@dtypegffffffֿ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_consistent2s*r/cCs\gd}tgd}tdd}t|d||}t||||dddfdkdS)NAABBTFTrbaseline_similarity_boundsr&rr')r arrayr rr(r r+r,r)r*r$r.rrr"test_predict_consistent_structured9s  *r9cCs4t|dtt}||jj||jksJdS)Nr&)rr(r)r*log_marginal_likelihoodkernel_thetar-rrrtest_lml_improvingBsr=cCs0t|dtt}t||jj|ddS)Nr&)rr(r)r*r r:r;r<r-rrrtest_lml_precomputedKsr?cCsJt|dtt}tj|jjjtj d}|j |ddt |jj|ddS)Nr&rF) clone_kernelr>) rr(r)r*r onesr;r<shapefloat64r:r )r$r. input_thetarrrtest_lml_without_cloning_kernelTsrEcCszt|dtt}||jjd\}}tt |dk|jj|jj dddfkB|jj|jj dddfkBs;JdS)Nr&T-C6?rr) rr(r)r*r:r;r<r allabsbounds)r$r.lml lml_gradientrrrtest_converged_to_local_maximum^s  rLcsJt|dtt|jd\}}t|jfddd}t||ddS)Nr&Tcs |dS)NF)r:)r<r.rrss z#test_lml_gradient..g|=)rr(r)r*r:r<rr )r$rJrKlml_gradient_approxrrMrtest_lml_gradientls rQc Csd\}}tj|}|||dd}t|jddtd|jdddk}tddtd g|d g|d }tj }t d D]#}t |||d  ||} | | j j} | |ttjjksdJ| }qCdS)N)rrr)axisrOrrr r)rFr!r)r$n_restarts_optimizer random_state)r random RandomStaterandnr sumCrinfrangerr(r:r;r<finfofloat32eps) global_random_seed n_samples n_featuresrngr)r*r$last_lmlrUgprJrrrtest_random_startsys( ,  rgcsFfdd}t||d}|tt||jj||jks!JdS)Nc stj}|||dd}}tdD].}t|td|dddftd|dddf}||dd}||krA||}}q||fS)NF) eval_gradientrrr)r rWrXr] atleast_1duniformmaximumminimum) obj_func initial_thetarIrd theta_optfunc_min_r<rrarr optimizers   2  z(test_custom_optimizer..optimizer)r$rt)rr(r)y_mcr:r;r<)r$rartr.rrsrtest_custom_optimizers     rvcCsPt|d}|tt|t}t|dd|t}t t |d|dS)Nr&r) rr(r)rur,X2r rZr+r r argmax)r$r.y_proby_predrrrtest_multi_classs    r{cCsPt|d}|ttt|dd}|tt|t}|t}t||dS)Nr&r)r$n_jobs)rr(r)rur,rwr )r$r.gpc_2ryy_prob_2rrrtest_multi_class_n_jobss      rc Cstddgd}t|d}d}tjt|d|ttWdn1s&wYtddgdtdd gd}t|d}t j d d G}t d |ttt |d ksXJt |djtsbJ|djjddksnJt |djtsxJ|djjddksJWdn1swYttd }tddgddgd}t|d}t j d d H}t d ||tt |d ksJt |djtsJ|djjddksJt |djtsJ|djjddksJWddS1swYdS)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)noise_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.)rrpytestwarnsrr(r)r*rwarningscatch_warnings simplefilterlen issubclasscategorymessageargsr tile) r$r.warning_message kernel_sumgpc_sumrX_tile kernel_dimsgpc_dimsrrrtest_warning_boundss@        "rzparams, error_type, err_msgz!kernel cannot be a CompoundKernelcCsNtdi|}tj||d|ttWddS1s wYdS)z0Check that expected error are raised during fit.rNr)rrraisesr(r)r*)params error_typeerr_msgr.rrrtest_gpc_fit_errors "rcCsTt|d}|tt|t\}}|jtjdfksJ|jtjdfks(JdS)>Checks that the latent mean and variance have the right shape.r&rN)rr(r)r*latent_mean_and_variancerB)r$r. latent_meanlatent_variancerrr'test_gpc_latent_mean_and_variance_shape s  rcCsVttd}|tttjtdd|tWddS1s$wYdS)rr&zdReturning the mean and variance of the latent function f is only supported for binary classificationrN) rrr(r)rurr ValueErrorrrMrrrAtest_gpc_latent_mean_and_variance_complain_on_more_than_2_classes,s   "rcCs@gd}tgd}tdd}t|d||}||dS)Nr0r4rr5r&)r r7r rr(rr8rrr9test_latent_mean_and_variance_works_on_structured_kernels:s  r):__doc__rnumpyr rscipy.optimizersklearn.exceptionsrsklearn.gaussian_processr sklearn.gaussian_process.kernelsrrrrr[4sklearn.gaussian_process.tests._mini_sequence_kernelr sklearn.utils._testingr r r atleast_2dlinspaceTr)rwr7ravelintr*fXemptyrBrur"kernelsnon_fixed_kernelsmark parametrizer/r9r=r?rErLrQrgrvr{rrrrrrrrrrrsz                         I