U e{j@sdZddlZddlZddlmZddlmZmZddl m Z m Z m Z ddl mZddlmZejdgfejdgfejdgfejdgfejdgfgZejdgfejd dgfejd dgfgZejejejejgZGd d d Z Gd d d Z!GdddZ"Gddde Z#GdddZ$GdddZ%Gddde Z&GdddZ'e(ddd d ggdd d ggde(dd d gdgd d gdgde)ddd d ggdd d ggde)ddd d ggd d gdgdej)ddddd d gdgd d gdgde)ddd d ggdd d gdgde)d dd d ggdd d gd d ggdejd d d!dd d ggd d gdejd d d!dd gd de*d d dd d"dd d ggd#d#gde*d d dd d"dd gd$dejd d d!ejd dd!@d d d d gd%dejd d d!ejd dd!@d d d d gd d ggd ddgfdej+,e-e.dd d&gj/e0dge-d'd'gd'd(gd'd)gd'd*gd'd+ggdgZ1ej23d,e1d-d.Z4d/d0Z5Gd1d2d2e Z6d3d4Z7d5d6Z8d7d8Z9ej23d,e*ed e:e;eje d;d<d=d>Z?ej23d,e@ed ejAe=d?d:gej2j>e d;d<d@dAZBej23d,e*ed Be@ed @e*ed e*ed e*ed e(ded Be)dBed BgdCdDZCej23d,e(de)de)dBgdEdFZDdS)GzD This module tests fitting and model evaluation with various inputs N)assert_allclose)fittingmodels)Fittable1DModel FittableModelModel) Parameter) HAS_SCIPYc@sDeZdZdZddZejdeddZ ejde ddZ d S) TestInputTypez This class tests that models accept numbers, lists and arrays. Add new models to one of the lists above to test for this. cCsNd|_d|_tddd|_tddd|_tjddddf\|_|_dS)Ng333333@g@r 皙?) xynparangex1y1mgridy2x2selfr`/var/www/html/bwave/python/venv/lib/python3.8/site-packages/astropy/modeling/tests/test_input.py setup_class,s zTestInputType.setup_class)modelparamscCs*||}||j||j||jdSN)rrrrrrmrrr test_input1D3s  zTestInputType.test_input1DcCs6||}||j|j||j|j||j|jdSr )rrrrrrr!rrr test_input2D:szTestInputType.test_input2DN) __name__ __module__ __qualname____doc__rpytestmark parametrizemodel1d_paramsr#model2d_paramsr$rrrrr %s    r c@seZdZdZddZddZddZdd Zd d Zd d Z ddZ ddZ e j je dde j deddZe j je dde j deddZe j je dde j deddZe j je dde j deddZdS) TestFittingz/Test various input options to fitting routines.cCs.td|_tjddddf\|_|_dS)Nr rrrrrrrrrrrEs zTestFitting.setup_classcCs^tddddg}td}ddddg|_||j}t}|||j|}t|j|dddS)z1 set 1D x, 1psetrr Hz>ZatolN) rarrayr Polynomial1D parametersrrLinearLSQFitterr)rexpectedp1rpfitrrrrtest_linear_fitter_1setIs  z#TestFitting.test_linear_fitter_1setcCstddgddgddgddgg}tjddd}ddddddd d g|_i}td D]}||g||j|<qRtjdd di|}||jd d }t }|||j|}t |j |dddS)z%1 set 1D x, 2 sets 1D y, 2 param_setsrr r r0n_models?@g@model_set_axisFrAr1r2N)r0) rr3rr4r5range param_namesrrr6r param_setsrr7r8rirr9rrrrtest_linear_fitter_NsetTs" z#TestFitting.test_linear_fitter_NsetcCsbtddddggj}td}ddddg|_||j}t}|||j|}t |j |d d d S) z71 pset, 1 set 1D x, 1 set 1D y, Chebyshev 1D polynomialg@g@gC@gq@r0rr r {Gz?r2N) rr3Tr Chebyshev1Dr5rrr6rrE)rr7Zch1rr9rrrrtest_linear_fitter_1dchebcs  z%TestFitting.test_linear_fitter_1dchebcCsbtddddggj}td}ddddg|_||j}t}|||j|}t |j |d d d S) zH 1 pset, 1 set 1D x, 1 set 1D y, Legendre 1D polynomial g@g@gm@gv@r0r r r@-q=r2N) rr3rJr Legendre1Dr5rrr6rrE)rr7Zleg1rr9rrrrtest_linear_fitter_1dlegendws  z'TestFitting.test_linear_fitter_1dlegendcCstd}ddddddg|_ddddddg}||j|j}t}|||j|j|}t|j|ddt||j|j|dddS) Nr rr r0r@rMr2)r Polynomial2Dr5rrrr6r)rp2r7zr9rrrrtest_linear_fitter_1set2ds z%TestFitting.test_linear_fitter_1set2dc CsZd}tjt|d>td}||j}tjddd}t}|||j|W5QRXdS)z A ValueError is raised if a 1 data set (1d x, 1d y) is fit with a model with multiple parameter sets. zHNumber of data sets .* is expected to equal the number of parameter setsmatchrPr r;N)r)raises ValueErrorrr4rrr6)rMESSAGEr8rr9rrrtest_wrong_numpsets  zTestFitting.test_wrong_numpsetc Cstddgddgddgddgddgddgg}tjddd}i}tdD]}d|g||j|<qHtjdd di|}||jd d }t}|||j|}t |j |d d dS)z;A case of 1 set of x and multiple sets of y and parameters.r rr r0r@rPr;rAFrBr1r2N)rP) rr3rr4rCrDrrr6rrErFrrrtest_wrong_psets$  zTestFitting.test_wrong_psetrequires scipyreasonfittercCsD|}tjdddd}||j}|||j|}t|jdddgdS)z.1 set 1D x, 1 set 1D y, 1 pset NonLinearFitterr r0皙?)meanstddevN)r Gaussian1Drrr5)rr`g1rrrrrtest_nonlinear_lsqt_1set_1ds  z'TestFitting.test_nonlinear_lsqt_1set_1dc Csb|}d}tjt|d@tjddgddgddgd d }||jd d }|||j|}W5QRXd S)z51 set 1D x, 1 set 1D y, 2 param_sets, NonLinearFitterz6Non-linear fitters can only fit one data set at a timerUgffffff$@r r0g @gq= ףp?rar )rbrcr<FrBN)r)rWrXrrdr)rr`rYrer_rrrtest_nonlinear_lsqt_Nset_1dsz'TestFitting.test_nonlinear_lsqt_Nset_1dcCsX|}tjddddddd}||j|j}|||j|j|}t|jddddddgdS) z.1 set 2d x, 1set 2D y, 1 pset, NonLinearFitterr r0r@333333?rar)Zx_meanZy_meanx_stddevy_stddevthetaN)r Gaussian2Drrrr5)rr`g2rSrrrrtest_nonlinear_lsqt_1set_2dsz'TestFitting.test_nonlinear_lsqt_1set_2dc Cs|}d}tjt|d^tjddgddgddgddgddgddgd d }||j|j}|||j|j|}W5QRXd S) z41 set 2d x, 1set 2D y, 2 param_sets, NonLinearFitterTInput argument .* does not have the correct dimensions in .* for a model set with .*rUr r0r@rirarr )rjrkrlr<N)r)rWrXrrmrflattenr)rr`rYrnrSrgrrrtest_nonlinear_lsqt_Nset_2ds z'TestFitting.test_nonlinear_lsqt_Nset_2dN)r%r&r'r(rr:rHrLrOrTrZr\r)r*skipifr r+fittersrfrhrorrrrrrr.Bs*         r.c@sXeZdZdZddZddZddZdd Zd d Zd d Z ddZ ddZ ddZ dS)TestEvaluationzs Test various input options to model evaluation TestFitting actually covers evaluation of polynomials cCs.td|_tjddddf\|_|_dS)Nr r/rrrrrs zTestEvaluation.setup_classcCshtjddgddgddgdd}||jdd}t|dd d f|d d d ftgksdtd S) zY This case covers: N param sets , 1 set 1D x --> N 1D y data r r0rar r;FrBrNr )rrdrrallZnonzeror3AssertionError)rrerrrrtest_non_linear_NYsetsz$TestEvaluation.test_non_linear_NYsetcCsbtjddgddgddgdd}t|j|jg}||}t|dddf|dddfd d dS) zU This case covers: N param sets , N sets 1D x --> N N sets 1D y data r r0rar r;Nrr rMr2)rrdrr3rr)rrexxrrrrtest_non_linear_NXYsetsz%TestEvaluation.test_non_linear_NXYsetcCs&td}||j}|jdks"tdS)z 1 data set, 1 pset, Polynomial1Dr@)rvN)rr4rshaperxrr8rrrrtest_p1_1set_1psets  z!TestEvaluation.test_p1_1set_1psetcCsbtjddd}|t|j|jgjdd}|jdks8tt|dddf|d ddfd d dS) z'N data sets, N param_sets, Polynomial1Dr@r r;rB)rvr rNr rMr2) rr4rr3rrJr|rxrr}rrrtest_p1_nset_npset"sz!TestEvaluation.test_p1_nset_npsetcCs*td}||j|j}|jdks&tdS)z#1 pset, 1 2D data set, Polynomial2DrP)r r N)rrQrrr|rx)rrRrSrrrtest_p2_1set_1pset*s z!TestEvaluation.test_p2_1set_1psetcCsNtjddd}t|j|jg}t|j|jg}|||}|jdksJtdS)z$N param_sets, N 2D data sets, Poly2drPr r;)r r r N)rrQrr3rrr|rx)rrRrzyyrSrrrtest_p2_nset_npset1s  z!TestEvaluation.test_p2_nset_npsetcCszt|j|jgj}d|d<d|d<d|d<tjdddgd d gdd }||d d }t|j|j||dd }t||dS)z Test model set with negative model_set_axis. In this case model_set_axis=-1 is identical to model_set_axis=1. d)rrr rc)r rrPr r r0r@)Zc0c1r<rrBN)rr3rrJrr4rr|)rrzr8rZyy1rrrtest_nset_domain:s  zTestEvaluation.test_nset_domaincCsPtddgddgg}tjddd|d}tjddddf\}}|||dS)Nr>g?r0g@g@)Z cov_matrixr )rr3rrmr)rZcovgrrrrrtest_evaluate_gauss2dKsz$TestEvaluation.test_evaluate_gauss2dN) r%r&r'r(rryr{r~rrrrrrrrrrus   ruc@s$eZdZeZeZeddZdS) TModel_1_1cCs |||Sr r)rr8rRrrrevaluateVszTModel_1_1.evaluateN)r%r&r'rr8rR staticmethodrrrrrrRsrc@sPeZdZdZddZddZddZdd Zd d Zd d Z ddZ ddZ dS)&TestSingleInputSingleOutputSingleModela A suite of tests to check various cases of parameter and input combinations on models with n_input = n_output = 1 on a toy model with n_models=1. Many of these tests mirror test cases in ``astropy.modeling.tests.test_parameters.TestParameterInitialization``, except that this tests how different parameter arrangements interact with different types of model inputs. cCsBtdd}|d}t|ts tt|dks2t|dks>tdS)O Scalar parameters with a scalar input should return a scalar. r r rroN)r isinstancefloatrxrndimrtrrrr#test_scalar_parameters_scalar_inputfs  zJTestSingleInputSingleOutputSingleModel.test_scalar_parameters_scalar_inputcCs^tdd}|tdd}t|tjs,tt|dks>tt|dddd d gksZtd S) Scalar parameters should broadcast with an array input to result in an array output of the same shape as the input. r r rPrrP r7N)rrrrndarrayrxr|rwrrrr%test_scalar_parameters_1d_array_inputqs  zLTestSingleInputSingleOutputSingleModel.test_scalar_parameters_1d_array_inputcCsltdd}|tdddd}t|tjs4tt|dksFtt|dd d gd d d ggkshtdS)rr r r[r r0rr r0rrrrrN rrrreshaperrrxr|rwrrrr%test_scalar_parameters_2d_array_input}s  zLTestSingleInputSingleOutputSingleModel.test_scalar_parameters_2d_array_inputcCstdd}|tddddd}t|tjs6tt|dksHtt|dd gd d gd d ggddgddgddgggkstdS)rr r r r0rr r0r rrrrrrc+WNrrrrr%test_scalar_parameters_3d_array_inputs zLTestSingleInputSingleOutputSingleModel.test_scalar_parameters_3d_array_inputcCsVtddgddg}|d}t|tjs*tt|dksd!d?gggd@dAdgdBdCdggdDdEdgdFdGd"ggggt|d#dHdIgd&dJdKggd)dLdMgd,dNdOgggdPd$dQgdRd'dSggdTd*dUgdVd-dWgggdXdYd%gdZd[d(ggd\d]d+gd^d_d.ggggd`S)arrIrrrrrrrrrrrr r r0rrrr rvrrrrrrrrrrrrrrgq= ף[@gq= ףk@gGzt@g(\[@gk@g)\t@gGz[@g(\k@g= ףpt@gfffff[@gQk@gRt@rrrrrrrrrrrrrrrrrrrrrrrrrgq= ףj@gGzs@gj@g)\s@g(\j@g= ףps@gQj@gRs@gq= ף@^@gGz0t@g(\B^@g)\0t@gGzD^@g= ףp1t@gfffffF^@gR1t@gQ``@gq= ףm@gzGa`@gm@g ףp=b`@g(\m@g33333c`@gQm@N)r rr|rxr)rrrrrrrrrrsb &$*zQTestSingleInputDoubleOutputSingleModel.test_mixed_array_parameters_1d_array_inputNrrrrrrs  (r)rinputsoutputs)rr rr0)r)rr r)rr r )Zc0_0g333333?g#>n?gQ>n?r r rPr=g?g?g?r>rcCsV|d|d}t|r*t||dn(tt|D]}t|||d|q6dS)Nrrr)rZisscalarrrClen)rresultrGrrrtest_mixed_inputs  r"cCsGdddt}|d}|ddgd}ddgddgdg}t||D]\}}t||qB|d}|dddg}ddddgg}t||D]\}}t||qdS)Ncs6eZdZdZdZdZeZddZfddZ Z S)ztest_more_outputs..MFr r0cSs||||||fSr r)rrrarrrrsz%test_more_outputs..M.evaluatecsVtj||\}}|j||j}dd|D}dd|D}|jt|ff||S)NcSsg|] }|jqSr)r|).0outrrr sz9test_more_outputs..M.__call__..cSsg|]}|dkrdn|qS)r rr)r$r|rrrr&s)superZprepare_inputsrr5Zprepare_outputstuple)rargskwargsrrgrZ output_shapes __class__rr__call__s z%test_more_outputs..M.__call__) r%r&r'Zstandard_broadcastingrrrr#rr. __classcell__rrr,rMs r0r r>r=r?)rzipr)r0cr!r7rerrrtest_more_outputss r5c@s(eZdZdZdZdZdZeddZdS)TInputFormatterz6 A toy model to test input/output formatting. r rrcCs||fSr rr7rrrr/szTInputFormatter.evaluateN) r%r&r'r(rrrrrrrrrr6&s r6cCs t}|dd}|dkstdS)Nr r r)r6rxrr!rrrtest_format_input_scalars4s r9cCs>t}|ddgddg}t|tddgtddgfdSNr r )r6rrr3r8rrrtest_format_input_arrays:sr;cCs@t}tddggjtddggjf}||}t||dSr:)r6rr3rJr)rinputr!rrr#test_format_input_arrays_transposed@s$r=)r@rP)Z lookup_tabler]r^c Cs|dd}t||dddt||dddd|_t||dddt||dddt||ddd}tjt|d|ddd W5QRXtjt|d|ddd d d W5QRXd }tjt|d|dW5QRXd S)W Test calling a model with positional, keywrd and a mixture of both arguments. r r r7)rr3rrz-Too many input arguments - expected 2, got .*rUr0rrr3+Missing input arguments - expected 2, got 1Nrrr)rWrXr positionalrYrrrtest_call_keyword_args_1Gs rFrc Cs|d}t||ddd|_t||ddd}tjt|d|dddW5QRXtjt|d|ddddd W5QRXd }tjt|d |W5QRXd S) r>r )r)r3.Too many input arguments - expected .*, got .*rUr r0rrAz+Missing input arguments - expected 1, got 0NrCrDrrrtest_call_keyword_args_2ks rHr'c Cs|dd}d|_t||dddt||dddd}tjt|d|dddW5QRXtjt|d|ddddd W5QRXd }tjt|d |W5QRXd S) r>r r r?r@rGrUr0rrAz+Missing input arguments - expected 2, got 0N)rrr)rWrXrDrrrtest_call_keyword_args_3s rIc Cs|dd}t||dddt||ddd|}d|_t||dddt||dddt||ddd}tjt|d|ddd W5QRXtjt|d|ddd d d W5QRXd }tjt|d|dW5QRXd S)r>r r )Zx0r)rr?r@rGrUr0rrArBN)rcopyrr)rWrXrDrrrtest_call_keyword_mappingss rK)Er(numpyrr)Z numpy.testingrZastropy.modelingrrZastropy.modeling.corerrrZastropy.modeling.parametersrZ"astropy.utils.compat.optional_depsr r4rNrKZShiftZScaler,rQZ Legendre2DZ Chebyshev2Dr-ZLevMarLSQFitterZ TRFLSQFitterZ LMLSQFitterZDogBoxLSQFitterrtr r.rurrrr rZIdentityMappingrmmathZ MultiplyUfuncr3ZlinspacerJrZbroadcast_modelsr*r+r"r5r6r9r;r=Z Rotation2DZ Pix2Sky_TANZ Tabular2DZonesrsrFrdZ Tabular1DrHrIrKrrrrs           ;V 0Q c                4