Investigate Array API usage in scanpy codebase

[ilan-gold]

What kind of feature would you like to request?

Additional function parameters / changed functionality / changed defaults?

Please describe your wishes

The array api is somewhat mature to the point of being using in scipy for interoperability between dense CPU/GPU functions (see https://docs.scipy.org/doc/scipy/dev/api-dev/array_api.html). It is coming with the new sparse format that may come into existence: pydata/sparse#618

So it would be good to start experimenting with this in scanpy