o tBh&@sRdZddlZddlZddlmZddlZddlmZddl m Z m Z m Z mZddlmZddlmZddlmZmZd d Zeedd d jZegd jZejeedkedZeeZ ej!ej"edZ#de#e dk<de#e dke dk@<de#e dk<e dddZ$e dde$e ddde dde dddgZ%dde%DZ&ej'(de%ddZ)dd 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%d/d0Z2ej'(de%d1d2Z3d3d4Z4ej'(d5de die5d6fgd7d8Z6dS)9z,Testing for Gaussian process classification N) approx_fprime)GaussianProcessClassifier)RBFCompoundKernelConstantKernel WhiteKernel) MiniSeqKernel)ConvergenceWarning)assert_almost_equalassert_array_equalcCs t|S)N)npsin)xr~/var/www/html/riverr-enterprise-integrations-main/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)N)AABB)TFTr)baseline_similarity_boundsr&rr')r arrayrrr(r r+r,)r)r*r$r.rrr"test_predict_consistent_structured9s  *r5cCs4t|dtt}||jj||jksJdS)Nr&)rr(r)r*log_marginal_likelihoodkernel_thetar-rrrtest_lml_improvingBsr9cCs0t|dtt}t||jj|ddS)Nr&)rr(r)r*r r6r7r8r-rrrtest_lml_precomputedKsr;cCsJt|dtt}tj|jjjtj d}|j |ddt |jj|ddS)Nr&rF) clone_kernelr:) rr(r)r*r onesr7r8shapefloat64r6r )r$r. input_thetarrrtest_lml_without_cloning_kernelTsrAcCszt|dtt}||jjd\}}tt |dk|jj|jj dddfkB|jj|jj dddfkBs;JdS)Nr&T-C6?rr) rr(r)r*r6r7r8r allabsbounds)r$r.lml lml_gradientrrrtest_converged_to_local_maximum^s  rHcsJt|dtt|jd\}}t|jfddd}t||ddS)Nr&Tcs |dS)NF)r6)r8r.rrss z#test_lml_gradient..g|=)rr(r)r*r6r8rr )r$rFrGlml_gradient_approxrrIrtest_lml_gradientls rMc Csd\}}tjd}|||dd}t|jddtd|jdddk}tddtd g|d g|d }tj }t d D]#}t ||dd  ||}| |j j} | |ttjjksdJ| }qCdS)N)rrrr)axisrKrr r)rBr!r)r$n_restarts_optimizer random_state)r random RandomStaterandnr sumCrinfrangerr(r6r7r8finfofloat32eps) n_samples n_featuresrngr)r*r$last_lmlrQgprFrrrtest_random_startsys$ ,  rbcCsBdd}t||d}|tt||jj||jksJdS)Nc Sstjd}|||dd}}tdD].}t|td|dddftd|dddf}||dd}||krA||}}q||fS)NrF) eval_gradientrr)r rSrTrY atleast_1duniformmaximumminimum) obj_func initial_thetarEr_ theta_optfunc_min_r8rrrr optimizers  2  z(test_custom_optimizer..optimizer)r$rn)rr(r)y_mcr6r7r8)r$rnr.rrrtest_custom_optimizers   rpcCsPt|d}|tt|t}t|dd|t}t t |d|dS)Nr&r) rr(r)ror,X2r rVr+r r argmax)r$r.y_proby_predrrrtest_multi_classs    rucCsPt|d}|ttt|dd}|tt|t}|t}t||dS)Nr&r)r$n_jobs)rr(r)ror,rqr )r$r.gpc_2rsy_prob_2rrrtest_multi_class_n_jobss      ryc 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.)rrpytestwarnsr r(r)r*rwarningscatch_warnings simplefilterlen issubclasscategorymessageargsr tile) r$r.warning_message kernel_sumgpc_sumr}X_tile kernel_dimsgpc_dimsrrrtest_warning_boundssP            "rzparams, error_type, err_msgz!kernel cannot be a CompoundKernelcCsNtdi|}tj||d|ttWddS1s wYdS)z0Check that expected error are raised during fit.rzNr)rrraisesr(r)r*)params error_typeerr_msgr.rrrtest_gpc_fit_errors "r)7__doc__rnumpyr scipy.optimizerrsklearn.gaussian_processr sklearn.gaussian_process.kernelsrrrrWr4sklearn.gaussian_process.tests._mini_sequence_kernelrsklearn.exceptionsr sklearn.utils._testingr r r atleast_2dlinspaceTr)rqr4ravelintr*fXemptyr>ror"kernelsnon_fixed_kernelsmark parametrizer/r5r9r;rArHrMrbrpruryr ValueErrorrrrrrsp                        M