Matrix: flatten resampler and test on refined geometry

Ports to flatten_chained_connection so that it can construct matrices for multiply refined meshes.

TODO:

• At the moment this is just in matrix.py. Should it be added to QBXLayerPotentialSource.resampler? Maybe hidden under a flatten=False argument, so it doesn't change the previous behavior.