PyPI - scipy - Versions diffs - 1.16.2__cp311-cp311-win_arm64.whl - Mend

scipy 1.16.2__cp311-cp311-win_arm64.whl

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.

Files changed (1530) hide show

scipy/__config__.py +161 -0
scipy/__init__.py +150 -0
scipy/_cyutility.cp311-win_arm64.lib +0 -0
scipy/_cyutility.cp311-win_arm64.pyd +0 -0
scipy/_distributor_init.py +18 -0
scipy/_lib/__init__.py +14 -0
scipy/_lib/_array_api.py +931 -0
scipy/_lib/_array_api_compat_vendor.py +9 -0
scipy/_lib/_array_api_no_0d.py +103 -0
scipy/_lib/_bunch.py +229 -0
scipy/_lib/_ccallback.py +251 -0
scipy/_lib/_ccallback_c.cp311-win_arm64.lib +0 -0
scipy/_lib/_ccallback_c.cp311-win_arm64.pyd +0 -0
scipy/_lib/_disjoint_set.py +254 -0
scipy/_lib/_docscrape.py +761 -0
scipy/_lib/_elementwise_iterative_method.py +346 -0
scipy/_lib/_fpumode.cp311-win_arm64.lib +0 -0
scipy/_lib/_fpumode.cp311-win_arm64.pyd +0 -0
scipy/_lib/_gcutils.py +105 -0
scipy/_lib/_pep440.py +487 -0
scipy/_lib/_sparse.py +41 -0
scipy/_lib/_test_ccallback.cp311-win_arm64.lib +0 -0
scipy/_lib/_test_ccallback.cp311-win_arm64.pyd +0 -0
scipy/_lib/_test_deprecation_call.cp311-win_arm64.lib +0 -0
scipy/_lib/_test_deprecation_call.cp311-win_arm64.pyd +0 -0
scipy/_lib/_test_deprecation_def.cp311-win_arm64.lib +0 -0
scipy/_lib/_test_deprecation_def.cp311-win_arm64.pyd +0 -0
scipy/_lib/_testutils.py +373 -0
scipy/_lib/_threadsafety.py +58 -0
scipy/_lib/_tmpdirs.py +86 -0
scipy/_lib/_uarray/LICENSE +29 -0
scipy/_lib/_uarray/__init__.py +116 -0
scipy/_lib/_uarray/_backend.py +707 -0
scipy/_lib/_uarray/_uarray.cp311-win_arm64.lib +0 -0
scipy/_lib/_uarray/_uarray.cp311-win_arm64.pyd +0 -0
scipy/_lib/_util.py +1283 -0
scipy/_lib/array_api_compat/__init__.py +22 -0
scipy/_lib/array_api_compat/_internal.py +59 -0
scipy/_lib/array_api_compat/common/__init__.py +1 -0
scipy/_lib/array_api_compat/common/_aliases.py +727 -0
scipy/_lib/array_api_compat/common/_fft.py +213 -0
scipy/_lib/array_api_compat/common/_helpers.py +1058 -0
scipy/_lib/array_api_compat/common/_linalg.py +232 -0
scipy/_lib/array_api_compat/common/_typing.py +192 -0
scipy/_lib/array_api_compat/cupy/__init__.py +13 -0
scipy/_lib/array_api_compat/cupy/_aliases.py +156 -0
scipy/_lib/array_api_compat/cupy/_info.py +336 -0
scipy/_lib/array_api_compat/cupy/_typing.py +31 -0
scipy/_lib/array_api_compat/cupy/fft.py +36 -0
scipy/_lib/array_api_compat/cupy/linalg.py +49 -0
scipy/_lib/array_api_compat/dask/__init__.py +0 -0
scipy/_lib/array_api_compat/dask/array/__init__.py +12 -0
scipy/_lib/array_api_compat/dask/array/_aliases.py +376 -0
scipy/_lib/array_api_compat/dask/array/_info.py +416 -0
scipy/_lib/array_api_compat/dask/array/fft.py +21 -0
scipy/_lib/array_api_compat/dask/array/linalg.py +72 -0
scipy/_lib/array_api_compat/numpy/__init__.py +28 -0
scipy/_lib/array_api_compat/numpy/_aliases.py +190 -0
scipy/_lib/array_api_compat/numpy/_info.py +366 -0
scipy/_lib/array_api_compat/numpy/_typing.py +30 -0
scipy/_lib/array_api_compat/numpy/fft.py +35 -0
scipy/_lib/array_api_compat/numpy/linalg.py +143 -0
scipy/_lib/array_api_compat/torch/__init__.py +22 -0
scipy/_lib/array_api_compat/torch/_aliases.py +855 -0
scipy/_lib/array_api_compat/torch/_info.py +369 -0
scipy/_lib/array_api_compat/torch/_typing.py +3 -0
scipy/_lib/array_api_compat/torch/fft.py +85 -0
scipy/_lib/array_api_compat/torch/linalg.py +121 -0
scipy/_lib/array_api_extra/__init__.py +38 -0
scipy/_lib/array_api_extra/_delegation.py +171 -0
scipy/_lib/array_api_extra/_lib/__init__.py +1 -0
scipy/_lib/array_api_extra/_lib/_at.py +463 -0
scipy/_lib/array_api_extra/_lib/_backends.py +46 -0
scipy/_lib/array_api_extra/_lib/_funcs.py +937 -0
scipy/_lib/array_api_extra/_lib/_lazy.py +357 -0
scipy/_lib/array_api_extra/_lib/_testing.py +278 -0
scipy/_lib/array_api_extra/_lib/_utils/__init__.py +1 -0
scipy/_lib/array_api_extra/_lib/_utils/_compat.py +74 -0
scipy/_lib/array_api_extra/_lib/_utils/_compat.pyi +45 -0
scipy/_lib/array_api_extra/_lib/_utils/_helpers.py +559 -0
scipy/_lib/array_api_extra/_lib/_utils/_typing.py +10 -0
scipy/_lib/array_api_extra/_lib/_utils/_typing.pyi +105 -0
scipy/_lib/array_api_extra/testing.py +359 -0
scipy/_lib/cobyqa/__init__.py +20 -0
scipy/_lib/cobyqa/framework.py +1240 -0
scipy/_lib/cobyqa/main.py +1506 -0
scipy/_lib/cobyqa/models.py +1529 -0
scipy/_lib/cobyqa/problem.py +1296 -0
scipy/_lib/cobyqa/settings.py +132 -0
scipy/_lib/cobyqa/subsolvers/__init__.py +14 -0
scipy/_lib/cobyqa/subsolvers/geometry.py +387 -0
scipy/_lib/cobyqa/subsolvers/optim.py +1203 -0
scipy/_lib/cobyqa/utils/__init__.py +18 -0
scipy/_lib/cobyqa/utils/exceptions.py +22 -0
scipy/_lib/cobyqa/utils/math.py +77 -0
scipy/_lib/cobyqa/utils/versions.py +67 -0
scipy/_lib/decorator.py +399 -0
scipy/_lib/deprecation.py +274 -0
scipy/_lib/doccer.py +366 -0
scipy/_lib/messagestream.cp311-win_arm64.lib +0 -0
scipy/_lib/messagestream.cp311-win_arm64.pyd +0 -0
scipy/_lib/pyprima/__init__.py +212 -0
scipy/_lib/pyprima/cobyla/__init__.py +0 -0
scipy/_lib/pyprima/cobyla/cobyla.py +559 -0
scipy/_lib/pyprima/cobyla/cobylb.py +714 -0
scipy/_lib/pyprima/cobyla/geometry.py +226 -0
scipy/_lib/pyprima/cobyla/initialize.py +215 -0
scipy/_lib/pyprima/cobyla/trustregion.py +492 -0
scipy/_lib/pyprima/cobyla/update.py +289 -0
scipy/_lib/pyprima/common/__init__.py +0 -0
scipy/_lib/pyprima/common/_bounds.py +34 -0
scipy/_lib/pyprima/common/_linear_constraints.py +46 -0
scipy/_lib/pyprima/common/_nonlinear_constraints.py +54 -0
scipy/_lib/pyprima/common/_project.py +173 -0
scipy/_lib/pyprima/common/checkbreak.py +93 -0
scipy/_lib/pyprima/common/consts.py +47 -0
scipy/_lib/pyprima/common/evaluate.py +99 -0
scipy/_lib/pyprima/common/history.py +38 -0
scipy/_lib/pyprima/common/infos.py +30 -0
scipy/_lib/pyprima/common/linalg.py +435 -0
scipy/_lib/pyprima/common/message.py +290 -0
scipy/_lib/pyprima/common/powalg.py +131 -0
scipy/_lib/pyprima/common/preproc.py +277 -0
scipy/_lib/pyprima/common/present.py +5 -0
scipy/_lib/pyprima/common/ratio.py +54 -0
scipy/_lib/pyprima/common/redrho.py +47 -0
scipy/_lib/pyprima/common/selectx.py +296 -0
scipy/_lib/tests/__init__.py +0 -0
scipy/_lib/tests/test__gcutils.py +110 -0
scipy/_lib/tests/test__pep440.py +67 -0
scipy/_lib/tests/test__testutils.py +32 -0
scipy/_lib/tests/test__threadsafety.py +51 -0
scipy/_lib/tests/test__util.py +641 -0
scipy/_lib/tests/test_array_api.py +322 -0
scipy/_lib/tests/test_bunch.py +169 -0
scipy/_lib/tests/test_ccallback.py +196 -0
scipy/_lib/tests/test_config.py +45 -0
scipy/_lib/tests/test_deprecation.py +10 -0
scipy/_lib/tests/test_doccer.py +143 -0
scipy/_lib/tests/test_import_cycles.py +18 -0
scipy/_lib/tests/test_public_api.py +482 -0
scipy/_lib/tests/test_scipy_version.py +28 -0
scipy/_lib/tests/test_tmpdirs.py +48 -0
scipy/_lib/tests/test_warnings.py +137 -0
scipy/_lib/uarray.py +31 -0
scipy/cluster/__init__.py +31 -0
scipy/cluster/_hierarchy.cp311-win_arm64.lib +0 -0
scipy/cluster/_hierarchy.cp311-win_arm64.pyd +0 -0
scipy/cluster/_optimal_leaf_ordering.cp311-win_arm64.lib +0 -0
scipy/cluster/_optimal_leaf_ordering.cp311-win_arm64.pyd +0 -0
scipy/cluster/_vq.cp311-win_arm64.lib +0 -0
scipy/cluster/_vq.cp311-win_arm64.pyd +0 -0
scipy/cluster/hierarchy.py +4348 -0
scipy/cluster/tests/__init__.py +0 -0
scipy/cluster/tests/hierarchy_test_data.py +145 -0
scipy/cluster/tests/test_disjoint_set.py +202 -0
scipy/cluster/tests/test_hierarchy.py +1238 -0
scipy/cluster/tests/test_vq.py +434 -0
scipy/cluster/vq.py +832 -0
scipy/conftest.py +683 -0
scipy/constants/__init__.py +358 -0
scipy/constants/_codata.py +2266 -0
scipy/constants/_constants.py +369 -0
scipy/constants/codata.py +21 -0
scipy/constants/constants.py +53 -0
scipy/constants/tests/__init__.py +0 -0
scipy/constants/tests/test_codata.py +78 -0
scipy/constants/tests/test_constants.py +83 -0
scipy/datasets/__init__.py +90 -0
scipy/datasets/_download_all.py +71 -0
scipy/datasets/_fetchers.py +225 -0
scipy/datasets/_registry.py +26 -0
scipy/datasets/_utils.py +81 -0
scipy/datasets/tests/__init__.py +0 -0
scipy/datasets/tests/test_data.py +128 -0
scipy/differentiate/__init__.py +27 -0
scipy/differentiate/_differentiate.py +1129 -0
scipy/differentiate/tests/__init__.py +0 -0
scipy/differentiate/tests/test_differentiate.py +694 -0
scipy/fft/__init__.py +114 -0
scipy/fft/_backend.py +196 -0
scipy/fft/_basic.py +1650 -0
scipy/fft/_basic_backend.py +197 -0
scipy/fft/_debug_backends.py +22 -0
scipy/fft/_fftlog.py +223 -0
scipy/fft/_fftlog_backend.py +200 -0
scipy/fft/_helper.py +348 -0
scipy/fft/_pocketfft/LICENSE.md +25 -0
scipy/fft/_pocketfft/__init__.py +9 -0
scipy/fft/_pocketfft/basic.py +251 -0
scipy/fft/_pocketfft/helper.py +249 -0
scipy/fft/_pocketfft/pypocketfft.cp311-win_arm64.lib +0 -0
scipy/fft/_pocketfft/pypocketfft.cp311-win_arm64.pyd +0 -0
scipy/fft/_pocketfft/realtransforms.py +109 -0
scipy/fft/_pocketfft/tests/__init__.py +0 -0
scipy/fft/_pocketfft/tests/test_basic.py +1011 -0
scipy/fft/_pocketfft/tests/test_real_transforms.py +505 -0
scipy/fft/_realtransforms.py +706 -0
scipy/fft/_realtransforms_backend.py +63 -0
scipy/fft/tests/__init__.py +0 -0
scipy/fft/tests/mock_backend.py +96 -0
scipy/fft/tests/test_backend.py +98 -0
scipy/fft/tests/test_basic.py +504 -0
scipy/fft/tests/test_fftlog.py +215 -0
scipy/fft/tests/test_helper.py +558 -0
scipy/fft/tests/test_multithreading.py +84 -0
scipy/fft/tests/test_real_transforms.py +247 -0
scipy/fftpack/__init__.py +103 -0
scipy/fftpack/_basic.py +428 -0
scipy/fftpack/_helper.py +115 -0
scipy/fftpack/_pseudo_diffs.py +554 -0
scipy/fftpack/_realtransforms.py +598 -0
scipy/fftpack/basic.py +20 -0
scipy/fftpack/convolve.cp311-win_arm64.lib +0 -0
scipy/fftpack/convolve.cp311-win_arm64.pyd +0 -0
scipy/fftpack/helper.py +19 -0
scipy/fftpack/pseudo_diffs.py +22 -0
scipy/fftpack/realtransforms.py +19 -0
scipy/fftpack/tests/__init__.py +0 -0
scipy/fftpack/tests/fftw_double_ref.npz +0 -0
scipy/fftpack/tests/fftw_longdouble_ref.npz +0 -0
scipy/fftpack/tests/fftw_single_ref.npz +0 -0
scipy/fftpack/tests/test.npz +0 -0
scipy/fftpack/tests/test_basic.py +877 -0
scipy/fftpack/tests/test_helper.py +54 -0
scipy/fftpack/tests/test_import.py +33 -0
scipy/fftpack/tests/test_pseudo_diffs.py +388 -0
scipy/fftpack/tests/test_real_transforms.py +836 -0
scipy/integrate/__init__.py +122 -0
scipy/integrate/_bvp.py +1160 -0
scipy/integrate/_cubature.py +729 -0
scipy/integrate/_dop.cp311-win_arm64.lib +0 -0
scipy/integrate/_dop.cp311-win_arm64.pyd +0 -0
scipy/integrate/_ivp/__init__.py +8 -0
scipy/integrate/_ivp/base.py +290 -0
scipy/integrate/_ivp/bdf.py +478 -0
scipy/integrate/_ivp/common.py +451 -0
scipy/integrate/_ivp/dop853_coefficients.py +193 -0
scipy/integrate/_ivp/ivp.py +755 -0
scipy/integrate/_ivp/lsoda.py +224 -0
scipy/integrate/_ivp/radau.py +572 -0
scipy/integrate/_ivp/rk.py +601 -0
scipy/integrate/_ivp/tests/__init__.py +0 -0
scipy/integrate/_ivp/tests/test_ivp.py +1287 -0
scipy/integrate/_ivp/tests/test_rk.py +37 -0
scipy/integrate/_lebedev.py +5450 -0
scipy/integrate/_lsoda.cp311-win_arm64.lib +0 -0
scipy/integrate/_lsoda.cp311-win_arm64.pyd +0 -0
scipy/integrate/_ode.py +1395 -0
scipy/integrate/_odepack.cp311-win_arm64.lib +0 -0
scipy/integrate/_odepack.cp311-win_arm64.pyd +0 -0
scipy/integrate/_odepack_py.py +273 -0
scipy/integrate/_quad_vec.py +674 -0
scipy/integrate/_quadpack.cp311-win_arm64.lib +0 -0
scipy/integrate/_quadpack.cp311-win_arm64.pyd +0 -0
scipy/integrate/_quadpack_py.py +1283 -0
scipy/integrate/_quadrature.py +1336 -0
scipy/integrate/_rules/__init__.py +12 -0
scipy/integrate/_rules/_base.py +518 -0
scipy/integrate/_rules/_gauss_kronrod.py +202 -0
scipy/integrate/_rules/_gauss_legendre.py +62 -0
scipy/integrate/_rules/_genz_malik.py +210 -0
scipy/integrate/_tanhsinh.py +1385 -0
scipy/integrate/_test_multivariate.cp311-win_arm64.lib +0 -0
scipy/integrate/_test_multivariate.cp311-win_arm64.pyd +0 -0
scipy/integrate/_test_odeint_banded.cp311-win_arm64.lib +0 -0
scipy/integrate/_test_odeint_banded.cp311-win_arm64.pyd +0 -0
scipy/integrate/_vode.cp311-win_arm64.lib +0 -0
scipy/integrate/_vode.cp311-win_arm64.pyd +0 -0
scipy/integrate/dop.py +15 -0
scipy/integrate/lsoda.py +15 -0
scipy/integrate/odepack.py +17 -0
scipy/integrate/quadpack.py +23 -0
scipy/integrate/tests/__init__.py +0 -0
scipy/integrate/tests/test__quad_vec.py +211 -0
scipy/integrate/tests/test_banded_ode_solvers.py +305 -0
scipy/integrate/tests/test_bvp.py +714 -0
scipy/integrate/tests/test_cubature.py +1375 -0
scipy/integrate/tests/test_integrate.py +840 -0
scipy/integrate/tests/test_odeint_jac.py +74 -0
scipy/integrate/tests/test_quadpack.py +680 -0
scipy/integrate/tests/test_quadrature.py +730 -0
scipy/integrate/tests/test_tanhsinh.py +1171 -0
scipy/integrate/vode.py +15 -0
scipy/interpolate/__init__.py +228 -0
scipy/interpolate/_bary_rational.py +715 -0
scipy/interpolate/_bsplines.py +2469 -0
scipy/interpolate/_cubic.py +973 -0
scipy/interpolate/_dfitpack.cp311-win_arm64.lib +0 -0
scipy/interpolate/_dfitpack.cp311-win_arm64.pyd +0 -0
scipy/interpolate/_dierckx.cp311-win_arm64.lib +0 -0
scipy/interpolate/_dierckx.cp311-win_arm64.pyd +0 -0
scipy/interpolate/_fitpack.cp311-win_arm64.lib +0 -0
scipy/interpolate/_fitpack.cp311-win_arm64.pyd +0 -0
scipy/interpolate/_fitpack2.py +2397 -0
scipy/interpolate/_fitpack_impl.py +811 -0
scipy/interpolate/_fitpack_py.py +898 -0
scipy/interpolate/_fitpack_repro.py +996 -0
scipy/interpolate/_interpnd.cp311-win_arm64.lib +0 -0
scipy/interpolate/_interpnd.cp311-win_arm64.pyd +0 -0
scipy/interpolate/_interpolate.py +2266 -0
scipy/interpolate/_ndbspline.py +415 -0
scipy/interpolate/_ndgriddata.py +329 -0
scipy/interpolate/_pade.py +67 -0
scipy/interpolate/_polyint.py +1025 -0
scipy/interpolate/_ppoly.cp311-win_arm64.lib +0 -0
scipy/interpolate/_ppoly.cp311-win_arm64.pyd +0 -0
scipy/interpolate/_rbf.py +290 -0
scipy/interpolate/_rbfinterp.py +550 -0
scipy/interpolate/_rbfinterp_pythran.cp311-win_arm64.lib +0 -0
scipy/interpolate/_rbfinterp_pythran.cp311-win_arm64.pyd +0 -0
scipy/interpolate/_rgi.py +764 -0
scipy/interpolate/_rgi_cython.cp311-win_arm64.lib +0 -0
scipy/interpolate/_rgi_cython.cp311-win_arm64.pyd +0 -0
scipy/interpolate/dfitpack.py +24 -0
scipy/interpolate/fitpack.py +31 -0
scipy/interpolate/fitpack2.py +29 -0
scipy/interpolate/interpnd.py +24 -0
scipy/interpolate/interpolate.py +30 -0
scipy/interpolate/ndgriddata.py +23 -0
scipy/interpolate/polyint.py +24 -0
scipy/interpolate/rbf.py +18 -0
scipy/interpolate/tests/__init__.py +0 -0
scipy/interpolate/tests/data/bug-1310.npz +0 -0
scipy/interpolate/tests/data/estimate_gradients_hang.npy +0 -0
scipy/interpolate/tests/data/gcvspl.npz +0 -0
scipy/interpolate/tests/test_bary_rational.py +368 -0
scipy/interpolate/tests/test_bsplines.py +3754 -0
scipy/interpolate/tests/test_fitpack.py +519 -0
scipy/interpolate/tests/test_fitpack2.py +1431 -0
scipy/interpolate/tests/test_gil.py +64 -0
scipy/interpolate/tests/test_interpnd.py +452 -0
scipy/interpolate/tests/test_interpolate.py +2630 -0
scipy/interpolate/tests/test_ndgriddata.py +308 -0
scipy/interpolate/tests/test_pade.py +107 -0
scipy/interpolate/tests/test_polyint.py +972 -0
scipy/interpolate/tests/test_rbf.py +246 -0
scipy/interpolate/tests/test_rbfinterp.py +534 -0
scipy/interpolate/tests/test_rgi.py +1151 -0
scipy/io/__init__.py +116 -0
scipy/io/_fast_matrix_market/__init__.py +600 -0
scipy/io/_fast_matrix_market/_fmm_core.cp311-win_arm64.lib +0 -0
scipy/io/_fast_matrix_market/_fmm_core.cp311-win_arm64.pyd +0 -0
scipy/io/_fortran.py +354 -0
scipy/io/_harwell_boeing/__init__.py +7 -0
scipy/io/_harwell_boeing/_fortran_format_parser.py +316 -0
scipy/io/_harwell_boeing/hb.py +571 -0
scipy/io/_harwell_boeing/tests/__init__.py +0 -0
scipy/io/_harwell_boeing/tests/test_fortran_format.py +74 -0
scipy/io/_harwell_boeing/tests/test_hb.py +70 -0
scipy/io/_idl.py +917 -0
scipy/io/_mmio.py +968 -0
scipy/io/_netcdf.py +1104 -0
scipy/io/_test_fortran.cp311-win_arm64.lib +0 -0
scipy/io/_test_fortran.cp311-win_arm64.pyd +0 -0
scipy/io/arff/__init__.py +28 -0
scipy/io/arff/_arffread.py +873 -0
scipy/io/arff/arffread.py +19 -0
scipy/io/arff/tests/__init__.py +0 -0
scipy/io/arff/tests/data/iris.arff +225 -0
scipy/io/arff/tests/data/missing.arff +8 -0
scipy/io/arff/tests/data/nodata.arff +11 -0
scipy/io/arff/tests/data/quoted_nominal.arff +13 -0
scipy/io/arff/tests/data/quoted_nominal_spaces.arff +13 -0
scipy/io/arff/tests/data/test1.arff +10 -0
scipy/io/arff/tests/data/test10.arff +8 -0
scipy/io/arff/tests/data/test11.arff +11 -0
scipy/io/arff/tests/data/test2.arff +15 -0
scipy/io/arff/tests/data/test3.arff +6 -0
scipy/io/arff/tests/data/test4.arff +11 -0
scipy/io/arff/tests/data/test5.arff +26 -0
scipy/io/arff/tests/data/test6.arff +12 -0
scipy/io/arff/tests/data/test7.arff +15 -0
scipy/io/arff/tests/data/test8.arff +12 -0
scipy/io/arff/tests/data/test9.arff +14 -0
scipy/io/arff/tests/test_arffread.py +421 -0
scipy/io/harwell_boeing.py +17 -0
scipy/io/idl.py +17 -0
scipy/io/matlab/__init__.py +66 -0
scipy/io/matlab/_byteordercodes.py +75 -0
scipy/io/matlab/_mio.py +375 -0
scipy/io/matlab/_mio4.py +632 -0
scipy/io/matlab/_mio5.py +901 -0
scipy/io/matlab/_mio5_params.py +281 -0
scipy/io/matlab/_mio5_utils.cp311-win_arm64.lib +0 -0
scipy/io/matlab/_mio5_utils.cp311-win_arm64.pyd +0 -0
scipy/io/matlab/_mio_utils.cp311-win_arm64.lib +0 -0
scipy/io/matlab/_mio_utils.cp311-win_arm64.pyd +0 -0
scipy/io/matlab/_miobase.py +435 -0
scipy/io/matlab/_streams.cp311-win_arm64.lib +0 -0
scipy/io/matlab/_streams.cp311-win_arm64.pyd +0 -0
scipy/io/matlab/byteordercodes.py +17 -0
scipy/io/matlab/mio.py +16 -0
scipy/io/matlab/mio4.py +17 -0
scipy/io/matlab/mio5.py +19 -0
scipy/io/matlab/mio5_params.py +18 -0
scipy/io/matlab/mio5_utils.py +17 -0
scipy/io/matlab/mio_utils.py +17 -0
scipy/io/matlab/miobase.py +16 -0
scipy/io/matlab/streams.py +16 -0
scipy/io/matlab/tests/__init__.py +0 -0
scipy/io/matlab/tests/data/bad_miuint32.mat +0 -0
scipy/io/matlab/tests/data/bad_miutf8_array_name.mat +0 -0
scipy/io/matlab/tests/data/big_endian.mat +0 -0
scipy/io/matlab/tests/data/broken_utf8.mat +0 -0
scipy/io/matlab/tests/data/corrupted_zlib_checksum.mat +0 -0
scipy/io/matlab/tests/data/corrupted_zlib_data.mat +0 -0
scipy/io/matlab/tests/data/debigged_m4.mat +0 -0
scipy/io/matlab/tests/data/japanese_utf8.txt +5 -0
scipy/io/matlab/tests/data/little_endian.mat +0 -0
scipy/io/matlab/tests/data/logical_sparse.mat +0 -0
scipy/io/matlab/tests/data/malformed1.mat +0 -0
scipy/io/matlab/tests/data/miuint32_for_miint32.mat +0 -0
scipy/io/matlab/tests/data/miutf8_array_name.mat +0 -0
scipy/io/matlab/tests/data/nasty_duplicate_fieldnames.mat +0 -0
scipy/io/matlab/tests/data/one_by_zero_char.mat +0 -0
scipy/io/matlab/tests/data/parabola.mat +0 -0
scipy/io/matlab/tests/data/single_empty_string.mat +0 -0
scipy/io/matlab/tests/data/some_functions.mat +0 -0
scipy/io/matlab/tests/data/sqr.mat +0 -0
scipy/io/matlab/tests/data/test3dmatrix_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/test3dmatrix_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/test3dmatrix_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/test3dmatrix_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/test_empty_struct.mat +0 -0
scipy/io/matlab/tests/data/test_mat4_le_floats.mat +0 -0
scipy/io/matlab/tests/data/test_skip_variable.mat +0 -0
scipy/io/matlab/tests/data/testbool_8_WIN64.mat +0 -0
scipy/io/matlab/tests/data/testcell_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testcell_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testcell_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testcell_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testcellnest_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testcellnest_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testcellnest_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testcellnest_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testcomplex_4.2c_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testcomplex_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testcomplex_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testcomplex_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testcomplex_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testdouble_4.2c_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testdouble_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testdouble_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testdouble_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testdouble_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testemptycell_5.3_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testemptycell_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testemptycell_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testemptycell_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testfunc_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testhdf5_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testmatrix_4.2c_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testmatrix_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testmatrix_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testmatrix_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testmatrix_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testminus_4.2c_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testminus_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testminus_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testminus_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testminus_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testmulti_4.2c_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testmulti_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testmulti_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testobject_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testobject_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testobject_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testobject_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testonechar_4.2c_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testonechar_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testonechar_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testonechar_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testonechar_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testscalarcell_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testsimplecell.mat +0 -0
scipy/io/matlab/tests/data/testsparse_4.2c_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testsparse_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testsparse_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testsparse_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testsparse_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testsparsecomplex_4.2c_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testsparsecomplex_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/testsparsecomplex_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testsparsecomplex_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testsparsecomplex_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testsparsefloat_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststring_4.2c_SOL2.mat +0 -0
scipy/io/matlab/tests/data/teststring_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/teststring_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststring_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststring_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststringarray_4.2c_SOL2.mat +0 -0
scipy/io/matlab/tests/data/teststringarray_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/teststringarray_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststringarray_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststringarray_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststruct_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/teststruct_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststruct_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststruct_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststructarr_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/teststructarr_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststructarr_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststructarr_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststructnest_6.1_SOL2.mat +0 -0
scipy/io/matlab/tests/data/teststructnest_6.5.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststructnest_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/teststructnest_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testunicode_7.1_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testunicode_7.4_GLNX86.mat +0 -0
scipy/io/matlab/tests/data/testvec_4_GLNX86.mat +0 -0
scipy/io/matlab/tests/test_byteordercodes.py +29 -0
scipy/io/matlab/tests/test_mio.py +1399 -0
scipy/io/matlab/tests/test_mio5_utils.py +179 -0
scipy/io/matlab/tests/test_mio_funcs.py +51 -0
scipy/io/matlab/tests/test_mio_utils.py +45 -0
scipy/io/matlab/tests/test_miobase.py +32 -0
scipy/io/matlab/tests/test_pathological.py +33 -0
scipy/io/matlab/tests/test_streams.py +241 -0
scipy/io/mmio.py +17 -0
scipy/io/netcdf.py +17 -0
scipy/io/tests/__init__.py +0 -0
scipy/io/tests/data/Transparent Busy.ani +0 -0
scipy/io/tests/data/array_float32_1d.sav +0 -0
scipy/io/tests/data/array_float32_2d.sav +0 -0
scipy/io/tests/data/array_float32_3d.sav +0 -0
scipy/io/tests/data/array_float32_4d.sav +0 -0
scipy/io/tests/data/array_float32_5d.sav +0 -0
scipy/io/tests/data/array_float32_6d.sav +0 -0
scipy/io/tests/data/array_float32_7d.sav +0 -0
scipy/io/tests/data/array_float32_8d.sav +0 -0
scipy/io/tests/data/array_float32_pointer_1d.sav +0 -0
scipy/io/tests/data/array_float32_pointer_2d.sav +0 -0
scipy/io/tests/data/array_float32_pointer_3d.sav +0 -0
scipy/io/tests/data/array_float32_pointer_4d.sav +0 -0
scipy/io/tests/data/array_float32_pointer_5d.sav +0 -0
scipy/io/tests/data/array_float32_pointer_6d.sav +0 -0
scipy/io/tests/data/array_float32_pointer_7d.sav +0 -0
scipy/io/tests/data/array_float32_pointer_8d.sav +0 -0
scipy/io/tests/data/example_1.nc +0 -0
scipy/io/tests/data/example_2.nc +0 -0
scipy/io/tests/data/example_3_maskedvals.nc +0 -0
scipy/io/tests/data/fortran-3x3d-2i.dat +0 -0
scipy/io/tests/data/fortran-mixed.dat +0 -0
scipy/io/tests/data/fortran-sf8-11x1x10.dat +0 -0
scipy/io/tests/data/fortran-sf8-15x10x22.dat +0 -0
scipy/io/tests/data/fortran-sf8-1x1x1.dat +0 -0
scipy/io/tests/data/fortran-sf8-1x1x5.dat +0 -0
scipy/io/tests/data/fortran-sf8-1x1x7.dat +0 -0
scipy/io/tests/data/fortran-sf8-1x3x5.dat +0 -0
scipy/io/tests/data/fortran-si4-11x1x10.dat +0 -0
scipy/io/tests/data/fortran-si4-15x10x22.dat +0 -0
scipy/io/tests/data/fortran-si4-1x1x1.dat +0 -0
scipy/io/tests/data/fortran-si4-1x1x5.dat +0 -0
scipy/io/tests/data/fortran-si4-1x1x7.dat +0 -0
scipy/io/tests/data/fortran-si4-1x3x5.dat +0 -0
scipy/io/tests/data/invalid_pointer.sav +0 -0
scipy/io/tests/data/null_pointer.sav +0 -0
scipy/io/tests/data/scalar_byte.sav +0 -0
scipy/io/tests/data/scalar_byte_descr.sav +0 -0
scipy/io/tests/data/scalar_complex32.sav +0 -0
scipy/io/tests/data/scalar_complex64.sav +0 -0
scipy/io/tests/data/scalar_float32.sav +0 -0
scipy/io/tests/data/scalar_float64.sav +0 -0
scipy/io/tests/data/scalar_heap_pointer.sav +0 -0
scipy/io/tests/data/scalar_int16.sav +0 -0
scipy/io/tests/data/scalar_int32.sav +0 -0
scipy/io/tests/data/scalar_int64.sav +0 -0
scipy/io/tests/data/scalar_string.sav +0 -0
scipy/io/tests/data/scalar_uint16.sav +0 -0
scipy/io/tests/data/scalar_uint32.sav +0 -0
scipy/io/tests/data/scalar_uint64.sav +0 -0
scipy/io/tests/data/struct_arrays.sav +0 -0
scipy/io/tests/data/struct_arrays_byte_idl80.sav +0 -0
scipy/io/tests/data/struct_arrays_replicated.sav +0 -0
scipy/io/tests/data/struct_arrays_replicated_3d.sav +0 -0
scipy/io/tests/data/struct_inherit.sav +0 -0
scipy/io/tests/data/struct_pointer_arrays.sav +0 -0
scipy/io/tests/data/struct_pointer_arrays_replicated.sav +0 -0
scipy/io/tests/data/struct_pointer_arrays_replicated_3d.sav +0 -0
scipy/io/tests/data/struct_pointers.sav +0 -0
scipy/io/tests/data/struct_pointers_replicated.sav +0 -0
scipy/io/tests/data/struct_pointers_replicated_3d.sav +0 -0
scipy/io/tests/data/struct_scalars.sav +0 -0
scipy/io/tests/data/struct_scalars_replicated.sav +0 -0
scipy/io/tests/data/struct_scalars_replicated_3d.sav +0 -0
scipy/io/tests/data/test-1234Hz-le-1ch-10S-20bit-extra.wav +0 -0
scipy/io/tests/data/test-44100Hz-2ch-32bit-float-be.wav +0 -0
scipy/io/tests/data/test-44100Hz-2ch-32bit-float-le.wav +0 -0
scipy/io/tests/data/test-44100Hz-be-1ch-4bytes.wav +0 -0
scipy/io/tests/data/test-44100Hz-le-1ch-4bytes-early-eof-no-data.wav +0 -0
scipy/io/tests/data/test-44100Hz-le-1ch-4bytes-early-eof.wav +0 -0
scipy/io/tests/data/test-44100Hz-le-1ch-4bytes-incomplete-chunk.wav +0 -0
scipy/io/tests/data/test-44100Hz-le-1ch-4bytes-rf64.wav +0 -0
scipy/io/tests/data/test-44100Hz-le-1ch-4bytes.wav +0 -0
scipy/io/tests/data/test-48000Hz-2ch-64bit-float-le-wavex.wav +0 -0
scipy/io/tests/data/test-8000Hz-be-3ch-5S-24bit.wav +0 -0
scipy/io/tests/data/test-8000Hz-le-1ch-1byte-ulaw.wav +0 -0
scipy/io/tests/data/test-8000Hz-le-2ch-1byteu.wav +0 -0
scipy/io/tests/data/test-8000Hz-le-3ch-5S-24bit-inconsistent.wav +0 -0
scipy/io/tests/data/test-8000Hz-le-3ch-5S-24bit-rf64.wav +0 -0
scipy/io/tests/data/test-8000Hz-le-3ch-5S-24bit.wav +0 -0
scipy/io/tests/data/test-8000Hz-le-3ch-5S-36bit.wav +0 -0
scipy/io/tests/data/test-8000Hz-le-3ch-5S-45bit.wav +0 -0
scipy/io/tests/data/test-8000Hz-le-3ch-5S-53bit.wav +0 -0
scipy/io/tests/data/test-8000Hz-le-3ch-5S-64bit.wav +0 -0
scipy/io/tests/data/test-8000Hz-le-4ch-9S-12bit.wav +0 -0
scipy/io/tests/data/test-8000Hz-le-5ch-9S-5bit.wav +0 -0
scipy/io/tests/data/various_compressed.sav +0 -0
scipy/io/tests/test_fortran.py +264 -0
scipy/io/tests/test_idl.py +483 -0
scipy/io/tests/test_mmio.py +831 -0
scipy/io/tests/test_netcdf.py +550 -0
scipy/io/tests/test_paths.py +93 -0
scipy/io/tests/test_wavfile.py +501 -0
scipy/io/wavfile.py +938 -0
scipy/linalg/__init__.pxd +1 -0
scipy/linalg/__init__.py +236 -0
scipy/linalg/_basic.py +2146 -0
scipy/linalg/_blas_subroutines.h +164 -0
scipy/linalg/_cythonized_array_utils.cp311-win_arm64.lib +0 -0
scipy/linalg/_cythonized_array_utils.cp311-win_arm64.pyd +0 -0
scipy/linalg/_cythonized_array_utils.pxd +40 -0
scipy/linalg/_cythonized_array_utils.pyi +16 -0
scipy/linalg/_decomp.py +1645 -0
scipy/linalg/_decomp_cholesky.py +413 -0
scipy/linalg/_decomp_cossin.py +236 -0
scipy/linalg/_decomp_interpolative.cp311-win_arm64.lib +0 -0
scipy/linalg/_decomp_interpolative.cp311-win_arm64.pyd +0 -0
scipy/linalg/_decomp_ldl.py +356 -0
scipy/linalg/_decomp_lu.py +401 -0
scipy/linalg/_decomp_lu_cython.cp311-win_arm64.lib +0 -0
scipy/linalg/_decomp_lu_cython.cp311-win_arm64.pyd +0 -0
scipy/linalg/_decomp_lu_cython.pyi +6 -0
scipy/linalg/_decomp_polar.py +113 -0
scipy/linalg/_decomp_qr.py +494 -0
scipy/linalg/_decomp_qz.py +452 -0
scipy/linalg/_decomp_schur.py +336 -0
scipy/linalg/_decomp_svd.py +545 -0
scipy/linalg/_decomp_update.cp311-win_arm64.lib +0 -0
scipy/linalg/_decomp_update.cp311-win_arm64.pyd +0 -0
scipy/linalg/_expm_frechet.py +417 -0
scipy/linalg/_fblas.cp311-win_arm64.lib +0 -0
scipy/linalg/_fblas.cp311-win_arm64.pyd +0 -0
scipy/linalg/_flapack.cp311-win_arm64.lib +0 -0
scipy/linalg/_flapack.cp311-win_arm64.pyd +0 -0
scipy/linalg/_lapack_subroutines.h +1521 -0
scipy/linalg/_linalg_pythran.cp311-win_arm64.lib +0 -0
scipy/linalg/_linalg_pythran.cp311-win_arm64.pyd +0 -0
scipy/linalg/_matfuncs.py +1050 -0
scipy/linalg/_matfuncs_expm.cp311-win_arm64.lib +0 -0
scipy/linalg/_matfuncs_expm.cp311-win_arm64.pyd +0 -0
scipy/linalg/_matfuncs_expm.pyi +6 -0
scipy/linalg/_matfuncs_inv_ssq.py +886 -0
scipy/linalg/_matfuncs_schur_sqrtm.cp311-win_arm64.lib +0 -0
scipy/linalg/_matfuncs_schur_sqrtm.cp311-win_arm64.pyd +0 -0
scipy/linalg/_matfuncs_sqrtm.py +107 -0
scipy/linalg/_matfuncs_sqrtm_triu.cp311-win_arm64.lib +0 -0
scipy/linalg/_matfuncs_sqrtm_triu.cp311-win_arm64.pyd +0 -0
scipy/linalg/_misc.py +191 -0
scipy/linalg/_procrustes.py +113 -0
scipy/linalg/_sketches.py +189 -0
scipy/linalg/_solve_toeplitz.cp311-win_arm64.lib +0 -0
scipy/linalg/_solve_toeplitz.cp311-win_arm64.pyd +0 -0
scipy/linalg/_solvers.py +862 -0
scipy/linalg/_special_matrices.py +1322 -0
scipy/linalg/_testutils.py +65 -0
scipy/linalg/basic.py +23 -0
scipy/linalg/blas.py +495 -0
scipy/linalg/cython_blas.cp311-win_arm64.lib +0 -0
scipy/linalg/cython_blas.cp311-win_arm64.pyd +0 -0
scipy/linalg/cython_blas.pxd +169 -0
scipy/linalg/cython_blas.pyx +1432 -0
scipy/linalg/cython_lapack.cp311-win_arm64.lib +0 -0
scipy/linalg/cython_lapack.cp311-win_arm64.pyd +0 -0
scipy/linalg/cython_lapack.pxd +1528 -0
scipy/linalg/cython_lapack.pyx +12045 -0
scipy/linalg/decomp.py +23 -0
scipy/linalg/decomp_cholesky.py +21 -0
scipy/linalg/decomp_lu.py +21 -0
scipy/linalg/decomp_qr.py +20 -0
scipy/linalg/decomp_schur.py +21 -0
scipy/linalg/decomp_svd.py +21 -0
scipy/linalg/interpolative.py +989 -0
scipy/linalg/lapack.py +1081 -0
scipy/linalg/matfuncs.py +23 -0
scipy/linalg/misc.py +21 -0
scipy/linalg/special_matrices.py +22 -0
scipy/linalg/tests/__init__.py +0 -0
scipy/linalg/tests/_cython_examples/extending.pyx +23 -0
scipy/linalg/tests/_cython_examples/meson.build +34 -0
scipy/linalg/tests/data/carex_15_data.npz +0 -0
scipy/linalg/tests/data/carex_18_data.npz +0 -0
scipy/linalg/tests/data/carex_19_data.npz +0 -0
scipy/linalg/tests/data/carex_20_data.npz +0 -0
scipy/linalg/tests/data/carex_6_data.npz +0 -0
scipy/linalg/tests/data/gendare_20170120_data.npz +0 -0
scipy/linalg/tests/test_basic.py +2074 -0
scipy/linalg/tests/test_batch.py +588 -0
scipy/linalg/tests/test_blas.py +1127 -0
scipy/linalg/tests/test_cython_blas.py +118 -0
scipy/linalg/tests/test_cython_lapack.py +22 -0
scipy/linalg/tests/test_cythonized_array_utils.py +130 -0
scipy/linalg/tests/test_decomp.py +3189 -0
scipy/linalg/tests/test_decomp_cholesky.py +268 -0
scipy/linalg/tests/test_decomp_cossin.py +314 -0
scipy/linalg/tests/test_decomp_ldl.py +137 -0
scipy/linalg/tests/test_decomp_lu.py +308 -0
scipy/linalg/tests/test_decomp_polar.py +110 -0
scipy/linalg/tests/test_decomp_update.py +1701 -0
scipy/linalg/tests/test_extending.py +46 -0
scipy/linalg/tests/test_fblas.py +607 -0
scipy/linalg/tests/test_interpolative.py +232 -0
scipy/linalg/tests/test_lapack.py +3620 -0
scipy/linalg/tests/test_matfuncs.py +1125 -0
scipy/linalg/tests/test_matmul_toeplitz.py +136 -0
scipy/linalg/tests/test_procrustes.py +214 -0
scipy/linalg/tests/test_sketches.py +118 -0
scipy/linalg/tests/test_solve_toeplitz.py +150 -0
scipy/linalg/tests/test_solvers.py +844 -0
scipy/linalg/tests/test_special_matrices.py +636 -0
scipy/misc/__init__.py +6 -0
scipy/misc/common.py +6 -0
scipy/misc/doccer.py +6 -0
scipy/ndimage/__init__.py +174 -0
scipy/ndimage/_ctest.cp311-win_arm64.lib +0 -0
scipy/ndimage/_ctest.cp311-win_arm64.pyd +0 -0
scipy/ndimage/_cytest.cp311-win_arm64.lib +0 -0
scipy/ndimage/_cytest.cp311-win_arm64.pyd +0 -0
scipy/ndimage/_delegators.py +303 -0
scipy/ndimage/_filters.py +2422 -0
scipy/ndimage/_fourier.py +306 -0
scipy/ndimage/_interpolation.py +1033 -0
scipy/ndimage/_measurements.py +1689 -0
scipy/ndimage/_morphology.py +2634 -0
scipy/ndimage/_nd_image.cp311-win_arm64.lib +0 -0
scipy/ndimage/_nd_image.cp311-win_arm64.pyd +0 -0
scipy/ndimage/_ndimage_api.py +16 -0
scipy/ndimage/_ni_docstrings.py +214 -0
scipy/ndimage/_ni_label.cp311-win_arm64.lib +0 -0
scipy/ndimage/_ni_label.cp311-win_arm64.pyd +0 -0
scipy/ndimage/_ni_support.py +139 -0
scipy/ndimage/_rank_filter_1d.cp311-win_arm64.lib +0 -0
scipy/ndimage/_rank_filter_1d.cp311-win_arm64.pyd +0 -0
scipy/ndimage/_support_alternative_backends.py +84 -0
scipy/ndimage/filters.py +27 -0
scipy/ndimage/fourier.py +21 -0
scipy/ndimage/interpolation.py +22 -0
scipy/ndimage/measurements.py +24 -0
scipy/ndimage/morphology.py +27 -0
scipy/ndimage/tests/__init__.py +12 -0
scipy/ndimage/tests/data/label_inputs.txt +21 -0
scipy/ndimage/tests/data/label_results.txt +294 -0
scipy/ndimage/tests/data/label_strels.txt +42 -0
scipy/ndimage/tests/dots.png +0 -0
scipy/ndimage/tests/test_c_api.py +102 -0
scipy/ndimage/tests/test_datatypes.py +67 -0
scipy/ndimage/tests/test_filters.py +3083 -0
scipy/ndimage/tests/test_fourier.py +187 -0
scipy/ndimage/tests/test_interpolation.py +1491 -0
scipy/ndimage/tests/test_measurements.py +1592 -0
scipy/ndimage/tests/test_morphology.py +2950 -0
scipy/ndimage/tests/test_ni_support.py +78 -0
scipy/ndimage/tests/test_splines.py +70 -0
scipy/odr/__init__.py +131 -0
scipy/odr/__odrpack.cp311-win_arm64.lib +0 -0
scipy/odr/__odrpack.cp311-win_arm64.pyd +0 -0
scipy/odr/_add_newdocs.py +34 -0
scipy/odr/_models.py +315 -0
scipy/odr/_odrpack.py +1154 -0
scipy/odr/models.py +20 -0
scipy/odr/odrpack.py +21 -0
scipy/odr/tests/__init__.py +0 -0
scipy/odr/tests/test_odr.py +607 -0
scipy/optimize/__init__.pxd +1 -0
scipy/optimize/__init__.py +460 -0
scipy/optimize/_basinhopping.py +741 -0
scipy/optimize/_bglu_dense.cp311-win_arm64.lib +0 -0
scipy/optimize/_bglu_dense.cp311-win_arm64.pyd +0 -0
scipy/optimize/_bracket.py +706 -0
scipy/optimize/_chandrupatla.py +551 -0
scipy/optimize/_cobyla_py.py +297 -0
scipy/optimize/_cobyqa_py.py +72 -0
scipy/optimize/_constraints.py +598 -0
scipy/optimize/_dcsrch.py +728 -0
scipy/optimize/_differentiable_functions.py +835 -0
scipy/optimize/_differentialevolution.py +1970 -0
scipy/optimize/_direct.cp311-win_arm64.lib +0 -0
scipy/optimize/_direct.cp311-win_arm64.pyd +0 -0
scipy/optimize/_direct_py.py +280 -0
scipy/optimize/_dual_annealing.py +732 -0
scipy/optimize/_elementwise.py +798 -0
scipy/optimize/_group_columns.cp311-win_arm64.lib +0 -0
scipy/optimize/_group_columns.cp311-win_arm64.pyd +0 -0
scipy/optimize/_hessian_update_strategy.py +479 -0
scipy/optimize/_highspy/__init__.py +0 -0
scipy/optimize/_highspy/_core.cp311-win_arm64.lib +0 -0
scipy/optimize/_highspy/_core.cp311-win_arm64.pyd +0 -0
scipy/optimize/_highspy/_highs_options.cp311-win_arm64.lib +0 -0
scipy/optimize/_highspy/_highs_options.cp311-win_arm64.pyd +0 -0
scipy/optimize/_highspy/_highs_wrapper.py +338 -0
scipy/optimize/_isotonic.py +157 -0
scipy/optimize/_lbfgsb.cp311-win_arm64.lib +0 -0
scipy/optimize/_lbfgsb.cp311-win_arm64.pyd +0 -0
scipy/optimize/_lbfgsb_py.py +634 -0
scipy/optimize/_linesearch.py +896 -0
scipy/optimize/_linprog.py +733 -0
scipy/optimize/_linprog_doc.py +1434 -0
scipy/optimize/_linprog_highs.py +422 -0
scipy/optimize/_linprog_ip.py +1141 -0
scipy/optimize/_linprog_rs.py +572 -0
scipy/optimize/_linprog_simplex.py +663 -0
scipy/optimize/_linprog_util.py +1521 -0
scipy/optimize/_lsap.cp311-win_arm64.lib +0 -0
scipy/optimize/_lsap.cp311-win_arm64.pyd +0 -0
scipy/optimize/_lsq/__init__.py +5 -0
scipy/optimize/_lsq/bvls.py +183 -0
scipy/optimize/_lsq/common.py +731 -0
scipy/optimize/_lsq/dogbox.py +345 -0
scipy/optimize/_lsq/givens_elimination.cp311-win_arm64.lib +0 -0
scipy/optimize/_lsq/givens_elimination.cp311-win_arm64.pyd +0 -0
scipy/optimize/_lsq/least_squares.py +1044 -0
scipy/optimize/_lsq/lsq_linear.py +361 -0
scipy/optimize/_lsq/trf.py +587 -0
scipy/optimize/_lsq/trf_linear.py +249 -0
scipy/optimize/_milp.py +394 -0
scipy/optimize/_minimize.py +1199 -0
scipy/optimize/_minpack.cp311-win_arm64.lib +0 -0
scipy/optimize/_minpack.cp311-win_arm64.pyd +0 -0
scipy/optimize/_minpack_py.py +1178 -0
scipy/optimize/_moduleTNC.cp311-win_arm64.lib +0 -0
scipy/optimize/_moduleTNC.cp311-win_arm64.pyd +0 -0
scipy/optimize/_nnls.py +96 -0
scipy/optimize/_nonlin.py +1634 -0
scipy/optimize/_numdiff.py +963 -0
scipy/optimize/_optimize.py +4169 -0
scipy/optimize/_pava_pybind.cp311-win_arm64.lib +0 -0
scipy/optimize/_pava_pybind.cp311-win_arm64.pyd +0 -0
scipy/optimize/_qap.py +760 -0
scipy/optimize/_remove_redundancy.py +522 -0
scipy/optimize/_root.py +732 -0
scipy/optimize/_root_scalar.py +538 -0
scipy/optimize/_shgo.py +1606 -0
scipy/optimize/_shgo_lib/__init__.py +0 -0
scipy/optimize/_shgo_lib/_complex.py +1225 -0
scipy/optimize/_shgo_lib/_vertex.py +460 -0
scipy/optimize/_slsqp_py.py +603 -0
scipy/optimize/_slsqplib.cp311-win_arm64.lib +0 -0
scipy/optimize/_slsqplib.cp311-win_arm64.pyd +0 -0
scipy/optimize/_spectral.py +260 -0
scipy/optimize/_tnc.py +438 -0
scipy/optimize/_trlib/__init__.py +12 -0
scipy/optimize/_trlib/_trlib.cp311-win_arm64.lib +0 -0
scipy/optimize/_trlib/_trlib.cp311-win_arm64.pyd +0 -0
scipy/optimize/_trustregion.py +318 -0
scipy/optimize/_trustregion_constr/__init__.py +6 -0
scipy/optimize/_trustregion_constr/canonical_constraint.py +390 -0
scipy/optimize/_trustregion_constr/equality_constrained_sqp.py +231 -0
scipy/optimize/_trustregion_constr/minimize_trustregion_constr.py +584 -0
scipy/optimize/_trustregion_constr/projections.py +411 -0
scipy/optimize/_trustregion_constr/qp_subproblem.py +637 -0
scipy/optimize/_trustregion_constr/report.py +49 -0
scipy/optimize/_trustregion_constr/tests/__init__.py +0 -0
scipy/optimize/_trustregion_constr/tests/test_canonical_constraint.py +296 -0
scipy/optimize/_trustregion_constr/tests/test_nested_minimize.py +39 -0
scipy/optimize/_trustregion_constr/tests/test_projections.py +214 -0
scipy/optimize/_trustregion_constr/tests/test_qp_subproblem.py +645 -0
scipy/optimize/_trustregion_constr/tests/test_report.py +34 -0
scipy/optimize/_trustregion_constr/tr_interior_point.py +361 -0
scipy/optimize/_trustregion_dogleg.py +122 -0
scipy/optimize/_trustregion_exact.py +437 -0
scipy/optimize/_trustregion_krylov.py +65 -0
scipy/optimize/_trustregion_ncg.py +126 -0
scipy/optimize/_tstutils.py +972 -0
scipy/optimize/_zeros.cp311-win_arm64.lib +0 -0
scipy/optimize/_zeros.cp311-win_arm64.pyd +0 -0
scipy/optimize/_zeros_py.py +1475 -0
scipy/optimize/cobyla.py +19 -0
scipy/optimize/cython_optimize/__init__.py +133 -0
scipy/optimize/cython_optimize/_zeros.cp311-win_arm64.lib +0 -0
scipy/optimize/cython_optimize/_zeros.cp311-win_arm64.pyd +0 -0
scipy/optimize/cython_optimize/_zeros.pxd +33 -0
scipy/optimize/cython_optimize/c_zeros.pxd +26 -0
scipy/optimize/cython_optimize.pxd +11 -0
scipy/optimize/elementwise.py +38 -0
scipy/optimize/lbfgsb.py +23 -0
scipy/optimize/linesearch.py +18 -0
scipy/optimize/minpack.py +27 -0
scipy/optimize/minpack2.py +17 -0
scipy/optimize/moduleTNC.py +19 -0
scipy/optimize/nonlin.py +29 -0
scipy/optimize/optimize.py +40 -0
scipy/optimize/slsqp.py +22 -0
scipy/optimize/tests/__init__.py +0 -0
scipy/optimize/tests/_cython_examples/extending.pyx +43 -0
scipy/optimize/tests/_cython_examples/meson.build +32 -0
scipy/optimize/tests/test__basinhopping.py +535 -0
scipy/optimize/tests/test__differential_evolution.py +1703 -0
scipy/optimize/tests/test__dual_annealing.py +416 -0
scipy/optimize/tests/test__linprog_clean_inputs.py +312 -0
scipy/optimize/tests/test__numdiff.py +885 -0
scipy/optimize/tests/test__remove_redundancy.py +228 -0
scipy/optimize/tests/test__root.py +124 -0
scipy/optimize/tests/test__shgo.py +1164 -0
scipy/optimize/tests/test__spectral.py +226 -0
scipy/optimize/tests/test_bracket.py +896 -0
scipy/optimize/tests/test_chandrupatla.py +982 -0
scipy/optimize/tests/test_cobyla.py +195 -0
scipy/optimize/tests/test_cobyqa.py +252 -0
scipy/optimize/tests/test_constraint_conversion.py +286 -0
scipy/optimize/tests/test_constraints.py +255 -0
scipy/optimize/tests/test_cython_optimize.py +92 -0
scipy/optimize/tests/test_differentiable_functions.py +1025 -0
scipy/optimize/tests/test_direct.py +321 -0
scipy/optimize/tests/test_extending.py +28 -0
scipy/optimize/tests/test_hessian_update_strategy.py +300 -0
scipy/optimize/tests/test_isotonic_regression.py +167 -0
scipy/optimize/tests/test_lbfgsb_hessinv.py +65 -0
scipy/optimize/tests/test_lbfgsb_setulb.py +122 -0
scipy/optimize/tests/test_least_squares.py +986 -0
scipy/optimize/tests/test_linear_assignment.py +116 -0
scipy/optimize/tests/test_linesearch.py +328 -0
scipy/optimize/tests/test_linprog.py +2577 -0
scipy/optimize/tests/test_lsq_common.py +297 -0
scipy/optimize/tests/test_lsq_linear.py +287 -0
scipy/optimize/tests/test_milp.py +459 -0
scipy/optimize/tests/test_minimize_constrained.py +845 -0
scipy/optimize/tests/test_minpack.py +1194 -0
scipy/optimize/tests/test_nnls.py +469 -0
scipy/optimize/tests/test_nonlin.py +572 -0
scipy/optimize/tests/test_optimize.py +3344 -0
scipy/optimize/tests/test_quadratic_assignment.py +455 -0
scipy/optimize/tests/test_regression.py +40 -0
scipy/optimize/tests/test_slsqp.py +645 -0
scipy/optimize/tests/test_tnc.py +345 -0
scipy/optimize/tests/test_trustregion.py +110 -0
scipy/optimize/tests/test_trustregion_exact.py +351 -0
scipy/optimize/tests/test_trustregion_krylov.py +170 -0
scipy/optimize/tests/test_zeros.py +998 -0
scipy/optimize/tnc.py +22 -0
scipy/optimize/zeros.py +26 -0
scipy/signal/__init__.py +316 -0
scipy/signal/_arraytools.py +264 -0
scipy/signal/_czt.py +575 -0
scipy/signal/_delegators.py +568 -0
scipy/signal/_filter_design.py +5893 -0
scipy/signal/_fir_filter_design.py +1458 -0
scipy/signal/_lti_conversion.py +534 -0
scipy/signal/_ltisys.py +3546 -0
scipy/signal/_max_len_seq.py +139 -0
scipy/signal/_max_len_seq_inner.cp311-win_arm64.lib +0 -0
scipy/signal/_max_len_seq_inner.cp311-win_arm64.pyd +0 -0
scipy/signal/_peak_finding.py +1310 -0
scipy/signal/_peak_finding_utils.cp311-win_arm64.lib +0 -0
scipy/signal/_peak_finding_utils.cp311-win_arm64.pyd +0 -0
scipy/signal/_polyutils.py +172 -0
scipy/signal/_savitzky_golay.py +357 -0
scipy/signal/_short_time_fft.py +2228 -0
scipy/signal/_signal_api.py +30 -0
scipy/signal/_signaltools.py +5309 -0
scipy/signal/_sigtools.cp311-win_arm64.lib +0 -0
scipy/signal/_sigtools.cp311-win_arm64.pyd +0 -0
scipy/signal/_sosfilt.cp311-win_arm64.lib +0 -0
scipy/signal/_sosfilt.cp311-win_arm64.pyd +0 -0
scipy/signal/_spectral_py.py +2471 -0
scipy/signal/_spline.cp311-win_arm64.lib +0 -0
scipy/signal/_spline.cp311-win_arm64.pyd +0 -0
scipy/signal/_spline.pyi +34 -0
scipy/signal/_spline_filters.py +848 -0
scipy/signal/_support_alternative_backends.py +73 -0
scipy/signal/_upfirdn.py +219 -0
scipy/signal/_upfirdn_apply.cp311-win_arm64.lib +0 -0
scipy/signal/_upfirdn_apply.cp311-win_arm64.pyd +0 -0
scipy/signal/_waveforms.py +687 -0
scipy/signal/_wavelets.py +29 -0
scipy/signal/bsplines.py +21 -0
scipy/signal/filter_design.py +28 -0
scipy/signal/fir_filter_design.py +21 -0
scipy/signal/lti_conversion.py +20 -0
scipy/signal/ltisys.py +25 -0
scipy/signal/signaltools.py +27 -0
scipy/signal/spectral.py +21 -0
scipy/signal/spline.py +18 -0
scipy/signal/tests/__init__.py +0 -0
scipy/signal/tests/_scipy_spectral_test_shim.py +311 -0
scipy/signal/tests/mpsig.py +122 -0
scipy/signal/tests/test_array_tools.py +111 -0
scipy/signal/tests/test_bsplines.py +365 -0
scipy/signal/tests/test_cont2discrete.py +424 -0
scipy/signal/tests/test_czt.py +221 -0
scipy/signal/tests/test_dltisys.py +599 -0
scipy/signal/tests/test_filter_design.py +4744 -0
scipy/signal/tests/test_fir_filter_design.py +851 -0
scipy/signal/tests/test_ltisys.py +1225 -0
scipy/signal/tests/test_max_len_seq.py +71 -0
scipy/signal/tests/test_peak_finding.py +915 -0
scipy/signal/tests/test_result_type.py +51 -0
scipy/signal/tests/test_savitzky_golay.py +363 -0
scipy/signal/tests/test_short_time_fft.py +1107 -0
scipy/signal/tests/test_signaltools.py +4735 -0
scipy/signal/tests/test_spectral.py +2141 -0
scipy/signal/tests/test_splines.py +427 -0
scipy/signal/tests/test_upfirdn.py +322 -0
scipy/signal/tests/test_waveforms.py +400 -0
scipy/signal/tests/test_wavelets.py +59 -0
scipy/signal/tests/test_windows.py +987 -0
scipy/signal/waveforms.py +20 -0
scipy/signal/wavelets.py +17 -0
scipy/signal/windows/__init__.py +52 -0
scipy/signal/windows/_windows.py +2513 -0
scipy/signal/windows/windows.py +23 -0
scipy/sparse/__init__.py +350 -0
scipy/sparse/_base.py +1613 -0
scipy/sparse/_bsr.py +880 -0
scipy/sparse/_compressed.py +1328 -0
scipy/sparse/_construct.py +1454 -0
scipy/sparse/_coo.py +1581 -0
scipy/sparse/_csc.py +367 -0
scipy/sparse/_csparsetools.cp311-win_arm64.lib +0 -0
scipy/sparse/_csparsetools.cp311-win_arm64.pyd +0 -0
scipy/sparse/_csr.py +558 -0
scipy/sparse/_data.py +569 -0
scipy/sparse/_dia.py +677 -0
scipy/sparse/_dok.py +669 -0
scipy/sparse/_extract.py +178 -0
scipy/sparse/_index.py +444 -0
scipy/sparse/_lil.py +632 -0
scipy/sparse/_matrix.py +169 -0
scipy/sparse/_matrix_io.py +167 -0
scipy/sparse/_sparsetools.cp311-win_arm64.lib +0 -0
scipy/sparse/_sparsetools.cp311-win_arm64.pyd +0 -0
scipy/sparse/_spfuncs.py +76 -0
scipy/sparse/_sputils.py +632 -0
scipy/sparse/base.py +24 -0
scipy/sparse/bsr.py +22 -0
scipy/sparse/compressed.py +20 -0
scipy/sparse/construct.py +38 -0
scipy/sparse/coo.py +23 -0
scipy/sparse/csc.py +22 -0
scipy/sparse/csgraph/__init__.py +210 -0
scipy/sparse/csgraph/_flow.cp311-win_arm64.lib +0 -0
scipy/sparse/csgraph/_flow.cp311-win_arm64.pyd +0 -0
scipy/sparse/csgraph/_laplacian.py +563 -0
scipy/sparse/csgraph/_matching.cp311-win_arm64.lib +0 -0
scipy/sparse/csgraph/_matching.cp311-win_arm64.pyd +0 -0
scipy/sparse/csgraph/_min_spanning_tree.cp311-win_arm64.lib +0 -0
scipy/sparse/csgraph/_min_spanning_tree.cp311-win_arm64.pyd +0 -0
scipy/sparse/csgraph/_reordering.cp311-win_arm64.lib +0 -0
scipy/sparse/csgraph/_reordering.cp311-win_arm64.pyd +0 -0
scipy/sparse/csgraph/_shortest_path.cp311-win_arm64.lib +0 -0
scipy/sparse/csgraph/_shortest_path.cp311-win_arm64.pyd +0 -0
scipy/sparse/csgraph/_tools.cp311-win_arm64.lib +0 -0
scipy/sparse/csgraph/_tools.cp311-win_arm64.pyd +0 -0
scipy/sparse/csgraph/_traversal.cp311-win_arm64.lib +0 -0
scipy/sparse/csgraph/_traversal.cp311-win_arm64.pyd +0 -0
scipy/sparse/csgraph/_validation.py +66 -0
scipy/sparse/csgraph/tests/__init__.py +0 -0
scipy/sparse/csgraph/tests/test_connected_components.py +119 -0
scipy/sparse/csgraph/tests/test_conversions.py +61 -0
scipy/sparse/csgraph/tests/test_flow.py +209 -0
scipy/sparse/csgraph/tests/test_graph_laplacian.py +368 -0
scipy/sparse/csgraph/tests/test_matching.py +307 -0
scipy/sparse/csgraph/tests/test_pydata_sparse.py +197 -0
scipy/sparse/csgraph/tests/test_reordering.py +70 -0
scipy/sparse/csgraph/tests/test_shortest_path.py +540 -0
scipy/sparse/csgraph/tests/test_spanning_tree.py +66 -0
scipy/sparse/csgraph/tests/test_traversal.py +148 -0
scipy/sparse/csr.py +22 -0
scipy/sparse/data.py +18 -0
scipy/sparse/dia.py +22 -0
scipy/sparse/dok.py +22 -0
scipy/sparse/extract.py +23 -0
scipy/sparse/lil.py +22 -0
scipy/sparse/linalg/__init__.py +148 -0
scipy/sparse/linalg/_dsolve/__init__.py +71 -0
scipy/sparse/linalg/_dsolve/_add_newdocs.py +147 -0
scipy/sparse/linalg/_dsolve/_superlu.cp311-win_arm64.lib +0 -0
scipy/sparse/linalg/_dsolve/_superlu.cp311-win_arm64.pyd +0 -0
scipy/sparse/linalg/_dsolve/linsolve.py +882 -0
scipy/sparse/linalg/_dsolve/tests/__init__.py +0 -0
scipy/sparse/linalg/_dsolve/tests/test_linsolve.py +928 -0
scipy/sparse/linalg/_eigen/__init__.py +22 -0
scipy/sparse/linalg/_eigen/_svds.py +540 -0
scipy/sparse/linalg/_eigen/_svds_doc.py +382 -0
scipy/sparse/linalg/_eigen/arpack/COPYING +45 -0
scipy/sparse/linalg/_eigen/arpack/__init__.py +20 -0
scipy/sparse/linalg/_eigen/arpack/_arpack.cp311-win_arm64.lib +0 -0
scipy/sparse/linalg/_eigen/arpack/_arpack.cp311-win_arm64.pyd +0 -0
scipy/sparse/linalg/_eigen/arpack/arpack.py +1706 -0
scipy/sparse/linalg/_eigen/arpack/tests/__init__.py +0 -0
scipy/sparse/linalg/_eigen/arpack/tests/test_arpack.py +717 -0
scipy/sparse/linalg/_eigen/lobpcg/__init__.py +16 -0
scipy/sparse/linalg/_eigen/lobpcg/lobpcg.py +1110 -0
scipy/sparse/linalg/_eigen/lobpcg/tests/__init__.py +0 -0
scipy/sparse/linalg/_eigen/lobpcg/tests/test_lobpcg.py +725 -0
scipy/sparse/linalg/_eigen/tests/__init__.py +0 -0
scipy/sparse/linalg/_eigen/tests/test_svds.py +886 -0
scipy/sparse/linalg/_expm_multiply.py +816 -0
scipy/sparse/linalg/_interface.py +920 -0
scipy/sparse/linalg/_isolve/__init__.py +20 -0
scipy/sparse/linalg/_isolve/_gcrotmk.py +503 -0
scipy/sparse/linalg/_isolve/iterative.py +1051 -0
scipy/sparse/linalg/_isolve/lgmres.py +230 -0
scipy/sparse/linalg/_isolve/lsmr.py +486 -0
scipy/sparse/linalg/_isolve/lsqr.py +589 -0
scipy/sparse/linalg/_isolve/minres.py +372 -0
scipy/sparse/linalg/_isolve/tests/__init__.py +0 -0
scipy/sparse/linalg/_isolve/tests/test_gcrotmk.py +183 -0
scipy/sparse/linalg/_isolve/tests/test_iterative.py +809 -0
scipy/sparse/linalg/_isolve/tests/test_lgmres.py +225 -0
scipy/sparse/linalg/_isolve/tests/test_lsmr.py +185 -0
scipy/sparse/linalg/_isolve/tests/test_lsqr.py +120 -0
scipy/sparse/linalg/_isolve/tests/test_minres.py +97 -0
scipy/sparse/linalg/_isolve/tests/test_utils.py +9 -0
scipy/sparse/linalg/_isolve/tfqmr.py +179 -0
scipy/sparse/linalg/_isolve/utils.py +121 -0
scipy/sparse/linalg/_matfuncs.py +940 -0
scipy/sparse/linalg/_norm.py +195 -0
scipy/sparse/linalg/_onenormest.py +467 -0
scipy/sparse/linalg/_propack/_cpropack.cp311-win_arm64.lib +0 -0
scipy/sparse/linalg/_propack/_cpropack.cp311-win_arm64.pyd +0 -0
scipy/sparse/linalg/_propack/_dpropack.cp311-win_arm64.lib +0 -0
scipy/sparse/linalg/_propack/_dpropack.cp311-win_arm64.pyd +0 -0
scipy/sparse/linalg/_propack/_spropack.cp311-win_arm64.lib +0 -0
scipy/sparse/linalg/_propack/_spropack.cp311-win_arm64.pyd +0 -0
scipy/sparse/linalg/_propack/_zpropack.cp311-win_arm64.lib +0 -0
scipy/sparse/linalg/_propack/_zpropack.cp311-win_arm64.pyd +0 -0
scipy/sparse/linalg/_special_sparse_arrays.py +949 -0
scipy/sparse/linalg/_svdp.py +309 -0
scipy/sparse/linalg/dsolve.py +22 -0
scipy/sparse/linalg/eigen.py +21 -0
scipy/sparse/linalg/interface.py +20 -0
scipy/sparse/linalg/isolve.py +22 -0
scipy/sparse/linalg/matfuncs.py +18 -0
scipy/sparse/linalg/tests/__init__.py +0 -0
scipy/sparse/linalg/tests/propack_test_data.npz +0 -0
scipy/sparse/linalg/tests/test_expm_multiply.py +367 -0
scipy/sparse/linalg/tests/test_interface.py +561 -0
scipy/sparse/linalg/tests/test_matfuncs.py +592 -0
scipy/sparse/linalg/tests/test_norm.py +154 -0
scipy/sparse/linalg/tests/test_onenormest.py +252 -0
scipy/sparse/linalg/tests/test_propack.py +165 -0
scipy/sparse/linalg/tests/test_pydata_sparse.py +272 -0
scipy/sparse/linalg/tests/test_special_sparse_arrays.py +337 -0
scipy/sparse/sparsetools.py +17 -0
scipy/sparse/spfuncs.py +17 -0
scipy/sparse/sputils.py +17 -0
scipy/sparse/tests/__init__.py +0 -0
scipy/sparse/tests/data/csc_py2.npz +0 -0
scipy/sparse/tests/data/csc_py3.npz +0 -0
scipy/sparse/tests/test_arithmetic1d.py +341 -0
scipy/sparse/tests/test_array_api.py +561 -0
scipy/sparse/tests/test_base.py +5870 -0
scipy/sparse/tests/test_common1d.py +447 -0
scipy/sparse/tests/test_construct.py +872 -0
scipy/sparse/tests/test_coo.py +1119 -0
scipy/sparse/tests/test_csc.py +98 -0
scipy/sparse/tests/test_csr.py +214 -0
scipy/sparse/tests/test_dok.py +209 -0
scipy/sparse/tests/test_extract.py +51 -0
scipy/sparse/tests/test_indexing1d.py +603 -0
scipy/sparse/tests/test_matrix_io.py +109 -0
scipy/sparse/tests/test_minmax1d.py +128 -0
scipy/sparse/tests/test_sparsetools.py +344 -0
scipy/sparse/tests/test_spfuncs.py +97 -0
scipy/sparse/tests/test_sputils.py +424 -0
scipy/spatial/__init__.py +129 -0
scipy/spatial/_ckdtree.cp311-win_arm64.lib +0 -0
scipy/spatial/_ckdtree.cp311-win_arm64.pyd +0 -0
scipy/spatial/_distance_pybind.cp311-win_arm64.lib +0 -0
scipy/spatial/_distance_pybind.cp311-win_arm64.pyd +0 -0
scipy/spatial/_distance_wrap.cp311-win_arm64.lib +0 -0
scipy/spatial/_distance_wrap.cp311-win_arm64.pyd +0 -0
scipy/spatial/_geometric_slerp.py +238 -0
scipy/spatial/_hausdorff.cp311-win_arm64.lib +0 -0
scipy/spatial/_hausdorff.cp311-win_arm64.pyd +0 -0
scipy/spatial/_kdtree.py +920 -0
scipy/spatial/_plotutils.py +274 -0
scipy/spatial/_procrustes.py +132 -0
scipy/spatial/_qhull.cp311-win_arm64.lib +0 -0
scipy/spatial/_qhull.cp311-win_arm64.pyd +0 -0
scipy/spatial/_qhull.pyi +213 -0
scipy/spatial/_spherical_voronoi.py +341 -0
scipy/spatial/_voronoi.cp311-win_arm64.lib +0 -0
scipy/spatial/_voronoi.cp311-win_arm64.pyd +0 -0
scipy/spatial/_voronoi.pyi +4 -0
scipy/spatial/ckdtree.py +18 -0
scipy/spatial/distance.py +3147 -0
scipy/spatial/distance.pyi +210 -0
scipy/spatial/kdtree.py +25 -0
scipy/spatial/qhull.py +25 -0
scipy/spatial/qhull_src/COPYING_QHULL.txt +39 -0
scipy/spatial/tests/__init__.py +0 -0
scipy/spatial/tests/data/cdist-X1.txt +10 -0
scipy/spatial/tests/data/cdist-X2.txt +20 -0
scipy/spatial/tests/data/degenerate_pointset.npz +0 -0
scipy/spatial/tests/data/iris.txt +150 -0
scipy/spatial/tests/data/pdist-boolean-inp.txt +20 -0
scipy/spatial/tests/data/pdist-chebyshev-ml-iris.txt +1 -0
scipy/spatial/tests/data/pdist-chebyshev-ml.txt +1 -0
scipy/spatial/tests/data/pdist-cityblock-ml-iris.txt +1 -0
scipy/spatial/tests/data/pdist-cityblock-ml.txt +1 -0
scipy/spatial/tests/data/pdist-correlation-ml-iris.txt +1 -0
scipy/spatial/tests/data/pdist-correlation-ml.txt +1 -0
scipy/spatial/tests/data/pdist-cosine-ml-iris.txt +1 -0
scipy/spatial/tests/data/pdist-cosine-ml.txt +1 -0
scipy/spatial/tests/data/pdist-double-inp.txt +20 -0
scipy/spatial/tests/data/pdist-euclidean-ml-iris.txt +1 -0
scipy/spatial/tests/data/pdist-euclidean-ml.txt +1 -0
scipy/spatial/tests/data/pdist-hamming-ml.txt +1 -0
scipy/spatial/tests/data/pdist-jaccard-ml.txt +1 -0
scipy/spatial/tests/data/pdist-jensenshannon-ml-iris.txt +1 -0
scipy/spatial/tests/data/pdist-jensenshannon-ml.txt +1 -0
scipy/spatial/tests/data/pdist-minkowski-3.2-ml-iris.txt +1 -0
scipy/spatial/tests/data/pdist-minkowski-3.2-ml.txt +1 -0
scipy/spatial/tests/data/pdist-minkowski-5.8-ml-iris.txt +1 -0
scipy/spatial/tests/data/pdist-seuclidean-ml-iris.txt +1 -0
scipy/spatial/tests/data/pdist-seuclidean-ml.txt +1 -0
scipy/spatial/tests/data/pdist-spearman-ml.txt +1 -0
scipy/spatial/tests/data/random-bool-data.txt +100 -0
scipy/spatial/tests/data/random-double-data.txt +100 -0
scipy/spatial/tests/data/random-int-data.txt +100 -0
scipy/spatial/tests/data/random-uint-data.txt +100 -0
scipy/spatial/tests/data/selfdual-4d-polytope.txt +27 -0
scipy/spatial/tests/test__plotutils.py +91 -0
scipy/spatial/tests/test__procrustes.py +116 -0
scipy/spatial/tests/test_distance.py +2389 -0
scipy/spatial/tests/test_hausdorff.py +199 -0
scipy/spatial/tests/test_kdtree.py +1536 -0
scipy/spatial/tests/test_qhull.py +1313 -0
scipy/spatial/tests/test_slerp.py +417 -0
scipy/spatial/tests/test_spherical_voronoi.py +358 -0
scipy/spatial/transform/__init__.py +31 -0
scipy/spatial/transform/_rigid_transform.cp311-win_arm64.lib +0 -0
scipy/spatial/transform/_rigid_transform.cp311-win_arm64.pyd +0 -0
scipy/spatial/transform/_rotation.cp311-win_arm64.lib +0 -0
scipy/spatial/transform/_rotation.cp311-win_arm64.pyd +0 -0
scipy/spatial/transform/_rotation_groups.py +140 -0
scipy/spatial/transform/_rotation_spline.py +460 -0
scipy/spatial/transform/rotation.py +21 -0
scipy/spatial/transform/tests/__init__.py +0 -0
scipy/spatial/transform/tests/test_rigid_transform.py +1221 -0
scipy/spatial/transform/tests/test_rotation.py +2569 -0
scipy/spatial/transform/tests/test_rotation_groups.py +169 -0
scipy/spatial/transform/tests/test_rotation_spline.py +183 -0
scipy/special/__init__.pxd +1 -0
scipy/special/__init__.py +841 -0
scipy/special/_add_newdocs.py +9961 -0
scipy/special/_basic.py +3576 -0
scipy/special/_comb.cp311-win_arm64.lib +0 -0
scipy/special/_comb.cp311-win_arm64.pyd +0 -0
scipy/special/_ellip_harm.py +214 -0
scipy/special/_ellip_harm_2.cp311-win_arm64.lib +0 -0
scipy/special/_ellip_harm_2.cp311-win_arm64.pyd +0 -0
scipy/special/_gufuncs.cp311-win_arm64.lib +0 -0
scipy/special/_gufuncs.cp311-win_arm64.pyd +0 -0
scipy/special/_input_validation.py +17 -0
scipy/special/_lambertw.py +149 -0
scipy/special/_logsumexp.py +426 -0
scipy/special/_mptestutils.py +453 -0
scipy/special/_multiufuncs.py +610 -0
scipy/special/_orthogonal.py +2592 -0
scipy/special/_orthogonal.pyi +330 -0
scipy/special/_precompute/__init__.py +0 -0
scipy/special/_precompute/cosine_cdf.py +17 -0
scipy/special/_precompute/expn_asy.py +54 -0
scipy/special/_precompute/gammainc_asy.py +116 -0
scipy/special/_precompute/gammainc_data.py +124 -0
scipy/special/_precompute/hyp2f1_data.py +484 -0
scipy/special/_precompute/lambertw.py +68 -0
scipy/special/_precompute/loggamma.py +43 -0
scipy/special/_precompute/struve_convergence.py +131 -0
scipy/special/_precompute/utils.py +38 -0
scipy/special/_precompute/wright_bessel.py +342 -0
scipy/special/_precompute/wright_bessel_data.py +152 -0
scipy/special/_precompute/wrightomega.py +41 -0
scipy/special/_precompute/zetac.py +27 -0
scipy/special/_sf_error.py +15 -0
scipy/special/_specfun.cp311-win_arm64.lib +0 -0
scipy/special/_specfun.cp311-win_arm64.pyd +0 -0
scipy/special/_special_ufuncs.cp311-win_arm64.lib +0 -0
scipy/special/_special_ufuncs.cp311-win_arm64.pyd +0 -0
scipy/special/_spfun_stats.py +106 -0
scipy/special/_spherical_bessel.py +397 -0
scipy/special/_support_alternative_backends.py +295 -0
scipy/special/_test_internal.cp311-win_arm64.lib +0 -0
scipy/special/_test_internal.cp311-win_arm64.pyd +0 -0
scipy/special/_test_internal.pyi +9 -0
scipy/special/_testutils.py +321 -0
scipy/special/_ufuncs.cp311-win_arm64.lib +0 -0
scipy/special/_ufuncs.cp311-win_arm64.pyd +0 -0
scipy/special/_ufuncs.pyi +522 -0
scipy/special/_ufuncs.pyx +13173 -0
scipy/special/_ufuncs_cxx.cp311-win_arm64.lib +0 -0
scipy/special/_ufuncs_cxx.cp311-win_arm64.pyd +0 -0
scipy/special/_ufuncs_cxx.pxd +142 -0
scipy/special/_ufuncs_cxx.pyx +427 -0
scipy/special/_ufuncs_cxx_defs.h +147 -0
scipy/special/_ufuncs_defs.h +57 -0
scipy/special/add_newdocs.py +15 -0
scipy/special/basic.py +87 -0
scipy/special/cython_special.cp311-win_arm64.lib +0 -0
scipy/special/cython_special.cp311-win_arm64.pyd +0 -0
scipy/special/cython_special.pxd +259 -0
scipy/special/cython_special.pyi +3 -0
scipy/special/orthogonal.py +45 -0
scipy/special/sf_error.py +20 -0
scipy/special/specfun.py +24 -0
scipy/special/spfun_stats.py +17 -0
scipy/special/tests/__init__.py +0 -0
scipy/special/tests/_cython_examples/extending.pyx +12 -0
scipy/special/tests/_cython_examples/meson.build +34 -0
scipy/special/tests/data/__init__.py +0 -0
scipy/special/tests/data/boost.npz +0 -0
scipy/special/tests/data/gsl.npz +0 -0
scipy/special/tests/data/local.npz +0 -0
scipy/special/tests/test_basic.py +4815 -0
scipy/special/tests/test_bdtr.py +112 -0
scipy/special/tests/test_boost_ufuncs.py +64 -0
scipy/special/tests/test_boxcox.py +125 -0
scipy/special/tests/test_cdflib.py +712 -0
scipy/special/tests/test_cdft_asymptotic.py +49 -0
scipy/special/tests/test_cephes_intp_cast.py +29 -0
scipy/special/tests/test_cosine_distr.py +83 -0
scipy/special/tests/test_cython_special.py +363 -0
scipy/special/tests/test_data.py +719 -0
scipy/special/tests/test_dd.py +42 -0
scipy/special/tests/test_digamma.py +45 -0
scipy/special/tests/test_ellip_harm.py +278 -0
scipy/special/tests/test_erfinv.py +89 -0
scipy/special/tests/test_exponential_integrals.py +118 -0
scipy/special/tests/test_extending.py +28 -0
scipy/special/tests/test_faddeeva.py +85 -0
scipy/special/tests/test_gamma.py +12 -0
scipy/special/tests/test_gammainc.py +152 -0
scipy/special/tests/test_hyp2f1.py +2566 -0
scipy/special/tests/test_hypergeometric.py +234 -0
scipy/special/tests/test_iv_ratio.py +249 -0
scipy/special/tests/test_kolmogorov.py +491 -0
scipy/special/tests/test_lambertw.py +109 -0
scipy/special/tests/test_legendre.py +1518 -0
scipy/special/tests/test_log1mexp.py +85 -0
scipy/special/tests/test_loggamma.py +70 -0
scipy/special/tests/test_logit.py +162 -0
scipy/special/tests/test_logsumexp.py +469 -0
scipy/special/tests/test_mpmath.py +2293 -0
scipy/special/tests/test_nan_inputs.py +65 -0
scipy/special/tests/test_ndtr.py +77 -0
scipy/special/tests/test_ndtri_exp.py +94 -0
scipy/special/tests/test_orthogonal.py +821 -0
scipy/special/tests/test_orthogonal_eval.py +275 -0
scipy/special/tests/test_owens_t.py +53 -0
scipy/special/tests/test_pcf.py +24 -0
scipy/special/tests/test_pdtr.py +48 -0
scipy/special/tests/test_powm1.py +65 -0
scipy/special/tests/test_precompute_expn_asy.py +24 -0
scipy/special/tests/test_precompute_gammainc.py +108 -0
scipy/special/tests/test_precompute_utils.py +36 -0
scipy/special/tests/test_round.py +18 -0
scipy/special/tests/test_sf_error.py +146 -0
scipy/special/tests/test_sici.py +36 -0
scipy/special/tests/test_specfun.py +48 -0
scipy/special/tests/test_spence.py +32 -0
scipy/special/tests/test_spfun_stats.py +61 -0
scipy/special/tests/test_sph_harm.py +85 -0
scipy/special/tests/test_spherical_bessel.py +400 -0
scipy/special/tests/test_support_alternative_backends.py +248 -0
scipy/special/tests/test_trig.py +72 -0
scipy/special/tests/test_ufunc_signatures.py +46 -0
scipy/special/tests/test_wright_bessel.py +205 -0
scipy/special/tests/test_wrightomega.py +117 -0
scipy/special/tests/test_zeta.py +301 -0
scipy/stats/__init__.py +670 -0
scipy/stats/_ansari_swilk_statistics.cp311-win_arm64.lib +0 -0
scipy/stats/_ansari_swilk_statistics.cp311-win_arm64.pyd +0 -0
scipy/stats/_axis_nan_policy.py +692 -0
scipy/stats/_biasedurn.cp311-win_arm64.lib +0 -0
scipy/stats/_biasedurn.cp311-win_arm64.pyd +0 -0
scipy/stats/_biasedurn.pxd +27 -0
scipy/stats/_binned_statistic.py +795 -0
scipy/stats/_binomtest.py +375 -0
scipy/stats/_bws_test.py +177 -0
scipy/stats/_censored_data.py +459 -0
scipy/stats/_common.py +5 -0
scipy/stats/_constants.py +42 -0
scipy/stats/_continued_fraction.py +387 -0
scipy/stats/_continuous_distns.py +12486 -0
scipy/stats/_correlation.py +210 -0
scipy/stats/_covariance.py +636 -0
scipy/stats/_crosstab.py +204 -0
scipy/stats/_discrete_distns.py +2098 -0
scipy/stats/_distn_infrastructure.py +4201 -0
scipy/stats/_distr_params.py +299 -0
scipy/stats/_distribution_infrastructure.py +5750 -0
scipy/stats/_entropy.py +428 -0
scipy/stats/_finite_differences.py +145 -0
scipy/stats/_fit.py +1351 -0
scipy/stats/_hypotests.py +2060 -0
scipy/stats/_kde.py +732 -0
scipy/stats/_ksstats.py +600 -0
scipy/stats/_levy_stable/__init__.py +1231 -0
scipy/stats/_levy_stable/levyst.cp311-win_arm64.lib +0 -0
scipy/stats/_levy_stable/levyst.cp311-win_arm64.pyd +0 -0
scipy/stats/_mannwhitneyu.py +492 -0
scipy/stats/_mgc.py +550 -0
scipy/stats/_morestats.py +4626 -0
scipy/stats/_mstats_basic.py +3658 -0
scipy/stats/_mstats_extras.py +521 -0
scipy/stats/_multicomp.py +449 -0
scipy/stats/_multivariate.py +7281 -0
scipy/stats/_new_distributions.py +452 -0
scipy/stats/_odds_ratio.py +466 -0
scipy/stats/_page_trend_test.py +486 -0
scipy/stats/_probability_distribution.py +1964 -0
scipy/stats/_qmc.py +2956 -0
scipy/stats/_qmc_cy.cp311-win_arm64.lib +0 -0
scipy/stats/_qmc_cy.cp311-win_arm64.pyd +0 -0
scipy/stats/_qmc_cy.pyi +54 -0
scipy/stats/_qmvnt.py +454 -0
scipy/stats/_qmvnt_cy.cp311-win_arm64.lib +0 -0
scipy/stats/_qmvnt_cy.cp311-win_arm64.pyd +0 -0
scipy/stats/_quantile.py +335 -0
scipy/stats/_rcont/__init__.py +4 -0
scipy/stats/_rcont/rcont.cp311-win_arm64.lib +0 -0
scipy/stats/_rcont/rcont.cp311-win_arm64.pyd +0 -0
scipy/stats/_relative_risk.py +263 -0
scipy/stats/_resampling.py +2352 -0
scipy/stats/_result_classes.py +40 -0
scipy/stats/_sampling.py +1314 -0
scipy/stats/_sensitivity_analysis.py +713 -0
scipy/stats/_sobol.cp311-win_arm64.lib +0 -0
scipy/stats/_sobol.cp311-win_arm64.pyd +0 -0
scipy/stats/_sobol.pyi +54 -0
scipy/stats/_sobol_direction_numbers.npz +0 -0
scipy/stats/_stats.cp311-win_arm64.lib +0 -0
scipy/stats/_stats.cp311-win_arm64.pyd +0 -0
scipy/stats/_stats.pxd +10 -0
scipy/stats/_stats_mstats_common.py +322 -0
scipy/stats/_stats_py.py +11089 -0
scipy/stats/_stats_pythran.cp311-win_arm64.lib +0 -0
scipy/stats/_stats_pythran.cp311-win_arm64.pyd +0 -0
scipy/stats/_survival.py +683 -0
scipy/stats/_tukeylambda_stats.py +199 -0
scipy/stats/_unuran/__init__.py +0 -0
scipy/stats/_unuran/unuran_wrapper.cp311-win_arm64.lib +0 -0
scipy/stats/_unuran/unuran_wrapper.cp311-win_arm64.pyd +0 -0
scipy/stats/_unuran/unuran_wrapper.pyi +179 -0
scipy/stats/_variation.py +126 -0
scipy/stats/_warnings_errors.py +38 -0
scipy/stats/_wilcoxon.py +265 -0
scipy/stats/biasedurn.py +16 -0
scipy/stats/contingency.py +521 -0
scipy/stats/distributions.py +24 -0
scipy/stats/kde.py +18 -0
scipy/stats/morestats.py +27 -0
scipy/stats/mstats.py +140 -0
scipy/stats/mstats_basic.py +42 -0
scipy/stats/mstats_extras.py +25 -0
scipy/stats/mvn.py +17 -0
scipy/stats/qmc.py +236 -0
scipy/stats/sampling.py +73 -0
scipy/stats/stats.py +41 -0
scipy/stats/tests/__init__.py +0 -0
scipy/stats/tests/common_tests.py +356 -0
scipy/stats/tests/data/_mvt.py +171 -0
scipy/stats/tests/data/fisher_exact_results_from_r.py +607 -0
scipy/stats/tests/data/jf_skew_t_gamlss_pdf_data.npy +0 -0
scipy/stats/tests/data/levy_stable/stable-Z1-cdf-sample-data.npy +0 -0
scipy/stats/tests/data/levy_stable/stable-Z1-pdf-sample-data.npy +0 -0
scipy/stats/tests/data/levy_stable/stable-loc-scale-sample-data.npy +0 -0
scipy/stats/tests/data/nist_anova/AtmWtAg.dat +108 -0
scipy/stats/tests/data/nist_anova/SiRstv.dat +85 -0
scipy/stats/tests/data/nist_anova/SmLs01.dat +249 -0
scipy/stats/tests/data/nist_anova/SmLs02.dat +1869 -0
scipy/stats/tests/data/nist_anova/SmLs03.dat +18069 -0
scipy/stats/tests/data/nist_anova/SmLs04.dat +249 -0
scipy/stats/tests/data/nist_anova/SmLs05.dat +1869 -0
scipy/stats/tests/data/nist_anova/SmLs06.dat +18069 -0
scipy/stats/tests/data/nist_anova/SmLs07.dat +249 -0
scipy/stats/tests/data/nist_anova/SmLs08.dat +1869 -0
scipy/stats/tests/data/nist_anova/SmLs09.dat +18069 -0
scipy/stats/tests/data/nist_linregress/Norris.dat +97 -0
scipy/stats/tests/data/rel_breitwigner_pdf_sample_data_ROOT.npy +0 -0
scipy/stats/tests/data/studentized_range_mpmath_ref.json +1499 -0
scipy/stats/tests/test_axis_nan_policy.py +1388 -0
scipy/stats/tests/test_binned_statistic.py +568 -0
scipy/stats/tests/test_censored_data.py +152 -0
scipy/stats/tests/test_contingency.py +294 -0
scipy/stats/tests/test_continued_fraction.py +173 -0
scipy/stats/tests/test_continuous.py +2198 -0
scipy/stats/tests/test_continuous_basic.py +1053 -0
scipy/stats/tests/test_continuous_fit_censored.py +683 -0
scipy/stats/tests/test_correlation.py +80 -0
scipy/stats/tests/test_crosstab.py +115 -0
scipy/stats/tests/test_discrete_basic.py +580 -0
scipy/stats/tests/test_discrete_distns.py +700 -0
scipy/stats/tests/test_distributions.py +10413 -0
scipy/stats/tests/test_entropy.py +322 -0
scipy/stats/tests/test_fast_gen_inversion.py +435 -0
scipy/stats/tests/test_fit.py +1090 -0
scipy/stats/tests/test_hypotests.py +1991 -0
scipy/stats/tests/test_kdeoth.py +676 -0
scipy/stats/tests/test_marray.py +289 -0
scipy/stats/tests/test_mgc.py +217 -0
scipy/stats/tests/test_morestats.py +3259 -0
scipy/stats/tests/test_mstats_basic.py +2071 -0
scipy/stats/tests/test_mstats_extras.py +172 -0
scipy/stats/tests/test_multicomp.py +405 -0
scipy/stats/tests/test_multivariate.py +4381 -0
scipy/stats/tests/test_odds_ratio.py +148 -0
scipy/stats/tests/test_qmc.py +1492 -0
scipy/stats/tests/test_quantile.py +199 -0
scipy/stats/tests/test_rank.py +345 -0
scipy/stats/tests/test_relative_risk.py +95 -0
scipy/stats/tests/test_resampling.py +2000 -0
scipy/stats/tests/test_sampling.py +1450 -0
scipy/stats/tests/test_sensitivity_analysis.py +310 -0
scipy/stats/tests/test_stats.py +9707 -0
scipy/stats/tests/test_survival.py +466 -0
scipy/stats/tests/test_tukeylambda_stats.py +85 -0
scipy/stats/tests/test_variation.py +216 -0
scipy/version.py +12 -0
scipy-1.16.2.dist-info/DELVEWHEEL +2 -0
scipy-1.16.2.dist-info/LICENSE.txt +912 -0
scipy-1.16.2.dist-info/METADATA +1061 -0
scipy-1.16.2.dist-info/RECORD +1530 -0
scipy-1.16.2.dist-info/WHEEL +4 -0
scipy.libs/msvcp140-5f1c5dd31916990d94181e07bc3afb32.dll +0 -0
scipy.libs/scipy_openblas-f3ac85b1f412f7e86514c923dc4058d1.dll +0 -0

scipy/special/tests/test_legendre.py ADDED Viewed

@@ -0,0 +1,1518 @@
+import math
+import numpy as np
+import pytest
+from numpy.testing import (assert_equal, assert_almost_equal, assert_array_almost_equal,
+    assert_allclose, suppress_warnings)
+from scipy import special
+from scipy.special import (legendre_p, legendre_p_all, assoc_legendre_p,
+    assoc_legendre_p_all, sph_legendre_p, sph_legendre_p_all)
+# The functions lpn, lpmn, clpmn, appearing below are
+# deprecated in favor of legendre_p_all, assoc_legendre_p_all, and
+# assoc_legendre_p_all (assoc_legendre_p_all covers lpmn and clpmn)
+# respectively. The deprecated functions listed above are implemented as
+# shims around their respective replacements. The replacements are tested
+# separately, but tests for the deprecated functions remain to verify the
+# correctness of the shims.
+# Base polynomials come from Abrahmowitz and Stegan
+class TestLegendre:
+    def test_legendre(self):
+        leg0 = special.legendre(0)
+        leg1 = special.legendre(1)
+        leg2 = special.legendre(2)
+        leg3 = special.legendre(3)
+        leg4 = special.legendre(4)
+        leg5 = special.legendre(5)
+        assert_equal(leg0.c, [1])
+        assert_equal(leg1.c, [1,0])
+        assert_almost_equal(leg2.c, np.array([3,0,-1])/2.0, decimal=13)
+        assert_almost_equal(leg3.c, np.array([5,0,-3,0])/2.0)
+        assert_almost_equal(leg4.c, np.array([35,0,-30,0,3])/8.0)
+        assert_almost_equal(leg5.c, np.array([63,0,-70,0,15,0])/8.0)
+    @pytest.mark.parametrize('n', [1, 2, 3, 4, 5])
+    @pytest.mark.parametrize('zr', [0.5241717, 12.80232, -9.699001,
+                                    0.5122437, 0.1714377])
+    @pytest.mark.parametrize('zi', [9.766818, 0.2999083, 8.24726, -22.84843,
+                                    -0.8792666])
+    def test_lpn_against_clpmn(self, n, zr, zi):
+        with suppress_warnings() as sup:
+            sup.filter(category=DeprecationWarning)
+            reslpn = special.lpn(n, zr + zi*1j)
+            resclpmn = special.clpmn(0, n, zr+zi*1j)
+        assert_allclose(reslpn[0], resclpmn[0][0])
+        assert_allclose(reslpn[1], resclpmn[1][0])
+class TestLegendreP:
+    @pytest.mark.parametrize("shape", [(10,), (4, 9), (3, 5, 7)])
+    def test_ode(self, shape):
+        rng = np.random.default_rng(1234)
+        n = rng.integers(0, 100, shape)
+        x = rng.uniform(-1, 1, shape)
+        p, p_jac, p_hess = legendre_p(n, x, diff_n=2)
+        assert p.shape == shape
+        assert p_jac.shape == p.shape
+        assert p_hess.shape == p_jac.shape
+        err = (1 - x * x) * p_hess - 2 * x * p_jac + n * (n + 1) * p
+        np.testing.assert_allclose(err, 0, atol=1e-10)
+    @pytest.mark.parametrize("n_max", [1, 2, 4, 8, 16, 32])
+    @pytest.mark.parametrize("x_shape", [(10,), (4, 9), (3, 5, 7)])
+    def test_all_ode(self, n_max, x_shape):
+        rng = np.random.default_rng(1234)
+        x = rng.uniform(-1, 1, x_shape)
+        p, p_jac, p_hess = legendre_p_all(n_max, x, diff_n=2)
+        n = np.arange(n_max + 1)
+        n = np.expand_dims(n, axis = tuple(range(1, x.ndim + 1)))
+        assert p.shape == (len(n),) + x.shape
+        assert p_jac.shape == p.shape
+        assert p_hess.shape == p_jac.shape
+        err = (1 - x * x) * p_hess - 2 * x * p_jac + n * (n + 1) * p
+        np.testing.assert_allclose(err, 0, atol=1e-10)
+    def test_legacy(self):
+        with suppress_warnings() as sup:
+            sup.filter(category=DeprecationWarning)
+            p, pd = special.lpn(2, 0.5)
+        assert_array_almost_equal(p, [1.00000, 0.50000, -0.12500], 4)
+        assert_array_almost_equal(pd, [0.00000, 1.00000, 1.50000], 4)
+class TestAssocLegendreP:
+    @pytest.mark.parametrize("shape", [(10,), (4, 9), (3, 5, 7, 10)])
+    @pytest.mark.parametrize("m_max", [5, 4])
+    @pytest.mark.parametrize("n_max", [7, 10])
+    def test_lpmn(self, shape, n_max, m_max):
+        rng = np.random.default_rng(1234)
+        x = rng.uniform(-0.99, 0.99, shape)
+        p_all, p_all_jac, p_all_hess = \
+            assoc_legendre_p_all(n_max, m_max, x, diff_n=2)
+        n = np.arange(n_max + 1)
+        n = np.expand_dims(n, axis = tuple(range(1, x.ndim + 2)))
+        m = np.concatenate([np.arange(m_max + 1), np.arange(-m_max, 0)])
+        m = np.expand_dims(m, axis = (0,) + tuple(range(2, x.ndim + 2)))
+        x = np.expand_dims(x, axis = (0, 1))
+        p, p_jac, p_hess = assoc_legendre_p(n, m, x, diff_n=2)
+        np.testing.assert_allclose(p, p_all)
+        np.testing.assert_allclose(p_jac, p_all_jac)
+        np.testing.assert_allclose(p_hess, p_all_hess)
+    @pytest.mark.parametrize("shape", [(10,), (4, 9), (3, 5, 7, 10)])
+    @pytest.mark.parametrize("norm", [True, False])
+    def test_ode(self, shape, norm):
+        rng = np.random.default_rng(1234)
+        n = rng.integers(0, 10, shape)
+        m = rng.integers(-10, 10, shape)
+        x = rng.uniform(-1, 1, shape)
+        p, p_jac, p_hess = assoc_legendre_p(n, m, x, norm=norm, diff_n=2)
+        assert p.shape == shape
+        assert p_jac.shape == p.shape
+        assert p_hess.shape == p_jac.shape
+        np.testing.assert_allclose((1 - x * x) * p_hess,
+            2 * x * p_jac - (n * (n + 1) - m * m / (1 - x * x)) * p,
+            rtol=1e-05, atol=1e-08)
+    @pytest.mark.parametrize("shape", [(10,), (4, 9), (3, 5, 7)])
+    def test_all(self, shape):
+        rng = np.random.default_rng(1234)
+        n_max = 20
+        m_max = 20
+        x = rng.uniform(-0.99, 0.99, shape)
+        p, p_jac, p_hess = assoc_legendre_p_all(n_max, m_max, x, diff_n=2)
+        m = np.concatenate([np.arange(m_max + 1), np.arange(-m_max, 0)])
+        n = np.arange(n_max + 1)
+        n = np.expand_dims(n, axis = tuple(range(1, x.ndim + 2)))
+        m = np.expand_dims(m, axis = (0,) + tuple(range(2, x.ndim + 2)))
+        np.testing.assert_allclose((1 - x * x) * p_hess,
+            2 * x * p_jac - (n * (n + 1) - m * m / (1 - x * x)) * p,
+            rtol=1e-05, atol=1e-08)
+    @pytest.mark.parametrize("shape", [(10,), (4, 9), (3, 5, 7)])
+    @pytest.mark.parametrize("norm", [True, False])
+    def test_specific(self, shape, norm):
+        rng = np.random.default_rng(1234)
+        x = rng.uniform(-0.99, 0.99, shape)
+        p, p_jac = assoc_legendre_p_all(4, 4, x, norm=norm, diff_n=1)
+        np.testing.assert_allclose(p[0, 0],
+            assoc_legendre_p_0_0(x, norm=norm))
+        np.testing.assert_allclose(p[0, 1], 0)
+        np.testing.assert_allclose(p[0, 2], 0)
+        np.testing.assert_allclose(p[0, 3], 0)
+        np.testing.assert_allclose(p[0, 4], 0)
+        np.testing.assert_allclose(p[0, -3], 0)
+        np.testing.assert_allclose(p[0, -2], 0)
+        np.testing.assert_allclose(p[0, -1], 0)
+        np.testing.assert_allclose(p[1, 0],
+            assoc_legendre_p_1_0(x, norm=norm))
+        np.testing.assert_allclose(p[1, 1],
+            assoc_legendre_p_1_1(x, norm=norm))
+        np.testing.assert_allclose(p[1, 2], 0)
+        np.testing.assert_allclose(p[1, 3], 0)
+        np.testing.assert_allclose(p[1, 4], 0)
+        np.testing.assert_allclose(p[1, -4], 0)
+        np.testing.assert_allclose(p[1, -3], 0)
+        np.testing.assert_allclose(p[1, -2], 0)
+        np.testing.assert_allclose(p[1, -1],
+            assoc_legendre_p_1_m1(x, norm=norm))
+        np.testing.assert_allclose(p[2, 0],
+            assoc_legendre_p_2_0(x, norm=norm))
+        np.testing.assert_allclose(p[2, 1],
+            assoc_legendre_p_2_1(x, norm=norm))
+        np.testing.assert_allclose(p[2, 2],
+            assoc_legendre_p_2_2(x, norm=norm))
+        np.testing.assert_allclose(p[2, 3], 0)
+        np.testing.assert_allclose(p[2, 4], 0)
+        np.testing.assert_allclose(p[2, -4], 0)
+        np.testing.assert_allclose(p[2, -3], 0)
+        np.testing.assert_allclose(p[2, -2],
+            assoc_legendre_p_2_m2(x, norm=norm))
+        np.testing.assert_allclose(p[2, -1],
+            assoc_legendre_p_2_m1(x, norm=norm))
+        np.testing.assert_allclose(p[3, 0],
+            assoc_legendre_p_3_0(x, norm=norm))
+        np.testing.assert_allclose(p[3, 1],
+            assoc_legendre_p_3_1(x, norm=norm))
+        np.testing.assert_allclose(p[3, 2],
+            assoc_legendre_p_3_2(x, norm=norm))
+        np.testing.assert_allclose(p[3, 3],
+            assoc_legendre_p_3_3(x, norm=norm))
+        np.testing.assert_allclose(p[3, 4], 0)
+        np.testing.assert_allclose(p[3, -4], 0)
+        np.testing.assert_allclose(p[3, -3],
+            assoc_legendre_p_3_m3(x, norm=norm))
+        np.testing.assert_allclose(p[3, -2],
+            assoc_legendre_p_3_m2(x, norm=norm))
+        np.testing.assert_allclose(p[3, -1],
+            assoc_legendre_p_3_m1(x, norm=norm))
+        np.testing.assert_allclose(p[4, 0],
+            assoc_legendre_p_4_0(x, norm=norm))
+        np.testing.assert_allclose(p[4, 1],
+            assoc_legendre_p_4_1(x, norm=norm))
+        np.testing.assert_allclose(p[4, 2],
+            assoc_legendre_p_4_2(x, norm=norm))
+        np.testing.assert_allclose(p[4, 3],
+            assoc_legendre_p_4_3(x, norm=norm))
+        np.testing.assert_allclose(p[4, 4],
+            assoc_legendre_p_4_4(x, norm=norm))
+        np.testing.assert_allclose(p[4, -4],
+            assoc_legendre_p_4_m4(x, norm=norm))
+        np.testing.assert_allclose(p[4, -3],
+            assoc_legendre_p_4_m3(x, norm=norm))
+        np.testing.assert_allclose(p[4, -2],
+            assoc_legendre_p_4_m2(x, norm=norm))
+        np.testing.assert_allclose(p[4, -1],
+            assoc_legendre_p_4_m1(x, norm=norm))
+        np.testing.assert_allclose(p_jac[0, 0],
+            assoc_legendre_p_0_0_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[0, 1], 0)
+        np.testing.assert_allclose(p_jac[0, 2], 0)
+        np.testing.assert_allclose(p_jac[0, 3], 0)
+        np.testing.assert_allclose(p_jac[0, 4], 0)
+        np.testing.assert_allclose(p_jac[0, -4], 0)
+        np.testing.assert_allclose(p_jac[0, -3], 0)
+        np.testing.assert_allclose(p_jac[0, -2], 0)
+        np.testing.assert_allclose(p_jac[0, -1], 0)
+        np.testing.assert_allclose(p_jac[1, 0],
+            assoc_legendre_p_1_0_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[1, 1],
+            assoc_legendre_p_1_1_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[1, 2], 0)
+        np.testing.assert_allclose(p_jac[1, 3], 0)
+        np.testing.assert_allclose(p_jac[1, 4], 0)
+        np.testing.assert_allclose(p_jac[1, -4], 0)
+        np.testing.assert_allclose(p_jac[1, -3], 0)
+        np.testing.assert_allclose(p_jac[1, -2], 0)
+        np.testing.assert_allclose(p_jac[1, -1],
+            assoc_legendre_p_1_m1_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[2, 0],
+            assoc_legendre_p_2_0_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[2, 1],
+            assoc_legendre_p_2_1_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[2, 2],
+            assoc_legendre_p_2_2_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[2, 3], 0)
+        np.testing.assert_allclose(p_jac[2, 4], 0)
+        np.testing.assert_allclose(p_jac[2, -4], 0)
+        np.testing.assert_allclose(p_jac[2, -3], 0)
+        np.testing.assert_allclose(p_jac[2, -2],
+            assoc_legendre_p_2_m2_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[2, -1],
+            assoc_legendre_p_2_m1_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[3, 0],
+            assoc_legendre_p_3_0_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[3, 1],
+            assoc_legendre_p_3_1_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[3, 2],
+            assoc_legendre_p_3_2_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[3, 3],
+            assoc_legendre_p_3_3_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[3, 4], 0)
+        np.testing.assert_allclose(p_jac[3, -4], 0)
+        np.testing.assert_allclose(p_jac[3, -3],
+            assoc_legendre_p_3_m3_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[3, -2],
+            assoc_legendre_p_3_m2_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[3, -1],
+            assoc_legendre_p_3_m1_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[4, 0],
+            assoc_legendre_p_4_0_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[4, 1],
+            assoc_legendre_p_4_1_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[4, 2],
+            assoc_legendre_p_4_2_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[4, 3],
+            assoc_legendre_p_4_3_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[4, 4],
+            assoc_legendre_p_4_4_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[4, -4],
+            assoc_legendre_p_4_m4_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[4, -3],
+            assoc_legendre_p_4_m3_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[4, -2],
+            assoc_legendre_p_4_m2_jac(x, norm=norm))
+        np.testing.assert_allclose(p_jac[4, -1],
+            assoc_legendre_p_4_m1_jac(x, norm=norm))
+    @pytest.mark.parametrize("m_max", [7])
+    @pytest.mark.parametrize("n_max", [10])
+    @pytest.mark.parametrize("x", [1, -1])
+    def test_all_limits(self, m_max, n_max, x):
+        p, p_jac = assoc_legendre_p_all(n_max, m_max, x, diff_n=1)
+        n = np.arange(n_max + 1)
+        np.testing.assert_allclose(p_jac[:, 0],
+            pow(x, n + 1) * n * (n + 1) / 2)
+        np.testing.assert_allclose(p_jac[:, 1],
+            np.where(n >= 1, pow(x, n) * np.inf, 0))
+        np.testing.assert_allclose(p_jac[:, 2],
+            np.where(n >= 2, -pow(x, n + 1) * (n + 2) * (n + 1) * n * (n - 1) / 4, 0))
+        np.testing.assert_allclose(p_jac[:, -2],
+            np.where(n >= 2, -pow(x, n + 1) / 4, 0))
+        np.testing.assert_allclose(p_jac[:, -1],
+            np.where(n >= 1, -pow(x, n) * np.inf, 0))
+        for m in range(3, m_max + 1):
+            np.testing.assert_allclose(p_jac[:, m], 0)
+            np.testing.assert_allclose(p_jac[:, -m], 0)
+    @pytest.mark.parametrize("m_max", [3, 5, 10])
+    @pytest.mark.parametrize("n_max", [10])
+    def test_legacy(self, m_max, n_max):
+        x = 0.5
+        p, p_jac = assoc_legendre_p_all(n_max, m_max, x, diff_n=1)
+        with suppress_warnings() as sup:
+            sup.filter(category=DeprecationWarning)
+            p_legacy, p_jac_legacy = special.lpmn(m_max, n_max, x)
+            for m in range(m_max + 1):
+                np.testing.assert_allclose(p_legacy[m], p[:, m])
+            p_legacy, p_jac_legacy = special.lpmn(-m_max, n_max, x)
+            for m in range(m_max + 1):
+                np.testing.assert_allclose(p_legacy[m], p[:, -m])
+class TestMultiAssocLegendreP:
+    @pytest.mark.parametrize("shape", [(1000,), (4, 9), (3, 5, 7)])
+    @pytest.mark.parametrize("branch_cut", [2, 3])
+    @pytest.mark.parametrize("z_min, z_max", [(-10 - 10j, 10 + 10j),
+        (-1, 1), (-10j, 10j)])
+    @pytest.mark.parametrize("norm", [True, False])
+    def test_specific(self, shape, branch_cut, z_min, z_max, norm):
+        rng = np.random.default_rng(1234)
+        z = rng.uniform(z_min.real, z_max.real, shape) + \
+            1j * rng.uniform(z_min.imag, z_max.imag, shape)
+        p, p_jac = assoc_legendre_p_all(4, 4,
+            z, branch_cut=branch_cut, norm=norm, diff_n=1)
+        np.testing.assert_allclose(p[0, 0],
+            assoc_legendre_p_0_0(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[0, 1], 0)
+        np.testing.assert_allclose(p[0, 2], 0)
+        np.testing.assert_allclose(p[0, 3], 0)
+        np.testing.assert_allclose(p[0, 4], 0)
+        np.testing.assert_allclose(p[0, -4], 0)
+        np.testing.assert_allclose(p[0, -3], 0)
+        np.testing.assert_allclose(p[0, -2], 0)
+        np.testing.assert_allclose(p[0, -1], 0)
+        np.testing.assert_allclose(p[1, 0],
+            assoc_legendre_p_1_0(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[1, 1],
+            assoc_legendre_p_1_1(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[1, 2], 0)
+        np.testing.assert_allclose(p[1, 3], 0)
+        np.testing.assert_allclose(p[1, 4], 0)
+        np.testing.assert_allclose(p[1, -4], 0)
+        np.testing.assert_allclose(p[1, -3], 0)
+        np.testing.assert_allclose(p[1, -2], 0)
+        np.testing.assert_allclose(p[1, -1],
+            assoc_legendre_p_1_m1(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[2, 0],
+            assoc_legendre_p_2_0(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[2, 1],
+            assoc_legendre_p_2_1(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[2, 2],
+            assoc_legendre_p_2_2(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[2, 3], 0)
+        np.testing.assert_allclose(p[2, 4], 0)
+        np.testing.assert_allclose(p[2, -4], 0)
+        np.testing.assert_allclose(p[2, -3], 0)
+        np.testing.assert_allclose(p[2, -2],
+            assoc_legendre_p_2_m2(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[2, -1],
+            assoc_legendre_p_2_m1(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[3, 0],
+            assoc_legendre_p_3_0(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[3, 1],
+            assoc_legendre_p_3_1(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[3, 2],
+            assoc_legendre_p_3_2(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[3, 3],
+            assoc_legendre_p_3_3(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[3, 4], 0)
+        np.testing.assert_allclose(p[3, -4], 0)
+        np.testing.assert_allclose(p[3, -3],
+            assoc_legendre_p_3_m3(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[3, -2],
+            assoc_legendre_p_3_m2(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[3, -1],
+            assoc_legendre_p_3_m1(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[4, 0],
+            assoc_legendre_p_4_0(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[4, 1],
+            assoc_legendre_p_4_1(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[4, 2],
+            assoc_legendre_p_4_2(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[4, 3],
+            assoc_legendre_p_4_3(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[4, 4],
+            assoc_legendre_p_4_4(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[4, -4],
+            assoc_legendre_p_4_m4(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[4, -3],
+            assoc_legendre_p_4_m3(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[4, -2],
+            assoc_legendre_p_4_m2(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p[4, -1],
+            assoc_legendre_p_4_m1(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[0, 0],
+            assoc_legendre_p_0_0_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[0, 1], 0)
+        np.testing.assert_allclose(p_jac[0, 2], 0)
+        np.testing.assert_allclose(p_jac[0, 3], 0)
+        np.testing.assert_allclose(p_jac[0, 4], 0)
+        np.testing.assert_allclose(p_jac[0, -4], 0)
+        np.testing.assert_allclose(p_jac[0, -3], 0)
+        np.testing.assert_allclose(p_jac[0, -2], 0)
+        np.testing.assert_allclose(p_jac[0, -1], 0)
+        np.testing.assert_allclose(p_jac[1, 0],
+            assoc_legendre_p_1_0_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[1, 1],
+            assoc_legendre_p_1_1_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[1, 2], 0)
+        np.testing.assert_allclose(p_jac[1, 3], 0)
+        np.testing.assert_allclose(p_jac[1, 4], 0)
+        np.testing.assert_allclose(p_jac[1, -4], 0)
+        np.testing.assert_allclose(p_jac[1, -3], 0)
+        np.testing.assert_allclose(p_jac[1, -2], 0)
+        np.testing.assert_allclose(p_jac[1, -1],
+            assoc_legendre_p_1_m1_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[2, 0],
+            assoc_legendre_p_2_0_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[2, 1],
+            assoc_legendre_p_2_1_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[2, 2],
+            assoc_legendre_p_2_2_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[2, 3], 0)
+        np.testing.assert_allclose(p_jac[2, 4], 0)
+        np.testing.assert_allclose(p_jac[2, -4], 0)
+        np.testing.assert_allclose(p_jac[2, -3], 0)
+        np.testing.assert_allclose(p_jac[2, -2],
+            assoc_legendre_p_2_m2_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[2, -1],
+            assoc_legendre_p_2_m1_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[3, 0],
+            assoc_legendre_p_3_0_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[3, 1],
+            assoc_legendre_p_3_1_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[3, 2],
+            assoc_legendre_p_3_2_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[3, 3],
+            assoc_legendre_p_3_3_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[3, 4], 0)
+        np.testing.assert_allclose(p_jac[3, -4], 0)
+        np.testing.assert_allclose(p_jac[3, -3],
+            assoc_legendre_p_3_m3_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[3, -2],
+            assoc_legendre_p_3_m2_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[3, -1],
+            assoc_legendre_p_3_m1_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[4, 0],
+            assoc_legendre_p_4_0_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[4, 1],
+            assoc_legendre_p_4_1_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[4, 2],
+            assoc_legendre_p_4_2_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[4, 3],
+            assoc_legendre_p_4_3_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[4, 4],
+            assoc_legendre_p_4_4_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[4, -4],
+            assoc_legendre_p_4_m4_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[4, -3],
+            assoc_legendre_p_4_m3_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[4, -2],
+            assoc_legendre_p_4_m2_jac(z, branch_cut=branch_cut, norm=norm))
+        np.testing.assert_allclose(p_jac[4, -1],
+            assoc_legendre_p_4_m1_jac(z, branch_cut=branch_cut, norm=norm))
+class TestSphLegendreP:
+    @pytest.mark.parametrize("shape", [(10,), (4, 9), (3, 5, 7)])
+    def test_specific(self, shape):
+        rng = np.random.default_rng(1234)
+        theta = rng.uniform(-np.pi, np.pi, shape)
+        p, p_jac = sph_legendre_p_all(4, 4, theta, diff_n=1)
+        np.testing.assert_allclose(p[0, 0],
+            sph_legendre_p_0_0(theta))
+        np.testing.assert_allclose(p[0, 1], 0)
+        np.testing.assert_allclose(p[0, 2], 0)
+        np.testing.assert_allclose(p[0, 3], 0)
+        np.testing.assert_allclose(p[0, 4], 0)
+        np.testing.assert_allclose(p[0, -3], 0)
+        np.testing.assert_allclose(p[0, -2], 0)
+        np.testing.assert_allclose(p[0, -1], 0)
+        np.testing.assert_allclose(p[1, 0],
+            sph_legendre_p_1_0(theta))
+        np.testing.assert_allclose(p[1, 1],
+            sph_legendre_p_1_1(theta))
+        np.testing.assert_allclose(p[1, 2], 0)
+        np.testing.assert_allclose(p[1, 3], 0)
+        np.testing.assert_allclose(p[1, 4], 0)
+        np.testing.assert_allclose(p[1, -4], 0)
+        np.testing.assert_allclose(p[1, -3], 0)
+        np.testing.assert_allclose(p[1, -2], 0)
+        np.testing.assert_allclose(p[1, -1],
+            sph_legendre_p_1_m1(theta))
+        np.testing.assert_allclose(p[2, 0],
+            sph_legendre_p_2_0(theta))
+        np.testing.assert_allclose(p[2, 1],
+            sph_legendre_p_2_1(theta))
+        np.testing.assert_allclose(p[2, 2],
+            sph_legendre_p_2_2(theta))
+        np.testing.assert_allclose(p[2, 3], 0)
+        np.testing.assert_allclose(p[2, 4], 0)
+        np.testing.assert_allclose(p[2, -4], 0)
+        np.testing.assert_allclose(p[2, -3], 0)
+        np.testing.assert_allclose(p[2, -2],
+            sph_legendre_p_2_m2(theta))
+        np.testing.assert_allclose(p[2, -1],
+            sph_legendre_p_2_m1(theta))
+        np.testing.assert_allclose(p[3, 0],
+            sph_legendre_p_3_0(theta))
+        np.testing.assert_allclose(p[3, 1],
+            sph_legendre_p_3_1(theta))
+        np.testing.assert_allclose(p[3, 2],
+            sph_legendre_p_3_2(theta))
+        np.testing.assert_allclose(p[3, 3],
+            sph_legendre_p_3_3(theta))
+        np.testing.assert_allclose(p[3, 4], 0)
+        np.testing.assert_allclose(p[3, -4], 0)
+        np.testing.assert_allclose(p[3, -3],
+            sph_legendre_p_3_m3(theta))
+        np.testing.assert_allclose(p[3, -2],
+            sph_legendre_p_3_m2(theta))
+        np.testing.assert_allclose(p[3, -1],
+            sph_legendre_p_3_m1(theta))
+        np.testing.assert_allclose(p[4, 0],
+            sph_legendre_p_4_0(theta))
+        np.testing.assert_allclose(p[4, 1],
+            sph_legendre_p_4_1(theta))
+        np.testing.assert_allclose(p[4, 2],
+            sph_legendre_p_4_2(theta))
+        np.testing.assert_allclose(p[4, 3],
+            sph_legendre_p_4_3(theta))
+        np.testing.assert_allclose(p[4, 4],
+            sph_legendre_p_4_4(theta))
+        np.testing.assert_allclose(p[4, -4],
+            sph_legendre_p_4_m4(theta))
+        np.testing.assert_allclose(p[4, -3],
+            sph_legendre_p_4_m3(theta))
+        np.testing.assert_allclose(p[4, -2],
+            sph_legendre_p_4_m2(theta))
+        np.testing.assert_allclose(p[4, -1],
+            sph_legendre_p_4_m1(theta))
+        np.testing.assert_allclose(p_jac[0, 0],
+            sph_legendre_p_0_0_jac(theta))
+        np.testing.assert_allclose(p_jac[0, 1], 0)
+        np.testing.assert_allclose(p_jac[0, 2], 0)
+        np.testing.assert_allclose(p_jac[0, 3], 0)
+        np.testing.assert_allclose(p_jac[0, 4], 0)
+        np.testing.assert_allclose(p_jac[0, -3], 0)
+        np.testing.assert_allclose(p_jac[0, -2], 0)
+        np.testing.assert_allclose(p_jac[0, -1], 0)
+        np.testing.assert_allclose(p_jac[1, 0],
+            sph_legendre_p_1_0_jac(theta))
+        np.testing.assert_allclose(p_jac[1, 1],
+            sph_legendre_p_1_1_jac(theta))
+        np.testing.assert_allclose(p_jac[1, 2], 0)
+        np.testing.assert_allclose(p_jac[1, 3], 0)
+        np.testing.assert_allclose(p_jac[1, 4], 0)
+        np.testing.assert_allclose(p_jac[1, -4], 0)
+        np.testing.assert_allclose(p_jac[1, -3], 0)
+        np.testing.assert_allclose(p_jac[1, -2], 0)
+        np.testing.assert_allclose(p_jac[1, -1],
+            sph_legendre_p_1_m1_jac(theta))
+        np.testing.assert_allclose(p_jac[2, 0],
+            sph_legendre_p_2_0_jac(theta))
+        np.testing.assert_allclose(p_jac[2, 1],
+            sph_legendre_p_2_1_jac(theta))
+        np.testing.assert_allclose(p_jac[2, 2],
+            sph_legendre_p_2_2_jac(theta))
+        np.testing.assert_allclose(p_jac[2, 3], 0)
+        np.testing.assert_allclose(p_jac[2, 4], 0)
+        np.testing.assert_allclose(p_jac[2, -4], 0)
+        np.testing.assert_allclose(p_jac[2, -3], 0)
+        np.testing.assert_allclose(p_jac[2, -2],
+            sph_legendre_p_2_m2_jac(theta))
+        np.testing.assert_allclose(p_jac[2, -1],
+            sph_legendre_p_2_m1_jac(theta))
+        np.testing.assert_allclose(p_jac[3, 0],
+            sph_legendre_p_3_0_jac(theta))
+        np.testing.assert_allclose(p_jac[3, 1],
+            sph_legendre_p_3_1_jac(theta))
+        np.testing.assert_allclose(p_jac[3, 2],
+            sph_legendre_p_3_2_jac(theta))
+        np.testing.assert_allclose(p_jac[3, 3],
+            sph_legendre_p_3_3_jac(theta))
+        np.testing.assert_allclose(p_jac[3, 4], 0)
+        np.testing.assert_allclose(p_jac[3, -4], 0)
+        np.testing.assert_allclose(p_jac[3, -3],
+            sph_legendre_p_3_m3_jac(theta))
+        np.testing.assert_allclose(p_jac[3, -2],
+            sph_legendre_p_3_m2_jac(theta))
+        np.testing.assert_allclose(p_jac[3, -1],
+            sph_legendre_p_3_m1_jac(theta))
+        np.testing.assert_allclose(p_jac[4, 0],
+            sph_legendre_p_4_0_jac(theta))
+        np.testing.assert_allclose(p_jac[4, 1],
+            sph_legendre_p_4_1_jac(theta))
+        np.testing.assert_allclose(p_jac[4, 2],
+            sph_legendre_p_4_2_jac(theta))
+        np.testing.assert_allclose(p_jac[4, 3],
+            sph_legendre_p_4_3_jac(theta))
+        np.testing.assert_allclose(p_jac[4, 4],
+            sph_legendre_p_4_4_jac(theta))
+        np.testing.assert_allclose(p_jac[4, -4],
+            sph_legendre_p_4_m4_jac(theta))
+        np.testing.assert_allclose(p_jac[4, -3],
+            sph_legendre_p_4_m3_jac(theta))
+        np.testing.assert_allclose(p_jac[4, -2],
+            sph_legendre_p_4_m2_jac(theta))
+        np.testing.assert_allclose(p_jac[4, -1],
+            sph_legendre_p_4_m1_jac(theta))
+    @pytest.mark.parametrize("shape", [(10,), (4, 9), (3, 5, 7, 10)])
+    def test_ode(self, shape):
+        rng = np.random.default_rng(1234)
+        n = rng.integers(0, 10, shape)
+        m = rng.integers(-10, 10, shape)
+        theta = rng.uniform(-np.pi, np.pi, shape)
+        p, p_jac, p_hess = sph_legendre_p(n, m, theta, diff_n=2)
+        assert p.shape == shape
+        assert p_jac.shape == p.shape
+        assert p_hess.shape == p_jac.shape
+        np.testing.assert_allclose(np.sin(theta) * p_hess, -np.cos(theta) * p_jac
+            - (n * (n + 1) * np.sin(theta) - m * m / np.sin(theta)) * p,
+            rtol=1e-05, atol=1e-08)
+class TestLegendreFunctions:
+    def test_clpmn(self):
+        z = 0.5+0.3j
+        with suppress_warnings() as sup:
+            sup.filter(category=DeprecationWarning)
+            clp = special.clpmn(2, 2, z, 3)
+        assert_array_almost_equal(clp,
+                   (np.array([[1.0000, z, 0.5*(3*z*z-1)],
+                           [0.0000, np.sqrt(z*z-1), 3*z*np.sqrt(z*z-1)],
+                           [0.0000, 0.0000, 3*(z*z-1)]]),
+                    np.array([[0.0000, 1.0000, 3*z],
+                           [0.0000, z/np.sqrt(z*z-1), 3*(2*z*z-1)/np.sqrt(z*z-1)],
+                           [0.0000, 0.0000, 6*z]])),
+                    7)
+    def test_clpmn_close_to_real_2(self):
+        eps = 1e-10
+        m = 1
+        n = 3
+        x = 0.5
+        with suppress_warnings() as sup:
+            sup.filter(category=DeprecationWarning)
+            clp_plus = special.clpmn(m, n, x+1j*eps, 2)[0][m, n]
+            clp_minus = special.clpmn(m, n, x-1j*eps, 2)[0][m, n]
+        assert_array_almost_equal(np.array([clp_plus, clp_minus]),
+                                  np.array([special.lpmv(m, n, x),
+                                         special.lpmv(m, n, x)]),
+                                  7)
+    def test_clpmn_close_to_real_3(self):
+        eps = 1e-10
+        m = 1
+        n = 3
+        x = 0.5
+        with suppress_warnings() as sup:
+            sup.filter(category=DeprecationWarning)
+            clp_plus = special.clpmn(m, n, x+1j*eps, 3)[0][m, n]
+            clp_minus = special.clpmn(m, n, x-1j*eps, 3)[0][m, n]
+        assert_array_almost_equal(np.array([clp_plus, clp_minus]),
+                                  np.array([special.lpmv(m, n, x)*np.exp(-0.5j*m*np.pi),
+                                         special.lpmv(m, n, x)*np.exp(0.5j*m*np.pi)]),
+                                  7)
+    def test_clpmn_across_unit_circle(self):
+        eps = 1e-7
+        m = 1
+        n = 1
+        x = 1j
+        with suppress_warnings() as sup:
+            sup.filter(category=DeprecationWarning)
+            for type in [2, 3]:
+                assert_almost_equal(special.clpmn(m, n, x+1j*eps, type)[0][m, n],
+                                special.clpmn(m, n, x-1j*eps, type)[0][m, n], 6)
+    def test_inf(self):
+        with suppress_warnings() as sup:
+            sup.filter(category=DeprecationWarning)
+            for z in (1, -1):
+                for n in range(4):
+                    for m in range(1, n):
+                        lp = special.clpmn(m, n, z)
+                        assert np.isinf(lp[1][1,1:]).all()
+                        lp = special.lpmn(m, n, z)
+                        assert np.isinf(lp[1][1,1:]).all()
+    def test_deriv_clpmn(self):
+        # data inside and outside of the unit circle
+        zvals = [0.5+0.5j, -0.5+0.5j, -0.5-0.5j, 0.5-0.5j,
+                 1+1j, -1+1j, -1-1j, 1-1j]
+        m = 2
+        n = 3
+        with suppress_warnings() as sup:
+            sup.filter(category=DeprecationWarning)
+            for type in [2, 3]:
+                for z in zvals:
+                    for h in [1e-3, 1e-3j]:
+                        approx_derivative = (special.clpmn(m, n, z+0.5*h, type)[0]
+                                            - special.clpmn(m, n, z-0.5*h, type)[0])/h
+                        assert_allclose(special.clpmn(m, n, z, type)[1],
+                                        approx_derivative,
+                                        rtol=1e-4)
+    """
+    @pytest.mark.parametrize("m_max", [3])
+    @pytest.mark.parametrize("n_max", [5])
+    @pytest.mark.parametrize("z", [-1])
+    def test_clpmn_all_limits(self, m_max, n_max, z):
+        rng = np.random.default_rng(1234)
+        type = 2
+        p, p_jac = special.clpmn_all(m_max, n_max, type, z, diff_n=1)
+        n = np.arange(n_max + 1)
+        np.testing.assert_allclose(p_jac[0], pow(z, n + 1) * n * (n + 1) / 2)
+        np.testing.assert_allclose(p_jac[1], np.where(n >= 1, pow(z, n) * np.inf, 0))
+        np.testing.assert_allclose(p_jac[2], np.where(n >= 2,
+            -pow(z, n + 1) * (n + 2) * (n + 1) * n * (n - 1) / 4, 0))
+        np.testing.assert_allclose(p_jac[-2], np.where(n >= 2, -pow(z, n + 1) / 4, 0))
+        np.testing.assert_allclose(p_jac[-1], np.where(n >= 1, -pow(z, n) * np.inf, 0))
+        for m in range(3, m_max + 1):
+            np.testing.assert_allclose(p_jac[m], 0)
+            np.testing.assert_allclose(p_jac[-m], 0)
+    """
+    def test_lpmv(self):
+        lp = special.lpmv(0,2,.5)
+        assert_almost_equal(lp,-0.125,7)
+        lp = special.lpmv(0,40,.001)
+        assert_almost_equal(lp,0.1252678976534484,7)
+        # XXX: this is outside the domain of the current implementation,
+        #      so ensure it returns a NaN rather than a wrong answer.
+        with np.errstate(all='ignore'):
+            lp = special.lpmv(-1,-1,.001)
+        assert lp != 0 or np.isnan(lp)
+    def test_lqmn(self):
+        lqmnf = special.lqmn(0,2,.5)
+        lqf = special.lqn(2,.5)
+        assert_array_almost_equal(lqmnf[0][0],lqf[0],4)
+        assert_array_almost_equal(lqmnf[1][0],lqf[1],4)
+    def test_lqmn_gt1(self):
+        """algorithm for real arguments changes at 1.0001
+           test against analytical result for m=2, n=1
+        """
+        x0 = 1.0001
+        delta = 0.00002
+        for x in (x0-delta, x0+delta):
+            lq = special.lqmn(2, 1, x)[0][-1, -1]
+            expected = 2/(x*x-1)
+            assert_almost_equal(lq, expected)
+    def test_lqmn_shape(self):
+        a, b = special.lqmn(4, 4, 1.1)
+        assert_equal(a.shape, (5, 5))
+        assert_equal(b.shape, (5, 5))
+        a, b = special.lqmn(4, 0, 1.1)
+        assert_equal(a.shape, (5, 1))
+        assert_equal(b.shape, (5, 1))
+    def test_lqn(self):
+        lqf = special.lqn(2,.5)
+        assert_array_almost_equal(lqf,(np.array([0.5493, -0.7253, -0.8187]),
+                                       np.array([1.3333, 1.216, -0.8427])),4)
+    @pytest.mark.parametrize("function", [special.lpn, special.lqn])
+    @pytest.mark.parametrize("n", [1, 2, 4, 8, 16, 32])
+    @pytest.mark.parametrize("z_complex", [False, True])
+    @pytest.mark.parametrize("z_inexact", [False, True])
+    @pytest.mark.parametrize(
+        "input_shape",
+        [
+            (), (1, ), (2, ), (2, 1), (1, 2), (2, 2), (2, 2, 1), (2, 2, 2)
+        ]
+    )
+    def test_array_inputs_lxn(self, function, n, z_complex, z_inexact, input_shape):
+        """Tests for correct output shapes."""
+        rng = np.random.default_rng(1234)
+        if z_inexact:
+            z = rng.integers(-3, 3, size=input_shape)
+        else:
+            z = rng.uniform(-1, 1, size=input_shape)
+        if z_complex:
+            z = 1j * z + 0.5j * z
+        with suppress_warnings() as sup:
+            sup.filter(category=DeprecationWarning)
+            P_z, P_d_z = function(n, z)
+        assert P_z.shape == (n + 1, ) + input_shape
+        assert P_d_z.shape == (n + 1, ) + input_shape
+    @pytest.mark.parametrize("function", [special.lqmn])
+    @pytest.mark.parametrize(
+        "m,n",
+        [(0, 1), (1, 2), (1, 4), (3, 8), (11, 16), (19, 32)]
+    )
+    @pytest.mark.parametrize("z_inexact", [False, True])
+    @pytest.mark.parametrize(
+        "input_shape", [
+            (), (1, ), (2, ), (2, 1), (1, 2), (2, 2), (2, 2, 1)
+        ]
+    )
+    def test_array_inputs_lxmn(self, function, m, n, z_inexact, input_shape):
+        """Tests for correct output shapes and dtypes."""
+        rng = np.random.default_rng(1234)
+        if z_inexact:
+            z = rng.integers(-3, 3, size=input_shape)
+        else:
+            z = rng.uniform(-1, 1, size=input_shape)
+        P_z, P_d_z = function(m, n, z)
+        assert P_z.shape == (m + 1, n + 1) + input_shape
+        assert P_d_z.shape == (m + 1, n + 1) + input_shape
+    @pytest.mark.parametrize("function", [special.clpmn, special.lqmn])
+    @pytest.mark.parametrize(
+        "m,n",
+        [(0, 1), (1, 2), (1, 4), (3, 8), (11, 16), (19, 32)]
+    )
+    @pytest.mark.parametrize(
+        "input_shape", [
+            (), (1, ), (2, ), (2, 1), (1, 2), (2, 2), (2, 2, 1)
+        ]
+    )
+    def test_array_inputs_clxmn(self, function, m, n, input_shape):
+        """Tests for correct output shapes and dtypes."""
+        rng = np.random.default_rng(1234)
+        z = rng.uniform(-1, 1, size=input_shape)
+        z = 1j * z + 0.5j * z
+        with suppress_warnings() as sup:
+            sup.filter(category=DeprecationWarning)
+            P_z, P_d_z = function(m, n, z)
+        assert P_z.shape == (m + 1, n + 1) + input_shape
+        assert P_d_z.shape == (m + 1, n + 1) + input_shape
+def assoc_legendre_factor(n, m, norm):
+    if norm:
+        return (math.sqrt((2 * n + 1) *
+            math.factorial(n - m) / (2 * math.factorial(n + m))))
+    return 1
+def assoc_legendre_p_0_0(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(0, 0, norm)
+    return np.full_like(z, fac)
+def assoc_legendre_p_1_0(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(1, 0, norm)
+    return fac * z
+def assoc_legendre_p_1_1(z, *, branch_cut=2, norm=False):
+    branch_sign = np.where(branch_cut == 3, np.where(np.signbit(np.real(z)), 1, -1), -1)
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(1, 1, norm)
+    w = np.sqrt(np.where(branch_cut == 3, z * z - 1, 1 - z * z))
+    return branch_cut_sign * branch_sign * fac * w
+def assoc_legendre_p_1_m1(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(1, -1, norm)
+    return (-branch_cut_sign * fac *
+        assoc_legendre_p_1_1(z, branch_cut=branch_cut) / 2)
+def assoc_legendre_p_2_0(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(2, 0, norm)
+    return fac * (3 * z * z - 1) / 2
+def assoc_legendre_p_2_1(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(2, 1, norm)
+    return (3 * fac * z *
+        assoc_legendre_p_1_1(z, branch_cut=branch_cut))
+def assoc_legendre_p_2_2(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(2, 2, norm)
+    return 3 * branch_cut_sign * fac * (1 - z * z)
+def assoc_legendre_p_2_m2(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(2, -2, norm)
+    return branch_cut_sign * fac * (1 - z * z) / 8
+def assoc_legendre_p_2_m1(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(2, -1, norm)
+    return (-branch_cut_sign * fac * z *
+        assoc_legendre_p_1_1(z, branch_cut=branch_cut) / 2)
+def assoc_legendre_p_3_0(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(3, 0, norm)
+    return fac * (5 * z * z - 3) * z / 2
+def assoc_legendre_p_3_1(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(3, 1, norm)
+    return (3 * fac * (5 * z * z - 1) *
+        assoc_legendre_p_1_1(z, branch_cut=branch_cut) / 2)
+def assoc_legendre_p_3_2(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(3, 2, norm)
+    return 15 * branch_cut_sign * fac * (1 - z * z) * z
+def assoc_legendre_p_3_3(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(3, 3, norm)
+    return (15 * branch_cut_sign * fac * (1 - z * z) *
+        assoc_legendre_p_1_1(z, branch_cut=branch_cut))
+def assoc_legendre_p_3_m3(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(3, -3, norm)
+    return (fac * (z * z - 1) *
+        assoc_legendre_p_1_1(z, branch_cut=branch_cut) / 48)
+def assoc_legendre_p_3_m2(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(3, -2, norm)
+    return branch_cut_sign * fac * (1 - z * z) * z / 8
+def assoc_legendre_p_3_m1(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(3, -1, norm)
+    return (branch_cut_sign * fac * (1 - 5 * z * z) *
+        assoc_legendre_p_1_1(z, branch_cut=branch_cut) / 8)
+def assoc_legendre_p_4_0(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(4, 0, norm)
+    return fac * ((35 * z * z - 30) * z * z + 3) / 8
+def assoc_legendre_p_4_1(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(4, 1, norm)
+    return (5 * fac * (7 * z * z - 3) * z *
+       assoc_legendre_p_1_1(z, branch_cut=branch_cut) / 2)
+def assoc_legendre_p_4_2(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(4, 2, norm)
+    return 15 * branch_cut_sign * fac * ((8 - 7 * z * z) * z * z - 1) / 2
+def assoc_legendre_p_4_3(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(4, 3, norm)
+    return (105 * branch_cut_sign * fac * (1 - z * z) * z *
+        assoc_legendre_p_1_1(z, branch_cut=branch_cut))
+def assoc_legendre_p_4_4(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(4, 4, norm)
+    return 105 * fac * np.square(z * z - 1)
+def assoc_legendre_p_4_m4(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(4, -4, norm)
+    return fac * np.square(z * z - 1) / 384
+def assoc_legendre_p_4_m3(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(4, -3, norm)
+    return (fac * (z * z - 1) * z *
+        assoc_legendre_p_1_1(z, branch_cut=branch_cut) / 48)
+def assoc_legendre_p_4_m2(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(4, -2, norm)
+    return branch_cut_sign * fac * ((8 - 7 * z * z) * z * z - 1) / 48
+def assoc_legendre_p_4_m1(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(4, -1, norm)
+    return (branch_cut_sign * fac * (3 - 7 * z * z) * z *
+        assoc_legendre_p_1_1(z, branch_cut=branch_cut) / 8)
+def assoc_legendre_p_1_1_jac_div_z(z, branch_cut=2):
+    branch_sign = np.where(branch_cut == 3, np.where(np.signbit(np.real(z)), 1, -1), -1)
+    out11_div_z = (-branch_sign /
+        np.sqrt(np.where(branch_cut == 3, z * z - 1, 1 - z * z)))
+    return out11_div_z
+def assoc_legendre_p_0_0_jac(z, *, branch_cut=2, norm=False):
+    return np.zeros_like(z)
+def assoc_legendre_p_1_0_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(1, 0, norm)
+    return np.full_like(z, fac)
+def assoc_legendre_p_1_1_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(1, 1, norm)
+    return (fac * z *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut))
+def assoc_legendre_p_1_m1_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(1, -1, norm)
+    return (-branch_cut_sign * fac * z *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut) / 2)
+def assoc_legendre_p_2_0_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(2, 0, norm)
+    return 3 * fac * z
+def assoc_legendre_p_2_1_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(2, 1, norm)
+    return (3 * fac * (2 * z * z - 1) *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut))
+def assoc_legendre_p_2_2_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(2, 2, norm)
+    return -6 * branch_cut_sign * fac * z
+def assoc_legendre_p_2_m1_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(2, -1, norm)
+    return (branch_cut_sign * fac * (1 - 2 * z * z) *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut) / 2)
+def assoc_legendre_p_2_m2_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(2, -2, norm)
+    return -branch_cut_sign * fac * z / 4
+def assoc_legendre_p_3_0_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(3, 0, norm)
+    return 3 * fac * (5 * z * z - 1) / 2
+def assoc_legendre_p_3_1_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(3, 1, norm)
+    return (3 * fac * (15 * z * z - 11) * z *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut) / 2)
+def assoc_legendre_p_3_2_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(3, 2, norm)
+    return 15 * branch_cut_sign * fac * (1 - 3 * z * z)
+def assoc_legendre_p_3_3_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(3, 3, norm)
+    return (45 * branch_cut_sign * fac * (1 - z * z) * z *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut))
+def assoc_legendre_p_3_m3_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(3, -3, norm)
+    return (fac * (z * z - 1) * z *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut) / 16)
+def assoc_legendre_p_3_m2_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(3, -2, norm)
+    return branch_cut_sign * fac * (1 - 3 * z * z) / 8
+def assoc_legendre_p_3_m1_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(3, -1, norm)
+    return (branch_cut_sign * fac * (11 - 15 * z * z) * z *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut) / 8)
+def assoc_legendre_p_4_0_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(4, 0, norm)
+    return 5 * fac * (7 * z * z - 3) * z / 2
+def assoc_legendre_p_4_1_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(4, 1, norm)
+    return (5 * fac * ((28 * z * z - 27) * z * z + 3) *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut) / 2)
+def assoc_legendre_p_4_2_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(4, 2, norm)
+    return 30 * branch_cut_sign * fac * (4 - 7 * z * z) * z
+def assoc_legendre_p_4_3_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(4, 3, norm)
+    return (105 * branch_cut_sign * fac * ((5 - 4 * z * z) * z * z - 1) *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut))
+def assoc_legendre_p_4_4_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(4, 4, norm)
+    return 420 * fac * (z * z - 1) * z
+def assoc_legendre_p_4_m4_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(4, -4, norm)
+    return fac * (z * z - 1) * z / 96
+def assoc_legendre_p_4_m3_jac(z, *, branch_cut=2, norm=False):
+    fac = assoc_legendre_factor(4, -3, norm)
+    return (fac * ((4 * z * z - 5) * z * z + 1) *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut) / 48)
+def assoc_legendre_p_4_m2_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(4, -2, norm)
+    return branch_cut_sign * fac * (4 - 7 * z * z) * z / 12
+def assoc_legendre_p_4_m1_jac(z, *, branch_cut=2, norm=False):
+    branch_cut_sign = np.where(branch_cut == 3, -1, 1)
+    fac = assoc_legendre_factor(4, -1, norm)
+    return (branch_cut_sign * fac * ((27 - 28 * z * z) * z * z - 3) *
+        assoc_legendre_p_1_1_jac_div_z(z, branch_cut=branch_cut) / 8)
+def sph_legendre_factor(n, m):
+    return assoc_legendre_factor(n, m, norm=True) / np.sqrt(2 * np.pi)
+def sph_legendre_p_0_0(theta):
+    fac = sph_legendre_factor(0, 0)
+    return np.full_like(theta, fac)
+def sph_legendre_p_1_0(theta):
+    fac = sph_legendre_factor(1, 0)
+    return fac * np.cos(theta)
+def sph_legendre_p_1_1(theta):
+    fac = sph_legendre_factor(1, 1)
+    return -fac * np.abs(np.sin(theta))
+def sph_legendre_p_1_m1(theta):
+    fac = sph_legendre_factor(1, -1)
+    return fac * np.abs(np.sin(theta)) / 2
+def sph_legendre_p_2_0(theta):
+    fac = sph_legendre_factor(2, 0)
+    return fac * (3 * np.square(np.cos(theta)) - 1) / 2
+def sph_legendre_p_2_1(theta):
+    fac = sph_legendre_factor(2, 1)
+    return -3 * fac * np.abs(np.sin(theta)) * np.cos(theta)
+def sph_legendre_p_2_2(theta):
+    fac = sph_legendre_factor(2, 2)
+    return 3 * fac * (1 - np.square(np.cos(theta)))
+def sph_legendre_p_2_m2(theta):
+    fac = sph_legendre_factor(2, -2)
+    return fac * (1 - np.square(np.cos(theta))) / 8
+def sph_legendre_p_2_m1(theta):
+    fac = sph_legendre_factor(2, -1)
+    return fac * np.cos(theta) * np.abs(np.sin(theta)) / 2
+def sph_legendre_p_3_0(theta):
+    fac = sph_legendre_factor(3, 0)
+    return (fac * (5 * np.square(np.cos(theta)) - 3) *
+        np.cos(theta) / 2)
+def sph_legendre_p_3_1(theta):
+    fac = sph_legendre_factor(3, 1)
+    return (-3 * fac * (5 * np.square(np.cos(theta)) - 1) *
+        np.abs(np.sin(theta)) / 2)
+def sph_legendre_p_3_2(theta):
+    fac = sph_legendre_factor(3, 2)
+    return (-15 * fac * (np.square(np.cos(theta)) - 1) *
+        np.cos(theta))
+def sph_legendre_p_3_3(theta):
+    fac = sph_legendre_factor(3, 3)
+    return -15 * fac * np.power(np.abs(np.sin(theta)), 3)
+def sph_legendre_p_3_m3(theta):
+    fac = sph_legendre_factor(3, -3)
+    return fac * np.power(np.abs(np.sin(theta)), 3) / 48
+def sph_legendre_p_3_m2(theta):
+    fac = sph_legendre_factor(3, -2)
+    return (-fac * (np.square(np.cos(theta)) - 1) *
+        np.cos(theta) / 8)
+def sph_legendre_p_3_m1(theta):
+    fac = sph_legendre_factor(3, -1)
+    return (fac * (5 * np.square(np.cos(theta)) - 1) *
+        np.abs(np.sin(theta)) / 8)
+def sph_legendre_p_4_0(theta):
+    fac = sph_legendre_factor(4, 0)
+    return (fac * (35 * np.square(np.square(np.cos(theta))) -
+        30 * np.square(np.cos(theta)) + 3) / 8)
+def sph_legendre_p_4_1(theta):
+    fac = sph_legendre_factor(4, 1)
+    return (-5 * fac * (7 * np.square(np.cos(theta)) - 3) *
+        np.cos(theta) * np.abs(np.sin(theta)) / 2)
+def sph_legendre_p_4_2(theta):
+    fac = sph_legendre_factor(4, 2)
+    return (-15 * fac * (7 * np.square(np.cos(theta)) - 1) *
+        (np.square(np.cos(theta)) - 1) / 2)
+def sph_legendre_p_4_3(theta):
+    fac = sph_legendre_factor(4, 3)
+    return -105 * fac * np.power(np.abs(np.sin(theta)), 3) * np.cos(theta)
+def sph_legendre_p_4_4(theta):
+    fac = sph_legendre_factor(4, 4)
+    return 105 * fac * np.square(np.square(np.cos(theta)) - 1)
+def sph_legendre_p_4_m4(theta):
+    fac = sph_legendre_factor(4, -4)
+    return fac * np.square(np.square(np.cos(theta)) - 1) / 384
+def sph_legendre_p_4_m3(theta):
+    fac = sph_legendre_factor(4, -3)
+    return (fac * np.power(np.abs(np.sin(theta)), 3) *
+        np.cos(theta) / 48)
+def sph_legendre_p_4_m2(theta):
+    fac = sph_legendre_factor(4, -2)
+    return (-fac * (7 * np.square(np.cos(theta)) - 1) *
+        (np.square(np.cos(theta)) - 1) / 48)
+def sph_legendre_p_4_m1(theta):
+    fac = sph_legendre_factor(4, -1)
+    return (fac * (7 * np.square(np.cos(theta)) - 3) *
+        np.cos(theta) * np.abs(np.sin(theta)) / 8)
+def sph_legendre_p_0_0_jac(theta):
+    return np.zeros_like(theta)
+def sph_legendre_p_1_0_jac(theta):
+    fac = sph_legendre_factor(1, 0)
+    return -fac * np.sin(theta)
+def sph_legendre_p_1_1_jac(theta):
+    fac = sph_legendre_factor(1, 1)
+    return -fac * np.cos(theta) * (2 * np.heaviside(np.sin(theta), 1) - 1)
+def sph_legendre_p_1_m1_jac(theta):
+    fac = sph_legendre_factor(1, -1)
+    return fac * np.cos(theta) * (2 * np.heaviside(np.sin(theta), 1) - 1) / 2
+def sph_legendre_p_2_0_jac(theta):
+    fac = sph_legendre_factor(2, 0)
+    return -3 * fac * np.cos(theta) * np.sin(theta)
+def sph_legendre_p_2_1_jac(theta):
+    fac = sph_legendre_factor(2, 1)
+    return (3 * fac * (-np.square(np.cos(theta)) *
+        (2 * np.heaviside(np.sin(theta), 1) - 1) +
+        np.abs(np.sin(theta)) * np.sin(theta)))
+def sph_legendre_p_2_2_jac(theta):
+    fac = sph_legendre_factor(2, 2)
+    return 6 * fac * np.sin(theta) * np.cos(theta)
+def sph_legendre_p_2_m2_jac(theta):
+    fac = sph_legendre_factor(2, -2)
+    return fac * np.sin(theta) * np.cos(theta) / 4
+def sph_legendre_p_2_m1_jac(theta):
+    fac = sph_legendre_factor(2, -1)
+    return (-fac * (-np.square(np.cos(theta)) *
+        (2 * np.heaviside(np.sin(theta), 1) - 1) +
+        np.abs(np.sin(theta)) * np.sin(theta)) / 2)
+def sph_legendre_p_3_0_jac(theta):
+    fac = sph_legendre_factor(3, 0)
+    return 3 * fac * (1 - 5 * np.square(np.cos(theta))) * np.sin(theta) / 2
+def sph_legendre_p_3_1_jac(theta):
+    fac = sph_legendre_factor(3, 1)
+    return (3 * fac * (11 - 15 * np.square(np.cos(theta))) * np.cos(theta) *
+        (2 * np.heaviside(np.sin(theta), 1) - 1) / 2)
+def sph_legendre_p_3_2_jac(theta):
+    fac = sph_legendre_factor(3, 2)
+    return 15 * fac * (3 * np.square(np.cos(theta)) - 1) * np.sin(theta)
+def sph_legendre_p_3_3_jac(theta):
+    fac = sph_legendre_factor(3, 3)
+    return -45 * fac * np.abs(np.sin(theta)) * np.sin(theta) * np.cos(theta)
+def sph_legendre_p_3_m3_jac(theta):
+    fac = sph_legendre_factor(3, -3)
+    return fac * np.abs(np.sin(theta)) * np.sin(theta) * np.cos(theta) / 16
+def sph_legendre_p_3_m2_jac(theta):
+    fac = sph_legendre_factor(3, -2)
+    return fac * (3 * np.square(np.cos(theta)) - 1) * np.sin(theta) / 8
+def sph_legendre_p_3_m1_jac(theta):
+    fac = sph_legendre_factor(3, -1)
+    return (-fac * (11 - 15 * np.square(np.cos(theta))) *
+        np.cos(theta) *
+        (2 * np.heaviside(np.sin(theta), 1) - 1) / 8)
+def sph_legendre_p_4_0_jac(theta):
+    fac = sph_legendre_factor(4, 0)
+    return (-5 * fac * (7 * np.square(np.cos(theta)) - 3) *
+        np.sin(theta) * np.cos(theta) / 2)
+def sph_legendre_p_4_1_jac(theta):
+    fac = sph_legendre_factor(4, 1)
+    return (5 * fac * (-3 + 27 * np.square(np.cos(theta)) -
+        28 * np.square(np.square(np.cos(theta)))) *
+        (2 * np.heaviside(np.sin(theta), 1) - 1) / 2)
+def sph_legendre_p_4_2_jac(theta):
+    fac = sph_legendre_factor(4, 2)
+    return (30 * fac * (7 * np.square(np.cos(theta)) - 4) *
+        np.sin(theta) * np.cos(theta))
+def sph_legendre_p_4_3_jac(theta):
+    fac = sph_legendre_factor(4, 3)
+    return (-105 * fac * (4 * np.square(np.cos(theta)) - 1) *
+        np.abs(np.sin(theta)) * np.sin(theta))
+def sph_legendre_p_4_4_jac(theta):
+    fac = sph_legendre_factor(4, 4)
+    return (-420 * fac * (np.square(np.cos(theta)) - 1) *
+        np.sin(theta) * np.cos(theta))
+def sph_legendre_p_4_m4_jac(theta):
+    fac = sph_legendre_factor(4, -4)
+    return (-fac * (np.square(np.cos(theta)) - 1) *
+        np.sin(theta) * np.cos(theta) / 96)
+def sph_legendre_p_4_m3_jac(theta):
+    fac = sph_legendre_factor(4, -3)
+    return (fac * (4 * np.square(np.cos(theta)) - 1) *
+        np.abs(np.sin(theta)) * np.sin(theta) / 48)
+def sph_legendre_p_4_m2_jac(theta):
+    fac = sph_legendre_factor(4, -2)
+    return (fac * (7 * np.square(np.cos(theta)) - 4) * np.sin(theta) *
+        np.cos(theta) / 12)
+def sph_legendre_p_4_m1_jac(theta):
+    fac = sph_legendre_factor(4, -1)
+    return (-fac * (-3 + 27 * np.square(np.cos(theta)) -
+        28 * np.square(np.square(np.cos(theta)))) *
+        (2 * np.heaviside(np.sin(theta), 1) - 1) / 8)