The default normal vectors associated with tesselated surfaces are computed for each triangle face. Some purposes (like nice rendering) require per-vertex normals instead. These are typically calculated by averaging the face normals. This function offers two methods, first, an average normal weighted by triangle area, and second the normal as an average of the surrounding triangle faces. The second option is faster (because it uses the pre-computed unit-normals), but the first method generates nicer surface-normals, particularly at regions with sharp edges or small triangles.
This should do essentially the same thing as Dirk-Jan Kroon's 'Patch Normals' submission #24330, and my code is probably slower because it is not compiled. However, it's a bit simpler to read and (hopefully) to use.
There is a condition where this produces the wrong result. If the same vertex index appears more than once in a column of the face matrix, then only the first face normal is added to that index. To illustrate, if we have: `a = 1:5; a([1,1]) = a([1,1]) + [5,5];` then index 1 of `a` will only be incremented by 5, not 10. This could be fixed by replacing the three `vn(f(:,1),:)` lines with: `vn = accumarray([f(:),ones(size(f(:)));f(:),2*ones(size(f(:))); f(:),3*ones(size(f(:)))], [repmat(fn(:,1),[3,1]);repmat(fn(:,2),[3,1]);repmat(fn(:,3),[3,1])]);` I'm sure there is a much easier way to do this, but I was trying to avoid a loop.