Commit 08594b48 authored by Christophe Geuzaine's avatar Christophe Geuzaine

getCurvature

parent 898ff28d
......@@ -25,7 +25,7 @@ for e in ent:
# singularities (edges/points)
normals = gmsh.model.getNormal(surf, param)
# get surface curvature
curvMax, curvMin, _, _ = gmsh.model.getPrincipalCurvatures(surf, param)
curv = gmsh.model.getCurvature(2, surf, param)
for i in range(0,len(coord),3):
nn.append(coord[i])
nn.append(coord[i+1])
......@@ -36,7 +36,7 @@ for e in ent:
cc.append(coord[i])
cc.append(coord[i+1])
cc.append(coord[i+2])
cc.append(curvMax[i/3])
cc.append(curv[i/3])
t = gmsh.view.add("normals")
gmsh.view.addListData(t, "VP", len(nn)/6, nn)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment