§
    µ³Ùg4:  ã                   óÀ   — d dl Z d dlZd dlZd dlZd dlmZ ddlmZ d dl	m
c mZ d dlmZ ddlmZ dgZd„ Z G d	„ d¦  «        Zd
„ Zej        fd„Zdej        dœd„ZdS )é    N)Úprodé   )Ú_bspl)Ú	csr_array)Ú_not_a_knotÚ	NdBSplinec                 óp   — t          j        | t           j        ¦  «        rt           j        S t           j        S )z>Return np.complex128 for complex dtypes, np.float64 otherwise.)ÚnpÚ
issubdtypeÚcomplexfloatingÚ
complex128Úfloat64©Údtypes    úR/granary/edev/quantum/lib/python3.11/site-packages/scipy/interpolate/_ndbspline.pyÚ
_get_dtyper      s)   € å	„}�U�BÔ.Ñ/Ô/ð ÝŒ}ÐåŒzÐó    c                   óp   — e Zd ZdZddœd„Zed„ ¦   «         Zed„ ¦   «         Zdddœd„Ze	dd
„¦   «         Z
dS )r   a©  Tensor product spline object.

    The value at point ``xp = (x1, x2, ..., xN)`` is evaluated as a linear
    combination of products of one-dimensional b-splines in each of the ``N``
    dimensions::

       c[i1, i2, ..., iN] * B(x1; i1, t1) * B(x2; i2, t2) * ... * B(xN; iN, tN)


    Here ``B(x; i, t)`` is the ``i``-th b-spline defined by the knot vector
    ``t`` evaluated at ``x``.

    Parameters
    ----------
    t : tuple of 1D ndarrays
        knot vectors in directions 1, 2, ... N,
        ``len(t[i]) == n[i] + k + 1``
    c : ndarray, shape (n1, n2, ..., nN, ...)
        b-spline coefficients
    k : int or length-d tuple of integers
        spline degrees.
        A single integer is interpreted as having this degree for
        all dimensions.
    extrapolate : bool, optional
        Whether to extrapolate out-of-bounds inputs, or return `nan`.
        Default is to extrapolate.

    Attributes
    ----------
    t : tuple of ndarrays
        Knots vectors.
    c : ndarray
        Coefficients of the tensor-product spline.
    k : tuple of integers
        Degrees for each dimension.
    extrapolate : bool, optional
        Whether to extrapolate or return nans for out-of-bounds inputs.
        Defaults to true.

    Methods
    -------
    __call__
    design_matrix

    See Also
    --------
    BSpline : a one-dimensional B-spline object
    NdPPoly : an N-dimensional piecewise tensor product polynomial

    N)Úextrapolatec                ó¼  — t          ||¦  «        \  | _        | _        \  | _        | _        |€d}t          |¦  «        | _        t          j        |¦  «        | _	        | j        j
        d         }| j	        j        |k     rt          d|› d�¦  «        ‚t          |¦  «        D ]�}| j        |         }| j        |         }|j
        d         |z
  dz
  }	| j	        j
        |         |	k    r<t          d|› d| j	        j
        |         › dt!          |¦  «        › d	|	› d
|› d�¦  «        ‚Œ‚t#          | j	        j        ¦  «        }
t          j        | j	        |
¬¦  «        | _	        d S )NTr   zCoefficients must be at least z-dimensional.r   z,Knots, coefficients and degree in dimension z are inconsistent: got z coefficients for z knots, need at least z for k=ú.r   )Ú_preprocess_inputsÚ_kÚ_indices_k1dÚ_tÚ_len_tÚboolr   r
   ÚasarrayÚcÚshapeÚndimÚ
ValueErrorÚrangeÚtÚkÚlenr   r   Úascontiguousarray)Úselfr$   r   r%   r   r!   ÚdÚtdÚkdÚnÚdts              r   Ú__init__zNdBSpline.__init__M   s„  € Ý=OÐPQÐSTÑ=UÔ=UÑ:ˆŒ�Ô"Ñ$: T¤W¨d¬kàÐØˆKÝ Ñ,Ô,ˆÔå”˜A‘”ˆŒàŒwŒ}˜QÔˆØŒ6Œ;˜ÒÐÝÐQ¸dÐQÐQÐQÑRÔRÐRå�t‘”ð 
	-ð 
	-ˆAØ”˜”ˆBØ”˜”ˆBØ”˜”˜bÑ  1Ñ$ˆAàŒvŒ|˜AŒ !Ò#Ð#Ý ð ",Ø%&ð",ð ",à)-¬¬°a¬ð",ð ",õ &)¨¡W¤Wð",ð ",ð EFð",ð ",ð ()ð	",ð ",ð ",ñ -ô -ð -ð $õ ˜œœÑ%Ô%ˆÝÔ% d¤f°BÐ7Ñ7Ô7ˆŒˆˆr   c                 ó*   — t          | j        ¦  «        S ©N)Útupler   ©r(   s    r   r%   zNdBSpline.ki   s   € å�T”W‰~Œ~Ðr   c                 ót   ‡ — t          ˆ fd„t          ‰ j        j        d         ¦  «        D ¦   «         ¦  «        S )Nc              3   óP   •K  — | ] }‰j         |d ‰j        |         …f         V — Œ!d S r0   )r   r   )Ú.0r)   r(   s     €r   ú	<genexpr>zNdBSpline.t.<locals>.<genexpr>p   s9   øè è € ÐRÐR°Q�T”W˜Q  ¤¨Q¤ Ð/Ô0ÐRÐRÐRÐRÐRÐRr   r   )r1   r#   r   r    r2   s   `r   r$   zNdBSpline.tm   s9   ø€ õ ÐRÐRÐRÐR½%ÀÄÄÈaÔ@PÑ:QÔ:QÐRÑRÔRÑRÔRÐRr   )Únur   c                óŽ  ‡— | j         j        d         }|€| j        }t          |¦  «        }|€"t	          j        |ft          j        ¬¦  «        }n‰t	          j        |t          j        ¬¦  «        }|j        dk    s|j        d         |k    r(t          d|›dt          | j        ¦  «        › d�¦  «        ‚t          |dk     ¦  «        rt          d|›�¦  «        ‚t	          j        |t          ¬¦  «        }|j        }|                     d	|d	         ¦  «        }t	          j        |¦  «        }|d	         |k    rt          d
|› d|› �¦  «        ‚| j        j        j        dk    }| j        }|r| j        j        |k    r| j        d         }|                     t          ¦  «        }|                     |j        d|…         dz   ¦  «        Š‰                     ¦   «         }t	          j        ˆfd„‰j        D ¦   «         t          j        ¬¦  «        }	‰j        d	         }
t	          j        |j        dd	…         |
fz   ‰j        ¬¦  «        }t1          j        || j         | j        | j        ||||
|	| j        |¦  «         |                     | j        j        ¦  «        }|                     |dd	…         | j        j        |d…         z   ¦  «        S )a@  Evaluate the tensor product b-spline at ``xi``.

        Parameters
        ----------
        xi : array_like, shape(..., ndim)
            The coordinates to evaluate the interpolator at.
            This can be a list or tuple of ndim-dimensional points
            or an array with the shape (num_points, ndim).
        nu : array_like, optional, shape (ndim,)
            Orders of derivatives to evaluate. Each must be non-negative.
            Defaults to the zeroth derivivative.
        extrapolate : bool, optional
            Whether to exrapolate based on first and last intervals in each
            dimension, or return `nan`. Default is to ``self.extrapolate``.

        Returns
        -------
        values : ndarray, shape ``xi.shape[:-1] + self.c.shape[ndim:]``
            Interpolated values at ``xi``
        r   Nr   r   z)invalid number of derivative orders nu = z for ndim = r   z'derivatives must be positive, got nu = éÿÿÿÿzShapes: xi.shape=z
 and ndim=r   ).N)r9   c                 ó.   •— g | ]}|‰j         j        z  ‘ŒS © )r   Úitemsize)r5   ÚsÚc1s     €r   ú
<listcomp>z&NdBSpline.__call__.<locals>.<listcomp>¯   s3   ø€ ð "7ð "7ð "7Ø&'ð #$ r¤xÔ'8Ñ"8ð "7ð "7ð "7r   )r   r    r   r   r
   ÚzerosÚintcr   r!   r"   r&   r$   ÚanyÚfloatÚreshaper'   r   r   ÚkindÚviewÚravelÚstridesÚintpÚemptyr   Úevaluate_ndbspliner   r   r   )r(   Úxir7   r   r!   Úxi_shapeÚwas_complexÚccÚc1rÚ_strides_c1Únum_c_trÚoutr>   s               @r   Ú__call__zNdBSpline.__call__r   sÒ  ø€ ð* ŒwŒ}˜QÔˆàÐØÔ*ˆKÝ˜;Ñ'Ô'ˆàˆ:Ý”˜4˜'­¬Ð1Ñ1Ô1ˆBˆBå”˜B¥b¤gÐ.Ñ.Ô.ˆBØŒw˜!Š|ˆ|˜rœx¨œ{¨dÒ2Ð2Ý ð-¸2ð -ð -Ý! $¤&™kœkð-ð -ð -ñ.ô .ð .õ �2˜’6‰{Œ{ð OÝ Ð!MÀbÐ!MÐ!MÑNÔNÐNõ ŒZ˜¥%Ð(Ñ(Ô(ˆØ”8ˆØ�ZŠZ˜˜H RœLÑ)Ô)ˆÝÔ! "Ñ%Ô%ˆà�BŒ<˜4ÒÐÝÐK°ÐKÐKÀTÐKÐKÑLÔLÐLð ”f”lÔ'¨3Ò.ˆØŒVˆØð 	#˜4œ6œ;¨$Ò.Ð.ð ”˜	Ô"ˆBØ�WŠW•U‰^Œ^ˆð �ZŠZ˜œ  $ œ¨%Ñ/Ñ0Ô0ˆØ�hŠh‰jŒjˆõ ”jð "7ð "7ð "7ð "7Ø+-¬:ð"7ñ "7ô "7Ý>@¼gðGñ Gô Gˆð ”8˜B”<ˆÝŒh�r”x   ”}¨ {Ñ2¸"¼(ÐCÑCÔCˆåÔ  Ø!%¤Ø!%¤Ø!%¤Ø!#Ø!,Ø!$Ø!)Ø!,Ø!%Ô!2Ø!$ñ
	'ô 
	'ð 
	'ð �hŠh�t”v”|Ñ$Ô$ˆØ�{Š{˜8 C R Cœ=¨4¬6¬<¸¸¸Ô+>Ñ>Ñ?Ô?Ð?r   Tc                 óB  ‡‡— t          j        |t          ¬¦  «        }|j        d         }t	          |¦  «        |k    r#t          dt	          |¦  «        › d|›d�¦  «        ‚t          ‰|¦  «        \  Š}\  }Št          ˆˆfd„t          |¦  «        D ¦   «         ¦  «        }|dd…         d	z   }	t          j	        |	ddd…         t           j
        ¬¦  «        ddd…                              ¦   «         }
t          j        ||‰‰||
¦  «        \  }}}t          |||f¦  «        S )
a  Construct the design matrix as a CSR format sparse array.

        Parameters
        ----------
        xvals :  ndarray, shape(npts, ndim)
            Data points. ``xvals[j, :]`` gives the ``j``-th data point as an
            ``ndim``-dimensional array.
        t : tuple of 1D ndarrays, length-ndim
            Knot vectors in directions 1, 2, ... ndim,
        k : int
            B-spline degree.
        extrapolate : bool, optional
            Whether to extrapolate out-of-bounds values of raise a `ValueError`

        Returns
        -------
        design_matrix : a CSR array
            Each row of the design matrix corresponds to a value in `xvals` and
            contains values of b-spline basis elements which are non-zero
            at this value.

        r   r9   z*Data and knots are inconsistent: len(t) = z for  ndim = r   c              3   ó@   •K  — | ]}‰|         ‰|         z
  d z
  V — ŒdS ©r   Nr;   )r5   r)   r%   Úlen_ts     €€r   r6   z*NdBSpline.design_matrix.<locals>.<genexpr>é   s4   øè è € ÐAÐA°˜˜aœ 1 Q¤4™¨!Ñ+ÐAÐAÐAÐAÐAÐAr   r   N)r   )r
   r   rC   r    r&   r"   r   r1   r#   ÚcumprodrI   Úcopyr   Ú
_colloc_ndr   )ÚclsÚxvalsr$   r%   r   r!   r   r   Úc_shapeÚcsÚcstridesÚdataÚindicesÚindptrrX   s      `          @r   Údesign_matrixzNdBSpline.design_matrixÃ   sF  øø€ õ0 ”
˜5­Ð.Ñ.Ô.ˆØŒ{˜2ŒˆÝˆq‰6Œ6�TŠ>ˆ>Ýð½SÀ¹V¼Vð ð Øðð ð ñô ð õ (:¸!¸QÑ'?Ô'?Ñ$ˆˆ<™˜"˜eõ
 ÐAÐAÐAÐAÐAµU¸4±[´[ÐAÑAÔAÑAÔAˆð �Q�R�RŒ[˜4ÑˆÝ”:˜b   2 œh­b¬gÐ6Ñ6Ô6°t°t¸°tÔ<×AÒAÑCÔCˆõ !&Ô 0°Ø02Ø05Ø01Ø0<Ø08ñ!:ô !:Ñˆˆg�võ ˜$ ¨Ð0Ñ1Ô1Ð1r   )T)Ú__name__Ú
__module__Ú__qualname__Ú__doc__r.   Úpropertyr%   r$   rT   Úclassmethodrd   r;   r   r   r   r      sÁ   € € € € € ð1ð 1ðd 04ð 8ð 8ð 8ð 8ð 8ð8 ðð ñ „Xðð ðSð Sñ „XðSð "&°4ð O@ð O@ð O@ð O@ð O@ðb ð32ð 32ð 32ñ „[ð32ð 32ð 32r   c           
      óþ  — t          |t          ¦  «        st          d|› d�¦  «        ‚t          |¦  «        }	 t          | ¦  «         n# t          $ r	 | f|z  } Y nw xY wt          j        d„ | D ¦   «         t
          j        ¬¦  «        } t          | ¦  «        |k    r0t          dt          |¦  «        › dt          | ¦  «        ›d�¦  «        ‚t          |¦  «        }t          |¦  «        D �]P}t          j        ||         ¦  «        }| |         }|j	        d         |z
  d	z
  }|dk     rt          d
|› d�¦  «        ‚|j
        d	k    rt          d|› d�¦  «        ‚||d	z   k     rt          dd|z  dz   › d|› d|› d�¦  «        ‚t          j        |¦  «        dk                          ¦   «         rt          d|› d�¦  «        ‚t          t          j        |||d	z   …         ¦  «        ¦  «        dk     rt          d|› d�¦  «        ‚t          j        |¦  «                             ¦   «         st          d|› d�¦  «        ‚�ŒRt          d„ | D ¦   «         ¦  «        }t          j        t          j        t%          |¦  «        ¦  «        |¦  «        }t          j        |t
          j        ¬¦  «        j                             ¦   «         }	t          |¦  «        }d„ |D ¦   «         }
t          j        |t/          |
¦  «        ft0          ¬¦  «        }|                     t
          j        ¦  «         t          |¦  «        D ]$}||         ||dt          ||         ¦  «        …f<   Œ%t          j        |
t
          j        ¬¦  «        }
| |	||
ffS )zÁHelpers: validate and preprocess NdBSpline inputs.

       Parameters
       ----------
       k : int or tuple
          Spline orders
       t_tpl : tuple or array-likes
          Knots.
    z-Expect `t` to be a tuple of array-likes. Got z	 instead.c                 ó6   — g | ]}t          j        |¦  «        ‘ŒS r;   )ÚoperatorÚindex)r5   Úkis     r   r?   z&_preprocess_inputs.<locals>.<listcomp>  s"   € Ð3Ð3Ð3¨2•H”N 2Ñ&Ô&Ð3Ð3Ð3r   r   z	len(t) = z != len(k) = r   r   r   zSpline degree in dimension z cannot be negative.zKnot vector in dimension z must be one-dimensional.zNeed at least é   z knots for degree z in dimension zKnots in dimension z# must be in a non-decreasing order.z.Need at least two internal knots in dimension z should not have nans or infs.c              3   ó    K  — | ]	}|d z   V — Œ
dS rW   r;   )r5   r+   s     r   r6   z%_preprocess_inputs.<locals>.<genexpr>2  s&   è è € Ð%Ð%˜R�"�q‘&Ð%Ð%Ð%Ð%Ð%Ð%r   c                 ó,   — g | ]}t          |¦  «        ‘ŒS r;   ©r&   )r5   Útis     r   r?   z&_preprocess_inputs.<locals>.<listcomp><  s   € Ð%Ð%Ð%˜�S�‰WŒWÐ%Ð%Ð%r   N)Ú
isinstancer1   r"   r&   Ú	TypeErrorr
   r   Úint32r#   r    r!   ÚdiffrB   ÚuniqueÚisfiniteÚallÚunravel_indexÚaranger   rI   ÚTrZ   rJ   ÚmaxrC   ÚfillÚnan)r%   Út_tplr!   r)   r*   r+   r,   r    rb   r   rX   r   s               r   r   r   ú   sò  € õ �e�UÑ#Ô#ð 
Ýð 1Ø %ð1ð 1ð 1ñ 
ô 
ð 	
õ
 ˆu‰:Œ:€DðÝˆA‰ŒˆˆøÝð ð ð àˆD�‰Iˆˆˆðøøøõ 	Œ
Ð3Ð3°Ð3Ñ3Ô3½2¼8ÐDÑDÔD€Aå
ˆ1�v„v�‚~€~ÝÐA¥S¨¡Z¤ZÐAÐAµS¸±V´VÐAÐAÐAÑBÔBÐBõ ˆu‰:Œ:€DÝ�4‰[Œ[ð 0ñ 0ˆÝŒZ˜˜aœÑ!Ô!ˆØˆqŒTˆØŒH�QŒK˜"Ñ˜qÑ ˆØ�Š6ˆ6Ýð +¸1ð +ð +ð +ñ ,ô ,ð ,àŒ7�aŠ<ˆ<Ýð 2¸ð 2ð 2ð 2ñ 3ô 3ð 3àˆr�A‰vŠ:ˆ:Ýð 8¨a°©d°Q©hð 8ð 8Ø!#ð8ð 8Ø34ð8ð 8ð 8ñ 9ô 9ð 9åŒG�B‰KŒK˜!ŠO× Ò Ñ"Ô"ð 	8Ýð 7°1ð 7ð 7ð 7ñ 8ô 8ð 8å�rŒy˜˜B˜q 1™u˜HœÑ&Ô&Ñ'Ô'¨!Ò+Ð+Ýð 0Ø+,ð0ð 0ð 0ñ 1ô 1ð 1åŒ{˜2‰Œ×"Ò"Ñ$Ô$ð 	0Ýð /°1ð /ð /ð /ñ 0ô 0ð 0ñ	0õ Ð%Ð% 1Ð%Ñ%Ô%Ñ%Ô%€EÝÔ�rœy­¨e©¬Ñ5Ô5°uÑ=Ô=€GÝ”:˜g­R¬WÐ5Ñ5Ô5Ô7×<Ò<Ñ>Ô>€Lõ ˆu‰:Œ:€DØ%Ð%˜uÐ%Ñ%Ô%€EÝ	Œ�4�˜U™œÐ$­EÐ	2Ñ	2Ô	2€BØ‡G‚G�BŒF�O„O€OÝ�4‰[Œ[ð )ð )ˆØ % a¤ˆˆ1ˆn�s�5˜”8‰}Œ}ˆnÐÑÐÝŒJ�u¥B¤HÐ-Ñ-Ô-€Eàˆl˜R ˜KÐ'Ð's   ¹A	 Á	AÁAc           
      ó  — t          j        |j        t           j        ¦  «        r0t	          | |j        |fi |¤Ž}t	          | |j        |fi |¤Ž}|d|z  z   S |j        dk    r�|j        d         dk    rpt          j	        |¦  «        }t          |j        d         ¦  «        D ]?} || |d d …|f         fi |¤Ž\  |d d …|f<   }|dk    rt          d|›d|›d|› d�¦  «        ‚Œ@|S  || |fi |¤Ž\  }}|dk    rt          d|›d	|›d�¦  «        ‚|S )
Ny              ð?rp   r   r   z	solver = z returns info =z for column r   z returns info = )r
   r   r   r   Ú_iter_solveÚrealÚimagr!   r    Ú
empty_liker#   r"   )	ÚaÚbÚsolverÚsolver_argsr…   r†   ÚresÚjÚinfos	            r   r„   r„   F  sg  € õ
 
„}�Q”W�bÔ0Ñ1Ô1ð Ý˜1˜aœf fÐ<Ð<°Ð<Ð<ˆÝ˜1˜aœf fÐ<Ð<°Ð<Ð<ˆØ�b˜‘g‰~Ðà„v�‚{€{�q”w˜q”z A’~�~ÝŒm˜AÑÔˆÝ�q”w˜q”zÑ"Ô"ð 	Rð 	RˆAØ$˜f Q¨¨!¨!¨!¨Q¨$¬Ð?Ð?°;Ð?Ð?‰OˆC����1�‰I�tØ�qŠyˆyÝ Ð!P FÐ!PÐ!P¸Ð!PÐ!PÈAÐ!PÐ!PÐ!PÑQÔQÐQð àˆ
à�F˜1˜aÐ/Ð/ ;Ð/Ð/‰	ˆˆTØ�1Š9ˆ9ÝÐ> Ð>Ð>°DÐ>Ð>Ð>Ñ?Ô?Ð?Øˆ
r   é   ©rŠ   c                óÜ  ‡ ‡— t          ‰ ¦  «        }t          d„ ‰ D ¦   «         ¦  «        }	 t          ‰¦  «         n# t          $ r	 ‰f|z  ŠY nw xY wt          ‰ ¦  «        D ]]\  }}t          t	          j        |¦  «        ¦  «        }	|	‰|         k    r+t          d|	› d|› d‰|         › d‰|         dz   › d�	¦  «        ‚Œ^t          ˆˆ fd„t          |¦  «        D ¦   «         ¦  «        }
t	          j        d	„ t          j
        ‰ Ž D ¦   «         t          ¬
¦  «        }t                               ||
‰¦  «        }|j        }t          |d|…         ¦  «        t          ||d…         ¦  «        f}|                     |¦  «        }|t"          j        k    r$t'          j        t*          |¬¦  «        }d|vrd|d<    |||fi |¤Ž}|                     |||d…         z   ¦  «        }t          |
|‰¦  «        S )a™  Construct an interpolating NdBspline.

    Parameters
    ----------
    points : tuple of ndarrays of float, with shapes (m1,), ... (mN,)
        The points defining the regular grid in N dimensions. The points in
        each dimension (i.e. every element of the `points` tuple) must be
        strictly ascending or descending.      
    values : ndarray of float, shape (m1, ..., mN, ...)
        The data on the regular grid in n dimensions.
    k : int, optional
        The spline degree. Must be odd. Default is cubic, k=3
    solver : a `scipy.sparse.linalg` solver (iterative or direct), optional.
        An iterative solver from `scipy.sparse.linalg` or a direct one,
        `sparse.sparse.linalg.spsolve`.
        Used to solve the sparse linear system
        ``design_matrix @ coefficients = rhs`` for the coefficients.
        Default is `scipy.sparse.linalg.gcrotmk`
    solver_args : dict, optional
        Additional arguments for the solver. The call signature is
        ``solver(csr_array, rhs_vector, **solver_args)``

    Returns
    -------
    spl : NdBSpline object

    Notes
    -----
    Boundary conditions are not-a-knot in all dimensions.
    c              3   ó4   K  — | ]}t          |¦  «        V — Œd S r0   rs   )r5   Úxs     r   r6   zmake_ndbspl.<locals>.<genexpr>~  s(   è è € Ð,Ð, •S˜‘V”VÐ,Ð,Ð,Ð,Ð,Ð,r   z
There are z points in dimension z, but order z requires at least  r   z points per dimension.c              3   ó‚   •K  — | ]9}t          t          j        ‰|         t          ¬ ¦  «        ‰|         ¦  «        V — Œ:dS )r   N)r   r
   r   rC   )r5   r)   r%   Úpointss     €€r   r6   zmake_ndbspl.<locals>.<genexpr>�  sX   øè è € ð $ð $Øõ �"œ* V¨A¤YµeÐ<Ñ<Ô<¸aÀ¼dÑCÔCð $ð $ð $ð $ð $ð $r   c                 ó   — g | ]}|‘ŒS r;   r;   )r5   Úxvs     r   r?   zmake_ndbspl.<locals>.<listcomp>�  s   € Ð@Ð@Ð@˜r˜Ð@Ð@Ð@r   r   Nr�   Úatolg�íµ ÷Æ°>)r&   r1   rv   Ú	enumerater
   Ú
atleast_1dr"   r#   r   Ú	itertoolsÚproductrC   r   rd   r    r   rD   ÚsslÚspsolveÚ	functoolsÚpartialr„   )r•   Úvaluesr%   rŠ   r‹   r!   rM   r)   ÚpointÚnumptsr$   r]   ÚmatrÚv_shapeÚ
vals_shapeÚvalsÚcoefs   ` `              r   Úmake_ndbsplr©   ^  sb  øø€ õ> ˆv‰;Œ;€DÝÐ,Ð, VÐ,Ñ,Ô,Ñ,Ô,€HðÝˆA‰ŒˆˆøÝð ð ð àˆD�‰Iˆˆˆðøøøõ ˜fÑ%Ô%ð Að A‰ˆˆ5Ý•R”] 5Ñ)Ô)Ñ*Ô*ˆØ�Q�q”TŠ>ˆ>Ýð @¨&ð @ð @Àqð @ð @Ø+,¨Q¬4ð@ð @à!" 1¤ a¡ð@ð @ð @ñ Aô Að Að õ
 	ð $ð $ð $ð $ð $Ý˜T‘{”{ð$ñ $ô $ñ 	$ô 	$€AåŒJÐ@Ð@¥YÔ%6¸Ð%?Ð@Ñ@Ô@ÍÐNÑNÔN€Eõ ×"Ò" 5¨!¨QÑ/Ô/€Dð
 Œl€GÝ�w˜u ˜u”~Ñ&Ô&­¨W°T°U°U¬^Ñ(<Ô(<Ð=€JØ�>Š>˜*Ñ%Ô%€Dà•”ÒÐÝÔ"¥;°vÐ>Ñ>Ô>ˆØ˜Ð$Ð$à"&ˆK˜Ñàˆ6�$˜Ð,Ð, Ð,Ð,€DØ�<Š<˜ 7¨4¨5¨5¤>Ñ1Ñ2Ô2€DÝ�Q˜˜aÑ Ô Ð s   ¬< ¼AÁA)r�   )r›   rŸ   rm   Únumpyr
   Úmathr   Ú r   Úscipy.sparse.linalgÚsparseÚlinalgr�   Úscipy.sparser   Ú	_bsplinesr   Ú__all__r   r   r   Úgcrotmkr„   r©   r;   r   r   ú<module>r´      s>  ðØ Ð Ð Ð Ø Ð Ð Ð Ø €€€Ø Ð Ð Ð à Ð Ð Ð Ð Ð à Ð Ð Ð Ð Ð à !Ð !Ð !Ð !Ð !Ð !Ð !Ð !Ð !Ø "Ð "Ð "Ð "Ð "Ð "à "Ð "Ð "Ð "Ð "Ð "àˆ-€ðð ð ð]2ð ]2ð ]2ð ]2ð ]2ñ ]2ô ]2ð ]2ð@I(ð I(ð I(ðX !œ[ð ð ð ð ð0E!¨s¬{ð E!ð E!ð E!ð E!ð E!ð E!ð E!r   