Hello,
It could be good to add vertex_normal computation directly when one reads a gifti mesh, so in the io.py module. It requires functions of trimesh to compute if from face_normal and avoid redundancies with functions in curvature.py.
Not very long, I can try to do that
Julien