+
    nDj                         ^ RI t^ RIHt ^ RIHt ^RIHt ^RIH	t
HtHt ]! R]P                  RR7      tR	 tR
 tR tR t	R tR tR# )    N)LinAlgError)get_lapack_funcs)_monomial_powers_impl)_build_system_build_evaluation_coefficients_polynomial_matrixsysv	preferred)dtypeilp64c           	          \        WW#WEV4      # N)&_pythran_build_evaluation_coefficients)xykernelepsilonpowersshiftscalexps   &&&&&&&&Y/data/cameron/venvs/s3viz/lib/python3.14/site-packages/scipy/interpolate/_rbfinterp_np.pyr   r      s     2	fve     c                     \        W4      # r   )_pythran_polynomial_matrix)r   r   r   s   &&&r   polynomial_matrixr      s    %a00r   c                     \        W4      p\        P                  ! V\        P                  R 7      p\	        V4      ^ 8X  d   VP                  ^ V 4      pV# ))r   )r   npasarrayint64lenreshape)ndimdegreer   outs   &&& r   _monomial_powersr&      s@    

-C
**S
)C
3x1}kk!T"Jr   c                     \        WW#WE4      # r   )_pythran_build_system)r   d	smoothingr   r   r   r   s   &&&&&&&r   r   r   $   s     y'JJr   c           	     j   \        WW#WEV4      w  rxr\        WxRRR7      w   rpV^ 8  d   \        RV)  R24      hV^ 8  dp   RpVP                  ^ ,          pV^ 8  dI   \	        W	,
          V
,          WV4      p\
        P                  P                  V4      pVV8  d
   RV RV R2p\        V4      hWV3# )	a  Build 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_bzThe z"-th argument had an illegal value.zSingular matrix.zqSingular matrix. The matrix of monomials evaluated at the data point coordinates does not have full column rank (/z).)	r   dsysv
ValueErrorshaper   r   linalgmatrix_rankr   )r   r)   r*   r   r   r   r   lhsrhsr   r   _coeffsinfomsgnmonospmatranks   &&&&&&&           r   _build_and_solve_systemr=   (   s    8 +	i"
Ce sTtLAq$ax4w&HIJJ	 aA:$ai%6CD99((.Df}!F!F82/  #r   c	           
      .    \        WW#WEWh4      p	W,          # r   )r   )
r   r   r   r   r   r   r   r7   r   vecs
   &&&&&&&&& r   compute_interpolationr@   [   s    
(	fveC <r   )numpyr   numpy.linalgr   scipy.linalg.lapackr   _rbfinterp_commonr   _rbfinterp_pythranr   r(   r   r   r   r   float64r/   r   r&   r=   r@    r   r   <module>rH      sO     $ 0 4  	rzzE1K1 fr   