geofile=fc_oogmsh.get_geo(3,3,'cylinderkey');
meshfile=fc_oogmsh.buildmesh3d(geofile,5);
Th=fc_simesh.siMesh(meshfile);
w={@(x,y,z) y.*cos(-(x.^2+y.^2)/10),@(x,y,z) -x.*cos(-(x.^2+y.^2)/10),@(x,y,z) z/5 };
W=Th.eval(w);
u=Th.eval(@(x,y,z) 3*x.^2-y.^3+z.^2+x.*y);
fc_tools.graphics.monitors.onGrid(1,2);
figure(1)
Th.plotmesh('d',1,'color','k', 'LineWidth',1.5)
hold on
Th.plotquiver(W,'LineWidth',1)
axis off;axis image
figure(2)
Th.plotmesh('d',1,'color','k', 'LineWidth',1.5)
hold on
Th.plotquiver(W,'colordata',u, 'scale',0.05)
axis off;axis image
colormap('jet');colorbar
fc_tools.graphics.monitors.autoGrid()
[fc-oogmsh] Input file : <fc-oogmsh>/geodir/3d/cylinderkey.geo
[fc-oogmsh] Mesh file <fc-oogmsh>/meshes/cylinderkey-5.msh [version 4.1] already exists.
-> Use "force" flag to rebuild if needed.
Mesh cylinderkey-5.msh is a 3-dimensional mesh
Force dimension to 3
