| 37 | In the most extreme case, as dx << 1 the function remains triangular, but as dx > 1, the function quickly drops to 0 (i.e. one cell removed from the origin). This is when the smoothing function would no longer serve its purpose (i.e. smoothing)... We could probably best avoid these issues by enforcing a larger kernel around the sink so that it most resembles an exponential, but then that would require changing the kernel of max refinement around a sink as well (which is currently set to 4 dx). We could also play with an exponential function that fell off faster. |