ENH: Make output of get_atlas_info conditional on system_info.verbosity
Need to INCREF dtype again since PyArray_NewFromDescr steals a reference
Remove excessive spacing
Merge pull request #4494 from embray/astropy-1987

BUG: Adds missing Py_INCREF in array_boolean_subscript
TST: add some more tests

test for gh-4494
test median returns array scalars and works with object arrays
BUG: Fix linspace for use with physical quantities

The fix for issue 3504 led to errors when using linspace with the quantities package. Multiplying with 1 instead of adding 0 prevents issues when using physical quantities.
Merge pull request #4503 from rproepp/3504-quantities

Fix linspace for use with physical quantities
Merge pull request #4507 from juliantaylor/test-add

TST: add some more tests
ENH: optimize srcstride0 copy with O3

enables autovectorization with gcc and improves performance by about 25% for cached data.
Merge pull request #4387 from abergeron/get_info

ENH: Make output of get_atlas_info conditional on system_info.verbosity
DOC: Add the change to disutils in the release notes.
Merge pull request #4514 from abergeron/get_info

Add the change to disutils in the release notes.
FIX: missing asarray import in numpy.libs.utils
TST: add unittest for byte_bounds
Merge pull request #4515 from ogrisel/undefined-asarray

FIX: missing asarray import in numpy.libs.utils
BUG: catch overflow of nominator in datetime conversions

closes gh-4486
Merge pull request #4518 from juliantaylor/datetime-overflow

BUG: catch overflow of nominator in datetime conversions
FIX: items/keys iterators in 3k, list wrapping
DOC: Link convolve with polymul

Matlab uses `conv` for both convolution and polynomial multiplication. Clarifying that numpy has functions for each.
ENH, MAINT: Check that covariance is positive-semidefinite.

numpy.random.multivariate_normal uses SVD to obtain the covariance
square root, but the result will be incorrect if the covariance is not
positive-semidefinite and no warning will be given. However, a check can
be made even while using the SVD by checking if the rows of u.T and v
have the same sign, where u and v are given by u, d, v = svd(a).

A deprecation warning is also fixed. When the default size was used an
empty list was passed to multiply.reduce, which resulted in a default
float 1. for a shape index. That is no longer legal. The problem is
fixed by using an array reshape, which also avoids potential overflow in
the original computation.

Closes #4489.