| 35 | | We see that in this case, smaller dx resembles a triangle, while dx~1 resembles a truer exponential. However, 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). |
| | 35 | We see that in this case, smaller dx resembles a triangle, while dx~1 resembles a truer exponential. However, 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). Or we could just leave it with an error message in the code when dx > 1 that this may not be the best smoothing function. |
