Normals for slightly complex surfaces MIGHT BE calculated wrongly
Working with testing new DualSPHysics beta capabilities, I have found that if one uses a slightly complex distribtution as:
The Cfg_Normals file shows wrongly calculated normals:
But svshapes command:
Does not show these artifacts - I am not sure what to trust, but going with svshapes for now. Also note that to actually perform the calculation it was necessary to use distanceh = 28 - since corner particles in this geometry is quite difficult to calculate for some reason.
Here I have not done anything crazy, simply just loaded a vtk file in, discretized using DualSPHysics normally and then done the normal calculation using xml directly.