Add the advection_matrix() function and a test for it.