本文共 1212 字，大约阅读时间需要 4 分钟。

the original tutor is like:

there are some points I wanna mention:

1. the discrete Gaussian curvature formula it listed is wrong. How? listen to me explain.

According to the paper "Discrete Differential-Geometry Operators for Triangulated 2-Manifolds" section 4.2, it says

In addition, the code also told us about that:

#include 
   
    #include 
    
     #include 
     
      #include 
      
       #include 
       
        #include 
        
         #include "tutorial_shared_path.h"int main(int argc, char *argv[]){ using namespace Eigen; using namespace std; MatrixXd V; MatrixXi F; igl::readOFF(TUTORIAL_SHARED_PATH "/bumpy.off",V,F); VectorXd K; // Compute integral of Gaussian curvature igl::gaussian_curvature(V,F,K); // Compute mass matrix SparseMatrix
         
           M,Minv; igl::massmatrix(V,F,igl::MASSMATRIX_TYPE_DEFAULT,M); igl::invert_diag(M,Minv); // Divide by area to get integral average, note this line K = (Minv*K).eval(); // Compute pseudocolor MatrixXd C; igl::jet(K,true,C); // Plot the mesh with pseudocolors igl::viewer::Viewer viewer; viewer.data.set_mesh(V, F); viewer.data.set_colors(C); viewer.launch();}