Looking through Schmidt et al. (2003) I found what I believe to be 3 small inaccuracies in curvature.m:
As far as I can tell the gradient function is not producing the desired fx and fy for the Evans method, using only 2 neighbours rather than the desired 6. Additionally, gradient calculates fx in the opposite direction.
It seems like + 2 * fx .* fy .* fxy should be - 2 * fx .* fy .* fxy in the tangential curvature.
It seems like + 2 * fx .* fy .* fxy should be - 2 * fx .* fy .* fxy in the planform curvature.
I am not an expert, so please let me know if I am mistaken.
Looking through Schmidt et al. (2003) I found what I believe to be 3 small inaccuracies in
curvature.m:gradientfunction is not producing the desiredfxandfyfor the Evans method, using only 2 neighbours rather than the desired 6. Additionally,gradientcalculatesfxin the opposite direction.+ 2 * fx .* fy .* fxyshould be- 2 * fx .* fy .* fxyin the tangential curvature.+ 2 * fx .* fy .* fxyshould be- 2 * fx .* fy .* fxyin the planform curvature.I am not an expert, so please let me know if I am mistaken.