\ci jddlZddlmZddlmZddlmZddlm Z m Z m ZdZ dZd Zd Z d Zd ZdS) N) LinAlgError)dgesv)_monomial_powers_impl) _build_system_build_evaluation_coefficients_polynomial_matrixc ,t|||||||SN)&_pythran_build_evaluation_coefficients)xykernelepsilonpowersshiftscalexps o/var/www/html/mdtn/previsions/meteo_cartes/venv/lib/python3.11/site-packages/scipy/interpolate/_rbfinterp_np.pyrrs% 2 1fgvue  c"t||Sr )_pythran_polynomial_matrix)r rrs rpolynomial_matrixrs %a 0 00rct||}tj|tj}t |dkr|d|}|S)N)dtyper)rnpasarrayint64lenreshape)ndimdegreerouts r_monomial_powersr$sN f - -C *S ) ) )C 3xx1}}kk!T"" Jrc*t||||||Sr )_pythran_build_system)rd smoothingrrrrs rrr!s Ay&'6 J JJrc t|||||||\}}} } t||dd\} } } } | dkrtd| d| dkrid}|jd}|dkrEt || z | z ||}t j|}||kr d|d|d }t|| | | fS) aBuild and solve the RBF interpolation system of equations. Parameters ---------- y : (P, N) float ndarray Data point coordinates. d : (P, S) float ndarray Data values at `y`. smoothing : (P,) float ndarray Smoothing parameter for each data point. kernel : str Name of the RBF. epsilon : float Shape parameter. powers : (R, N) int ndarray The exponents for each monomial in the polynomial. Returns ------- coeffs : (P + R, S) float ndarray Coefficients for each RBF and monomial. shift : (N,) float ndarray Domain shift used to create the polynomial matrix. scale : (N,) float ndarray Domain scaling used to create the polynomial matrix. T) overwrite_a overwrite_brzThe z"-th argument had an illegal value.zSingular matrix.zqSingular matrix. The matrix of monomials evaluated at the data point coordinates does not have full column rank (/z).) rr ValueErrorshaperrlinalg matrix_rankr)rr'r(rrrrlhsrhsrr_coeffsinfomsgnmonospmatranks r_build_and_solve_systemr:%s8+ 1i&"  CeUsCTtLLLAq&$ axxIIIIJJJ  a A::$a%i%6CCD9((..Df}}/!//$*/// # % rc 8t||||||||} | |zSr )r) r rrrrrrr4rvecs rcompute_interpolationr=Xs. ( 1fgvueR  C <r)numpyr numpy.linalgrscipy.linalg.lapackr_rbfinterp_commonr_rbfinterp_pythranrr&rr r rrr$r:r=rrrDs$$$$$$%%%%%%444444111KKK1 1 1 fr