geofile=fc_simesh.get_geo(3,2,'demisphere5'); meshfile=fc_oogmsh.gmsh.buildmesh3ds(geofile,10); 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 }; W=Th.eval(w); u=Th.eval(@(x,y,z) 3*x.^2-y.^3+z.^2+x.*y); fc_tools.graphics.monitors.onGrid(2,2,'figures',1:3); figure(1) fc_siplt.plotmesh(Th,'d',1,'color','k', 'LineWidth',1.5) hold on;axis off;axis image fc_siplt.plotquiver(Th,W,'LineWidth',1) figure(2) fc_siplt.plotmesh(Th,'d',1,'color','k', 'LineWidth',1.5) hold on;axis off;axis image fc_siplt.plotquiver(Th,W,'LineWidth',1,'merge',false,'inlegend',true) legend('show') figure(3) fc_siplt.plotmesh(Th,'d',1,'color','k', 'LineWidth',1.5) hold on;axis off;axis image fc_siplt.plotquiver(Th,W,'colordata',u, 'scale',0.1) colormap('jet');colorbar
[fc-oogmsh] Input file : <fc-oogmsh>/geodir/3ds/demisphere5.geo [fc-oogmsh] Mesh file <fc-oogmsh>/meshes/demisphere5-10.msh [version 4] already exists. -> Use "force" flag to rebuild if needed. Mesh demisphere5-10.msh is a 3-dimensional mesh Force dimension to 3


