Feature: Add matrix testing for different Python versions to CI workflow

[yardasol]

Background and motivation

As an improvement to our automated CI workflow, it would be nice to make our tests run in a matrix of different Python versions.

Implementation details

We use GH Actions to do this

Potential snags

This would increase the energy usage of our CI workflow :(