SciPy 1.17.0 Release Notes

Note: SciPy 1.17.0 is not released yet!

SciPy 1.17.0 is the culmination of 6 months of hard work. It contains
many new features, numerous bug-fixes, improved test coverage and better
documentation. There have been a number of deprecations and API changes
in this release, which are documented below. All users are encouraged to
upgrade to this release, as there are a large number of bug-fixes and
optimizations. Before upgrading, we recommend that users check that
their own code does not use deprecated SciPy functionality (to do so,
run your code with python -Wd and check for DeprecationWarning s).
Our development attention will now shift to bug-fix releases on the
1.17.x branch, and on adding new features on the main branch.

This release requires Python 3.11-3.14 and NumPy 1.26.4 or greater.

Highlights of this release

• Many SciPy functions have gained native support for batching of N-dimensional
  array input and additional support for the array API standard. An overall
  summary of the latter is now available in a set of tables.
• In scipy.sparse, coo_array now has full support for indexing across
  dimensions without needing to convert between sparse formats. ARPACK
  and PROPACK rewrites from Fortran77 to C now empower the use of external
  pseudorandom number generators.
• In scipy.spatial, transform.Rotation and transform.RigidTransform
  have been extended to support N-D arrays. geometric_slerp now has support
  for extrapolation.
• scipy.stats has gained the matrix t and logistic distributions and many
  performance and accuracy improvements.
• Initial support for 64-bit integer (ILP64) BLAS and LAPACK libraries has
  been added, including for MKL, Apple Accelerate and OpenBLAS. Please
  report any issues with ILP64 you encounter.

New features

scipy.integrate improvements

• The integration routines dopri5, dopri853, LSODA, vode, and
  zvode have been ported from Fortran77 to C.
• scipy.integrate.quad now has a fast path for returning 0 when the integration
  interval is empty.

scipy.cluster improvements

• scipy.cluster.hierarchy.is_isomorphic has improved performance and array
  API support.

scipy.interpolate improvements

• A new bc_type argument has been added to scipy.interpolate.make_splrep
  and scipy.interpolate.make_splprep to control the boundary conditions for
  spline fitting. Allowed values are "not-a-knot" (default) and
  "periodic".
• A new derivative method has been added to the
  scipy.interpolate.NdBSpline class, to construct a new spline representing a
  partial derivative of the given spline. This method is similar to the
  BSpline.derivative method of 1-D spline objects.
• Performance of "cubic" and "quintic" modes of
  scipy.interpolate.RegularGridInterpolator has been improved.
• Numerical stability of scipy.interpolate.AAA has been improved.
• scipy.interpolate.FloaterHormannInterpolator added support for
  multidimensional, batched inputs and gained a new axis parameter to
  select the interpolation axis.

scipy.linalg improvements

• scipy.linalg.inv routine has been improved:

  • it now attempts to detect the structure of its argument and selects an
    appropriate low-level matrix inversion routine. A new assume_a keyword
    allows to bypass the structure detection if the structure is known. For
    batched inputs, the detection is run for each 2D slice, unless an explicit
    value for assume_a is provided (in which case, the structure is
    assumed to be the same for all 2-D slices of the batch);
  • the new lower={True,False} keyword argument has been added to help
    select the upper or lower triangle of the input matrix for symmetric
    inputs; refer to the docstring of scipy.linalg.inv for details;
  • the routine emits a LinAlgWarning if it detects an ill-conditioned
    input;
  • performance for batched inputs has been improved.

• scipy.linalg.fiedler has gained native support for batched inputs.

• performance has improved for scipy.linalg.solve with batched inputs
  for certain matrix structures.

scipy.optimize improvements

• optimize.minimize(method="trust-exact") now accepts a
  solver-specific "subproblem_maxiter" option. This option can be used to
  assure that the algorithm converges for functions with an ill-conditioned
  Hessian.
• Callback functions used by optimize.minimize(method="slsqp") can
  opt into the new callback interface by accepting a single keyword argument
  intermediate_result.

scipy.signal improvements

• scipy.signal.abcd_normalize gained more informative error messages and the
  documentation was improved.
• scipy.signal.get_window now accepts the suffixes '_periodic' and
  '_symmetric' to distinguish between periodic and symmetric windows
  (overriding the fftbin parameter). This benefits the functions
  coherence, csd, periodogram, welch, spectrogram,
  stft, istft, resample, resample_poly, firwin,
  firwin2, firwin_2d, check_COLA and check_NOLA, which utilize
  get_window but do not expose the fftbin parameter.
• scipy.signal.hilbert2 gained the new keyword axes for specifying the
  axes along which the two-dimensional analytic signal should be calculated.
  Furthermore, the documentation of scipy.signal.hilbert and
  scipy.signal.hilbert2 was significantly improved.

scipy.sparse improvements

• coo_array now supports indexing. This includes slices, arrays,
  np.newaxis, Ellipsis, in 1D, 2D and the new nD. So COO format now
  has full support for nD and COO now allows indexing without converting
  formats.
• Additional sparse construction functions include expand_dims,
  swapaxes, permute_dims, and nD support for the kron function.
• ARPACK Fortran77 library is ported to C. Among many changes, it is now
  possible to use external random generators including NumPy PRNGs for
  reproducible runs. Previously this was not the case due to internal seeding
  behavior of the original ARPACK code.
• Similarly, PROPACK Fortran77 library is also ported to C with the same PRNG
  enhancements and other improvements.
• scipy.sparse.dok_array now supports an update method which can be
  used to update the sparse array using a dict, dict.items()-like iterable,
  or another dok_array matrix. It performs additional validation that keys
  are valid index tuples.
• scipy.sparse.dia_array.tocsr is approximately three times faster and
  some unneccesary copy operations have been removed from sparse format
  interconversions more broadly.
• Added scipy.sparse.linalg.funm_multiply_krylov, a restarted Krylov method
  for evaluating y = f(tA) b.

scipy.spatial improvements

• The spatial.transform module has gained an array API standard compatible
  backend.

• transform.Rotation and transform.RigidTransform have been extended
  from 0D single values and 1D arrays to N-D arrays, with standard indexing and
  broadcasting rules. Both now have the following additions:

  • A shape property.
  • A shape argument to their identity() constructors, which should be
    preferred over the existing num argument. This has also been added as an
    argument for Rotation.random() (RigidTransform does not currently
    have a random constructor).
  • An axis argument to their mean() functions.

• The resulting shapes for transform.Rotation.from_euler /
  from_davenport have changed to make them consistent with broadcasting
  rules. Angle inputs to Euler angles must now strictly match the number of
  provided axes in the last dimension. The resulting Rotation has the shape
  np.atleast_1d(angles).shape[:-1]. Angle inputs to Davenport angles must
  also match the number of axes in the last dimension. The resulting Rotation
  has the shape np.broadcast_shapes(np.atleast_2d(axes).shape[:-2], np.atleast_1d(angles).shape[:-1]).

• Rotation.from_matrix has gained an assume_valid argument that allows for
  performance improvements when users can guarantee valid matrix inputs.
  from_matrix is now also faster in cases where a known orthogonal matrix
  is used.

• The scipy.spatial.geometric_slerp function can now extrapolate. When given a
  value outside the range [0, 1], geometric_slerp() will continue with
  the same rotation outside this range. For example, if spherically
  interpolating with start being a point on the equator, and end
  being a point at the north pole, then a value of t=-1 would give you a
  point at the south pole.

• Rotation.as_euler and Rotation.as_davenport methods have gained a
  suppress_warnings parameter to enable suppression of gimbal lock warnings.

scipy.special improvements

• The following functions for statistical applications have significantly
  improved parameter ranges and reduced error rates: btdtria, btdtrib,
  chdtriv, chndtr, chndtrix, chndtridf, chndtrinc, fdtr,
  fdtrc, fdtri, gdtria, gdtrix, pdtrik, stdtr and
  stdtrit.
• The incomplete beta functions betainc, betaincc, betaincinv and
  betainccinv are improved for...

SciPy 1.16.3 Release Notes

SciPy 1.16.3 is a bug-fix release with no new features compared to 1.16.2.

Authors

• Name (commits)
• ChrisAB (1) +
• Lucas Colley (1)
• Ralf Gommers (3)
• Matt Haberland (8)
• Nick ODell (2)
• Ilhan Polat (1)
• Tyler Reddy (28)
• Lucas Roberts (2)

A total of 8 people contributed to this release.
People with a "+" by their names contributed a patch for the first time.
This list of names is automatically generated, and may not be fully complete.

The full issue and pull request lists, and the release asset hashes are available
in the associated README.txt file.

SciPy 1.16.2 Release Notes

SciPy 1.16.2 is a bug-fix release with no new features
compared to 1.16.1. This is the first stable release of
SciPy to provide Windows on ARM wheels on PyPI.

Authors

• Name (commits)
• Dietrich Brunn (1)
• Ralf Gommers (6)
• Adam Jones (1)
• Gleb Khmyznikov (1) +
• Jost Migenda (1) +
• newyork_loki (1)
• Nick ODell (3)
• Dimitri Papadopoulos Orfanos (1)
• Ilhan Polat (2)
• Tyler Reddy (26)
• Mugunthan Selvanayagam (1) +
• Shuhei Watanabe (1) +

A total of 12 people contributed to this release.
People with a "+" by their names contributed a patch for the first time.
This list of names is automatically generated, and may not be fully complete.

The full issue and pull request lists, and the release asset hashes are available
in the associated README.txt file.

SciPy 1.16.1 Release Notes

SciPy 1.16.1 is a bug-fix release that adds support for Python 3.14.0rc1,
including PyPI wheels.

Authors

• Name (commits)
• Evgeni Burovski (1)
• Rob Falck (1)
• Ralf Gommers (7)
• Geoffrey Gunter (1) +
• Matt Haberland (2)
• Joren Hammudoglu (1)
• Andrew Nelson (2)
• newyork_loki (1) +
• Ilhan Polat (1)
• Tyler Reddy (25)
• Daniel Schmitz (1)
• Dan Schult (2)

A total of 12 people contributed to this release.
People with a "+" by their names contributed a patch for the first time.
This list of names is automatically generated, and may not be fully complete.

The full issue and pull request lists, and the release asset hashes are available
in the associated README.txt file.

SciPy 1.16.0 Release Notes

SciPy 1.16.0 is the culmination of 6 months of hard work. It contains
many new features, numerous bug-fixes, improved test coverage and better
documentation. There have been a number of deprecations and API changes
in this release, which are documented below. All users are encouraged to
upgrade to this release, as there are a large number of bug-fixes and
optimizations. Before upgrading, we recommend that users check that
their own code does not use deprecated SciPy functionality (to do so,
run your code with python -Wd and check for DeprecationWarning s).
Our development attention will now shift to bug-fix releases on the
1.16.x branch, and on adding new features on the main branch.

This release requires Python 3.11-3.13 and NumPy 1.25.2 or greater.

Highlights of this release

• Improved experimental support for the Python array API standard, including
  new support in scipy.signal, and additional support in scipy.stats and
  scipy.special. Improved support for JAX and Dask backends has been added,
  with notable support in scipy.cluster.hierarchy, many functions in
  scipy.special, and many of the trimmed statistics functions.
• scipy.optimize now uses the new Python implementation from the
  PRIMA package for COBYLA. The PRIMA implementation fixes many bugs
  in the old Fortran 77 implementation with a better performance on average.
• scipy.sparse.coo_array now supports n-D arrays with reshaping, arithmetic and
  reduction operations like sum/mean/min/max. No n-D indexing or
  scipy.sparse.random_array support yet.
• Updated guide and tools for migration from sparse matrices to sparse arrays.
• Nearly all functions in the scipy.linalg namespace that accept array
  arguments now support N-dimensional arrays to be processed as a batch.
• Two new scipy.signal functions, scipy.signal.firwin_2d and
  scipy.signal.closest_STFT_dual_window, for creation of a 2-D FIR filter and
  scipy.signal.ShortTimeFFT dual window calculation, respectively.
• A new class, scipy.spatial.transform.RigidTransform, provides functionality
  to convert between different representations of rigid transforms in 3-D
  space.
• A new function scipy.ndimage.vectorized_filter for generic filters that
  take advantage of a vectorized Python callable was added.

New features

scipy.io improvements

• scipy.io.savemat now provides informative warnings for invalid field names.
• scipy.io.mmread now provides a clearer error message when provided with
  a source file path that does not exist.
• scipy.io.wavfile.read can now read non-seekable files.

scipy.integrate improvements

• The error estimate of scipy.integrate.tanhsinh was improved.

scipy.interpolate improvements

• Batch support was added to scipy.interpolate.make_smoothing_spline.

scipy.linalg improvements

• Nearly all functions in the scipy.linalg namespace that accept array
  arguments now support N-dimensional arrays to be processed as a batch.
  See linalg_batch for details.
• scipy.linalg.sqrtm is rewritten in C and its performance is improved. It
  also tries harder to return real-valued results for real-valued inputs if
  possible. See the function docstring for more details. In this version the
  input argument disp and the optional output argument errest are
  deprecated and will be removed four versions later. Similarly, after
  changing the underlying algorithm to recursion, the blocksize keyword
  argument has no effect and will be removed two versions later.
• Wrappers for ?stevd, ?langb, ?sytri, ?hetri and
  ?gbcon were added to scipy.linalg.lapack.
• The default driver of scipy.linalg.eigh_tridiagonal was improved.
• scipy.linalg.solve can now estimate the reciprocal condition number and
  the matrix norm calculation is more efficient.

scipy.ndimage improvements

• A new function scipy.ndimage.vectorized_filter for generic filters that
  take advantage of a vectorized Python callable was added.
• scipy.ndimage.rotate has improved performance, especially on ARM platforms.

scipy.optimize improvements

• COBYLA was updated to use the new Python implementation from the
  PRIMApackage.
  The PRIMA implementation fixes many bugs
  in the old Fortran 77 implementation. In addition, it results in fewer function evaluations on average
  but it depends on the problem and for some
  problems it can result in more function evaluations or a less optimal
  result. For those cases the user can try modifying the initial and final
  trust region radii given by rhobeg and tol respectively. A larger
  rhobeg can help the algorithm take bigger steps initially, while a
  smaller tol can help it continue and find a better solution.
  For more information, see the PRIMA documentation.
• Several of the scipy.optimize.minimize methods, and the
  scipy.optimize.least_squares function, have been given a workers
  keyword. This allows parallelization of some calculations via a map-like
  callable, such as multiprocessing.Pool. These parallelization
  opportunities typically occur during numerical differentiation. This can
  greatly speed up minimization when the objective function is expensive to
  calculate.
• The lm method of scipy.optimize.least_squares can now accept
  3-point and cs for the jac keyword.
• The SLSQP Fortran 77 code was ported to C. When this method is used now the
  constraint multipliers are exposed to the user through the multiplier
  keyword of the returned scipy.optimize.OptimizeResult object.
• NNLS code has been corrected and rewritten in C to address the performance
  regression introduced in 1.15.x
• scipy.optimize.root now warns for invalid inner parameters when using the
  newton_krylov method
• The return value of minimization with method='L-BFGS-B' now has
  a faster hess_inv.todense() implementation. Time complexity has improved
  from cubic to quadratic.
• scipy.optimize.least_squares has a new callback argument that is applicable
  to the trf and dogbox methods. callback may be used to track
  optimization results at each step or to provide custom conditions for
  stopping.

scipy.signal improvements

• A new function scipy.signal.firwin_2d for the creation of a 2-D FIR Filter
  using the 1-D window method was added.
• scipy.signal.cspline1d_eval and scipy.signal.qspline1d_eval now provide
  an informative error on empty input rather than hitting the recursion limit.
• A new function scipy.signal.closest_STFT_dual_window to calculate the
  scipy.signal.ShortTimeFFT dual window of a given window closest to a
  desired dual window.
• A new classmethod scipy.signal.ShortTimeFFT.from_win_equals_dual to
  create a scipy.signal.ShortTimeFFT instance where the window and its dual
  are equal up to a scaling factor. It allows to create short-time Fourier
  transforms which are unitary mappings.
• The performance of scipy.signal.convolve2d was improved.

scipy.sparse improvements

• scipy.sparse.coo_array now supports n-D arrays using binary and reduction
  operations.
• Faster operations between two DIA arrays/matrices for: add, sub, multiply,
  matmul.
• scipy.sparse.csgraph.dijkstra shortest_path is more efficient.
• scipy.sparse.csgraph.yen has performance improvements.
• Support for lazy loading of sparse.csgraph and sparse.linalg was
  added.

scipy.spatial improvements

• A new class, scipy.spatial.transform.RigidTransform, provides functionality
  to convert between different representations of rigid transforms in 3-D
  space, its application to vectors and transform composition.
  It follows the same design approach as scipy.spatial.transform.Rotation.
• scipy.spatial.transform.Rotation now has an appropriate __repr__ method,
  and improved performance for its scipy.spatial.transform.Rotation.apply
  method.

scipy.stats improvements

• A new function scipy.stats.quantile, an array API compatible function for
  quantile estimation, was added.
• scipy.stats.make_distribution was extended to work with existing discrete
  distributions and to facilitate the creation of custom distributions in the
  new random variable infrastructure.
• A new distribution, scipy.stats.Binomial, was added.
• An equal_var keyword was added to scipy.stats.tukey_hsd (enables the
  Games-Howell test) and scipy.stats.f_oneway (enables Welch ANOVA).
• The moment calculation for scipy.stats.gennorm was improved.
• The scipy.stats.mode implementation was vectorized, for faster batch
  calculation.
• Support for axis, nan_policy, and keepdims keywords was added to
  scipy.stats.power_divergence, scipy.stats.chisquare,
  scipy.stats.pointbiserialr, scipy.stats.kendalltau,
  scipy.stats.weightedtau, scipy.stats.theilslopes,
  scipy.stats.siegelslopes, scipy.stats.boxcox_llf, and
  scipy.stats.linregress.
• Support for keepdims and nan_policy keywords was added to
  `scipy....

SciPy 1.16.0 Release Notes

Note: SciPy 1.16.0 is not released yet!

SciPy 1.16.0 is the culmination of 6 months of hard work. It contains
many new features, numerous bug-fixes, improved test coverage and better
documentation. There have been a number of deprecations and API changes
in this release, which are documented below. All users are encouraged to
upgrade to this release, as there are a large number of bug-fixes and
optimizations. Before upgrading, we recommend that users check that
their own code does not use deprecated SciPy functionality (to do so,
run your code with python -Wd and check for DeprecationWarning s).
Our development attention will now shift to bug-fix releases on the
1.16.x branch, and on adding new features on the main branch.

This release requires Python 3.11-3.13 and NumPy 1.25.2 or greater.

Highlights of this release

• Improved experimental support for the Python array API standard, including
  new support in scipy.signal, and additional support in scipy.stats and
  scipy.special. Improved support for JAX and Dask backends has been added,
  with notable support in scipy.cluster.hierarchy, many functions in
  scipy.special, and many of the trimmed statistics functions.
• scipy.optimize now uses the new Python implementation from the PRIMA package for COBYLA.
  The PRIMA implementation fixes many bugs in the old Fortran 77 implementation with
  a better performance on average.
• scipy.sparse.coo_array now supports n-D arrays with reshaping, arithmetic and
  reduction operations like sum/mean/min/max. No n-D indexing or
  scipy.sparse.random_array support yet.
• Updated guide and tools for migration from sparse matrices to sparse arrays.
• Nearly all functions in the scipy.linalg namespace that accept array
  arguments now support N-dimensional arrays to be processed as a batch.
• Two new scipy.signal functions, scipy.signal.firwin_2d and
  scipy.signal.closest_STFT_dual_window, for creation of a 2-D FIR filter and
  scipy.signal.ShortTimeFFT dual window calculation, respectively.
• A new class, scipy.spatial.transform.RigidTransform, provides functionality
  to convert between different representations of rigid transforms in 3-D
  space.

New features

scipy.io improvements

• scipy.io.savemat now provides informative warnings for invalid field names.
• scipy.io.mmread now provides a clearer error message when provided with
  a source file path that does not exist.
• scipy.io.wavfile.read can now read non-seekable files.

scipy.integrate improvements

• The error estimate of scipy.integrate.tanhsinh was improved.

scipy.interpolate improvements

• Batch support was added to scipy.interpolate.make_smoothing_spline.

scipy.linalg improvements

• Nearly all functions in the scipy.linalg namespace that accept array
  arguments now support N-dimensional arrays to be processed as a batch.
  See :ref:linalg_batch for details.
• scipy.linalg.sqrtm is rewritten in C and its performance is improved. It
  also tries harder to return real-valued results for real-valued inputs if
  possible. See the function docstring for more details. In this version the
  input argument disp and the optional output argument errest are
  deprecated and will be removed four versions later. Similarly, after
  changing the underlying algorithm to recursion, the blocksize keyword
  argument has no effect and will be removed two versions later.
• Wrappers for ?stevd, ?langb, ?sytri, ?hetri and
  ?gbcon were added to scipy.linalg.lapack.
• The default driver of scipy.linalg.eigh_tridiagonal was improved.
• scipy.linalg.solve can now estimate the reciprocal condition number and
  the matrix norm calculation is more efficient.

scipy.ndimage improvements

• A new function scipy.ndimage.vectorized_filter for generic filters that
  take advantage of a vectorized Python callable was added.
• scipy.ndimage.rotate has improved performance, especially on ARM platforms.

scipy.optimize improvements

• COBYLA was updated to use the new Python implementation from the
  PRIMA package.
  The PRIMA implementation fixes many bugs
  in the old Fortran 77 implementation. In addition, it results in fewer function evaluations on average,
  but it depends on the problem and for some
  problems it can result in more function evaluations or a less optimal
  result. For those cases the user can try modifying the initial and final
  trust region radii given by rhobeg and tol respectively. A larger
  rhobeg can help the algorithm take bigger steps initially, while a
  smaller tol can help it continue and find a better solution.
  For more information, see the PRIMA documentation.
• Several of the scipy.optimize.minimize methods, and the
  scipy.optimize.least_squares function, have been given a workers
  keyword. This allows parallelization of some calculations via a map-like
  callable, such as multiprocessing.Pool. These parallelization
  opportunities typically occur during numerical differentiation. This can
  greatly speed up minimization when the objective function is expensive to
  calculate.
• The lm method of scipy.optimize.least_squares can now accept
  3-point and cs for the jac keyword.
• The SLSQP Fortran 77 code was ported to C. When this method is used now the
  constraint multipliers are exposed to the user through the multiplier
  keyword of the returned scipy.optimize.OptimizeResult object.
• NNLS code has been corrected and rewritten in C to address the performance
  regression introduced in 1.15.x
• scipy.optimize.root now warns for invalid inner parameters when using the
  newton_krylov method
• The return value of minimization with method='L-BFGS-B' now has
  a faster hess_inv.todense() implementation. Time complexity has improved
  from cubic to quadratic.
• scipy.optimize.least_squares has a new callback argument that is applicable
  to the trf and dogbox methods. callback may be used to track
  optimization results at each step or to provide custom conditions for
  stopping.

scipy.signal improvements

• A new function scipy.signal.firwin_2d for the creation of a 2-D FIR Filter
  using the 1-D window method was added.
• scipy.signal.cspline1d_eval and scipy.signal.qspline1d_eval now provide
  an informative error on empty input rather than hitting the recursion limit.
• A new function scipy.signal.closest_STFT_dual_window to calculate the
  scipy.signal.ShortTimeFFT dual window of a given window closest to a
  desired dual window.
• A new classmethod scipy.signal.ShortTimeFFT.from_win_equals_dual to
  create a scipy.signal.ShortTimeFFT instance where the window and its dual
  are equal up to a scaling factor. It allows to create short-time Fourier
  transforms which are unitary mappings.
• The performance of scipy.signal.convolve2d was improved.

scipy.sparse improvements

• scipy.sparse.coo_array now supports n-D arrays using binary and reduction
  operations.
• Faster operations between two DIA arrays/matrices for: add, sub, multiply,
  matmul.
• scipy.sparse.csgraph.dijkstra shortest_path is more efficient.
• scipy.sparse.csgraph.yen has performance improvements.
• Support for lazy loading of sparse.csgraph and sparse.linalg was
  added.

scipy.spatial improvements

• A new class, scipy.spatial.transform.RigidTransform, provides functionality
  to convert between different representations of rigid transforms in 3-D
  space, its application to vectors and transform composition.
  It follows the same design approach as scipy.spatial.transform.Rotation.
• scipy.spatial.transform.Rotation now has an appropriate __repr__ method,
  and improved performance for its scipy.spatial.transform.Rotation.apply
  method.

scipy.stats improvements

• A new function scipy.stats.quantile, an array API compatible function for
  quantile estimation, was added.
• scipy.stats.make_distribution was extended to work with existing discrete
  distributions and to facilitate the creation of custom distributions in the
  new random variable infrastructure.
• A new distribution, scipy.stats.Binomial, was added.
• An equal_var keyword was added to scipy.stats.tukey_hsd (enables the
  Games-Howell test) and scipy.stats.f_oneway (enables Welch ANOVA).
• The moment calculation for scipy.stats.gennorm was improved.
• The scipy.stats.mode implementation was vectorized, for faster batch
  calculation.
• Support for axis, nan_policy, and keepdims keywords was added to
  scipy.stats.power_divergence, scipy.stats.chisquare,
  scipy.stats.pointbiserialr, scipy.stats.kendalltau,
  scipy.stats.weightedtau, scipy.stats.theilslopes,
  scipy.stats.siegelslopes, scipy.stats.boxcox_llf, and
  scipy.stats.linregress.
• Support for keepdims and nan_policy keywords was added to
  scipy.stats.gstd.
• The performance of scipy.stats.special_ortho_group and scipy.stats.pearsonr
  was improved.
• Support for an ``rng`...

SciPy 1.16.0 Release Notes

Note: SciPy 1.16.0 is not released yet!

SciPy 1.16.0 is the culmination of 6 months of hard work. It contains
many new features, numerous bug-fixes, improved test coverage and better
documentation. There have been a number of deprecations and API changes
in this release, which are documented below. All users are encouraged to
upgrade to this release, as there are a large number of bug-fixes and
optimizations. Before upgrading, we recommend that users check that
their own code does not use deprecated SciPy functionality (to do so,
run your code with python -Wd and check for DeprecationWarning s).
Our development attention will now shift to bug-fix releases on the
1.16.x branch, and on adding new features on the main branch.

This release requires Python 3.11-3.13 and NumPy 1.25.2 or greater.

Highlights of this release

• Improved experimental support for the Python array API standard, including
  new support in scipy.signal, and additional support in scipy.stats and
  scipy.special. Improved support for JAX and Dask backends has been added,
  with notable support in scipy.cluster.hierarchy, many functions in
  scipy.special, and many of the trimmed statistics functions.
• scipy.optimize now uses the new Python implementation from the
  PRIMA package for COBYLA. The PRIMA implementation fixes many bugs
  in the old Fortran 77 implementation with a better performance on average.
• scipy.sparse.coo_array now supports n-D arrays with reshaping, arithmetic and
  reduction operations like sum/mean/min/max. No n-D indexing or
  scipy.sparse.random_array support yet.
• Updated guide and tools for migration from sparse matrices to sparse arrays.
• All functions in the scipy.linalg namespace that accept array arguments
  now support N-dimensional arrays to be processed as a batch.
• Two new scipy.signal functions, scipy.signal.firwin_2d and
  scipy.signal.closest_STFT_dual_window, for creation of a 2-D FIR filter and
  scipy.signal.ShortTimeFFT dual window calculation, respectively.
• A new class, scipy.spatial.transform.RigidTransform, provides functionality
  to convert between different representations of rigid transforms in 3-D
  space.

New features

scipy.io improvements

• scipy.io.savemat now provides informative warnings for invalid field names.
• scipy.io.mmread now provides a clearer error message when provided with
  a source file path that does not exist.
• scipy.io.wavfile.read can now read non-seekable files.

scipy.integrate improvements

• The error estimate of scipy.integrate.tanhsinh was improved.

scipy.interpolate improvements

• Batch support was added to scipy.interpolate.make_smoothing_spline.

scipy.linalg improvements

• All functions in the scipy.linalg namespace that accept array arguments
  now support N-dimensional arrays to be processed as a batch.
  See :ref:linalg_batch for details.
• scipy.linalg.sqrtm is rewritten in C and its performance is improved. It
  also tries harder to return real-valued results for real-valued inputs if
  possible. See the function docstring for more details. In this version the
  input argument disp and the optional output argument errest are
  deprecated and will be removed four versions later. Similarly, after
  changing the underlying algorithm to recursion, the blocksize keyword
  argument has no effect and will be removed two versions later.
• Wrappers for ?stevd, ?langb, ?sytri, ?hetri and
  ?gbcon were added to scipy.linalg.lapack.
• The default driver of scipy.linalg.eigh_tridiagonal was improved.
• scipy.linalg.solve can now estimate the reciprocal condition number and
  the matrix norm calculation is more efficient.

scipy.ndimage improvements

• A new function scipy.ndimage.vectorized_filter for generic filters that
  take advantage of a vectorized Python callable was added.
• scipy.ndimage.rotate has improved performance, especially on ARM platforms.

scipy.optimize improvements

• COBYLA was updated to use the new Python implementation from the PRIMA package.
  The PRIMA implementation fixes many bugs
  in the old Fortran 77 implementation. In addition, it results in fewer function evaluations on average,
  but it depends on the problem and for some
  problems it can result in more function evaluations or a less optimal
  result. For those cases the user can try modifying the initial and final
  trust region radii given by rhobeg and tol respectively. A larger
  rhobeg can help the algorithm take bigger steps initially, while a
  smaller tol can help it continue and find a better solution.
  For more information, see the PRIMA documentation.
• Several of the scipy.optimize.minimize methods, and the
  scipy.optimize.least_squares function, have been given a workers
  keyword. This allows parallelization of some calculations via a map-like
  callable, such as multiprocessing.Pool. These parallelization
  opportunities typically occur during numerical differentiation. This can
  greatly speed up minimization when the objective function is expensive to
  calculate.
• The lm method of scipy.optimize.least_squares can now accept
  3-point and cs for the jac keyword.
• The SLSQP Fortran 77 code was ported to C. When this method is used now the
  constraint multipliers are exposed to the user through the multiplier
  keyword of the returned scipy.optimize.OptimizeResult object.
• NNLS code has been corrected and rewritten in C to address the performance
  regression introduced in 1.15.x
• scipy.optimize.root now warns for invalid inner parameters when using the
  newton_krylov method
• The return value of minimization with method='L-BFGS-B' now has
  a faster hess_inv.todense() implementation. Time complexity has improved
  from cubic to quadratic.
• scipy.optimize.least_squares has a new callback argument that is applicable
  to the trf and dogbox methods. callback may be used to track
  optimization results at each step or to provide custom conditions for
  stopping.

scipy.signal improvements

• A new function scipy.signal.firwin_2d for the creation of a 2-D FIR Filter
  using the 1-D window method was added.
• scipy.signal.cspline1d_eval and scipy.signal.qspline1d_eval now provide
  an informative error on empty input rather than hitting the recursion limit.
• A new function scipy.signal.closest_STFT_dual_window to calculate the
  ~scipy.signal.ShortTimeFFT dual window of a given window closest to a
  desired dual window.
• A new classmethod scipy.signal.ShortTimeFFT.from_win_equals_dual to
  create a ~scipy.signal.ShortTimeFFT instance where the window and its dual
  are equal up to a scaling factor. It allows to create short-time Fourier
  transforms which are unitary mappings.
• The performance of scipy.signal.convolve2d was improved.

scipy.sparse improvements

• scipy.sparse.coo_array now supports n-D arrays using binary and reduction
  operations.
• Faster operations between two DIA arrays/matrices for: add, sub, multiply,
  matmul.
• scipy.sparse.csgraph.dijkstra shortest_path is more efficient.
• scipy.sparse.csgraph.yen has performance improvements.
• Support for lazy loading of sparse.csgraph and sparse.linalg was
  added.

scipy.spatial improvements

• A new class, scipy.spatial.transform.RigidTransform, provides functionality
  to convert between different representations of rigid transforms in 3-D
  space, its application to vectors and transform composition.
  It follows the same design approach as scipy.spatial.transform.Rotation.
• scipy.spatial.transform.Rotation now has an appropriate __repr__ method,
  and improved performance for its scipy.spatial.transform.Rotation.apply
  method.

scipy.stats improvements

• A new function scipy.stats.quantile, an array API compatible function for
  quantile estimation, was added.
• scipy.stats.make_distribution was extended to work with existing discrete
  distributions and to facilitate the creation of custom distributions in the
  new random variable infrastructure.
• A new distribution, scipy.stats.Binomial, was added.
• An equal_var keyword was added to scipy.stats.tukey_hsd (enables the
  Games-Howell test) and scipy.stats.f_oneway (enables Welch ANOVA).
• The moment calculation for scipy.stats.gennorm was improved.
• The scipy.stats.mode implementation was vectorized, for faster batch
  calculation.
• Support for axis, nan_policy, and keepdims keywords was added to
  scipy.stats.power_divergence, scipy.stats.chisquare,
  scipy.stats.pointbiserialr, scipy.stats.kendalltau,
  scipy.stats.weightedtau, scipy.stats.theilslopes,
  scipy.stats.siegelslopes, and scipy.stats.boxcox_llf.
• The performance of scipy.stats.special_ortho_group and scipy.stats.pearsonr
  was improved.

Array API Standard Support

Experimental support for array libraries other than NumPy has been added to
multiple submodules in recent ...

SciPy 1.15.3 Release Notes

SciPy 1.15.3 is a bug-fix release with no new features
compared to 1.15.2.

For the complete issue and PR lists see the raw release notes.

Authors

• Name (commits)
• aiudirog (1) +
• Nickolai Belakovski (1)
• Florian Bourgey (1) +
• Richard Strong Bowen (2) +
• Jake Bowhay (1)
• Dietrich Brunn (2)
• Evgeni Burovski (1)
• Lucas Colley (1)
• Ralf Gommers (1)
• Saarthak Gupta (1) +
• Matt Haberland (4)
• Chengyu Han (1) +
• Lukas Huber (1) +
• Nick ODell (2)
• Ilhan Polat (4)
• Tyler Reddy (52)
• Neil Schemenauer (1) +
• Dan Schult (1)
• sildater (1) +
• Gagandeep Singh (4)
• Albert Steppi (2)
• Matthias Urlichs (1) +
• David Varela (1) +
• ਗਗਨਦੀਪ ਸਿੰਘ (Gagandeep Singh) (3)

A total of 24 people contributed to this release.
People with a "+" by their names contributed a patch for the first time.
This list of names is automatically generated, and may not be fully complete.

SciPy 1.15.2 Release Notes

SciPy 1.15.2 is a bug-fix release with no new features
compared to 1.15.1. Free-threaded Python 3.13 wheels
for Linux ARM platform are available on PyPI starting with
this release.

Authors

• Name (commits)
• Peter Bell (1)
• Charles Bousseau (1) +
• Jake Bowhay (3)
• Matthew Brett (1)
• Ralf Gommers (3)
• Rohit Goswami (1)
• Matt Haberland (4)
• Parth Nobel (1) +
• Tyler Reddy (33)
• Daniel Schmitz (2)
• Dan Schult (5)
• Scott Shambaugh (2)
• Edgar Andrés Margffoy Tuay (1)
• Warren Weckesser (4)

A total of 14 people contributed to this release.
People with a "+" by their names contributed a patch for the first time.
This list of names is automatically generated, and may not be fully complete.

SciPy 1.15.1 Release Notes

SciPy 1.15.1 is a bug-fix release with no new features
compared to 1.15.0. Importantly, an issue with the
import of scipy.optimize breaking other packages
has been fixed.

Authors

• Name (commits)
• Ralf Gommers (3)
• Rohit Goswami (1)
• Matt Haberland (2)
• Tyler Reddy (7)
• Daniel Schmitz (1)

A total of 5 people contributed to this release.
People with a "+" by their names contributed a patch for the first time.
This list of names is automatically generated, and may not be fully complete.