Ei hSSKrSSKJr SSKJr SSKJr SSKJ r J r J r Sr SrS rS r S rS rg) N) LinAlgError)dgesv)_monomial_powers_impl) _build_system_build_evaluation_coefficients_polynomial_matrixc [XX#XEU5$N)&_pythran_build_evaluation_coefficients)xykernelepsilonpowersshiftscalexps \/home/mande/repo/quber/.venv/lib/python3.13/site-packages/scipy/interpolate/_rbfinterp_np.pyrrs 2 fve c[X5$r )_pythran_polynomial_matrix)r rrs rpolynomial_matrixrs %a 00rc[X5n[R"U[RS9n[ U5S:XaUR SU5nU$)N)dtyper)rnpasarrayint64lenreshape)ndimdegreerouts r_monomial_powersr$s@  -C **S )C 3x1}kk!T" Jrc[XX#XE5$r )_pythran_build_system)rd smoothingrrrrs rrr!s y' JJrc D[XX#XEU5upxp[XxSSS9u pn U S:a[SU *S35eU S:aaSnURSnUS:a?[ X - U - XV5n[ R RU5nUU:a SUSUS 3n[U5eXU 4$) 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+ i" CesTtLAq$ ax4w&HIJJ  a A:$ai%6CD99((.Df}!F!F82/ #  rc &[XX#XEXh5n X-$r )r) r rrrrrrr4rvecs rcompute_interpolationr=Xs ( fve C <r)numpyr numpy.linalgrscipy.linalg.lapackr_rbfinterp_commonr_rbfinterp_pythranrr&rr r rrr$r:r=rrrDs;$%41K1 fr