Struggling with my measure of diffusivity. Need to build a unit test for tridiagonal matrix operations on multi-dimensional tensors (e.g., finite differencing across one-dimension).