+ LV-j!fR.t^RIt^RIHt^RIt^RIHt]'d^RIHt Rt RRRllt R#)geometric_slerpN) TYPE_CHECKING) euclideanc8\P!W.4p\PPVP4wrE^\P !V4^8,^, pVPVPR\P 3,,pVPVPR\P 3,,p\P!W4p\PPV4p\P!W4p Vwr\P!W),4p\P!W),4pWR\P 3,,WR\P 3,,,#):NNN) npvstacklinalgqrTdiagnewaxisdotdetarctan2sincos) startendtbasisQRsignscsomegas &&& o/Volumes/fast/ai/experiments/ui-tars-smoke/.venv/lib/python3.14/site-packages/scipy/spatial/_geometric_slerp.py_geometric_slerpr s IIul #E 99<< DA q !A %E eggam$$A eggam$$A uA aA JJq EJE qyA qyA Q ]# #cam,<&< <<c PV^8dQhRRRRRRR\R\P/#)rrz npt.ArrayLikerrtolreturn)floatrndarray)formats"r __annotate__r&!sE_#_# _# _#_#  _# ZZ _#rc\\P!V\PR7p\P!V\PR7p\P!V4pVP^8d \ R4hVP^8wgVP^8wd \ R4hVP VP 8wd \ R4hVP ^8gVP ^8d \ R4h\P !W4'd"\P!WVP 4#W3FKp\P!\PPV4RR^R7'dKB\ R 4h \V\4'g \ R 4h\P!V4p\W4p\P!VR ^VR7'd\P !R ^R 7\P!V\PR7pVP ^8Xd#\P"!^VP 34#VP^8Xd0\%VV\P&!V44P)4#\%VVV4#)a Geometric spherical linear interpolation. The interpolation occurs along a unit-radius great circle arc in arbitrary dimensional space. Parameters ---------- start : (n_dimensions, ) array-like Single n-dimensional input coordinate in a 1-D array-like object. `n` must be greater than 1. end : (n_dimensions, ) array-like Single n-dimensional input coordinate in a 1-D array-like object. `n` must be greater than 1. t : float or (n_points,) 1D array-like A float or 1D array-like of doubles representing interpolation parameters. A common approach is to generate the array with ``np.linspace(0, 1, n_pts)`` for linearly spaced points. Ascending, descending, and scrambled orders are permitted. .. versionchanged:: 1.17.0 Extrapolation is permitted, allowing values below `0` and above `1`. tol : float The absolute tolerance for determining if the start and end coordinates are antipodes. Returns ------- result : (t.size, D) An array of doubles containing the interpolated spherical path and including start and end when 0 and 1 t are used. The interpolated values should correspond to the same sort order provided in the t array. The result may be 1-dimensional if ``t`` is a float. Raises ------ ValueError If ``start`` or ``end`` are not on the unit n-sphere, or for a variety of degenerate conditions. See Also -------- scipy.spatial.transform.Slerp : 3-D Slerp that works with quaternions Notes ----- The implementation is based on the mathematical formula provided in [1]_, and the first known presentation of this algorithm, derived from study of 4-D geometry, is credited to Glenn Davis in a footnote of the original quaternion Slerp publication by Ken Shoemake [2]_. .. versionadded:: 1.5.0 References ---------- .. [1] https://en.wikipedia.org/wiki/Slerp#Geometric_Slerp .. [2] Ken Shoemake (1985) Animating rotation with quaternion curves. ACM SIGGRAPH Computer Graphics, 19(3): 245-254. Examples -------- Interpolate four linearly-spaced values on the circumference of a circle spanning 90 degrees: >>> import numpy as np >>> from scipy.spatial import geometric_slerp >>> import matplotlib.pyplot as plt >>> fig = plt.figure() >>> ax = fig.add_subplot(111) >>> start = np.array([1, 0]) >>> end = np.array([0, 1]) >>> t_vals = np.linspace(0, 1, 4) >>> result = geometric_slerp(start, ... end, ... t_vals) The interpolated results should be at 30 degree intervals recognizable on the unit circle: >>> ax.scatter(result[...,0], result[...,1], c='k') >>> circle = plt.Circle((0, 0), 1, color='grey') >>> ax.add_artist(circle) >>> ax.set_aspect('equal') >>> plt.show() Attempting to interpolate between antipodes on a circle is ambiguous because there are two possible paths, and on a sphere there are infinite possible paths on the geodesic surface. Nonetheless, one of the ambiguous paths is returned along with a warning: >>> import warnings >>> opposite_pole = np.array([-1, 0]) >>> with warnings.catch_warnings(): ... warnings.simplefilter("ignore", UserWarning) ... geometric_slerp(start, ... opposite_pole, ... t_vals) array([[ 1.00000000e+00, 0.00000000e+00], [ 5.00000000e-01, 8.66025404e-01], [-5.00000000e-01, 8.66025404e-01], [-1.00000000e+00, 1.22464680e-16]]) Extend the original example to a sphere and plot interpolation points in 3D: >>> from mpl_toolkits.mplot3d import proj3d >>> fig = plt.figure() >>> ax = fig.add_subplot(111, projection='3d') Plot the unit sphere for reference (optional): >>> u = np.linspace(0, 2 * np.pi, 100) >>> v = np.linspace(0, np.pi, 100) >>> x = np.outer(np.cos(u), np.sin(v)) >>> y = np.outer(np.sin(u), np.sin(v)) >>> z = np.outer(np.ones(np.size(u)), np.cos(v)) >>> ax.plot_surface(x, y, z, color='y', alpha=0.1) Interpolating over a larger number of points may provide the appearance of a smooth curve on the surface of the sphere, which is also useful for discretized integration calculations on a sphere surface: >>> start = np.array([1, 0, 0]) >>> end = np.array([0, 0, 1]) >>> t_vals = np.linspace(0, 1, 200) >>> result = geometric_slerp(start, ... end, ... t_vals) >>> ax.plot(result[...,0], ... result[...,1], ... result[...,2], ... c='k') >>> plt.show() It is also possible to perform extrapolations outside the interpolation interval by using interpolation parameter values below 0 or above 1. For example, the above example may be adjusted to extrapolate to an antipodal position: >>> fig = plt.figure() >>> ax = fig.add_subplot(111, projection='3d') >>> ax.plot_surface(x, y, z, color='y', alpha=0.1) >>> start = np.array([1, 0, 0]) >>> end = np.array([0, 0, 1]) >>> t_vals = np.linspace(0, 2, 400) >>> result = geometric_slerp(start, ... end, ... t_vals) >>> ax.plot(result[...,0], result[...,1], result[...,2], c='k') >>> plt.show() )dtypez:The interpolation parameter value must be one dimensional.z1Start and end coordinates must be one-dimensionalz;The dimensions of start and end must match (have same size)zLThe start and end coordinates must both be in at least two-dimensional spaceg?g& .>)rtolatolz(start and end are not on a unit n-sphereztol must be a floatg@z_start and end are antipodes using the specified tolerance; this may cause ambiguous slerp paths) stacklevel)rasarrayfloat64ndim ValueErrorsize array_equallinspaceallcloser norm isinstancer#fabsrwarningswarnemptyr atleast_1dravel)rrrr!coord coord_dists&&&& rrr!sH JJuBJJ /E **S +C 1 Avvz:; ; zzQ#((a-34 4 zzSXX;< < zzA~A!" " ~~e!!{{500{{299>>%0# $ !##34 4  c5 ! !.//ggcl5&J {{:s55 ="# $ 1BJJ'Avv{xxEJJ((vv{ # " a 0227%' :  # !# #r)gHz>) __all__r7typingrnumpyrscipy.spatial.distancer numpy.typingnptrrrrrEs0   ,=(_#_#r