Curvature Estimationl On triangle mesh

If you use it for academic purposes please cite :
1. "Estimating Curvatures and Their Derivatives on Triangle Meshes" by Szymon Rusinkiewicz (2004)
2. Y. Ben Shabat, A. Fischer, "Design of Adaptive Porous Micro-structures using curvature analysis for Additive Manufacturing" the 25th CIRP Design conference. 2015, Haifa, Israel.
This code was implemented according to [1].

GetCurvatures computes the curvature tensor and the principal curvatures at
%each vertex of a mesh given in a face vertex data structure
%INPUT:
-FV -struct - Triangle mesh face vertex data structure (containing FV.face and
FV.Vertices)
-toggleDerivatives - scalar 1 or 0 indicating whether or not to calculate curvature derivatives
%OUTPUT:
-PrincipalCurvatures - 2XN matrix (where N is the number of vertices containing the principal curvatures k1 and k2 at each vertex
-PrincipalDir1 - NX3 matrix containing the direction of the k1 principal
curvature.
-PrincipalDir2 - NX3 matrix containing the direction of the k2 principal curvature.
-FaceCMatrix - 4XM matrix (where M is the number of faces) containing the 4
%coefficients of the curvature tensr of each face
VertexCMatrix- 4XN matrix (where M is the number of faces) containing the 4 coefficients of the curvature tensor of each tensor.
-Cmagnitude - NX1 matrix containing the square sum of the curvature tensor coefficients at each vertex (invariant scalar indicating the change of curvature)