test_helix.m

来自「有限元编程做网格划分的」· M 代码 · 共 7 行

[fens,gcells] = t4block(2.5*pi,0.95, 0.425, 50, 4, 2);
fens = transform_apply(fens,@(x, data) (x+ [0, 2.5, 0]), []);
climbPerRevolution= 2.3;
fens = transform_2_helix(fens,climbPerRevolution);
bg=mesh_bdry(gcells);
drawmesh({fens,bg},'gcells','facecolor','red')