ENH: add `convolve`

This is not in the spec and is unlikely to get into it.
Is however available in numpy, cupy and jax; for pytorch it can be emulated as in https://github.com/pytorch/pytorch/blob/main/torch/_numpy/_funcs_impl.py#L612
An "array agnostic" implementation is unlikely to be usable however. 

