This code (as well as other analogical codes that I saw on fileexchange) will take the triangulation of my surface as an input parameter (i.e., the mesh). But I was wondering maybe there is a simpler way to calculate the curvature of the surface since I have a smooth surface that has been already fitted into my triangulated point cloud?
It is more complicated for 3d surface. Indeed you would need to calculate normals on triangular mesh (i.e., first, for triangles, then for vertices). I do not think surfnorm works for triangular mesh (isosurface). Then the problem will consist of calculating pricncipal curvatures (minimum and maximum curvatures out of all directions) from the second fundamental form.