+ LV-jfW<^RIt^RIt^RIt^RIt^RIt^RIt^RIt^RIt^RIt^RI t ^RI H t H t HtHtHtHt^RIHuHt^RIHt^RIt^RIHtHtHtHtHtHtH t H!t!H"t"H#t#H$t$H%t%H&t&H't'H(t(H)t)H*t*^RI+H,t-^RI.H/uH,t0^RI1H2t2H3t3H4t4H5t5H6t6^RI7H8t8H9t9^RIH:t:H;t;HuH?t@^RIAHBtB^R ICHDtD^R IEHFtF^R IHGtH^R IHItJ^R IHKtK]PPtM]!]4!RR44tN]!]4!RR44tORtPRtQRtRRtSRtT] 3RltURPRltV^#^] 3RltWRtX!RR4tY]!]4!RR44tZRt[RQR lt\]PPR!R"R#.4t^]!]4!R$R%44t_!R&R'4t`RRR(ltaR)tbR*tc!R+R,4tdR-te]!])4!R.R/44tfR0tgR1thR2ti!R3R44tj]!](4!R5R644tk!R7R84tl!R9R:4tmR;tnR<to]!]:4!R=R>44tpR?tq!R@RA4tr!RBRC4ts]!];4!RDRE]s44tt]!];4!RFRG]s44tu]!]<4!RHRI44tv]!]<4!RJRK44tw!RLRM4tx!RNRO4tyR#)SN)xp_assert_equalxp_assert_closexp_default_dtype concat_1dmake_xp_test_casexp_ravel)raises)BSplineBPolyPPolymake_interp_splinemake_lsq_splinesplevsplrepsplprepsplder splantidersprootsplintinsert CubicSpline NdBSplinemake_smoothing_splineRegularGridInterpolator) _not_a_knot_augknt_woodbury_algorithm_periodic_knots_make_interp_per_full_matr) Fperiodic root_rati)generate_knots make_splrep make_splprep) AxisError)_run_concurrent_barrier) make_ndbspl) _dfitpack) _bsplines)_dierckxc(a]tRt^1toRtRtRtRtRtRt ] !RRR 7] !R R R 7] !R R R 7R444t Rt Rt RtRtRtRtRtRtRtRtRtRtRtRtRtRtRtR tR!tR"tR#t R$t!R%t"]#PHPKR&]&!R<^44R'4t'R(t(R)t)]#PHPKR*.R=O4]#PHPKR+]&!^44R,44t*R-t+R.t,R/t-]#PHPKR0.R>O4R14t.]#PHPKR0.R>O4R24t/R3t0R4t1R5t2]#PHPg]4PjR68HR7]6R87R94t7R:t8Vt9R;#)? TestBSplinec \\\3\3/\ ^R.R.^R7B\ P !RR7;_uu_4\\\3/\ ^\ P.R.^R7BRRR4\\\3/\ ^\ P.R.^R7B\\\3/\ ^R .R.^R7B\\\3/\ ^.^..R.^R7B\\\3/\ .R O^.^R7B\\\3/\ .R ORR.^R7B\\\3/\ .R O.RORR7B\\\3/\ .R O.RORR7B\\\3/\ .RO.RO^R7B^ ^r2VPW#,^,VPR 7pVP\ PPV44p\WEV4p\WFP4\WVP4W6P 8XgQhR# +'giEL;i)??tckignore)invalidNcubic@dtyperr.rr.r=)rAr0@@@)r0r0r0)rArr.r.r=r?)r0r.r.) assert_raises TypeError ValueErrorr dictnperrstatenaninfarangefloat64asarrayrandomrr2r3r4)selfxpnr4r2r3bs&& v/Volumes/fast/ai/experiments/ui-tars-smoke/.venv/lib/python3.14/site-packages/scipy/interpolate/tests/test_bsplines.py test_ctorTestBSpline.test_ctor4sy*-w 1!SbTQ/ 1 [[ * * *g RBFF tq1Q R+j'NTQKB41-MNj'JTQGtq-IJj'MTaS1#J2$!-LM j'KTI!-JKj' <RH: < i O1\WM Oi K1\SI K j' B, a@ B11 IIac!e2::I . JJryy''* + A! 3333CCxx;+ * *s 3I I/ cr\4pVPp\VPV^,RRR7\VPV^,RRR7VP V^,8XgQh\ P!\4;_uu_4RVnRRR4R# +'giR#;i)rV瞯<atolrtolfooN) _make_random_splinetckrr2r3r4pytestr AttributeError)rQrTr_s& rUtest_tckTestBSpline.test_tckWsz  !eeSV%e<SV%e<ssc!f}}]]> * *AE+ * * *s B%% B6 c\VP.RO4VPR^.4^R7pVP^^^ 4pVPVR8VPRVPR7R4p\ V!V4V4R#)rr0rCr1r9rDN)rr0r=)r rOlinspacewhererNr)rQrRrTxxexpecteds&& rUtest_call_namespaceTestBSpline.test_call_namespacecso bjj, B70Cq I [[Ar "88BGRZZ"**Z%EsK"x(c,\P!^^^ 4p\^^.R.^R7p\V!V4\P!V4R,4\.RO^^.^R7p\V!V4\P !VR8RR44R#)rrCr1ffffff?rDN)rrmr.)rIrer r ones_likerf)rQrgrTs& rU test_degree_0TestBSpline.test_degree_0lsu [[Ar " q!f *"r||B/#56 lq!f 2"rxxT 3<=rkc >VP.RO4pVP.R O4p^p\W#V4pVPRR^2VPR7p\ V!V4V^,\ WaR7,V^,\ V^, VR7,,V^,\ V^, VR7,,RR7\ \PWbV34wrxp \WxW34p \ V!V4VPV 4RR7R#) rr0rCr9rR+=r[Nr>)r0r=r?) rOr rerNrB_012maprIr) rQrRr2r3r4rTxx_npt_npc_np splev_results && rU test_degree_1TestBSpline.test_degree_1us JJ ' JJ{ #  A!  KKS"BJJK 7 aD aDq 1Q4acb(9#9 9AaDqsrARqQ [[b" %240RT2 @rkz dask.arrayz"_naive_eval is not dask-compatiblereasonz jax.numpyztoo slow; XXX a slow-if marker?torchz OOB on CIc \VR7pVPwr4pVPW5,W5)^, ,^24pV!V4pTPVUu.uFp\ WWEVR7NK up4p \ WyRR7TPVUu.uFp\ WWEVR7NK up4p \ WzRR7R#uupiuupi)rrrsrtN)r^r_restack _naive_evalr _naive_eval_2) rQrRrTr2r3r4rgy_brwy_ny_n2s && rUtest_rndm_naive_eval TestBSpline.test_rndm_naive_evals 2 &%%a [[qAw +ehhC1 A!26CDu-xx2F2aqQb92FG. DGs CCc \4pVPwr#p\P!W$,W$)^, ,^24p\ V!V4\ WRW434RR7R#r.rsrtN)r^r_rIrerrrQrTr2r3r4rgs& rUtest_rndm_splevTestBSpline.test_rndm_splevsL  !%%a [[qAw +"uRQ3%@rkc\PPR4p\P!VP^44pVP^4p\ W#4p\ V!pVP VPrv\P!Wg,Wg)^, ,^P4p\V!V4\W4RR7R#)rsrtN) rIrP RandomStatesortrr r2r4rerr) rQrngrwyr_rTr2r4rgs & rUtest_rndm_splrepTestBSpline.test_rndm_splrepsii##D) GGCJJrN # JJrNQl SMssACC1 [[qAw +"uR~E:rkcZ\VR7pVPVP4VnVPVPVP ,VPVP )^, ,^dVP R7p\V!V4VPV44R#)rrr9N)r^rnr3rer2r4rNr)rQrRrTrgs&& rUtest_rndm_unityTestBSpline.test_rndm_unityso 2 &ll133 [[QSS133tAv;2::[ F"r||B/0rkc\PPR4p^^rC\P!VPV44pVPV^^3R7p\ VP WV34wrV\ WVV4pWT,WT)^, ,rWV, VP VPR44,,p V!V 4PR8XgQhR#)rsizeNr?r@)r?r@r)rIrPrrrvrOr shape) rQrRrrSr4r2r3rTtmtprgs && rUtest_vectorizationTestBSpline.test_vectorizationsii##D)11 GGCJJqM " JJQ1IJ &2::v& A! qAwB 7bjjI)>?? ?u{{o---rkc \PPR4p^!^r2\P!VPW#,^,44pVPV4p\PWQPV^,43,p\ WEV4\ WFV4rVR,V^,, p \P !V^,V , VR,V ,^24p \V!V 4V!V 4RR7\V!V 4\WWS34RR7\V!V 4\WWc34RR7R#)rrsrtNr;) rIrPrrr_r rerr) rQrrSr4r2r3c_padrTb_paddtrgs & rU test_len_cTestBSpline.test_len_csii##D)11 GGCJJqs1u% & JJqMaAaC()1#WQq%95 rUQqT\ [[1AbEBJ 3"uRyu5"uRQ3%@"uRU7eDrkc \4pVPwr4pW5,W5)^, ,rvRF1p\V!Wg.V4V!VR,VR, .V4RRR7K3 R#)r.绽|=& .>Hz>rZNTFr^r_r) rQnum_parallel_threadsrTr2_r4rrextraps && rUtest_endpointsTestBSpline.test_endpointssb  !%%aqAwB#F Arh/rEz2:6?dQU W$rkc \4pVPwr4p\V!W5^,V)^, R, 4V!W5^,V)^, R,4RR7R#)r.rrrtNr)rQrrTr2rr4s&& rUtest_continuityTestBSpline.test_continuitysV  !%%a!aC1+-.!aC1+2E0F rkc \VR7pVPwr4pVR,V^,, pVPW5,V, W5)^, ,V,^24pW5,V8WsV)^, ,8,p\V!Wx,RR7V!Wx,RR74\ \ P WsV34wrp VP \WW3^R74p \V!VRR7V 4R#)rrTrF)extNr;)r^r_rerrvrIrOr) rQrRrTr2r3r4rrgmaskxx_npryrzr{s && rU test_extrapTestBSpline.test_extraps 2 &%%a rUQqT\ [[AbdGbL" 5r baRT7l+ "(5"(6 8  RAJ7Tzz%dA"FG "$/>rkc\4pVPwr#pV^,^, VR,^,.pV!V4p\P!\P!V44'dQhR#)rNr;)r^r_rIallisnan)rQrTr2rr4rgyys& rUtest_default_extrapTestBSpline.test_default_extrapsZ  !%%adQh" " rU66"((2,'''''rkc \PPR4p\P!VP^44pVP^4p\ VP W434wr4^p\ W4VRR7pVP^,V^,, pVR ,V^,, pVPW5,V, W7,V,^24p W5,WV,, W7,W5,, ,,p \ \P WV34wrp \V4VP8XdRMRp\V!V 4VP \WW344VR7VP .R O4p W5,WV,, W7,W5,, ,,p \V!V RR7V!V RR7\V4VP8XdRMRR7R #) rperiodicr-q=gH׊>rtTrsgƠ>Nr;)r;r?r.) rIrPrrrvrOr rrerrNrr)rQrRrr2r3r4rTrSrrgxyxy_npryrzr[s&& rUtest_periodic_extrap TestBSpline.test_periodic_extrapsyii##D) GGCJJqM " JJqM2::v&  A! 4 GGAJ!a%  rUQqT\ [[AD2Ir 2 TRA$Y14!$;/ / RAJ7T(, :u bE2::eE$?;<4 ZZ ( TRA$Y14!$;/ / bj ) bd #*2."**<$ rkc\4pVPwr#p\P!W#V34p\P !W$,W$),^d4p\ V!V4V!V4RRR7R#)drsrZN)r^r_r from_splinerIrer)rQrTr2r3r4pprgs& rU test_ppolyTestBSpline.test_ppolysY  !%%a   ay ) [[quc *"r"vE>rkc \4pVPwr#p\P!V^,VR,^24p\PWR3,p\ ^V^,4F&p\ WRW43VR7p\Wq!WVR7RR7K( \V!WT^,R7\P!V4RR7R#)rdernursrtNr;) r^r_rIrerrangerr zeros_like)rQrTr2r3r4rgryds& rUtest_derivative_rndm TestBSpline.test_derivative_rndm#s  !%%a [[1qub ) UU25\AaC=Crq9#.B B" E :! "1 r}}R'8uErkc\^p.ROp\PPR4p\P^^VP^4^^3,p\ W$V4p\P !.RO4p\ V!Wf^8g,R, 4V!Wf^8g,R,44\P!V!R4V!R44'dQh\P !^^.4p\ V!VR, ^R7V!VR,^R74\P !^^.4p\P!V!VR, ^R7V!VR,^R74'dQh\P!V!VR, ^R7V!VR,^R74'dQhR#) r=rrrN) r;r;rr.r.r?r@rrrrr)r.r?r@rg2H@gη@)rIrPrrr rOrallclose) rQr4r2rr3rTrwx0x1s & rUtest_derivative_jumps!TestBSpline.test_derivative_jumps0sS  2ii##D) EE!Q 1 q!+ , A!  JJ| $!F)e+,!F)e+, .;;q{AgJ7777ZZA "u*+"u*+ - ZZA ;;qe2Ab5jQ4GHHHH;;qUq11QY13EFFFFFrkc VPR^^4p\P!VP.R O4R7p\ \ PW#P VP34wrEpVP\WEWcP344p\V!V4VRR7\V4VP8XdRMRp\V!V4VP\V4VPR7VR7\P!VP.R O4R7pVP^^^ VPR7p\V!V4VPV^8W",RV, ^,4RR7R#) r.)r2rsrtrr9rBNr;rr.r=r?rr.r.r=)rer basis_elementrOrvrIr2r3rr4rrrNB_0123rf) rQrRrgrTrryrzr{r[s && rUtest_basis_element_quadratic(TestBSpline.test_basis_element_quadraticJs [[Q #  ! !BJJ|$< = RaccN;Tzz%dCC/@"AB "|%8&r*bjj8Ud"rzz&*BJJzGdS  ! !BJJ|$< = [[Ar[ 4"aB 3% Arkc \4pVPwr#p\P!W$,W$)^, ,^4p\ V!V4\ WRW44RR7R#r)r^r_rIrer_sum_basis_elementsrs& rUtest_basis_element_rndm#TestBSpline.test_basis_element_rndmZsK  !%%a [[qAw +"22!?eLrkc\4pVPwr#pVR,p\W%V4p\W!PPV4p\W!PP V4p\ P!W$,W$)^, ,^4p\V!V4PV!V4RR7\V!V4P V!V4RR7R#)r0rsrtNy?@) r^r_r r3realimagrIrer) rQrTr2r3r4ccb_reb_imrgs & rU test_cmplxTestBSpline.test_cmplx`s  !%%a (^ A1 q##((A&q##((A& [[qAw +" DH59" DH59rkc\P!VP.RO44pVPV!VP44'gQhR#)rNr)r rrOrrK)rQrRrTs&& rUtest_nanTestBSpline.test_nanms8  ! !"**\": ;xx"&& """"rkc \^VR7pVPwr4p\W4V4pVPW5,W5)^, ,^4p\ ^V4F.pVP 4p\ V!Wx4V!V4RRR7K0 R#)r)r4rRrrZN)r^r_r rer derivativer) rQrRrTr2r3r4b0rgjs && rUtest_derivative_method"TestBSpline.test_derivative_methodrst ! +%%a Q1  [[qAw +q!A A BrIqu5u Erkc\VR7pVPwr4pVPW5,W5)^, ,^4p\VP 4P 4!V4V!V4RRR7VP WDV3^R7pVP WD3^R7p\W4V4p\VP 4P 4!V4V!V4RRR7R#)rrrsrZaxisN)r^r_rerantiderivativerrr )rQrRrTr2r3r4rgs&& rUtest_antiderivative_method&TestBSpline.test_antiderivative_method{s 2 &%%a [[qAw +((*557;"E 7 HHaAYQH ' HHaV!H $ A! ((*557;"E 7rkc 4 \P!VP.R O44p\P!VP ^^4RRR7'gQh\P!VP ^^4RRR7'gQh\P!VP ^^4RRR7'gQh\P!VP ^^4RRR7'gQh\P!VP ^^4RRR7'gQh\P!VP ^^4RRR7'gQh\P!VP R ^4RRR7'gQh\P!VP R ^RR7RRR7'gQh\P!VP R ^RR7RRR7'gQh\P!VP ^R RR7RRR7'gQh\P!VP ^R RR7\ P!^R VP4RR7'gQhRVn VP4pVPV!^4V!^4, VPR 7p\P!VP ^^4V4'gQh\P!VP ^^4R V,4'gQh\P!VP RR4V4'gQh\P!VP RR4^V,4'gQh\VP RR 4VPV!R 4V!R4, 44\VP R ^4VPV!^4V!^4, V!^4,V!R 4, 44\VP R^4VPV!^4V!^4, V!^4,V!R 4, 44\VP R ^4VPV!^4V!^4, V!^4,V!R 4, ^V,,44\VP ^R 4VPV!^4V!^4, 44\VP RR4VPV!^4V!^4, 44\VP ^R4VPV!^4V!^4, ^V,, 44R #)rrrsabs_tolrATrFrr9?Nr<r;iiig+@i) r rrOmathisclose integrate_implrr_rr rNr)rQrRrTi period_ints&& rU test_integralTestBSpline.test_integrals)  ! !"**Y"7 8||AKK1-sEBBBB||AKK1-xGGGG||AKK1-tUCCCC||AKK1-sEBBBB||AKK1-xGGGG||AKK1-tUCCCC||AKKA.UCCCC||AKKA4K@#uUUUU||AKKA5KA3PUVVVV||AKK25KA4QVWWWW||AKK25KA!LLB6GG GG#   ZZ!qt 2::Z> ||AKK1-z::::||AKK1-rJ????||AKKB/<<<<||AKKB/Z@@@@ C- 1S6AcF?3 5 C+ 1Q4!A$;1#5##>? A Hf5 1Q4!A$;1#5##>? A C0 1Q4!A$;1#5##>Z#OP R  Ar*BJJqtad{,CD B,bjj1!.EF Ar* 1Q4!A$;Z#?@ Brkc.ROp\W4pRVn\P!V4pRF/wrE\ VP WE4VP WE44K1 R#)rrNr>))r)rr)r )r rr rrr)rQrwrTprrs& rUtest_integrate_ppoly TestBSpline.test_integrate_ppolysS  q $"   a 5FB AKK/KK/ 16rkc\P!.RO4pRFLpVP^^VR7p\V\P 4'gQhVP ^8XdKLQh R#)rrNr<r)r rr isinstancerIndarrayndim)rQrTrress& rUtest_integrate_0D_always$TestBSpline.test_integrate_0D_alwayssW  ! !) ,(K++a +.BiN)__name__ __module__ __qualname____firstlineno____static_attributes__r.rkrUBr-s rkr4N)rr.r=r=)r r __class__rr )rQr4rTs& rUtest_subclassingTestBSpline.test_subclassingsg   OOL ){{a||~''1,,,!++q000rkr c ^^rCVP^^W4,^,4p.ROpV^,pVPWs4\V4p\PP R4pVP VP VR74p \WYWAR7p V PPWg,3VRV,Wg^,R,8XgQhVP R4pV !V4PVRVVP,Wg^,R,8XgQhV P)^, V P3F#p \\\3/\WYWKR7BK% \WYWAR7P4\WYWAR7P^4\WYWAR7P4\WYWAR7P^43F p V P V P 8XdK Qh R#)rrr N)r2r3r4r rrr)rertuplerIrPrrOr r3rr'rEr%rHrr r ) rQr rRrSr4r2shpos_axisrr3rTaxb1s &&& rU test_axisTestBSpline.test_axiss11 KK1aeai ( !8 ( 2Yii##D) JJszzrz* + A! 'ssyyR\Ob(m;b!oMMMM ZZ "u{{b(mbhh6QJKHHHHFF7Q;'B )W : 1Q8 :( 1.99;1.99!<1.==?1.==a@BB77aff$ $$ Brkc $^pVP.RO4pVP.RO.RO.4p\W4VRR7p\W4R,V4p\W4R ,V4p\V!R4VPV!R4V!R4.44R#) r=r r8N)rr.r=r?r@rrr;)r;r=rr;)r=rr.)rNNN)r.rE)rOr rr)rQrRr4r2r3splspl0spl1s&& rU test_neg_axisTestBSpline.test_neg_axiss{  JJ, - JJ 6 7aAB'qD'1%qD'1%C"((DItCy+A"BCrkc Rp^ p^pRF pV!W#V4K \^^^4F pV!W#R4K ^p\^^4F pV!W#R4K R#)a Splines with different boundary conditions are built on different types of vectors of knots. As far as design matrix depends only on vector of knots, `k` and `x` it is useful to make tests for different boundary conditions (and as following different vectors of knots). c\PPR4p\P!VP V4^(,^, 4pVP V4^(,^, pVR8Xd VR,V^&\ WEWR7p\P !\VP4V, ^, 4p\VPWq4!V4p\P!VVPV4P4p \WP,VRR7\WRR7R#)zA To avoid repetition of code the following function is provided. rrr4bc_typersrtNr;)rIrPrr random_sampler eyelenr2r design_matrixtoarrayrr3) rSr4rNrrwrrr3 des_matr_def des_matr_csrs &&& rUrun_design_matrix_testsHTestBSpline.test_design_matrix_bc_types..run_design_matrix_testss))''-C))!,r1B67A!!!$r)B.A*$u!%aaADs466{Q*+A"466103L"00151244;GI  L66115 A LU Crk not-a-knotrNclampednatural)r)rQrVrSr4bcs& rUtest_design_matrix_bc_types'TestBSpline.test_design_matrix_bc_typessg D*  (B #A" -)q!QA #A, 7  q!A #A* 5rkrdegreec  X\PPR4pVP^ V^,,4p\P!V4\P !V4reTp\P \P!V^, V^, V4\P!WV^V^,,4\P!V^,V^,V43,p\P!\V4V, ^, 4p \WWq4p \V !V4\P!WHWq4P44\P!V^ , V^, VR,V^ ,.4pV'gE\P !\"4;_uu_4\P!WHWq4RRR4R#\V !V4\P!WHWq4P44R# +'giR#;i)z;Test that design_matrix(x) is equivalent to BSpline(..)(x).rrN)rIrPrrOaminamaxrrerPrQr rrRrSarrayr`r rG) rQrr_rrwxminxmaxr4r2r3bsplines &&& rU'test_design_matrix_same_as_BSpline_call3TestBSpline.test_design_matrix_same_as_BSpline_call%s~ii##D)   bFQJ/ 0WWQZd  EE"++dQhq&9++d!vz*:;++dQhq&9: ; FF3q6A:> "!/ AJ--aACKKM HHdRi4#:tbyA Bz**%%aA;+*  %%aA;CCE +**s ?H H) c\PPR4p^ p^p\P!VP V4^(,^, 4pVP V4^(,^, p\ WEVR7p\ ^^4F[pVRVpVRVp \P!VVPV4P4p \WP,V RR7K] R#)rr4Nrsrt) rIrPrrrOr rr rRr2rSrr3) rQrrSr4rwrrrxcycrUs & rUtest_design_matrix_x_shapes'TestBSpline.test_design_matrix_x_shapesAsii##D)   GGC%%a(2-2 3   a 2 % *!!!,q!A2AB2AB"00151244;GI  L6612E B rkc|.ROp\P!RV^4P4p\V.RO.RR7R#)r0rBrsrtN)r0r0r0rBrCrDrDrD)g?gmਪ?gK}\UU?rA)r rRrSr)rQr2des_matrs& rUtest_design_matrix_t_shapes'TestBSpline.test_design_matrix_t_shapesRs7 ,((Q2::<;<" $rkc p\PPR4p^ p^p\P!VP V4^(,^, 4pVP V4^(,^, p\ WEVR7p\ \4;_uu_4\P!WFPRRR1,V4RRR4^p.ROp.ROp\ \4;_uu_4\P!WGV4RRR4R# +'giLS;i +'giR#;i)rrjNr;)rAr0rBrCrD@r) rIrPrrrOr rErGr rRr2)rQrrSr4rwrrr2s& rUtest_design_matrix_asserts&TestBSpline.test_design_matrix_assertsZsii##D)   GGC%%a(2-2 3   a 2 % *!!!,: & &  ! !!VVDbD\1 5'  $  : & &  ! !! *' & ' & ' & &s,D/D$ D! $ D5 rNc\PPR4p\P!VP^44pVP^4pVR8Xd V^,VR&\ W4VR7p\ P !WQR7p\P!^^^4p\V!V4V!V4RR7\W4VR7p\VPVPRR7R#)rrrNrYrtNr;) rIrPrrrr from_power_basisrerr r3) rQrNrrwrcbrrgbspl_news && rUtest_from_power_basis!TestBSpline.test_from_power_basislsii##D) GGCJJrN # JJrN j aDAbE w /''< [[Ar "2Ru5%aG< 7rkc(\PPR4p\P!VP^44pVP^4VP^4R,,pVR8Xd V^,VR&\ W4VR7p\ P !WQR7p\W4PVR7p\W4PVR7p\VPVPRVP,,RR7R#)rr/rrxrYrtNr;) rIrPrrrr ryr rrrr3) rQrNrrwrrzr bspl_new_real bspl_new_imags && rUtest_from_power_basis_complex)TestBSpline.test_from_power_basis_complex|sii##D) GGCJJrN # JJrNSZZ^b0 0 j aDAbE w /''<*1ffgF *1ffgF  "}2F FUSrkc \P!.RO4p\P!.RO4p\P!\ WRR7RR7p\ VP .RORR7R#) a5 For x = [0, 1, 2, 3, 4] and y = [1, 1, 1, 1, 1] the coefficients of Cubic Spline in the power basis: $[[0, 0, 0, 0, 0],\$ $[0, 0, 0, 0, 0],\$ $[0, 0, 0, 0, 0],\$ $[1, 1, 1, 1, 1]]$ It could be shown explicitly that coefficients of the interpolating function in B-spline basis are c = [1, 1, 1, 1, 1, 1, 1] r[rxrYrtNr>)r.r.r.r.r.)r0r.r.r.r.r.r.)rIrcr ryrrr3)rQrwrrs& rUtest_from_power_basis_exmp&TestBSpline.test_from_power_basis_exmpsM HH_ % HH_ %'' A)(L09; 7eDrkct\P!^^.4p\P!R.4pVPRR7VPRR7\P!^^^ 4pVPRR7\ W^R7p\ V!V4\P !V4R,4R#)rrCFwriter1N)rIrcsetflagsrer rrn)rQr2r3rgrTs& rUtest_read_onlyTestBSpline.test_read_onlys HHaV  HHcUO     [[Ar " %  a ""r||B/#56rkcDa\SR7pV3Rlp\^ W24R#)rrc<VPwr pSPW#,W#)^, ,R4pV!V4R#)r.i'N)r_re)rrTr2r4rgrRs&& rU worker_fn/TestBSpline.test_concurrency..worker_fns3eeGA!QT1RT7E2B bErkN)r^r&)rQrRrTrs&f rUtest_concurrencyTestBSpline.test_concurrencys 2 &  I1rkcygwinz'threading.get_native_id not implemented)rr c t\4p\P!^^^ 4pV!V4p\P!4p\P !\ VPRV R244RVPPVPPR7pVPVR&\P !\ VPRV R244RVPPVPPR7pVPVR&WbnWrn \V!V4V4R#)rr2z.datzw+)moder:rrEr3N) r^rIre threading get_native_idmemmapstrjoinr2r:rr3r)rQtmpdirrTrgrhtidt_mmc_mms&& rU test_memmapTestBSpline.test_memmaps ! [[Ar "R5%%'yyV[[1SE78t ssyy ;##QyyV[[1SE78t ssyy ;##Q"x(rkr.Nr)FTr)r[rZrrX):r/r0r1r2rVrbriror|rskip_xp_backendsrrrrrrrrrrrrrrrrrrrrrr"r)r6r`mark parametrizerrArIr]rgrmrqrur|rrrrxfailsysplatformrarr3__classdictcell__ __classdict__s@rUr,r,1s!F )>D @l+OPk*KLgk2 /3MQ /A ;1 .E&W? ( 6? FG4A M :# F 7+BZ 1!1 [[VU2q\2%3%:D(6T [[],EF [[XuQx01G4C"$+$ [[Y)CD 8D 8 [[Y)CD TD TE& 7 2 [[  8 )  )rkr,cRa]tRtRto]P P R.RO4R4t]P P R.RO4R4tRt ]P P R.RO4R4t ]P P R R R .4R 4t R t Rt RtVtR #) TestInsertixvalc VP^VPR7pVPV4^,p\W4^R7pVPVP VP 3p\\W4!pVPV4p\VPVPVP4RR7\VPVPVPRVP )^, 4RR7WR ,8wdTMVRR p VPV RVR,VRR ,,34p \V!V 4V!V 4RR7VPV4^,p \W:^R7p \W2P!WJ3^R7^R7p V PV4p \V PVPRR7\V PVP!VPV4PV PV4P3^R7RR7WR ,8wdTMVRR p VPV RVR,VRR ,,34p \V !V 4V !V 4RR7R#) r;r9rjrYrtNrr.NNr r;)rMrNsinr _t_cr4r r insert_knotrr2rOr3concatcosr)rQrrRrwrrFr_spl_1fspl_1rgy1spl_y1spl_yyspl_yy1s&&& rU test_insertTestInsert.test_inserts IIarzzI * FF1IqL +vvsvvsuu%&+,%FHH!5EBFHHYsuufQh,?!@uMB%-QQsV YYC23B01 2Br7VVAY\#AQ/#Axxax'@AF$$T* 5777 II HHcood+--v/A/A$/G/I/IJQRH S B%-QQsV YYC23B01 2r GBKe!DE(FHH!5EBFHHYsuufQh,?!@uMB%-QQsV YYC23B01 2Br7rkcz\PPR4p^ ^rCVP\P!VP W4,^,R744pVPVP V^^3R74p\ WVV4pVPVP WT^,,WT)^, ,R74pVPV4p VPVP WT^,,WT)^, ,^!R74p \V!V 4V !V 4RR7R#)90r)lowhighrrrrYrtN) rIrP default_rngrOruniformr rr) rQrRrrSr4r2r3rFxvrrgs && rUtest_insert_randomTestInsert.test_insert_randomsii##E*11 JJrwws{{A{67 8 JJs{{Aq {2 3aA ZZ A#Qr!tW = ># ZZ A#Qr!tW2 F GBr7rkrc VP^VPR7pVPV4^,p\W4^R7wrVp\ VP WV34wrV\ WVVRR7pVPV4p \WPRR7wrp\V PVP V 4RR7\V PR V)^, VP V R V)^, 4RR7\PPR 4P!^^^)R 7p VP V 4p \V !V 4VP \#WW344RR7R #) r;r9rjrrTperrYrtNrr)rMrNrrrvrOr rrr_rr2r3rIrPrrr)rQrrRrwrr2r3r4rFrtfcfrrgs&&& rUtest_insert_periodicTestInsert.test_insert_periodics IIarzzI * FF1IqL#a2::v&aA:6#2wwD1 Be<!A 2er!t9(=EJ %%d+333K ZZ b 2::eE;.G#HuUrkrNrcVP^VPR7^,\P,pVP V4VP V4rT\ W4RV,,^R7pWn\ W4^R7pWn\ W5^R7pWnRp VPV 4p VPV 4p VPV 4p \V PV PRR7\V PV PRR7\V PV PRV P,,RR7R#)r;r9r/rj @rYrtN) rMrNrIpirrr rrrr2r3) rQrrRrwy_rey_imrFspl_respl_imrrspl_1respl_1ims &&& rU test_complexTestInsert.test_complex0s IIarzzI *Q . 6VVAYq d 2d7Na8%#Aq1(#Aq1( #$$R($$R(77R \!9FrkcRp^p\P!^.V^,,.RO,^.V^,,,4p\P!\V4V, ^, 4p\ W4VRR7p\ \ 4;_uu_4\WWB3RR7RRR4\ \ 4;_uu_4VPV4RRR4R# +'giLC;i +'giR#;i)rrrTrNr=r?r@r) rIrconesrQr rErGrr)rQrr4r2r3rFs& rU+test_insert_periodic_too_few_internal_knots6TestInsert.test_insert_periodic_too_few_internal_knotsGs  HHaS!A#Y-QqS 9 : GGCFQJN #aA:6 : & & 21yd +': & & OOB ' &' &' & &sC& C9& C6 9 D c^p\P!^.V^,,.RO,^.V^,,,4p\P!\V4V, ^, 4p\ W#V4p\ \ 4;_uu_4VPR4RRR4\ \ 4;_uu_4VP^4RRR4\ \ 4;_uu_4VP^^R7RRR4R# +'giLv;i +'giLX;i +'giR#;i)r?Nrrr;)rIrcrrQr rErGr)rQr4r2r3rFs& rUtest_insert_no_extrap TestInsert.test_insert_no_extrap[s  HHaS!A#Y-QqS 9 : GGCFQJN #aA : & & OOB ': & & OOA ': & & OOAO #' & ' &' &' & &s$DD*9D= D' * D: = E r.)rAr0r8r@@@))rAr=)r0r?)rr)r@r=)rr=) r皙?rBrD@@@gffffff@r)r/r0r1r2r`rrrrrrrrrr3rrs@rUrrs [[V%AB!=C!=H [[C 8 8 8 [[T$8VV" [[]T:,>?G@G, ( $ $rkrc RRlpRFYp\VR7p\\V44F2wr4V!WC4\^V^,4FpV!WCVRR4K K4 K[ R#)rc VVPwrVp\P!V4p\PV^,R, RVR,VR,,,VR ,R,3,p\ \ WWg3V4V!W4W4RV RVP 2R7R#) rrrr:Nr.Nzder = z k = )r[r\err_msgNr;)r_rIuniquerrrr4) rTrrr[r\r2r3r4rws &&&&& rU check_splev,test_knots_multiplicity..check_splevos%%a IIaL EE!A$s(C22/2s: ;aQC0!A)se6!##.G IrkrjrN)rrsrsr.r=r?r@r)r^ enumerate_make_multiplesr)rr4rTrr@rs rUtest_knots_multiplicityrks^I ! $q12EA  Q!}B3u5%3rkctV^8Xd(W2,Tu;8:dW2^,,8dR#R#W2V,,W2,8XdRpMHWV,, W2V,,W2,, , \W^, W#4,pW2V,^,,W2^,,8Xd RpWE,#W2V,^,,V, W2V,^,,W2^,,, , \W^, V^,V4,pWE,#)zk Naive way to compute B-spline basis functions. Useful only for testing! computes B(x; t[i],..., t[i+k+1]) r0rA_naive_Brwr4rr2c1c2s&&&& rUrrs  Avda(!aC&(s1(c11v~ A$hQ3!$ '(1c1*@ @1Qx1qS6  G!AhlQs1uXA#. /(1c1Q32J J GrkcVaaaaaSSS,8XdSoMVPSS4^, oSS,Su;8:dSS^,,8:gQhQhSS8d SSP^,S, 8gQh\VVVVV3Rl\^S^,444#)z5 Naive B-spline evaluation. Useful only for testing! c3x<"TF/pSSV, ,\SSSV, S4,xK1 R#5iNr).0rr3rr4r2rws& rU _naive_eval..s0F 1q1vAqsA... s7:) searchsortedrsumr)rwr2r3r4rRrsffff$@rUrrs AaDy  OOAq !A % Q41 !A#    6a!''!*q.(( ( FFa1 F FFrkc&aaaaSP^,S^,, pVS^,8gQhSP^,V8gQhSS,Su;8:dSV,8:gQhQh\VVVV3Rl\V444#)z'Naive B-spline evaluation, another way.c3\<"TF!pSV,\SSVS4,xK# R#5irrrrr3r4r2rws& rUr _naive_eval_2..s(=Hqqthq!Q***H),)rrr)rwr2r3r4rRrSsffff$ rUrrsr  acA !8O8 771:?? Q41 !    =E!H= ==rkc`\V4V^,, pWC^,8gQhVP^,V8gQhRp\V4FZp\P!WWc,^,RR7!V4pWRV,\ P !V4,, pK\ V#)r.rAFr)rQrrr rrI nan_to_num)rwr2r3r4rSsrrTs&&&& rUrrs A!A#A !8O8 771:?? A 1X  ! !!ac!e*% @ C qTBMM!$ $$ Hrkc\P!V4p\P!W^8V^8,V^8V^8,V^8V^8*,.RRR.4pVPV4#)z*A linear B-spline function B(x | 0, 1, 2).cR#rAr.rws&rUB_012..sbrkcV#rr.rs&rUrrsArkcRV, #rBr.rs&rUrrsArk)rI atleast_1d piecewiserO)rwrRresults&& rUrurusp aA \\!1uQ/!VA.!VQ/1)+~FHF ::f rkc\P!V4pV^8V^8V^8,V^8.pV^8Xd RRR.pMV^8Xd RRR.pM\RV 24h\P!WV4pV#)z0A quadratic B-spline function B(x | 0, 1, 2, 3).c W,R, #rr.rs&rUrB_0123..s 13r6rkc0RVR, ^,, #)rCg?rr.rs&rUrrs41T6A+-rkc0RV, ^,^, #)rCr.rs&rUrrsBqD19q=rkcR#r0r.rs&rUrr2rkcR#)rBgr.rs&rUrrs3rkcR#rr.rs&rUrrrrkznever be here: der=)rIrrGr)rwrcondsfuncspiecess&& rUrrs aA UQUq1u%q1u -E ax!-(* .se455 \\!E *F Mrkc6\PP^{4p\P!VPW,^,44pVPV4pVP V4VP V4rT\ P !WEV4#){)rIrPrrrOr construct_fast)rSr4rRrr2r3s&&& rUr^r^sj ))   $C  13q5!"A 1 A ::a="**Q-q  ! !! **rkc#"VPVPr!VPP4pV^,VR&V^,V^&\ W1V4xVPP4pV^,VRV^,%\ W1V4xVPP4pVR,W2)^, R%\ W1V4xR#5i)zIncrease knot multiplicity.:NNr;)r3r4r2copyr )rTr3r4t1s& rUrrs 33q B2BuI VBrF "  B!uBt!H "  B2Br!tuI " sCCcba]tRtRtoRtRtRtRtRtRt Rt R t R t R t R tR tVtR#) TestInteropic\P!^R\P,^)4p\P!V4p\ W4pVP VP VP3VnWVuVn Vn Vn \P!^R\P,^4Vn \PVP VP VP 3,p\P!WD34Vn\!VP VPVP4VnR#)rrDN)rIrerrr r2r3r4r_rgrrTxnewc_dstackrr b2)rQrgrrTrs& rU setup_methodTestInterop.setup_methods [[BruuHb ) VVBZ r &CCacc?#%1 $&KK2bee8R0 UU133QSS= !))RH%!##tww,rkc .VPVPVPr2p\\ W4V!V4RRR7\\ WP 4V!V4RRR7\\ P!VUu.uFp\ WB4NK up4V!V4RRR7\\RR7;_uu_4\ W4RRR4\\^VPP44R,pVPPV4pVPWcP 3p\\ P!\ W44V!V4PV4RRR7R#uupi +'giL;i)rYrZzCalling splev.. with BSplinematchNr)r/rTr2rrr_rIrOrErGr<rr3r' transposer2r4)rQr/rTr2rwr=rr_s& rU test_splevTestInterop.test_splevs"ii d$e% 9dEE*$e% 9 #>AE!K#>?$e% 9 :-K L L $OM5BDDII& '$ . TT^^B ttR 5#344**2.U H$? M LsE?  F F cVPVPr!\W4p\P!W4wrEp\ V^,VRR7\ V^,VRR7V^,V8XgQh\WRR7wpp\ V^,VRR7\ V^,VRR7V^,V8XgQh\ W4p \ W)RR7\ V!p \ W*!V4RR7R#)rrYrtT) full_outputN)rgrrrrrr ) rQrwrr_r2r3r4tck_frrrTs & rU test_splrepTestInterop.test_splrepsww1Ql,,q$aA.A.1v{{ $7q!Qa!%0a!%0Qx1}}1]E* SM1Q4e,rkcVPVPr!\PW"3,p\ \ 4;_uu_4\ W4RRR4\ \ 4;_uu_4\P !W4RRR4\ \RR7;_uu_4\ VR,VR,4RRR4\ \RR7;_uu_4\P !VR,VR,4RRR4R# +'giL;i +'giL;i +'giL;i +'giR#;i)Nm > k must holdr6:Nr?N) rgrrIr0rErGrrrF)rQrwry2s& rUtest_splrep_errorsTestInterop.test_splrep_errors/sww1 UU14[ : & & 1M' : & & LL '9,= > > 1R5!B% ? 9,= > > LL2" &? >' & & &? > > > >s0 D2D)*D<&&E D& ) D9 < E  E c\P!^\PR7PR4p\ V4wr#\ P!V4wrE\ W5RR7\ \P!\W244VRR7\ \P!\W444VRR7\ V^RR7wwrgp\ W7RR7\ \P!\Wv44VRR7R#)r9rYrtT)r r=Nr?r) rIrMrNrrrrrOr) rQrwrTur_u1b_fu_frs & rU test_splprepTestInterop.test_splprep?s IIb + 3 3F ;qz--" E* 5;/? 5=115A&a1$? Aq!U+ 5?3QUCrkcp\P!^<4PR4p\\RR7;_uu_4\ V4RRR4\\RR7;_uu_4\ P !V4RRR4\P!^^(^R7p\\RR7;_uu_4\ V.4RRR4\\RR7;_uu_4\ P !V.4RRR4.ROp\\RR7;_uu_4\ V.4RRR4\\RR7;_uu_4\ P !V.4RRR4.R Op.R Op\\RR7;_uu_4\ V.RV.!RRR4R# +'giELm;i +'giELI;i +'giEL;i +'giL;i +'giL;i +'giL;i +'giR#;i) r?ztoo many values to unpackr6N)numrBzInvalid inputsr) >IrQ >KrR)r.r?r=r@)rg333333?g?r.) rIrMrrErGrrrerF)rQrwrIs& rUtest_splprep_errorsTestInterop.test_splprep_errorsNsd IIe  $ $Y / :-H I I AJJ :-H I I MM! J KK21 % 9,= > > QCL? 9,= > > MM1# ? E :-= > > QCL? :-= > > MM1# ?   :-= > > qc4^ $? >/J I I I I I ? > > > > ? > > > ? > >sS F/,G= G+G+( G>HH$/ G  G  G( + G; > H  H! $ H5 cVPVPr!\P!.R O4\P,p\ \ V4VRRR7\ \ VPVPVP34VRRR7\\RR7;_uu_4\ V^2R7RRR4VPP^^^4p\P!\ VPWBP3^2R74pVPR 8XgQh\ WS, \P!V4RR7R# +'giL;i) rrrZzCalling sproot.. with BSpliner6)mestNrrt)rrr8r)r?r=r@)rTr2rIrcrrrr2r3r4rErGr9rOrr)rQrTr2rootsc2rrrs& rU test_sprootTestInterop.test_sprootks2-.ruu4q 5t$?QSS!##/TM:-L M M 2B NddnnQ1% ZZc440r: ;xx9$$$ BMM"$5EBN Ms /E E& c lVPVPr!\\^^V4\^^VP4RRR7\\^^V4VP ^^4RRR7\ \RR7;_uu_4\^^V4RRR4VPP^^^4p\P!\^^VPW2P344pVPR8XgQh\V\^^V4RRR7R# +'giL;i) rrsF)r[check_0dzCalling splint.. with BSpliner6N)r[ check_shape)r?r=)rTr2rrr_rrErGr3r9rIrOr2r4r)rQrTr2rXintegrs& rU test_splintTestInterop.test_splint|s2q!Qq!QUU+%% Iq!Q Aq) G:-L M M 1a NddnnQ1%F1a"$$TT):;<||v%%%q!Qe HN Ms D## D3 c VPVP3EFpp\VP4\VP4, pVPP 4pV^8dG\ PV\ P!V3VPR,,43,pRFp\V4p\P!VPP 4W1P34p\VPV^,RR7\VPV^,RR7VPV^,8XgQh\V\4'gQh\V\ 4'dKQh EKs R#rrrYrtNr.r=r?)rTr2rQr2r3r*rIrzerosrrrr4rr%r r<rQrTctb_crSbdtck_ds& rU test_splderTestInterop.test_splders &&$''"AQSSCH$B##((*CAveeC2%#))B-*?!@@AAY acchhj#ss%;<eAhU;eAhU;ttuQx'''!"g....!%////#rkc VPVP3EFpp\VP4\VP4, pVPP 4pV^8dG\ PV\ P!V3VPR,,43,pRFp\V4p\P!VPP 4W1P34p\VPV^,RR7\VPV^,RR7VPV^,8XgQh\V\4'gQh\V\ 4'dKQh EKs R#rc)rTr2rQr2r3r*rIrrerrrr4rr%r r<rfs& rUtest_splantiderTestInterop.test_splantiders&&$''"AQSSCH$B##((*CAveeC2%#))B-*?!@@A]((!##((*c33)?@eAhU;eAhU;ttuQx'''!"g....!%////#rkcVPVPVPr2pVPP^,pRVPV,VPV^,,,,p\ WQ4\ WQPVP VP34rv\\W64\W74RR7\V\4'gQh\V\4'gQh\\VP P44pVP PVR,R,4p \ WRPWP34p \ WR4p \\ P"!\W:44P^^^4V !V4RR7\V \4'gQh\V \4'gQhR#)r=rrYrtrNr8)rTr2rgr2rrr3r4rrr%r r<rr'r9rIrO) rQrTr2rgrtnbntck_nr=r0tck_n2bn2s & rUrTestInterop.test_insertsXFFDGGTWWr CCHHM !##a&133qs8# $2M6"ssACCo#>Eb b(u 6"g&&&&%''''5# $ TT^^BrFTM *TT2tt,-Rn  5#45??1aHBe -#w''''&%((((rk)rTr2rr_r/rgrN)r/r0r1r2r3r:r?rDrMrSrZr`rkrnrr3rrs@rUr-r-sE -H6-0' D%:C"H(0"0"))rkr-ca]tRtRto]P !RR]P ,4t]P!]4t Rt Rt Rt Rt ]PP!R.R*O4R 4t]PP!R.R*O4R 4tR tR t]PP!R.R+O4R 4tRtRtRt]PP!R.R,O4R4tRtRtRtRtRt]PP?RR7R4t Rt!Rt"Rt#Rt$Rt%Rt&R t'R!t(]PP!R.R-O4R"4t)R#t*R$t+R%t,R&t-R't.R(t/Vt0R)#). TestInterpirArBcnVPVP4VPVP43#r)rOrgr)rQrRs&&rU_get_xyTestInterp._get_xys'zz$''"BJJtww$777rkc\\4;_uu_4\VPVPRR7RRR4R# +'giR#;i)r8rjN)rErFr rgrrQs&rUtest_non_int_orderTestInterp.test_non_int_orders0 9 % % tww3 7& % % %s #A A cVPV4wr#\W#^R7p\V!V4VRRR7\W#^RR7p\V!V4VRRR7R#)rrjrsrZr4r Nr;rzr rrQrRrgrrTs&& rU test_order_0TestInterp.test_order_0Pb! r +"rE: r 4"rE:rkcVPV4wr#\W#^R7p\V!V4VRRR7\W#^RR7p\V!V4VRRR7R#)r.rjrsrZrNr;rrs&& rU test_linearTestInterp.test_linearrrkr4c.ROp.ROp\\RR7;_uu_4\W#VR7RRR4R# +'giR#;i)rz Shapes of xr6rjNrr.r=r?r@r)rr.r=r?r@rrrrErGr rQr4rwrs&& rUtest_incompatible_x_y TestInterp.test_incompatible_x_ys2  $ :] ; ; qq )< ; ; ;s : A c.ROp.ROp\\RR7;_uu_4\W#VR7RRR4.ROp\\RR7;_uu_4\W#VR7RRR4.ROp\P!V4P R 4p\\RR7;_uu_4\W#VR7RRR4R# +'giL;i +'giL};i +'giR#;i) rzx to not have duplicatesr6rjNzExpect x to be a 1D strictly)rr.r.r=r?r@r)rr=r.r?r@rr.r;)rErGr rIrOrrs&& rU test_broken_xTestInterp.test_broken_xs   :-G H H qq )I  :-K L L qq )M  JJqM ! !' * :-K L L qq )M LI HM L M L Ls#CC-C+ C  C( + C< cVPV4wr#RF$p\W#V4p\V!V4VRRR7K& R#)r=rsrZNr=r?r@rrrr)rQrRrgrr4rTs&& rUtest_not_a_knotTestInterp.test_not_a_knots:b!#A"21-A AbE2E >$rkc VPV4wr#\W#^RR7p\V!V4VRRR7\^^4F/p\V!V^,VR7V!VR ,VR7RR7K1 \W#^RR R7p\V!V4VRRR7\^^4F/p\V!V^,VR7V!VR ,VR7RR7K1 R #) rrrMrsrZrdy=rtr4rNr Nr;)rzr rr)rQrRrgrrTrs&& rU test_periodicTestInterp.test_periodicsb! rJ ?"rE:q!A AbeNAbfO% H rJR H"rE:q!A AbeNAbfO% Hrkcn^p\PPR4p\P!VP V4^ ,4pVP V4^d,pVR,V^&VP V4VP V4re\ WVVRR7p\V!V4VRR7R#)rrrrMrsrtNr;)rIrPrrrOrOr r)rQr4rRrSrrwrrTs&&& rUtest_periodic_randomTestInterp.test_periodic_randoms ii##D) GGC%%a(2- .   a 3 &u!zz!}bjjm1 qq* =!ae,rkcVPP^,p\PP R4pVP V4^,\P ,p\P!V4pRV^&^\P ,VR &\P!^V34p\P!V4V^&\P!V4V^&VPV4VPV4rT\WE^R^R7p\V4F&p\V!WG,4VRV3,RR7K( \V!V^,4V!VR ,4RR7R#) rrrArrrErsrtNr;)rgrrIrPrrOrrrerrrOr rr)rQrRrSrrwrrTrs&& rUtest_periodic_axisTestInterp.test_periodic_axis$s  GGMM! ii##D)   a 1 $ruu , GGAJ!BEE " HHaV vvay!vvay!zz!}bjjm1 qq*1 EqA AadGQq!tW5 9!A$1R56rkc b\PPR4p^p^p\P!VP V44pVP V4pVR,^, V^&\ \ 4;_uu_4\WEVRR7RRR4R# +'giR#;i)rrrMNr;)rIrPrrrOrErGr )rQrr4rSrwrs& rUtest_periodic_points_exception)TestInterp.test_periodic_points_exception5sii##D)   GGC%%a( )   a uqy! : & & qq* =' & & &s B B. c \PPR4p^p^p\P!VP V44pVP V4p\P !V^V,,4p\ \4;_uu_4\WEW&R4RRR4R# +'giR#;i)rrN) rIrPrrrOrerErGr )rQrr4rSrwrr2s& rUtest_periodic_knots_exception(TestInterp.test_periodic_knots_exception@sii##D)   GGC%%a( )   a  HHQQY  : & & qQ: 6' & & &s B-- B> c \VPVPVRR7p\VPVPRVR7p\ VPV4p\ WB!VP4RR7\ ^V4F9p\ VPW5R7p\ WB!VPVR7RR7K; R #) rrMT)rr4rsrtrrrN)r rgrrrrr)rQr4rTr_rFrs&& rUtest_periodic_splevTestInterp.test_periodic_splevKs tww1j ITWWdgg415DGGS!QtwwZe4q!A,C C477q!1 >rkc8\VPVP^RR7p\VPVPRR7p\ V!VP4V!VP4RR7\ P PR4p^p\ P!VPV4^ ,4pVPV4^d,pVR,V^&\WV^RR7p\WVRR7p\ V!V4V!V4RR7R#) r?rrMrxrsrtrNr;) r rgrrrrIrPrrrO)rQrTcubrrSrwrs& rUtest_periodic_cubicTestInterp.test_periodic_cubicXs tww1j I$''477J?$'' CLu=ii##D)  GGC%%a(2- .   a 3 &u! qq* =! 3!c!f51rkcZaaa^o\VPVPSRR7p\VPS4o\ VPVPSS4o\ P !VVV3Rl4p\V!VP4V!VP4RR7R#)r?rrMc,<\VSSS\R7#)rr)rrI)rwr3r4r2s&rUr6TestInterp.test_periodic_full_matrix..osK1ar$BrkrsrtN)r rgrrrrI vectorizer)rQrTr@r3r4r2s& @@@rUtest_periodic_full_matrix$TestInterp.test_periodic_full_matrixhsu  tww1j I DGGQ ' &twwA > \\B C$'' BtwwKe ?': & & qW 5': & & qgY 7': & & qR 0' 1 : & & qaV 4' &' &' &' &' & ' & &s;C1 D:D(D*D=1 D  D  D' * D: = E c\VP^4pV^,pR.R.rT\\RR7;_uu_4\W#WE3R7RRR4R.R.rT\\RR7;_uu_4\W#WE3R7RRR4R# +'giLG;i +'giR#;i) rzBad boundary conditions at 0.r6rxNzBad boundary conditions at 6.)rrr.r)ir)rMrErGr )rQrRrwrrrs&& rUtest_deriv_order_too_large%TestInterp.test_deriv_order_too_larges IIaL qDx&1 :-L M M qaV 4Nx'1 :-L M M qaV 4N M N M N M MsB.B B  B+ c^pVPV4wr4VRV,,pR.R .re\W4W%V3R7p\V!V4VRRR7\P!V!V^,^4V^,^,RR7'gQh\P!V!VR ,^4V^,^,RR7'gQhR F%p\W4VR7p\V!V4VRRR7K' R#) r?r/rxrsrZrrjN)r.y@)r.y@@r;rr.)rzr rcmathr)rQrRr4rgrrrrTs&& rUrTestInterp.test_complexs b! #b&[!zK=u rq%. A"rE:}}Qr!ua[%(1+uEEEE}}Qr"vq\58A;FFFFA"2Q/A AbE2E >rkcVP^ VPR7pVP^ VPR7pRFp\W#VR7pV!V4K R#) r9rjNr)rMint32r )rQrRrwrr4rTs&& rU test_int_xyTestInterp.test_int_xysM IIbI ) IIbI )A"11-A aDrkc|VPR^^d4pVR,pVR,pRFp\W4VR7K R#)r.:NNrrjNr;r)rer )rQrRrgrwrr4s&& rUtest_sliced_inputTestInterp.test_sliced_inputs: [[Q $ sG sGA qq )rkcVP^ VPR7pV^,pVPVPVP)3F>p\P !VR4P V4p\\\W#4K@ R#)rr9Nr;) rMrNrKrLxpxatsetrErGr )rQrRrwrzs&& rUtest_check_finiteTestInterp.test_check_finite"sg IIb I + qD&&"&&266'*Aq" !!!$A *&8! ?+rkc\\^ 44pVUu.uF q3^,NK pp\W$VR7R#uupi)rrjN)listrr )rQr4rwars&& rUtest_list_inputTestInterp.test_list_input+s8 rO 1aTT1 11% s;cVPV4wr#VPW#3^R7p^RR.3.p^RR.3.p\W#^WE3R7p\V!V4VRRR7\V!V^,^4VP V^,^,VP R 7RRR7\V!VR ,^4VP V^,^,VP R 7RRR7R #) r.r r0rBrCrDrMrsrZr9Nr;)rzrr rrOrN)rQrRrgrrrrTs&& rUtest_multiple_rhsTestInterp.test_multiple_rhs2sb! XXrhQX 'b"Xb"X rUN C"rE: beQK JJuQx{"**J 5E   bfaL JJuQx{"**J 5E rkc\PPR4p^^r2\P!VPVR74pVPV^^^3R7p\ WEV4pVP P V^^^38XgQh^VPR43.p^VPR43.p\ WEW'V3R7pVP P W2,^, ^^^38XgQhR#)rrrxNrrr)rIrPrrr r3r) rQrr4rSrwrrTd_ld_rs & rU test_shapesTestInterp.test_shapesCsii##D)"1 GGCJJAJ& ' JJQ1aLJ ) qQ 'ssyyQ1aL(((3::i()*3::i()* qQc ;ssyyQUQY1a0000rkc VPV4wr#VPV4p\W#^RR7p\W#^R.R.3R7p\VPVPRR7\W#^RR7p\W#^R.R.3R7p\VPVPRR7\W#^RR7p\W#^RR.3R7p\VPVPRR7\W#^RR7p\W#^RR7p\VPVPRR7\ \ 4;_uu_4\W#^R R7RRR4VPVPV4VPV43^R 7p^RR.3.p^RR.3.p\W#^Wg3R7p\W#^RR7p\VPVPRR7\PPR 4p^^r\P!VPV R 74p VPV ^^^3R 7p VPV 4VPV 4r^VPR43.p ^VPR43.p\WWV3R 7p\WV RR 7p\VPVPRR7R# +'giELy;i)r?r[rMrYrtrZNrArXtypor rrrxr)r[rZr)NrZrrYr)rzrr rr3rErGrrrIrPrrrOre)rQrRrgrr@r2rrrr4rSrwrrrs&& rUtest_string_aliasesTestInterp.test_string_aliasesRstb! VVBZ !Y ? !vh5I Jbdd/ !(>@ !*0F8(<>bdd/ !5F G !dXJ5G Hbdd/ !\ B !T :bdd/: & & rF ;'XXrvvbz266":.QX 7b"Xb"X !e^ D !(>@bdd/ii##D)"1 GGCJJAJ& ' JJQ1aLJ )zz!}bjjm1288I&'(288I&'( as < a ;bdd/1' & &s ?J>> K c\PPR4p^^rC\P!VPVR74pVPVR7p\ WS4p\ WVWs4pVP V4p\VP WVV34wrp \WW;4p \V PVRRR7R#)rrrsrZN) rIrPrrrmake_interp_full_matrrOrvr rr3) rQrRrr4rSrxy_npryrrwrr2rTs && rUtest_full_matrixTestInterp.test_full_matrixsii##D)!1wwszzqz)*zzqz!4# "4t 7 ZZ^bjj4t"45a qQ *Re%8rkc  D\PPR4p^p\^^ ^4EFp\ V^, ^, 4p\P !VP^V344p\^V^,4FpVRV)1VR13;;,\P !VP^W&, 344, uu&WVR1RV)13;;,\P !VP^W&, 344, uu&K VPWD34pWuRV1V)R13&VPWD34pWV)R1RV13&\P !W234p \\WD)^, R44FHwrjV ^8d \P!WZR7WRV 13&K+\P!WZR7WV R13&KJ VPV4p \\WWV4\PPW[4RR7EK R#)z Random elements in diagonal matrix with blocks in the left lower and right upper corners checking the implementation of Woodbury algorithm. rN)offsetrsrtr;) rIrPrrintdiagflatrerdiagonalrrlinalgsolve) rQrrSr4rrrurlldrrTs & rU test_woodburyTestInterp.test_woodburys ii##D) q"aA!a%1%F CJJ1v./A1fqj)#A2#qr' bkk#**aZ*@AA "crc' bkk#**aZ*@AA *V,-B#%gvgwx V,-B#%vgh ! A!%! R"@Aq5!{{17A!eH!{{17AeH B  1 A /ra@IIOOA1 ? ?#!rkr.Nrrr)r.r=r?r)1r/r0r1r2rIrerrgrrrzr~rrr`rrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr3rrs@rUrxrxs RBEE "B B88;; [[S,/*0*  [[S,/ *0 *? I  [[S"45 -6 -7" > 7 [[S,/ ?0 ?2 = ,;, T [[j) E* E>& ;5. 5?"*@ [[S,/&0&  " 130j 9??rkrxcVPVP8XgQhVPVPV,^,8XgQhVPp\P!WD3\PR7p\ V4FdpW,pWrV,8XdTpM\P !W'4^, p\ P!W#Wx4p WWhV, V^,13&Kf \P!WQ4p V #)zAssemble an spline order k with knots t to interpolate y(x) using full matrices. Not-a-knot BC only. This routine is here for testing only (even though it's functional). r9) rrIrerNrrr*evaluate_all_bsplslr ) rwrr2r4rSArrleftbbr3s &&&& rUrrs 66QVV   66QVVaZ!^ ## # A !rzz*A 1Xt Q4<D??1+a/D ' 'd 9 !!VDF]  A Hrkcz\\PWV34wrpVPpVPV, ^, p\P!WE3\P R7p\ V4FdpW,pWV,8XdTp M\P!W(4^, p \P!W#W4p WWyV, V ^,13&Kf \P!VPV4p \P!VPV4p \P!W4p WV 33#)z,Make the least-square spline, full matrices.r9)rvrIrOrrerNrrr*rdotTrr )rwrr2r4rrSrrrrrr4Yr3s&&&& rUmake_lsq_full_matrixrs"**qQi(GA! A  QA !rzz*A 1Xt Q4<D??1+a/D ' 'd 9 !!VDF]  qssAA qssAA A !f9rkmethodnorm-eqqrc ,a]tRtRto]P P R4t^ ^utt ]P!]P ]44t ]P ]4t ] !]P!] ^,] R,^4] 4t]R4t]R4tRt]R4t]R4tRt]R 4tR t]R 4t]R 4t]R 4t]R4t]R4t]P>PAR]!!]"!^^444R4t#Rt$Rt%Vt&R#)TestLSQircVPVPVPVP3wr#rE\ W#WE4wrg\ W#WEVR7p\ VPV4VPPVPV, ^, 38XgQhVwr\PPWRR7wp p \ VPV 4R#)r)rcondNr;) rwrr2r4rrrr3rrrIr lstsq) rQrrwrr2r4c0AYrTaarrrs && rU test_lstsqTestLSQ.test_lstsqsVVTVVTVVTVV3 a%aA1 A!v 6R ssyyQVVaZ!^----iioob2o6 Aq!R rkc .\VPVPVPVP34OVP N5wr4rVVP V4p\W4WVVR7p\W4WVWqR7p \VPV PRR7\VPV PRR7VP V P 8XgQhR#)r!wrrsrtN) rvrOrwrr2r4rnrrr3) rQrrRrwrr2r4r+rTb_ws &&& rU test_weightsTestLSQ.test_weightssHc"**tvvtvvtvv&>?GG a LLO A!v 6aAA=SUU/SUU/sscee||rkc b.\VPVPVPVP34OVP N5wr#rE\ PPR4PVP^,R7pVPV4p\W#WEVRR7p\W#WEVRR7p\W#WERR7p \VPVPRR7VPVPV PVP, 4R84'dQhR#) rrrr*rr!rsrtN)rvrOrwrr2r4rIrPrrrrrr3rabs) rQrRrwrr2r4r+b_neb_qrb_no_ws && rUtest_weights_sameTestLSQ.test_weights_samesGc"**tvvtvvtvv&>?GG a II ! !$ ' / /QWWQZ / @ JJqMqQQyAqQQt< qD9U366"&&DFF!23e;<<<<" $5;7;aQ%'';78::ss/ $7s + D*%D/c~.\VPVPVP34OVPNVP N5wr#rE^p\ PPWV3R7pVPV4p\W'W4RR7p\W'W4RR7p \VPV PRR7R#)r?rrr!rrYrtN) rvrOrwr2r4rSrIrPrrr3) rQrRrwr2r4rSr9rr2b_neqs && rUtest_multiple_rhs_3TestLSQ.test_multiple_rhs_3:sGc"**tvvtvv&67GGG a II  1)  , JJqMqQ$7a9=e4rkc.\VPVPVP34OVPN5wr4pVPVP R,4p\ W6WEVR7p\ W2PV4WEVR7p\ W2PV4WEVR7p \V!V4V!V4RV !V4,,RRR7R#)r0r!r/rYrZN?@) rvrOrwr2r4rrrrr) rQrrRrwr2r4rlrTrrs &&& rUrTestLSQ.test_complexCs=3rzzDFFDFF#34! HHUOQXXe_ahhuo6 eafn # WR[u=rkc\P!R^^d4pVR,pVR,p\V^4p\W4V^VR7R#)r.:NNr?rHNr;)rIrerr)rQrrgrwrr2s&& rUrTestLSQ.test_sliced_inputzsA[[Q $ sG sG AqMa1V4rkc &\P!^ 4P\4pV^,p\ V^4p\P \P \P )3FpWSR&\\\W#WAR7K R#) r!Nr;) rIrMrIrNrrKrLrErGr)rQrrwrr2rs&& rUtest_checkfiniteTestLSQ.test_checkfinitesd IIbM  ' qD AqM&&"&&266'*AbE *oqQ N+rkcVPVPVPrCpVPRR7VPRR7VPRR7\ W#WAR7R#)Fr)rwrr2rN)rwrr2rrrJs&& rUrTestLSQ.test_read_onlysO&&$&&$&&a      !A5rkr4c4VPVPr2\\P!V^,VR ,^4V4p\ W#WARR7p\ W#WARR7pVR,VRR ,R, p\ V!V4V!V4RR7R#) rrrHrrNrBrYrtr;)rwrrrIrerr)rQr4rwrr2 spl_norm_eqspl_qrrgs&& rUtest_qr_vs_norm_eqTestLSQ.test_qr_vs_norm_eqs|vvtvv1 BKK!aeQ/ 3%aA9E  qd;eafn # B%@rkc\P!VP^4p\P!VP^4p\ VPVPVP ^RR7p\ WVP ^RR7pVR,VRR,R, p\ V!V4V!V4RR7R#) r=rrHrNrBrYrtr;)rIrepeatrwrrr2r)rQrwrrspl_2rgs& rUtest_duplicatesTestLSQ.test_duplicatess IIdffa  IIdffa !DIdff$?eafn #b 5959rkr.Nr;)'r/r0r1r2rIrPrrrSr4rrwrrrer2parametrize_lsq_methodsr'r-r4rr;r?rrErrRrrXrr`rrrrr_rdr3rrs@rUrrs ))   %C qDAq  1 A 1 A AaD!B%+Q/A ! !   =:: 0 05 M MM,55 > >55OO66 [[S$uQ{"34A5A::rkrc@a]tRtRtoRtRt]R4tRtRt Vt R#) PackedMatrixia_A simplified CSR format for when non-zeros in each row are consecutive. Assuming that each row of an `(m, nc)` matrix 1) only has `nz` non-zeros, and 2) these non-zeros are consecutive, we only store an `(m, nz)` matrix of non-zeros and a 1D array of row offsets. This way, a row `i` of the original matrix A is ``A[i, offset[i]: offset[i] + nz]``. cWnW nW0nVP^8XgQhVP^8XgQhVP^,VP^,8XgQhR#)r=N)rrncr'r)rQrrrjs&&&&rU__init__PackedMatrix.__init__sQ vv{{{{awwqzV\\!_,,,rkcTVPP^,VP3#r8)rrrjr}s&rUrPackedMatrix.shapesvv||A''rkc\P!VP4pVPP^,p\ VP^,4Fwp\ VP VPV,, V4pVPVRV13,WVPV,VPV,V,13&Ky V#)r.N)rIrerrrminrjr)rQoutnelemrnels& rUtodensePackedMatrix.todenseshhtzz" Qsyy|$Adgg A.6C:>&&DSD/C4;;q>$++a.3"666 7% rk)rrjrN) r/r0r1r2__doc__rkpropertyrrtr3rrs@rUrhrhs--((rkrhc ^RIHpVPpVPpVPpVP wrxVP ^,V8XgQhVP 4p VP 4p \W'4Fp W[,p \W4Fp V \W48dK,V!W^3,W^3,4wrpVW^3&\^V4F2p\WWV3,WV3,4wWV3&WV^, 3&K4 RWR3&\V P ^,4F+p\WWV3,WV3,4wWV3&WV3&K- K K \\V P ^,44p\V \P!V\PR7V4pVV 3#)zbThis is a python counterpart of the `_qr_reduce` routine, declared in interpolate/src/__fitpack.h )dlartgrAr9r;)scipy.linalg.lapackryrrrjrr*rrpfprotarrhrIrcint64)a_prstartrowryrrrjrnzRrroirr3r rroffsR_ps&&& rU _qr_reduce_pyrs+ A ZZF B GGEA 771:?? A B 8  YrACJQ!tWa1g.GA!AdG1b\%+A!qD'1T7%C"Q$ac6"AeH288A;'%+A"T(B!tH%E"a4"T(( * aggaj! "D q"((4rxx8" =C 7Nrkc`W,W,,pV)V,W,,pWE3#)z@Givens rotate [a, b]. [aa] = [ c s] @ [a] [bb] [-s c] [b] r.)r3r rrTr&rs&&&& rUr{r{s) #)B AB 6MrkcNVPpVPwr4VPpVP^,VP^,8XgQh\P!VRV4pW^, ,W%^, ^3,, We^, R3&\ V^, RR4Frp\ WEV, 4pW'^V1R3,Wg^,Wx,1R3,,P^R7p W,V , W'^3,, WgR3&Kt V#)zxBacksubsitution solve upper triangular banded `R @ c = y.` `R` is in the "packed" format: `R[i, :]` is `a[i, i:i+k+1]` N.r r;)rrrjrIrrrpr) rrrrrrjr3rrssumms && rUfpbackrs A GGEA B 771: ## # afAT7Q!tQwZ'AdCiL 2a4R "dmQsUD.!Ac!%in$55:::BdTka1g-S& ! HrkcDa]tRtRtoRtRtRtRtRtRt Rt Vt R #) TestGivensQRic^p\P!V\R7pV^,^^V,, ,p\W24pW4WR3#)r?r9)rIrMrNr)rQrSr4rwrr2s&& rU_get_xytTestGivensQR._get_xyts@  IIau % qD1ac7N  Qzrkc ^ pVPV4wr#rE\P!W$V4p\P!VP 44wrxVP V,p VPwrWP^,V, ^, 8XgQhVPRRV^,1,p \P!V \PR7p VPPW^,4p \WV 4pVR,p\P !WW4\#\P$!VP 4V,VP 4V, 4\P&!V4RR7\#\P$!\)WR,, 4\)WR,,44\P&!V 4RR7\P*!W4p\P,!VP.VP0W/WE\P2!V4V4wpp\#VVR,RR7R#) rNr9rYrtgvIh%L=-a=rEN)rEr)rr rRrrrtrrindicesrIascontiguousarrayr|datarrhr* qr_reducerminimumrr0r rrrjrn)rQrSrwrr2r4a_csrqrqTyrrjrrry_c_fullc_bandedrs& rU test_vs_fullTestGivensQR.test_vs_full#s ]]1% a%%aA.uuU]]_%ccAg WWQZ!^a''''w!A#w'%%fBHH= JJ  qA# & B ' wZ1b-  199;?#$99;?457]]15EE S 3sX~#6#&sX~#689;s9KRW Y!!!addA1aRTU!QU;rkc^ pVPV4wr#rE\P!W$V4pVPwrxWP^,V, ^, 8XgQhVPRRV^,1,p \ P !V \ PR7p VPPWu^,4p \WV4p VR,p \W4wr\P!WW4\V PV PRR7\!V P"V P"RR7V P$V P$8XgQh\WRR7R#)rNr9rYrtF) check_dtyper)rr rRrrrIrr|rrrhrr*rrrrrrj)rQrSrwrr2r4rrrjrrrrRRrs& rUtest_py_vs_compiled TestGivensQR.test_py_vs_compiledFs ]]1% a%%aA. WWQZ!^a''''w!A#w'%%fBHH= JJ  qA# & B ' wZq%1b.acc. 188?uu}}U+rkct^ pVPV4wr#rE\P!^V^,\R7p\P !W$WV4wrxp VP ^,p \P!W$V4p WR,,P4p V PPW^,34p V PRRV^,1,P\P4p\W}RR7\!W4WP ^,V, ^, 8XgQhR#)rr9NrYrtr)rrIrMrNr* data_matrixrr rRtocsrrrrrIr|rr)rQrSrwrr2r4r+rrrjrra_wA_offset_s& rUtest_data_matrixTestGivensQR.test_data_matrixas ]]1% a IIa1E * ,,Q18 2 GGAJ%%aA.z!((* XX  qA#h '++g1g&--bhh7E*(WWQZ!^a''''rkc ^ pVPV4wr#rE\PW3^,3,p\P!W$V\P !V44wrgp\ WgV4p \P!WgW4\W4p \P!WhW#WE\P !V4V4wp p \WRR7R#)rrsrtN) rrIr0r*rrnrhrrr) rQrSrwrr2r4rrrjrr3rrs & rU test_fpbackTestGivensQR.test_fpbackrs ]]1% a EE!T'N ,,Q1bll1oF 2 B '1b, 1L??1!bll1oqIAqE*rkc 0^ pVPV4wr#rE\P!V^,3\R7pVFUp\ \ P !WEWqV^,4V4\ \ P !WEWq^V,4V4KW R#)rr9N)rrIrerNrr*r)rQrSrwrr2r4 zero_arrayrs& rUtest_evaluate_all_bspl#TestGivensQR.test_evaluate_all_bspls{ ]]1% aXXq1uhe4 D **1!a%@* N **1!A#>  Lrkr.N) r/r0r1r2rrrrrrr3rrs@rUrrs+! ? ? !! '@ ?+' & & & & & & & &' & & ' &@ ? ?sH=H2/I"II.6 J J2 I  I  I+ . I?  J  J& c \P!\R44;_uu_4pVR,pVR,pVR,pRRR4\XX4!V4p\ VXRRRR7R# +'giL4;i) a} Data is generated in the following way: >>> np.random.seed(1234) >>> n = 100 >>> x = np.sort(np.random.random_sample(n) * 4 - 2) >>> y = np.sin(x) + np.random.normal(scale=.5, size=n) >>> np.savetxt('x.csv', x) >>> np.savetxt('y.csv', y) We obtain the result of performing the GCV smoothing splines package (by Woltring, gcvspl) on the sample data points using its version for Octave (https://github.com/srkuberski/gcvspl). In order to use this implementation, one should clone the repository and open the folder in Octave. In Octave, we load up ``x`` and ``y`` (generated from Python code above): >>> x = csvread('x.csv'); >>> y = csvread('y.csv'); Then, in order to access the implementation, we compile gcvspl files in Octave: >>> mex gcvsplmex.c gcvspl.c >>> mex spldermex.c gcvspl.c The first function computes the vector of unknowns from the dataset (x, y) while the second one evaluates the spline in certain points with known vector of coefficients. >>> c = gcvsplmex( x, y, 2 ); >>> y0 = spldermex( x, c, 2, x, 0 ); If we want to compare the results of the gcvspl code, we can save ``y0`` in csv file: >>> csvwrite('y0.csv', y0); z gcvspl.npzrwry_GCVSPLN-C6?F)r[r\r)rIloadrrr)rQrrwrry_comprs& rUtest_compare_with_GCVSPL,TestSmoothingSpline.test_compare_with_GCVSPLsjRWWY|, - -S AS AJ'H . (1-a0 4US. -s A00 B c L\PPR4p^dp\P!VP V4^,^, 4pV^,\P !^V,4,V^,,VP RRV4,pVPV4VPV4rT\WERR7p\WE^RR7pVPV^,VR ,^V,4p\V!V4V!V4RR7R #) z~ In case the regularization parameter is 0, the resulting spline is an interpolation spline with natural boundary conditions. rrAr)lamr[rxrYrtNr;) rIrPrrrOrrrOrr rer) rQrRrrSrwr spline_GCV spline_interpgrids && rUtest_non_regularized_case-TestSmoothingSpline.test_non_regularized_cases ii##D)  GGC%%a(1,q0 1 qD266!a%= 1a4 '#**Ra*@ @zz!}bjjm1*1R8 *1IF {{1Q42A. 4(%d+" $rkc8\PPR4p^dp\P!VP V4^,^, 4pV^,\P !^V,4,V^,,VP RRV4,p\VPWE34wrE\WE4pVP\^d4^ R7FpVPV4p\P!V\V44P!R4\WEV4p \#V!WG,4WW,, 4p \#V !WG,4WW,, 4p W8gK\%RV R RV R 24h R#) rrArrg>@zJSpline with weights should be closer to the points than the original one: z.4z < N)rIrPrrrOrrrvrOrchoicerrrrrrr0rG) rQrRrrSrwrrFindr+spl_worigweighteds && rUtest_weighted_smoothing_spline2TestSmoothingSpline.test_weighted_smoothing_spline sCii##D)  GGC%%a(1,q0 1 qD266!a%= 1a4 '#**Ra*@ @2::v&#A)::eCjr:2C A FF1c#h  # #C ()!2Es16{QV+,D5=1612H $CCG)3$,R="2333rkr.N) r/r0r1r2rrrr`r fail_slowrr3rrs@rUrrs9 (>6Tp$* [[133rkrc <aaaaaa a Vwo o Vwoo\S4S, ^, pVS^,8gQh\S4S, ^, oSS^,8gQh\VVVVVV V 3Rl\V444p\P!V4#)z-A naive 2D tensort product spline evaluation.c 3<"TFGp\S4F5pSW3,\SSVS4,\S SVS4,xK7 KI R#5ir)rr4) rixiyr3r4nytxtyrwrs & rUrbspline2..' sO7!"E"Ibi!Aq"b/)AaBO;;,5<!sAA)rQrrrIrO) rr2r3r4nxr(rrrrwrs &&ff @@@@@rUbspline2r s DAq FB R1qB !A#II R1qB !A#II 77b 7 7C ::c?rkctV^8Xd(W2,Tu;8:dW2^,,8dR#R#W2V,,W2,8XdRpMHWV,, W2V,,W2,, , \W^, W#4,pW2V,^,,W2^,,8Xd RpWE,#W2V,^,,V, W2V,^,,W2^,,, , \W^, V^,V4,pWE,#)rr0rAr4rs&&&& rUr4r4, sAvda(!aC&(s1(c11v~ A$hQ3!$ '!AsA/ 91Qx1qS6  7N!AhlQs1uXA#. /!AsAaC2C C 7Nrkcaaaa\S4S, ^, pVS^,8d\S4V8gQh\VVVV3Rl\V444#)r.c3\<"TF!pSV,\SSVS4,xK# R#5irrrs& rUrbspline..= s'6Xqta1am##Xr)rQrr)rwr2r3r4rSsffff rUrfrf: sD A QA 1H3q6Q;' ' 6U1X6 66rkc0a]tRtRtoRRltRtRtVtR#) NdBSpline0i@ c r\V4pV\VP48:gQh\V4\;QJd.RV4FNK 5M !RV44Vn\;QJd.RV4FNK 5M !RV44VnW nR# \dT3T,pLi;i)aWTensor product spline object. c[i1, i2, ..., id] * B(x1, i1) * B(x2, i2) * ... * B(xd, id) Parameters ---------- c : ndarray, shape (n1, n2, ..., nd, ...) b-spline coefficients t : tuple of 1D ndarrays knot vectors in directions 1, 2, ... d ``len(t[i]) == n[i] + k + 1`` k : int or length-d tuple of integers spline degrees. c3N"TFp\P!V4xK R#5ir)operatorindex)rkis& rUr&NdBSpline0.__init__..Y s6Abx~~b))As#%c3Z"TF!p\P!V\R7xK# R#5i)r9N)rIrOrN)rtis& rUrrZ s?Qrrzz"E22Qs)+N)rQrrFr<r4r2r3)rQr2r3r4r's&&&& rUrkNdBSpline0.__init__A s1vs177|###  F 6A66A66?Q??Q??  T A s BB65B6cF\VP4p\V4V8XgQhR.V,p\V4FpVPV,W,reVPV,pWeV,8XdWsV&M\P !WV4^, W4&WSV,,Tu;8:dWSV,^,,8:gQhQhW4,V8d W4,\V4V, 8dKQh \ V4p^p\V4Uu.uF:p\W4,VPV,, W4,^,4NK< p p\P!V !Fp VPV ,\P!\V4Uu.uF?p\W,VPV,W,VPV,4NKA up4,p W, pK \P!V4#uupiuupi)none) rQr2rr4rIrr< itertoolsproductr3prodr4rO) rQrwr'rr tdxdr4ritersidxterms && rU__call__NdBSpline0.__call__] s466{1v~~JtOtAVVAYq AU{!r.2d8r/R!QZ/ // //419B! !3 33 !H a+,AD$&&)SVTVVAY*O3>*@"AAD NF-zz&!! I*@sAH%AH )r3r4r2Nr?)r/r0r1r2rkrr3rrs@rUrr@ s8""rkrca]tRtRtoRt]3Rlt]3RltRtRt Rt Rt R t R t R tR t]3R ltRtRtRtRt]P,P/RRR.4]P,P/RRR.4R44t]P,P/RR R!.4R4tRtRtRtRtRtRtRt Vt!R#)" TestNdBSplinei} c`\PPR4p^ ^rC^p\P!VP W4,^,R74pVP W53R7pVP V4pVP V4p\ WgV4p\V3Wt4p VP ^R7p VP V 4p \V !V R,4V!V 4RR7V !V R,4PV P^,VP^,38XgQhR#)rrNrsrtr) rIrPrrrrOr rrr) rQrRrrSr4n_trr2r3rTnbxis && rUtest_1DTestNdBSpline.test_1D sii##E*11 GGCKKQUQYK/ 0 KKaYK ' JJqM JJqM A!  tQ " [[b[ ! ZZ^ 2g;"E +"W+$$!aggaj(AAAArkc&VP^4pV^,p\W#^R7pV^,^V,,p\W%^R7pVPVP3pVPR,VPR,,pWx^3#)rrjrNrErMr r2r3 rQrRrwrrFy_1rt2rs && rU make_2d_caseTestNdBSpline.make_2d_case sx IIaL qD +dQqSj"1Q/eeUWW  UU7^eggg. .qyrkcVP^4pV^,p\W#^R7pVP^VPR7R,pV^,^V,,p\W%^R7pVPVP3pVPR,VPR,,pWxVP VP 3#)rrjr9rrr )rMr rNr2r3r4r s && rU make_2d_mixedTestNdBSpline.make_2d_mixed s IIaL qD + IIarzzI *S 0dQqSj"1Q/eeUWW  UU7^eggg. .suuegg%%rkc `.R OpVPVR7wr4pVUUu.uF*wrgV^,V^,^V,,,NK, pppVU U u.uFIp \T VU u.uFp \P!V 4NK up \P!V4V4NKK p p p \P!V \PR7p \ VPV 4VPWPR7RRR7\ W4^R7p V !V4P\V438XgQh\ V !V4VPWPR7RR7\P!V4p\PVR &VPV4p\VPV !V44VP.R O44\PPR4p V PRR 7^,pVPV4pV !V4pVPR8XgQhVPVR4P pVR,VR,rv\ \#WR7V^,V^,^V,,,RR7R #uuppiuup iuup p i)rrrr9Frsr^r[rjrtrrNrr8)r8r.)rrr)TFF)r@r?r=)r@r?)r;r=)r.)r..)rrrIrOrNrrrrQrKrrrPrrrrr)rQrRrr rr4rwrtargetrrr2bspl2rrrrrs&& rUtest_2D_separableTestNdBSpline.test_2D_separable s, /%%%, 467B&1!Q$!Q$1*%%B7  B (*+1A+ 2    ZZ"** - 2 6 <$)" $ "A&Ry3r7+---b  6 <5 JZZ^664 ZZ^r+RZZ8L-MNii##E* [[i[ (1 , ZZ^r||v%%%jjW%''6{CK1/11qs +% 9M8 , s0JJ*&J%!J*%J*c^p.R OpVUUu.uF*wrEV^,V^,^V,,,NK, pppVPVR7wrxp VPWW3^R7p R p \Wz^R7p V !V 4p \WxV 4!V 4pV PR 8XgQh\ V VPV.^,4RR7V !V4P\ P!V4RRV P PVR,8XgQh\ V !V4VPWaPR7R,RR R 7VPV R4p\W^R7pV!V 4p V PR8XgQhVPRVPR7V,p\ V VRR7V!V4P\ P!V4RRVP PVR,8XgQh\ V!V4VPWaPR7R,RR R 7R#uuppi)r=rrr rjrsrtNr9FgvIh%,=rrrr@r;rrrr=r=)r=r=rENN) rrrrrrIr3rOrNrr)rQrRr'rrwrrr rr4c2_4rbspl2_4r val_singlec2_22bspl2_22 target2_22s&& rUtest_2D_separable_2!TestNdBSpline.test_2D_separable_2 s /467B&1!Q$!Q$1*%%B7%%%, xx(qx1 B*rq)"- ||t###*!12 @ r{  BHHRL"$5 8N$NNNN  6 8+s 0C1C7c \PPR4p^p\P^^^^\P!VP ^R74^,^^^^3 ,p\P^^^^\P!VP ^R74^,^^^^3 ,pVP VP V, ^, VP V, ^, 3R7p\W43WRR7pRp\V!V4\WsV3WR4RR7\P.RO.RO3,p\V!V4VUu.uFp\WV3WR4NK upRR7R#uupi) rrrjrsrtN)r0r0r.rr=皙?g?@) rIrPrrrrrrrrr0) rQrr4rrr3rFrrs & rUtest_2D_randomTestNdBSpline.test_2D_random s9ii##E*  UU1aArwws{{{':;a?Aq!K L UU1aArwws{{{':;a?Aq!K L KKbggaik27719Q;7K 8!) B "Xq45 BUU;"#$B@BC""2h5C" $CsE; c:VP4wrr4.ROpVUUu.uF*wrgV^,V^,^V,,,NK, ppp\WW43R7p V !V4P\V438XgQh\ V !V4VRR7R#uuppi)ffffff?rjrsrtN)r3r)r8g333333@)rr)rrrrQr) rQr rkxkyrrwrrrs & rU test_2D_mixedTestNdBSpline.test_2D_mixed0 s++- 1467B&1!Q$!Q$1*%%B7"RH-Ry3r7+---b U ,8s0Bc&VPVR7wr#rE.ROp\W#WE3R7pRRRRRR^V,^3R^V,^3R^^V,3R^^V,3R^V,^V,3R /pVP4FvwrTPVU U u.uF wrV !W4NK up p VPR 7p V!WiR 7p\ WR R 7VP V 4!V4p\ WR R 7Kx RF^p\\4;_uu_4V!VVR 7RRR4\\4;_uu_4VP V4RRR4K` R#uup p i +'giLM;i +'giK;i)rrrjcZ^V^,,V^,^V,,,#rr.rwrs&&rUr2TestNdBSpline.test_2D_derivative..@ sQTQTAaCZ!8rkcL^V^,,^V,^,,#rr.r;s&&rUrr<A sQTQqS1W!5rkcLV^,V^,^V,,,#rr.r;s&&rUrr<B sAA!!4rkc^#r8r.r;s&&rUrr<C Arkc^#r8r.r;s&&rUrr<D r@rkc^#r8r.r;s&&rUrr<E r@rkc^#r8r.r;s&&rUrr<F r@rkc^#r8r.r;s&&rUrr<G sqrkr9rrsrtNr4rr.r.)rr))r;rr;rr.) rritemsrOrNrrrErG)rQrRr rr5r6rr test_casesr expected_fnrwr expected_valsdirect via_methodbad_nus&& rUtest_2D_derivative TestNdBSpline.test_2D_derivative9 sk++r+2 1"RH- 8 5 4 rT1I~ rT1I~ "I~ "I~ rT1R4L.   *//1OBJJ/12rtqQ"r2"**'M 2%F F >))"-b1J JE B 2#Fz**bV$+z**  (+* #3+****sE&  E,E?, E< ? F c \PPR4p^^r2\P^^^^\P!VP ^R74^,^^^^3 ,p\P^^^^\P!VP ^R74^,^^^^3 ,pVP VP V, ^, VP V, ^, 3R7p\P.RO.RO3,p\WE3WbV3R7p\WE3WbV3R7p \V!V4VU u.uF q!V 4NK up RR7R#uup i)rrrjrsrtNr,r-) rIrPrrrrrr0rrr rQrr5r6rrr3rrbspl2_0rRs & rUtest_2D_mixed_random"TestNdBSpline.test_2D_mixed_random_ s(ii##E*AB UU1aArwws{{{':;a?Aq!K L UU1aArwws{{{':;a?Aq!K L KKbgglQ." q0@AK B UU;"#$2(Ab2bXqH5b /12rr2 @2s E% cZ\P!^4p\P!^4R,p\W^,^R7p\W"^,^V,,^R7pVPR,VPR,,p\ VP VP 3WSP VP 34pV^,R,V^,^V,,R,,p\W3V4p\P!W4U U u.uFwrW3NK p p p V!V 4p \P!V 4P4'dQh\W!V 4RR7\V PVP4VRR7R#uup p i)rrrjNrsrtrr )rIrMr r3rr2r4rrrranyrrr) rQrwrspl_xspl_yrrvaluesrgirrTrbxis & rUtest_tx_neq_tyTestNdBSpline.test_tx_neq_tyo s5 IIaL IIaL3 "1da0"1dQqSjA6 WWW  0 0%''577+R''5771CDQ$AqD1Q3J#88%qff5!*!2!21!8 9!8qf!8 92h88C=$$&&&&SW51 FLL16F :s&F'cVP^4pV^,p\W#^R7pV^,^V,,p\W%^R7pV^,^V,,^,p\W'^R7pVPVPVP3p VPR,VPR,,VPR,,p W^3#)rrjr)NrEN)NNrEr ) rQrRrwrrFr ry_2rcr rs && rU make_3d_caseTestNdBSpline.make_3d_case s IIaL qD +dQqSj"1Q/dQqSj1n"1Q/eeUWWegg &eeM"ggm$%ggm$%qyrkc\PPR4pVPRR7^,wr#pV^,V^,^V,,,V^,^V,,^,,pVP 4wrgp\ Wg^R7p \ W#V4U u.uFqNK p p V !V 4p V PR8XgQh\WRR7R#uup i)rrrjrsrtNr? rd) rIrPrrr`rziprr) rQrrwrrrt3c3r4bspl3rrrs & rUtest_3D_separableTestNdBSpline.test_3D_separable sii##E*++7++a/aAA!$1qs Q7%%' "A&Q1 &Aa &r||u$$$U3's2 C(c nVPVR7wr#p\W#^R7p\PP R4pVP RR7^,wrxp \ WxV 4U u.uFqNK p p VPV 4p RRRRRRRRRR/p V P4F{wrV UUUu.uFwpppV!VVV4NK ppppVPVVPR 7p\V!WR 7VR R 7\VPV4!V 4VR R 7K} R #uup iuupppi)rrrjrrc^V^,,V^,^V,,,V^,^V,,^,,#rr.rwrrs&&&rUr2TestNdBSpline.test_3D_derivative.. s2q1a4x1a4!A#:'>!Q$1*q.'Qrkc^V,V^,^V,,,V^,^V,,^,,#rr.rns&&&rUrro s.q1u1qs ';q!taczA~'Nrkc^V,^V^,,^,,V^,^V,,^,,#rqr.rns&&&rUrro s.q1u!Q$ ';q!taczA~'NrkcZ^V,^V^,,^,,^,#rqr.rns&&&rUrro sq1u!Q$ ';a'?rkcR#r r.rns&&&rUrro ssrkr9rrsrtNrc)r.rr)r=rr)r=r.r)r=r.r?)r=r.r@) r`rrIrPrrrfrOrGrNrr)rQrRrgrhr4rirrwrrrxi_nprrHrrIxi_yi_zi_rJs&& rUtest_3D_derivative TestNdBSpline.test_3D_derivative s)%%%, "A&ii##E*++7++a/aaL)LqL) ZZ  Q N N ? *   *//1OBMRSUMCc[c37UMSJJ}BJJJGM E",m% H E,,R04m% P 2*Ts , D+5D0 c \PPR4p^p\P^^^^\P!VP ^R74^,^^^^3 ,p\P^^^^\P!VP ^R74^,^^^^3 ,p\P^^^^\P!VP ^R74^,^^^^3 ,pVP VP V, ^, VP V, ^, VP V, ^, 3R7p\W4V3WbR7p\W4V3WbR7pRp \V!V 4V!V 4RR7\P.RO.RO.R O3,p \V!V 4V U u.uF q!V 4NK up RR7R#uup i) rrrjrsrtN)r0r0r.r,r-g?r3gffffff?) rIrPrrrrrrrrr0) rQrr4rrtzr3rFspl_0rrRs & rUtest_3D_randomTestNdBSpline.test_3D_random sii##E*  UU1aArwws{{{':;a?Aq!K L UU1aArwws{{{':;a?Aq!K L UU1aArwws{{{':;a?Aq!K L KKbggaik27719Q; ! DK E a-BB<0 Br7 UU;""#$ Bb!9b%)b!9F!9s G% c \PPR4p^p\P^^^^\P!VP ^R74^,^^^^3 ,p\P^^^^\P!VP ^R74^,^^^^3 ,p\P^^^^\P!VP ^R74^,^^^^3 ,pVP VP V, ^, VP V, ^, VP V, ^, 3R7VP VP V, ^, VP V, ^, VP V, ^, 3R7R,,p\W4V3WbR7p\W4V3VPVR7p\W4V3VPVR7p \P.RO.RO.R O3,p \V!V 4V!V 4RV !V 4,,RR7R#) rrr/rjrsrtNr,r-r|) rIrPrrrrrrrrr0r) rQrr4rrr}r3rFrrrs & rUtest_3D_random_complex$TestNdBSpline.test_3D_random_complex sii##E*  UU1aArwws{{{':;a?Aq!K L UU1aArwws{{{':;a?Aq!K L UU1aArwws{{{':;a?Aq!K L [[rwwqy{BGGAIaK1E[ F [[rwwqy{BGGAIaK1E[ Fr IJ a-BB<15BB<15 UU;""#$ Br Rr ]2 @rk cls_extrapNT call_extrapcVP4wr4p\W4^VR7p.RO.RO.ROrp\\PWxV 34wrxp \ WxV 4U u.uFqNK p p V^,V^,^V,,,V ^,^V ,,^,,p V!WR7p \ WRR7R#uup i r?)r4rrrsrtN)r;r)rDrr)r;g@r`rrvrIrOrfr)rQrrrgrhr4rirwrrrrrrs&&& rUtest_extrapolate_3D_separable+TestNdBSpline.test_extrapolate_3D_separable s%%' "A:>abjj1),aQ1 &Aa &AA!$1qs Q7r3U3 's CrcVP4wr#pVwrV\W#^VR7p.RO.RO.ROrp\\PWV 34wrp \ WV 4U u.uFqNK p p V^,V ^,^V ,,,V ^,^V ,,^,,p V!WR7p\ WRR7R#uup irr)rQrrgrhr4rrrirwrrrrrrs&& rUtest_extrapolate_3D_separable_2-TestNdBSpline.test_extrapolate_3D_separable_2 s%%' "( "A:>abjj1),aQ1 &Aa &AA!$1qs Q7r3U3 's C c8VP4wrp\W^R7p.RO.RO.R Orvp\\PWVV34wrVp\ WVV4Uu.uFqNK p pV^,V^,^V,,,V^,^V,,^,,p V!V RR7p \P !V ^,4'gQh\P !V R ,4'gQh\V ^R V ^R RR7R#uupi) r?rjFrrsrtN)rr.r)rDrrr;)r;rr)r`rrvrIrOrfrr) rQrgrhr4rirwrrrrrrs & rU#test_extrapolate_false_3D_separable1TestNdBSpline.test_extrapolate_false_3D_separable s%%' "A&nnabjj1),aQ1 &Aa &AA!$1qs Q7ru-xxq """"xxr ####q fQrl? 's Dc |VP4wrp\W^R7p\P!R^\P^^^\P.4p\P!R R^\P^RR.4p\P!R R^^\PRR.4p\ WVV4Uu.uFqNK p pV^,V^,^V,,,V^,^V,,^,,p \P !V4\P !V4,\P !V4,p \PW&V!V 4p \P !W,4P4'gQh\WRR7R#uupi) r?rjrrrrsrtNrrDr;) r`rrIrOrKrfrrr) rQrgrhr4rirwrrrrrrrs & rU test_x_nan_3DTestNdBSpline.test_x_nan_3D s<%%' "A& JJArvvq!Q7 8 JJCBFFAsC8 9 JJCArvvsC8 9Q1 &Aa &AA!$1qs Q7xx{RXXa[(288A;6vv rxx %))++++U3's> F9c d\PPR4p^^r2\P!VP ^^^R74p\P V^,3V,WDR,3V,3,p\P!VP ^^^R74p\P V^,3V,WUR,3V,3,pVR,P P'dQhVR,P P'dQhVP VP^,V, ^, VP^,V, ^, 3R7pVPpVP P'dQh\P.R O.R O3,p\VR,VR,3WbV3R7p\VR,VR,3WbV3R7p \V!V4VU u.uF q!V 4NK up RR7R#uup i) rrNNr=rrjrsrtNr;r,r-)rIrPrrrrflags c_contiguousrrr0rrrrQs & rUtest_non_c_contiguous#TestNdBSpline.test_non_c_contiguous$ sii##E*AB WWS[[QQR[8 9 UUBqE8B;VIbL0 1 WWS[[QQR[8 9 UUBqE8B;VIbL0 1c7==-----c7==----- KKbggqj2o1277A:?Q3FGK H CC77''''' UU;"#$2c7BsG,a8<bgr#w/2h?b /12rr2 @2sH- cVP4wrp\W^R7p\^4FpRW,PnK RVPn\W^R7pV!R4V!R48XgQhR#)r?rjFNrd)r`rrr writeable)rQrgrhr4rirbspl3_s& rU test_readonlyTestNdBSpline.test_readonly= sl%%' "A&qA$)BEKK !"2Q'Y6)#4444rkcVP4wrp\P!.RO.RO.4p\WV4P WAV4p\P!WAW3V.4pVP ^,VP ^,8XgQh\ VP4VP4RR7\\4;_uu_4\P!.ROW.^,4RRR4\\RR7;_uu_4\P!^^..W.^,4RRR4R# +'giLV;i +'giR#;i)r.gؗҜ ? ?  # #aVHb#a% 8@ ?' &@ ? ?s"E#E E  E% c \PPR4p^p\P^^^^\P!VP ^R74^,^^^^3 ,p\P^^^^\P!VP ^R74^,^^^^3 ,p\P^^^^\P!VP ^R74^,^^^^3 ,pVP VP V, ^, VP V, ^, VP V, ^, 3R7p\W4V3WbR7pRp\^ W4R#)rrrjcV\P.RO.RO.RO3,pV!V4R#)r.Nr,r-r|)rIr0)rrFrs&& rUr1TestNdBSpline.test_concurrency..worker_fnb s&{&&'(B GrkN) rIrPrrrrrrr&) rQrr4rrr}r3rFrs & rUrTestNdBSpline.test_concurrencyX s%ii##E*  UU1aArwws{{{':;a?Aq!K L UU1aArwws{{{':;a?Aq!K L UU1aArwws{{{':;a?Aq!K L KKbggaik27719Q; ! DK E a-  I3rkr.)FT)TN)"r/r0r1r2rrIrrrr&r)r0r7rNrSr\r`rjryrrr`rrrrrrrrrrr3rrs@rUrr} sB,!  " &*9X)$V>&$&,$)L@ G(!& 4Q0G&@& [[\D$<8 [[]T4L9 4:9 4 [[X |'DE4F4 @ 4$@2 5 944rkrcPa]tRtRtoRtRt]PPR.RO4R4t Rt Rt Rt R t ]PPR ]P]P .4R 4tR t]PPR .RO4R4tRtRtVtR#) TestMakeNDik c\P!^4p\P!^4R,pVR ,^,V^,^V,,R ,,p\P!W4UUu.uFwrEWE3NK ppp\ W3V^R7p\ V!V4VP 4RR7\W^,^R7p\W"^,^V,,^R7p VPR ,V PR ,,p \ WPR^R7^RI H p V !W3VR R 7p \ V !V4V!V4R R7R#uuppi)rrNrjrYrtrrZ)rlinearr!rsrr ) rIrMrrr'rravelr r3scipy.interpolater) rQrwrrYrrTrrrWrXrRGIrZs & rUtest_2D_separable_simple#TestMakeND.test_2D_separable_simplel s IIaL IIaL3 7Q!Q$1*g!66!*!2!21!8 9!8qf!8 9A66Q/R&,,.u=#1da0"1dQqSjA6 WWW  0 0FFQ7 E1&&2Bb6:sE)cL\P!^4p\P!^4p\P!W4UUu.uFwr4W43NK pppVR,^,V^,^V,,R,,p\P!WfWf34p\ W3V^\ PR7pV!V4p \P!WfWf34P\4p V PR8XgQh\V P^^^4V RR7VPR4p \ W3V ^\ PR7pV!V4p V PR 8XgQh\V P^^^^4V PR4RR7R#uuppi) rNr4solverrsrtrr )$r@r)rr=r=) rIrMrrr1r'sslspsolverIrNrrr) rQrwrrrTrrYvalues4rrrvalues22s & rUtest_2D_separable_trailing_dims*TestMakeND.test_2D_separable_trailing_dims sR IIaL IIaL!*!2!21!8 9!8qf!8 97Q!Q$1*g!66))VV<=A67a DbFF;<CCEJ||w&&&q!Q/U ,??<0A68qEb||z)))q!Q2|45 B':sF r4c\P!^4p\P!^4R,p\P!W#4UUu.uFwrEWE3NK pppV^,R,V^,^V,,R,,p\ W#3Wq\ P R7p\V!V4VP4RR7R#uuppi)rrNrrYrtrr ) rIrMrrr'rrrr) rQr4rwrrrTrrYrs && rUr7TestMakeND.test_2D_mixed s IIaL IIaL3 !*!2!21!8 9!8qf!8 9Q$AqD1Q3J#88A66s{{CR&,,.u= :s Cc~\P!^4p\P!^4R,p\PVR&VR,^,V^,^V,,R,,p\\4;_uu_4\ W3V^R7RRR4R# +'giR#;i)rrNrjr;rr )rIrMrKrErGr')rQrwrrYs& rU test_2D_nansTestMakeND.test_2D_nans sy IIaL IIaL3 "7Q!Q$1*g!66 : & & ! ,' & & &s B++ B< c \P!.RO4p\P!.RO4p\P!.RO.RO.RO.RO.RO.RO.4pWV3#)r)rrBrCrDrr)r.r=r.r=r.r.)r.r=r?r=r.r.)r.r=r=r=r.r.)rIrc)rQrwrrs& rU_get_sample_2d_dataTestMakeND._get_sample_2d_data sP HH. / HH. / HH""""""   QwrkcVP4wrp\W3V^R7p\W3VRR7p\P!.RO.RO.4P p\ V!V4V!V4RR7R#) r.rjrr!rsrtNr.gffffff@g333333@rffffff @333333?r?r.rrrDrtr0r?rr'rrIrcrrrQrwrrrrZrs& rUtest_2D_vs_RGI_linear TestMakeND.test_2D_vs_RGI_linear sk**,aA61*%qfaA XX66899:  R#b'6rkcVP4wrp\W3V^\PR7p\ W3VRR7p\ P !.RO.RO.4Pp\V!V4V!V4RR7R#) r?r cubic_legacyr!rsrtNrr rr'rrrrIrcrrrs& rUtest_2D_vs_RGI_cubicTestMakeND.test_2D_vs_RGI_cubic sq**,aA61#++>%qfaG XX66899:  R#b'6rkrcVP4wr#p\W#3V^VRR7p\W#3VRR7p\P!.R O.R O.4P p\ V!V4V!V4RRR7R#) r?ư>)r4rr\rr!rsrrZNrrr)rQrrwrrrrZrs&& rUtest_2D_vs_RGI_cubic_iterative)TestMakeND.test_2D_vs_RGI_cubic_iterative ss **,aA61&tD%qfaG XX66899:  R#b'DArkcVP4wrp\W3V^\PR7p\ W3VRR7p\ P !.RO.RO.4Pp\V!V4V!V4RR7R#) rrquintic_legacyr!rsrtNrrrrs& rUtest_2D_vs_RGI_quintic!TestMakeND.test_2D_vs_RGI_quintic sr**,aA61#++>%qfa8HI XX66899:  R#b'6rkzk, methc\PPR4p\P!VP ^R74p\P!VP ^R74p\P!VP ^R74pVP R R7p\ WEV3Wq\ PR7p\WEV3WrR7p \PP RRR R7p \V!V 4V !V 4RR7R #) i@rrr!gffffff?r/rrsrtNr:)rdr?) rIrPrcumsumrr'rrrr) rQr4methrndmrwrrrYrrZrs &&& rUtest_3D_random_vs_RGI TestMakeND.test_3D_random_vs_RGI syy$$V, IIdlll* + IIdlll* + IIdlll* +9-A!9f#++F%qQiE YY  3Sw  ?R#b'6rkc tVP4wrpR^/p\\RR7;_uu_4\W3V3R^/VBRRR4\\RR7;_uu_4\W3\P !W3343R^/VBRRR4R# +'giLY;i +'giR#;i)maxiterrr6r4N)rrErGr'rIr1)rQrwrr solver_argss& rUtest_solver_err_not_converged(TestMakeND.test_solver_err_not_converged s**,a !n :X 6 6  6Q 6+ 67:X 6 6  1& 1 FQ F+ F7 67 67 6 6sB#&B& B# & B7 r.N))r?r?rE)r?r.)r.r?rH))r.r)r?r)rr)r/r0r1r2rrr`rrr7rrrrrgmresgcrotmkrrrrr3rrs@rUrrk s7(B4 [[S"JK>L>- 77 [[X 3;;'?@ BA B7 [[N 7 7GGrkrc\a]tRtRtoRtRtRtRtRtRt Rt R t R t R t R tVtR #) TestFpcheci c^p\P!^ 4P^^4p\P!^ 4p\P!\ RR7;_uu_4\ P!W2V4RRR4\P!\ RR7;_uu_4\ P!W#V4RRR4R# +'giLV;i +'giR#;i)r.z 1D sequencer6N)rIrMrr`r rG_bfpcheckrQr4r2rws& rU test_1D_x_tTestFpchec.test_1D_x_t s  IIbM ! !!Q ' IIbM ]]:] ; ; JJqQ <]]:] ; ; JJqQ < ;< ;< ; ;s#C 'C C  C- c^p^V^,,^, pV^ ,p\P!V4p\P!V4p\P!WTV4^ 8XgQh\P !\ RR7;_uu_4\P!WTV4RRR4^V^,,^,pW!, ^, p\P!V4p\P!V4p\P!WTV4^ 8XgQh\P !\ RR7;_uu_4\P!WTV4RRR4R# +'giL;i +'giR#;i)r?z Need k+1*r6N) rIrMdfitpackfpchecr`r rGrr)rQr4rSrr2rws& rUtest_condition_1TestFpchec.test_condition_1 s  AY] F IIaL IIaLqQ'2--- ]]:[ 9 9 JJqQ : qsGaK EAI IIaL IIaLqQ'2--- ]]:[ 9 9 JJqQ : 9: 9: 9 9sE'E:' E7 : F c(^p^.V^,,^.,^.V^,,,p.ROp\P!W2V4^8XgQh\P!W2V4eQhVP 4pV^,VR&\P!W4V4^8XgQh\ P !\RR7;_uu_4\P!W4V4RRR4VP 4pVR,V^&\P!W4V4^8XgQh\ P !\RR7;_uu_4\P!W4V4RRR4R# +'giL;i +'giR#;i)r?Nz Last k knots*r6zFirst k knots*)r.r=r?r@rr;rrrrr*r`r rG)rQr4r2rwtts& rUtest_condition_2TestFpchec.test_condition_2( s'  C1IOqc1Q3i ' qQ'1,,,zz!"*** VVXA2qa(B... ]]:_ = = JJqa >VVX21qa(B... ]]:-= > > JJqa ? > > = ? > >s E- F- E= F c0^p^.V^,,^^.,^.V^,,,p.ROp\P!W2V4^8XgQh\P!W2V4eQh^.V^,,^^.,^.V^,,,p\P!W2V4^8XgQh\P !\ RR7;_uu_4\P!W2V4RRR4R# +'giR#;i)r?NzInternal knots*r6r.r=r?rr@r)rrrrr`r rGrs& rUtest_condition_3TestFpchec.test_condition_3> s  C1IA !ac * "qQ'1,,,zz!"*** C1IA !ac *qQ'2--- ]]:-> ? ? JJqQ @ ? ? ?s "D D c^p^.V^,,^.V^,,,p.ROp\P!W2V4^8XgQh\P!W2V4eQhVP 4pV^,V^&\P!WBV4^8XgQh\P!W2V4eQhVP 4pV^,^, V^&\P!WBV4^(8XgQh\ P !\RR7;_uu_4\P!WBV4RRR4VP 4pVR,^,VR&\P!WBV4^(8XgQh\ P !\RR7;_uu_4\P!WBV4RRR4R# +'giL;i +'giR#;i)r?NzOut of bounds*r6rr;r)rQr4r2rwrgs& rUtest_condition_4TestFpchec.test_condition_4K sr  C1IQqS ! "qQ'1,,,zz!"*** VVX!1ra(A---zz!"*** VVX!q1ra(B... ]]:-= > > JJra ?VVX22ra(B... ]]:-= > > JJra ? > ? > ? > >s!G&G G  G, c^p.ROp.ROp\P!W2V4^28XgQh\P!\RR7;_uu_4\ P !W2V4RRR4.ROp\P!W2V4^28XgQh\P!\RR7;_uu_4\ P !W2V4RRR4R# +'giLx;i +'giR#;i)r.Schoenberg-Whitney*r6N)rrr.r=r=)r.r.r.)rrrrrr`r rGrrrs& rUtest_condition_5_x1xm TestFpchec.test_condition_5_x1xmj s   qQ'2--- ]]:-B C C JJqQ D qQ'2--- ]]:-B C C JJqQ D C D C D C Cs C3C( C% ( C9 c^p.ROpRR.p\P!W2V4^8XgQh\P!W2V4eQhR#)r.rg333333?N)rrr.r.)rrrrrs& rUtest_condition_5_k1TestFpchec.test_condition_5_k1x sE   #JqQ'1,,,zz!"***rkc^p^.V^,,^.,^.V^,,,p^.^,p\P!W2V4^28XgQh\P!\RR7;_uu_4\ P !W2V4RRR4^.V^,,^.,^.V^,,,p^.^,p\P!W2V4^28XgQh\P!\RR7;_uu_4\ P !W2V4RRR4R# +'giL;i +'giR#;i)r?rr6Nrrs& rUtest_condition_5_1TestFpchec.test_condition_5_1 s  C1IOqc1Q3i ' CEqQ'2--- ]]:-B C C JJqQ DC1IOqc1Q3i ' CEqQ'2--- ]]:-B C C JJqQ D C D C D C Cs>D;E; E  E c^p^.V^,,^^.,^.V^,,,pR.^,^.,p\P!W2V4^28XgQh\P!\RR7;_uu_4\ P !W2V4RRR4R.^,^^.,p\P!W2V4^8XgQh\ P !W2V4eQhR# +'giL^;i)r?r.rr6Nrrs& rUtest_condition_5_2TestFpchec.test_condition_5_2 s  C1IA !ac * E!GqcMqQ'2--- ]]:-B C C JJqQ DE!Gq!f qQ'1,,,zz!"*** D Cs C55 D c^p.ROp.ROp\P!W2V4^28XgQh\P!\RR7;_uu_4\ P !W2V4RRR4R# +'giR#;i)r.rr6N) rrr=r?r@rrrr)r.r.r.@rrrrrs& rUtest_condition_5_3TestFpchec.test_condition_5_3 sZ  ' )qQ'2--- ]]:-B C C JJqQ D C C Cs A// B r.N)r/r0r1r2rrrrrrrr r rr3rrs@rUrr s<   *!,  !>  +   +  rkrc\P!WW")4p\\V4^, 4Uu.uF(qSWE,WE^,,P 4NK* ppW4^R,p\\V44FApWu,^, pWe;;,V, uu&We^,;;,V,uu&KC VR;;,VR,, uu&\ \ V4\ V4RR7Wd3#uupi)z)Split the knot interval into "runs". rYrtr;)rIrrrQrr) rwr2r4 residualsrrfpartscarriescarrys &&&& rU_splitr s aG $B6;CGAI6F G6Fb1g&**,6FF G1R!G 3w<  Q U s u !  2J)B-JCKYe< :Hs.D c\WW#4wrERpRp\\V44F?pWX^,,WX,, ^8gK&WH,V8gK5TpWH,pKA VR8Xd \R4hWV,WV^,,,^,^,p W ,p \P !W4p \P VRV WV R3,p V #)z Insert a new knot given reduals.z5Internal error, please report it to SciPy developers.Nig}Ô%I)rrrQrGrIrr) rwr2r4rrridx_max fpart_maxr idx_newknotnew_knotidx_tt_news &&&& rU _add_knotr  sa+JFGI 3v;  c7RU?Q 69y#8G I $PQQ;AI.2q8K~H OOA (E EE!FU)Xy0 1E Lrkc.a]tRtRtoRt]P PR.RO4R4tRt Rt Rt Rt R t R tR t]P PR .RO4]P PR .RO4R44tRtRtRtVtR#)TestGenerateKnotsi c8\P!^\R7pV^,R^V,, ,p^p\P!R.V^,,R.V^,,,4p\ WW4R7pV!V4V, ^,p^RIHpVPWW64p\WW64p \WRR7\ WW8R7p V !V4V, ^,p VPWW;4p \WW;4p \WRR7R #) r;r9r0rAr)r4r2_fitpack_reprorYrtN) rIrMrNrcrrr%add_knotr r)rQrwrr4r2rFr_frnew_tnew_t_pyspl2 residuals2new_t2 new_t2_pys& rUtest_split_add_knot%TestGenerateKnots.test_split_add_knot s IIau % qD2q1u:   HHbT1Q3Z2$!*, -aa-VaZ!O ; Q10Q10e4qq21gkA% a6a6 6rkr4c HVP^VPR7pVPW2P,^, 4p\ \ W4V^R74R,p\ W4V^R7^,pVPWbPR7p\WVRR7R#)r;r9r4r rYrtNr;) rMrNrrrr"rrOr)rQr4rRrwrr2rs&&& rUtest_s0TestGenerateKnots.test_s0 sy IIarzzI * FF1UU719  Q/ 0 4 AA #A & ZZ**Z -E*rkc  ^ pVPW!PR7pV^,p\\W4^^R74pVP \ \ P!V4^44p\VR,VRR7R#)rr9r1rYrtNr;)rMrNrr"rOrrIrrQrRrSrwrknotsrhs&& rU test_s0_1TestGenerateKnots.test_s0_1 sd  IIazzI * qD^AA34::k"**Q-;<b 8%8rkc ^pVPV4pV^,p\\W4^^R74pVP\ \ P!V4^44p\ VR,VRR7R#)r1rYrtNr;)rMrr"rOrrIrr5s&& rU test_s0_n20TestGenerateKnots.test_s0_n20 s]  IIaL qD^AA34::k"**Q-;<b 8%8rkc \P!^ 4pV^,p\\4;_uu_4\ \ W^^^ R74RRR4R# +'giR#;i)rr4r nestN)rIrMrErGrr"rQrwrs& rU test_s0_nestTestGenerateKnots.test_s0_nest sE IIbM qD : & & QR8 9' & & &s A A+ c  0VP^VPR7pVPV\P,^, 4p^p\ \ W#VRR74p.RO.RO.R O.R O.R O.pVUu.uFqqPWqPR7NK pp\V4\V48XgQh\WV4Fwr\WRR7K \W#VRR7wpp \VR ,VPV4RR7R#uupi) aY To generate the `wanted` list below apply the following diff and rerun the test. The stdout will contain successive iterations of the `t` array. $ git diff scipy/interpolate/fitpack/fpcurf.f diff --git a/scipy/interpolate/fitpack/fpcurf.f b/scipy/interpolate/fitpack/fpcurf.f index 1afb1900f1..d817e51ad8 100644 --- a/scipy/interpolate/fitpack/fpcurf.f +++ b/scipy/interpolate/fitpack/fpcurf.f @@ -216,6 +216,9 @@ c t(j+k) <= x(i) <= t(j+k+1) and store it in fpint(j),j=1,2,...nrint. do 190 l=1,nplus c add a new knot. call fpknot(x,m,t,n,fpint,nrdata,nrint,nest,1) + print*, l, nest, ': ', t + print*, "n, nmax = ", n, nmax + c if n=nmax we locate the knots as for interpolation. if(n.eq.nmax) go to 10 c test whether we cannot further increase the number of knots. r9rr1rYrtN)rArArArArrrr) rArArArArDrrrr rArArArArBrDrrrr) rArArArArBrDrrrrr) rArArArArBrCrDrtrrrrr;) rMrNrrIrrr"rOrQrfrr) rQrRrwrr4r6wantedwantr2rrs && rU test_s_switchTestGenerateKnots.test_s_switch s. IIarzzI * FF1RUU719  ^AA6726:>A   BHH**T*4H5zS[((('EA A .(d+1ab 2::a=u=Is9$Dc`\\^44p\WR^R7p\V4R#)r;r)r r4N)rrr"next)rQrwgens& rUr!TestGenerateKnots.test_list_input> s$ qNQSA. S rkc VP^VPR7pVPW!P,^, 4pRp\ \ W#^V^ R74p\ VR,VP.R OVPR7RR7\\4;_uu_4\ \ W#^^R74RRR4R# +'giR#;i) r;r9rr>rYrt)r4r?Nr;rD) rMrNrrrr"rrOrErG)rQrRrwrr r6s&& rU test_nestTestGenerateKnots.test_nestD s IIarzzI * FF1UU719  ^AA<= "I JJ?rzzJ R : & & 2 3' & & &s )C  C c \P!^4p\P!V\P,^, 4p\ \ 4;_uu_4\ \W\P!^ 4R74RRR4\ \ 4;_uu_4\ \W\P!^4)R74RRR4R# +'giL];i +'giR#;i)r;r+N) rIrMrrrErGrr"rr@s& rUr-TestGenerateKnots.test_weightsU s IIaL FF1RUU719  : & &  " 6 7': & &  4 5' &' &' & &s +C *,C3 C0 3 D nptsr c \PPR4p^ \P!VP VR74,p\P !V\P ,^ , 4\P!V^, ^,)4,p^p\WEWaR7^,p\\WEWaR74R,p\WRR7R#)rrr1rYrtNr;) rIrPrrrrrexprrr"r) rQr rSrrwrr4r2rs &&& rUtest_vs_splrep TestGenerateKnots.test_vs_splrep_ syy$$U+ rwwt|||./ / FF1RUU72: !A# !2 2  11 "1 % .0 1" 5E*rkc f^p\P!V4pV^,p\\W#^RR74p\P !\ 4;_uu_4p\W#^RR7pRRR4\X4^8XgQh\VR,X^,4R# +'giL>;i)Jz5r1Nr;) rIrMrr"r`warnsRuntimeWarningrrQr)rQrSrwrr6rr_s& rUtest_s_too_small"TestGenerateKnots.test_s_too_smallr s  IIaL qD^AA78 \\. ) )Qe,C*1v{{b 3q6** )s B  B0 c\.RO.RO.RO^R7p\P!\RR7;_uu_4\ V4RRR4R# +'giR#;i)rA)r+r zweights are zeror6N)rAr0rBrC)rDrtrr)rArArArA)r"r`r rGr)rQrKs& rUtest_zero_weights#TestGenerateKnots.test_zero_weights s;]M]aP ]]:-? @ @ IA @ @ @s A A! r.Nr)2r)r{Gz?r)r/r0r1r2r.r`rrr2r7r;rArGrrNr-rVr]r`r3rrs@rUr"r" s70 [[S/2+3+99:*>X 4"6 [[V]3 [[S.1+24+" +rkr"c*VP^,pWV, ^, ,W,, pV^V,, ^, pW^,W!, ^, p\P!V^4pVR;;,R, uu&VR;;,R,uu&\V\P!W!, ^, 4V4!WaR7p\P !V^, VP^,3\ R7p\^VP^,^4F.p WyR3,Wy^,R3,, W^,R3&K0 WV, V,\P!V4, ,pV#)zStraitforward way to compute the discontinuity matrix. For testing ONLY. This routine returns a dense matrix, while `_fitpack_repro.disc` returns a packed one. rr:r.Nr=rr9rE) rrIrbr rPemptyrNrr factorial) r2r4rSdeltanrintrtiirmatrrs && rU disc_naiverl s  A !eaiL14 E !GaKE Q3qs1uB ))B CHHII266!%!)$a(3A 88U1Waggaj) 7D 1aggaj! $Q$!aCF)+T1W % 5[1 t~~a0 00D Krkc4a]tRtRtoRtRRltRtRtVtR#)F_densei zdThe r.h.s. of ``f(p) = s``, an analog of _fitpack_repro.F Uses full matrices, so is for tests only. NcWnW nW0nW@nVf\P !V\ R7MTVnWpnVPP^8XgQh\V\P!VP^,V, ^, 4WGR7!V4pWPR,,Vn ^RIHp \!V P#W44!P%4VnVP^8XgQhW P,p \P(V \P*!VP&P^,43,VnWPnR#)Nr9rr$r)rwrr2r4rIrnrNr+rr'r rPra_denserr%rhdiscrtb_denserrerr ) rQrwrr2r4r r+rrpr'rs &&&&&&&& rUrkF_dense.__init__ s12au-&vv{{a!RVVAGGAJNQ$67TUVW0 ;#SXXa^4<<> vv{{ Z%%BHHT\\%7%7%:;;<rkcV\P!VPVPV, 34p^RIHpHpV!VRR7wrVVPVP,pVP^,pV!VRV1RV13,VRV4p \VPWPVPR7p \P!VP^,V !VP 4VP", ^,,4p WnWP&, #)r)rr economic)rNr)rIvstackrprr scipy.linalgrr rrrr r2r4rrr+rwrrFr ) rQr!abrr rrqyrjr3rFfps && rUrF_dense.__call__ s YY dllQ&67 8 +":& SS477] WWQZ !CRC"H+r#2w 'dffaT5E5EF VVDFFAITVVtvv!5 99 :FF{rk) rprrrr4r rFr2r+rwrr)NT) r/r0r1r2rvrkrr3rrs@rUrnrn s.rkrnca]tRtRtoRt]3RltRtRt] PPR.RO4R4t Rt ] PPR.RO4R 4t] PPR 4t] PPR.RO4R 4tR t] PPR RR.4] PPR.RO4R44t] PPR RR.4R4t] PPR RR.4] PPR.RO4R44t] PPR.RO4] PPR.RO4R44t] PPR.RO4] PPR.RO4R44tRtVtR#)_TestMakeSplrepBasei Nc VPR8Xd@VP^^\P,^ 4pVP V4pRpW#V3#VP^^^ 4pVP VR,^, 4^,pRpW#V3#)rgq@H&?gQ @)rNrerIrr)rQrRrwrr s&& rU _get_xykt_TestMakeSplrepBase._get_xykt s~ <<: % Aqw+AAA7N Aq"%A$ #Q&AA7Nrkc ` \P!^^ ^ 4p\P!^^ ^ 4p\\4;_uu_4\ WVP R7RRR4\\4;_uu_4\ ^^RVP R7RRR4\\4;_uu_4\P !VP^^^34p\ WRVP R7RRR4\P !^ 4p\\4;_uu_4\ W^,VRVP R7RRR4\P !^ 4)p\\4;_uu_4\ W^,VRVP R7RRR4\P !VP^,^34p\\4;_uu_4\ W^,VRVP R7RRR4\\4;_uu_4\ VRRR 1,V^,RVP R7RRR4\\4;_uu_4\ W^,RRVP R7RRR4\\4;_uu_4\ W^,R VP R7RRR4\\4;_uu_4\ W^,^^RVP R7RRR4\\4;_uu_4\ W^,^^ VP R 7RRR4\\4;_uu_4\ \P!^4\P!^ 4RVP R7RRR4R# +'giELi;i +'giELD;i +'giEL;i +'giEL;i +'giELw;i +'giEL";i +'giEL;i +'giEL;i +'giEL;i +'giELh;i +'giEL<;i +'giR#;i) rrxNrr rN)r+r rNr8r4r rN)r4r?r rN)r r?rNr;) rIrerErGr#rNrrrrFrMrQrwrr+s& rUtest_input_errors%_TestMakeSplrepBase.test_input_errors s KK2r " KK2r " : & & dll 3': & & 1T\\ :': & &Aq)*A T\\ :' GGBK : & & a41T\\ B'WWR[L : & & a41T\\ B' GGQWWQZO $ : & & a41T\\ B': & & $B$Adll C'9 % % a43#t|| D&: & & a42t|| <': & & a411T\\ J': & & a412t|| D': & &  ! biilc4<< P' &a' & &' & &' & & ' & & ' & & ' & &' & && % %' & &' & &' & &' & &sO?O8=O(*!O<!P*!P$*+P84!Q 4 Q 3"Q44!R4AR O  O% ( O9 < P  P! $ P5 8 Q Q Q1 4 R  R  R- c x\\WW4VPR74R,p\WW4VPR7p\WWSW@PR7p\ VP VP RR7\ VP VP RR7VPVP8XgQhR#)r)r2r4r rNrYrtNr;)rr"rNr#rr2r3r4)rQrwrr4r r2spl_autospl_ts&&&&& rU_test_with_knots$_TestMakeSplrepBase._test_with_knotss  E Fr Jqqt||DAAaF EGG%8 EGG%8zzUWW$$$rkr4cPVP4wr#pVPW#W4R#r)rr)rQr4rwrr s&& rUtest_with_knots#_TestMakeSplrepBase.test_with_knots!s".."a aA)rkc T\WW0PR7p\WW0PR7p\\ WW0PR74R,p\ VP VP RR7\ VPVRR7\ VPVRR7R#)rMrYrtNr;)r#rNr rr"rr3r2)rQrwrr4rFspl_ir2s&&&& rU_test_default_s#_TestMakeSplrepBase._test_default_s'sp!!\\:"11llC <<@ A" EuwwU3qu-/rkcPVP4wr#pVPW#V4R#r)rr)rQr4rwrrs&& rUtest_default_s"_TestMakeSplrepBase.test_default_s0s".."a Q1%rkc RpVPR8XdE\P!^^\P,^4p\P!V4pM\P !^4pV^,p\ P!4;_uu_4\ P!R\4\W#^WPR8HR7pRRR4\ P!4;_uu_4\ P!R\4\W#^WPR7pRRR4\XPX^,4\\PVP^.VP ^,,3,V^,RR7R# +'giL;i +'giL;i)gKH9rr5r4r rNrrrt)rNrIrerrrMwarningscatch_warnings simplefilterr\rr#rr2rr3r4)rQr rwrr_rFs& rUr]$_TestMakeSplrepBase.test_s_too_small5s*  <<: % Aqw+Aq A " A1A  $ $ & &  ! ! allj.HJC ' $ $ & &  ! ! aa1llCC ' s1v&ceeaS#%%']23AU ,' &' &s7F04F- F* - F= c\^ pVPR8XdE\P!^^\P,V4p\P!V4pM\P !V4pV^,p\ W4WPR7p\ W4VRVPR7pVPP^8XgQhVPP^8XgQh\ W4^^V,, ,VRVPR7pVPP^8XgQhR#)rrrMh㈵>rN) rNrIrerrrMr#r3r')rQr4rSrwrrFrrcs&& rU test_shape_TestMakeSplrepBase.test_shapeRs  <<: % Aqw*Aq A ! A1A!!\\:AAt||DuuzzQww||q   A1ac7{a4Nww||q   rkc &^ p^^\P,r2\P!W#V^,4p\P!V4p\ \ 4;_uu_4\ WERRR7RRR4R# +'giR#;i)r:0yE>nonsenserN)rIrrerUrErGr#rQNrrTrwrs& rUtest_error_on_invalid_bc_type1_TestMakeSplrepBase.test_error_on_invalid_bc_typegs] !BEE'1 KKa!e $ FF1I : & & j 9' & & &s &A?? B rNrc \P!^^\P,^ 4p\P!V4p\P!^^\ V44p\ W4WRRVR8HR7p\ W4VRW!R7p\VPV^,4\\PVP^.VP^,,3,V^,RR7R#)rrrr+r4r rr+r r4rNrtN) rIrerrrQrr#rr2rr3r4)rQrNr4rwrr+r_rFs&&& rU%test_make_splrep_with_unequal_weights9_TestMakeSplrepBase.test_make_splrep_with_unequal_weightsps KK1RUU7B ' FF1I KK1c!f %QQt'Z2GI!!tqBs1v&ceeaS#%%']23AT +rkc Rp\P!^^\P,^ 4p\P!V4p\PW43,P wrVVP PRJgQhVP PRJgQh\P!^^\V44p\PWw3,P wrVP PRJgQh\W4V^RVR8HR7p \WTVR^VR7p V!W4\W6VR^VR7p V!W4\W4VR^VR7p V!W4\WVVR^VR7p V!W4R#)c\VPV^,4\\PVP^.VP ^,,3,V^,RR7R#)rrrtN)rr2rIrr3r4)rFr_s&&rUcheckO_TestMakeSplrepBase.test_make_splrep_with_non_c_contiguous_input..checksG CEE3q6 * BEE#%%!ceeAg"67F /rkFrrrrN) rIrerrr0rrrrQrr#) rQrNrrwrrrr+w1rr_rFs && rU,test_make_splrep_with_non_c_contiguous_input@_TestMakeSplrepBase.test_make_splrep_with_non_c_contiguous_inputsB / KK1RUU7B ' FF1Iqtxx$$---xx$$--- KK1c!f %ad xx$$---QQ!t'Z2GI"17C c!17C c!"7C c"B$!")+ crkc $\P!^^^ 4p\P!^\P,V,4pV^,VR,re^pRp\P!V.V^,,V.V^,,,4p\ W4WVVWxRVR7 p \ V \4'gQhV PP^,^V^,,8XgQhV PV8XgQh\V PRV^,\P!V.V^,,44\V PV^,)R\P!V.V^,,44R#)rrN)xbxer4r r2r?rNr;) rIrerrrcr#r%r r2rr4rrO) rQrNr4rwrrrr r2rFs &&& rU%test_make_splrep_impl_no_optimization9_TestMakeSplrepBase.test_make_splrep_impl_no_optimizations* KK1b ! FF1ruu9q= !1quB   HHbTQU^rda!en4 5!2tW>#w''''uu{{1~a!e,,,uuzzdqs RZZA%?@!fg B41q5>(BCrkrSr c \PP^{4p\P^\P!VP ^^\P ,V^, R74^\P ,3,p\P!V4\P!V4,p\WEW PR8HR7wrgp\WEW PR7p V^ 8XdV^8XdVPR8Xg>\V PVRR7\V PVRV)^, RR7R#R#)r%rrr rrrYrtN)rIrPrrrrrrrrrNr#rr2r3 rQrSr rrwrr2r3r4rFs &&& rU(test_make_splrep_matches_splrep_periodic<_TestMakeSplrepBase.test_make_splrep_matches_splrep_periodicsii##C( EE!RWWS[[AIAE[BCQYN O FF1Iq !)CEa!!\\:RAFt||z'A CEE15 1 CEE1Wqb1f:E :(Brkc 8\PP^{4p\P^\P!VP ^^\P ,V^, R74^\P ,3,p\P!V4\P!V4,p\WEW PR8HR7wrgp\WEW PVR7p \V PVRV)^, RR7R#)r%rrr)r rNr2NrYrt)rIrPrrrrrrrrrNr#rr3rs &&& rU"test_make_splrep_with_splrep_knots6_TestMakeSplrepBase.test_make_splrep_with_splrep_knotssii##C( EE!RWWS[[AIAE[BCQYN O FF1Iq !)CEa!!\\Q?q1"q&z6rkr.r.r=r?r@rrrdr)r3rGrd)rr;rr.rd)r/r0r1r2rNrIrrrr`rrrrr thread_unsafer]rrrrrrrr3rrs@rUr}r} sG 5Qn% [[S"45*6* 0 [[S"45&6& [[,,8 [[S),!-!(: [[YT(:; [[S/2 +3< + [[YT(:;+<+\ [[YT(:; [[S/2D3<D, [[S"34 [[S"56 ;75 ; [[S"34 [[S"567757rkr}ca]tRtRto]P P R.RO4R4t]P P R.RO4R4tRt Rt Rt Rt R t R tR tR tR tRtVtR#)TestMakeSplrepir4cB^RIHpVP4wr4p\P!^.V^,,RR.,^.V^,,,4pV!W4R,WaV4p\ W4WaV4pRFp \ V!V 4V!V 4RR7K R#) rFr8rDNrYrtrr.rr) scipy.interpolate._fitpack_reprorrrIrcrnr) rQr4rrwrr r2ff_dr!s && rUtest_fitpack_FTestMakeSplrep.test_fitpack_Fs 7.."a HHaS!A#Y#s+qc1Q3i7 8 a7Q1 %aA!$A AaD#a&u 5rkc ^RIHpVP4wr4p\P!^.V^,,RR.,^.V^,,,4p\P !VP ^,\R7pV!W4R ,WaWWR7p\W4WaWWR7p \W4WaV4p R FJp \V!V 4V !V 4RR7\P!V !V 4V !V 4RR7'gKJQh R#) rrr8rDr9NrQrYrtrr) rrrrIrcrMrrNrnrr) rQr4rrwrr r2r+fwfw_drr!s && rUtest_fitpack_F_with_weights*TestMakeSplrep.test_fitpack_F_with_weightss 7.."a HHaS!A#Y#s+qc1Q3i7 8 IIaggaj . qG*aA +qQ1*aA!$A BqE47 7{{3q647?? ??rkc ^RIHuHp\PP R4p\P ^^^^\P!VP^R74^,^^^^3 ,p\V4^rT\VPW54!P4p\W54pVP^,V^V,, ^, 8XgQh\WgRR7R#)rNrrrYrt)r interpolater%rIrPrrrrrQrhrqrtrlrr)rQr'rr2rSr4DD_denses& rUtest_disc_matrixTestMakeSplrep.test_disc_matrixs 76ii##E* EE!Q1bggckkqk&9:1*CEE>rkc b^pVP4wr#p\\W#WR74R,p\W#WR7p\W#WQVR7p\ VP VP RR7\ VP VP RR7VPVP8XgQhR#)r?r1)r2r4r rYrtNr;)rrr"r#rr2r3r4)rQr4rwrr r2rrs& rUrTestMakeSplrep.test_with_knotss .."a / 0 4qq.AAa0 EGG%8 EGG%8zzUWW$$$rkc^ pVPW!PR7pV^,p^p\W4V^R7pVPP^,^V^,,8XgQhR#rr9r1N)rMrNr#r2r)rQrRrSrwrr4rFs&& rUtest_no_internal_knots%TestMakeSplrep.test_no_internal_knots*sY  IIazzI * qD !!q)uu{{1~AaC(((rkc^ pVPW!PR7pV^,p\W4^R7p\W4^R7p\ VP VP RR7R#)rr9rjrYrtN)rMrNr#r rr3)rQrRrSrwrrFrs&& rUrTestMakeSplrep.test_default_s3sP  IIazzI * qD!!$"11-uwwU3rkc ^p\P!V4pV^,p\P!\4;_uu_4p\ W#^RR7p\ W#^RR7p\VPV^,4\\PVP^.VP^,,3,V^,RR7RRR4\X4^8XgQhR# +'giL$;i)rYrZr1rrtN)rIrMr`r[r\rr#rr2rrr3r4rQ)rQrSrwrrr_rFs& rUr]TestMakeSplrep.test_s_too_small<s  IIaL qD \\. ) )Qe,Caa51C CEE3q6 * BEE#%%!ceeAg"67F 0 * 1v{{ * )s B C)) C9 c \P!.RO\PR7p\P!.RO\PR7p\P!R.VP^,,\PR7p\ \ 4;_uu_4\ WV^^ R7RRR4R# +'giR#;i)4@r9g5{2?)r+r4r N)rgR9c@g ףp=e@gףp= Oi@gQm@gQo@g)\(p@gfffffq@gq= ףq@g(\Yr@g)\r@g(\r@gQs@gQ3s@g33333Gs@皙as@rghs@)gJ@gH@gLH@gffffffG@gE@g)\(D@gC@gB@g333333@@g<@g3333339@g5@g2@g,@g%@g@g @rA)rIrOrNrrErGr#rs& rUtest_issue_22704TestMakeSplrep.test_issue_22704Js JJ%-/JJ 8 JJ&(ZZ 1 JJy1771:-RZZ @ : & & Q" -' & & &s 'C C c^ ^r2VPW!PR7pV^,p\WEVR7p\WEVRR7pVPP^8XgQhVPP^8XgQh\WE^^V,, ,VRR7pVPP^8XgQhR#)rr9rjrr1N)rMrNr#r3r') rQrRrSr4rwrrFrrcs && rUrTestMakeSplrep.test_shapeXs11 IIazzI * qD!!$AA.uuzzQww||q   A1ac7{a48ww||q   rkc^ ^r2VPW!PR7pV^,p\WE^^R7p\WE^^R7pVPP^8XgQhVPP^8XgQhVP P ^,W#,^,8XgQhVP P ^,^V^,,8XgQhR#r)rMrNr#r3r'r2r)rQrRrSr4rwrr~rs&& rUtest_s0_vs_notTestMakeSplrep.test_s0_vs_noths11 IIazzI * qDAA+AA+ww||q   ww||q   ww}}Q1519,,,ww}}Q1A;...rkr.Nr)r/r0r1r2r`rrrrrrrrrr]rrrr3rrs@rUrrs [[S"45 66 6 [[S"45@6@ 0 ? %)4  .! / /rkrca]tRtRtoRt]P PR.RO4R4t]P PR.RO4R4t Rt Rt ]P PR^R .4R 4t ]P PR .RO4]P PRR .4R 44t RtVtR#)TestMakeSplrepPeriodicixrr4cVP^^ ^ 4pVPR4p\W4V^VPR7pVPP ^,^V^,,8XgQhR#)rrNre)rerr#rNr2r)rQr4rRrwrrFs&&& rUr-TestMakeSplrepPeriodic.test_no_internal_knots}sY KK2r " GGEN!!q$,,?uu{{1~AaC(((rkc^ p\P!^^\P,V4p\P!V4\P!V4,p\ W4V^VP R7p\ W4V^VP R7pVPP^8XgQhVPP^8XgQhVPP^,V^V,,8XgQhR#)rrN) rIrerrrr#rNr3r'r2r)rQr4rSrwrr~rs&& rUr%TestMakeSplrepPeriodic.test_s0_vs_nots  KK1RUU7A & FF1Iq !AADLLAAADLLAww||q   ww||q   ww}}Q1q1u9,,,rkcv^ p^^VP,rCVPW4V^,VPR7pVPV4p\ WVRVP R7p\ \WW4VRRR7VPV4VPV4,p\ WVRVP R7p\ \WW4VRRR7^VPV4,VPV4^,,p\ WVRVP R7p\ \WW4VRRR7R #) rr9rrrrrZrrN) rrerNrr#rNrrr)rQrRrrrTrwrrFs&& rU test_periodic_with_periodic_data7TestMakeSplrepPeriodic.test_periodic_with_periodic_datas !BEE'1 KKa!e2::K 6 FF1I!$ =a qt$? FF1Iq !!%>a qt$? bffQiK"&&)A+ %!$ =a qt$?rkc :^ p^^\P,r2\P!W#V^,4p\P!V4p\ \ 4;_uu_4\ WERVPR7RRR4R# +'giR#;i)rrrN)rIrrerUrErGr#rNrs& rU$test_periodic_with_non_periodic_data;TestMakeSplrepPeriodic.test_periodic_with_non_periodic_datasa !BEE'1 KKa!e $ FF1I : & & dll ;' & & &s &B  B r rZc \P!RR.4p\P!RR.4pV^8d\P!RR.4MRp\W#V^RVR7p\W#V^^VR7pV^8d\ VP V^,4\ \P VP^.VP^,,3,V^,4R#)rAr0rtNr)r+r4rNr )r+r4rr ) rIrcr#rrr2rr3r4)rQr rwrr+rFr_s&& rU'test_make_splrep_periodic_m_eq_2_k_eq_1>TestMakeSplrepPeriodic.test_make_splrep_periodic_m_eq_2_k_eq_1s HHc3Z  HHc3Z $%EBHHc3Z t!!q*BQQ!a0 q5 CEE3q6 *ceeaS#%%']23A rkk_fprc ^ p\P!^^V4p\P!^\P,V,4pVwrg\ WEWb^R7pRp \ WER ,V^,WbR7p \ P!4;_uu_4\ P!R\4\V R^W, 3\PWr, 33VR,4p RRR4V Pp \P!RR .4p V !V 4p\V^,V^,4R# +'giLX;i) rrrNr1r5rBMbP?rAr0r)rIrerrrr rrrr\r!rLrFrcr)rQrr rSrwrr4rzr_fp0splinerbsx_checky_checks&&& rUtest_fperiodic_basic_fit/TestMakeSplrepPeriodic.test_fperiodic_basic_fits  KK1a  FF1ruu9q= !QQ+1j#a&A;  $ $ & &  ! !(N ;&#CG rvvrv6F'GUSA' ZZ((C:&W+ GAJ/' &s AD;; E r.Nr))r.-C6)r=r )r?gAȘ)r/r0r1r2rNr`rrrrrrrrr3rrs@rUrrxsG [[S"45)6) [[S"45 -6 -@"< [[S1e*-  .   [[V%QR [[S4&)0*S0rkrca]tRtRto^ ^]3Rlt]PPR.R O4R4t ]PPR.R O4R4t Rt Rt Rt R tVtR #) TestMakeSplprepicVPWPR7VP,V, pVPV4VP V4.pWEV3#)r9)rMrNrrr)rQrr4rRrwrs&&&& rU_get_xykTestMakeSplprep._get_xyksG IIazzI *RUU 2Q 6 VVAYq "Qwrkr c P^ ^r2\P!V4\P,V, p\P!V4\P!V4.p^^R^R^ R^ /p\ WQR7wwrxr9\ WQR7wr\WRR7\V PVRR7\V4Wa,8XgQh\P!V4Pp \V PV RR7\V !V 4\WxV^R7!V 4RR7R#) rrrrr rYrtr N)rIrMrrrrr$rr2rQrOrr3r ) rQr rr4rwr num_knotsr2r3u_rFrIrs && rUtest_simple_vs_splprep&TestMakeSplprep.test_simple_vs_splpreps11 IIaL255 1 $ VVAYq "CD%u=   qa% E* qu-1v%%%ZZ]__r. Aaa03% ArkcVP4wr#p\V\4'gQh\P!V4^,^8XgQh\ W1R7wrE\ \P !V4VR7wrg\WWVR7\V^,V^,RR7\V^,4\V^,48XgQh\V^,V^,RR7V^,V^,8XgQh\P!\WT44\P!V48XgQh\W1R7wr\WWRR7\VPV^,RR7\VPPV^,RR7VPV^,8XgQhV!V4P\P!V48XgQh\\P !V4VR7wr\WWRR7\VPV^,RR7\VPPV^,RR7VPV^,8XgQhV!V4P\P!V48XgQh\\ 4;_uu_4\\P !V4PVR7RRR4R# +'giR#;i)rrrtrYN)r r%rrIrrrOrrQrr$r2r3rr4rErG) rQr rrr_rItck_au_arFs && rUtest_array_not_list#TestMakeSplprep.test_array_not_lists--/a!T""""xx{1~"""RZZ]a0 Q'Aau53q6{c%(m+++Aau51vq!!!xxa &"((1+555a%U+uQxe4q6uua   1v||rxx{***bjjmq1U+uQxe4q6uua   1v||rxx{*** : & & AA .' & & &s ,K66 L cVP^ VR7wr#p\V4wrV\V!V4VPV4RR7R#)rrrRrYrtNr r$rrrQrRrwrr4rFrIs&& rUtest_default_s_is_zero&TestMakeSplprep.test_default_s_is_zeros;--"-,aaA %8rkcVP^ VR7wr#p\V^R7wrV\VRR7wrx\WhRR7\V!V4VPV4RR7\V!V4VPV4RR7VPVP8XgQhVP P VP P 8XgQhR#)rrrrYrtrN)r r$rrr r3r) rQrRrwrr4ru_ispl_nu_ns && rUtest_s_zero_vs_near_zero(TestMakeSplprep.test_s_zero_vs_near_zero%s--"-,a!!q) !!u- u-c BHHQKe<c BHHQKd;zzUZZ'''ww}} ---rkcV\P!^\R7p\\4;_uu_4\ V4RRR4\\4;_uu_4\ V^R7RRR4\\4;_uu_4\ VRR7RRR4\ V.RR7wr#\ V.RR7wrEV!V4PR8XgQh\V!V4V.RR7R# +'giL;i +'giL;i +'giL;i) r;r9NrrrrYrtr.r;) rIrMrNrErGrr$rr)rQrwr_rrFrIs& rUrTestMakeSplprep.test_1D2s IIau % : & & AJ': & & a ': & & c "'1#&qcT*1v||v%%%A%0' &' &' &s# C2D D2 D  D  D( r.N)rrrrr/r0r1r2rIr r`rrrrrr$rr3rrs@rUr r stqR  [[S"67A8A: [[S"67/8/B9 .11rkr ca]tRtRto^ ^]3Rlt]PPR.R O4R4t ]PPR.R O4R4t Rt Rt Rt R tVtR #) TestMakeSplprepPeriodiciDcVP^^VP,WPR7pVPV4VP V4.pWEV3#)rr9)rerrNrr)rQrSr4rRrwrs&&&& rUr  TestMakeSplprepPeriodic._get_xykGsE KK1RUU7AZZK 8 VVAYq "Qwrkr c ^ p\P!^^\P,V4p\P!V4\P!V4.p^^R^R^R^/p\ WA^R7wwrgr\ WARR7wr\WRR7\V P4WQ,8XgQh\V !V 4\WgV^R 7!V 4RRR 7R #) rrrrrrrrYrtr rZN) rIrerrrrr$rrQr2r ) rQr rSrwrrr2r3r4rrFrIs && rUr.TestMakeSplprepPeriodic.test_simple_vs_splprepLs  KK1RUU7A & VVAYq "D"dBb9  A. qaj9 E*355zY\))) Aaa 8 ;" 0rkc&VP4wr#p\V\4'gQh\P!V4^,^8XgQh\ W1^R7wrE\ \P !V4V^R7wrg\WWVR7\V^,V^,RR7\V^,4\V^,48XgQh\V^,V^,4Fwr\WRR7K V^,V^,8XgQh\P!\WT44\P!V48XgQh\W1RR7wr\WWRR7V PV^,8XgQhV !V4P\P!V48XgQh\\P !V4VRR7wr\WWRR7V PV^,8XgQhV !V4P\P!V48XgQh\\4;_uu_4\\P !V4PVRR7RRR4R# +'giR#;i)rrrtrYrrN)r r%rrIrrrOrrQrfrr$r4rErGr) rQr rrr_rIrrrrrFs && rUr+TestMakeSplprepPeriodic.test_array_not_listes--/a!T""""xx{1~"""Q'RZZ]aQ7 Q'Aau53q6{c%(m+++#a&%(+FB B /,1vq!!!xxa &"((1+555aj9U+uua   1v||rxx{***bjjmq*EU+uua   1v||rxx{*** : & & AAz B' & & &s -I?? J cVP^ VR7wr#p\VRR7wrV\V!V4VPV4RR7R#)rrSrRrrxrYrtNrrs&& rUr.TestMakeSplprepPeriodic.test_default_s_is_zeros>--"-,aa4A %8rkc:VP^ VR7wr#p\V^RR7wrV\VRRR7wrx\WhRR7VPV4p \V!V4V RR7\V!V4V RRR 7VPVP8XgQhR #) rr3rrrrYrtrrrZN)r r$rrr ) rQrRrwrr4rr!r"r#y_arrs && rUr$0TestMakeSplprepPeriodic.test_s_zero_vs_near_zeros--"-,a!!q*= !!ujA u- c E6c E4@zzUZZ'''rkc\P!^^\P,^4p\P!V4p\ \ 4;_uu_4\ V^R7RRR4\ \ 4;_uu_4\V^RR7RRR4\ \ 4;_uu_4\VRRR7RRR4\V.RRR7wr#V!V4PR8XgQh\V!V4V.RR7R# +'giL;i +'giL;i +'giLw;i) rrNrrrrYrtr') rIrerrrErGrr$rr)rQrwrFrIs& rUrTestMakeSplprepPeriodic.test_1Ds KK1RUU7A & FF1I : & & A1 ': & & a 4': & & c: 6'qcUJ?1v||v%%%A%0' &' &' &s$DD!6D4 D ! D1 4 E r.N)rrrrr)rs@rUr+r+DsuqR  [[S"780900 [[S"78C9C<9 (11rkr+cBa]tRtRtoRtRtR RltRtRtRt Vt R#) BatchSplineic \P!W#R4pVPRRVnVP RVPR,4Uu.uF qd!W3/VBNK upVnW0nR#uupi)r.Nr;)rImoveaxisr _batch_shaper_splines_axis)rQrwrr rkwargsyis&&&&$, rUrkBatchSpline.__init__se KK $GGCRL;<99RQS;UV;UR00;UV  WsA2cVPUu.uF q"!V4NK pp\P!W0PVP,4pVP'd#\P !VRVP 4#T#uupir)r?rIrr>rr=r@)rQrwrrs&& rUrBatchSpline.__call__sd%)]] 3]6VAY] 3 JJq++agg5 612r{{1b$**-?a? 4sBNcVPUu.uFqDPWV4NK pp\P!WPP4#uupir)r?rrIrr>)rQrrTrrrs&&&& rUrBatchSpline.integrates@?C}} M}V  aK 0} Mzz!..// NsA c\P!V4pVPUu.uFq3PV4NK upVnV#uupir)r*deepcopyr?rrQrr(rs&& rUrBatchSpline.derivatives>mmD!mmD!@C M f--b1 M  NrL)r@r>r?r) r/r0r1r2rkrrrr r3rrs@rUr;r;s$ @ 0 rkr;c a]tRtRto]P P R]/3]/3]RR/3] RR/3] RR/3.4]P P R.R O4]P P R .RO4R 444t R t Vt R #) TestBatchizmake_spline, kwargsrr0rrr eval_shaper c\PPR4p^ p^^^V3pRp\P!.VOVN5!p \P!VPV4RV4p V\ 8XdK^p V ^,3V ^,,V R,3V ^,,,p V\ WR7,pV!W3RV/VBp \W3RVRV/VBpVP!VRV/p \PPV !V 4V!V 44V P^4VP^4r\PPV !V 4V!V 44V P^4VP^4r\PPV !V 4V!V 44\PPV P!V!VP!V!4R#)lUU{)r2r4r rrN)rrr;)rIrPrrer=rrHr;rtestingassert_allcloser rr)rQ make_splinerAr rQrrSrdomainrwrr4r2r(refs&&&&& rU test_batchTestBatch.test_batchsii##$45 Aq!  KK # # # KK 5)2t 4 / )A11q5!QrUHA$66AdQn,F!4T4V4!HTH+HH KK 1j 1 ""3q63q62%%a(#*<*>!$cnnQ&7S ""3q63q62 ""3==]==&;QRrkr.N)r.r.rrF) r/r0r1r2r`rrr rrrXr3rrs@rUrPrPs [[2 b !  $ % . Hi0 1 Hd+ ,    [[\+;< [[VZ0S1= SrkrPr8rrZ)zrrrrrrr*rrnumpyrIscipy._lib._array_apirrrrrrscipy._lib.array_api_extra_libarray_api_extrarr`r rErr r r r rrrrrrrrrrrrrrwr rscipy.sparse.linalgsparserscipy.interpolate._bsplinesrrrrrrr r!r"r#r$scipy.interpolate._fitpack_implr _fitpack_implrscipy._lib._utilr%scipy._lib._testutilsr&scipy.interpolate._ndbspliner'r(rr)rr*rrr,rrrrrrrurr^rr-rxrrrrfrrhrr{rrrrrr4rfrrrrrr r"rlrnr}rrr r+r;rPr.rkrUrhs1    )(*  !!EEBGG//&953-&;;//7_ )_ )_ )D7S$S$S$l6,$ G> $b+$_)_)D%&`?`?'`?F 88!++11(Y)X  (nLnLb* ()J3J3*J3^  7 :":"z9j4j4j4ZWGWGta a L(.>"ll#l^2,,^Q7Q7h;R/(R/ R/j;\00\0 \0~< h1h1!h1V< b1b1!b1J8$S$Srk