direct-solver: use stage2_density_discr in tests.

This required modifying quite a few things:
* require getting expansion centers for stage2
* expansion radii for stage2
* tweaks to resampling in matrix.py