meshfile=fc_oogmsh.gmsh.buildmesh3ds('demisphere5',10);
  Th=fc_simesh.siMesh(meshfile);

  u=Th.eval(@(x,y,z) 3*x.^2-y.^3+z.^2+x.*y);
  fc_tools.graphics.monitors.onGrid(1,2);
  figure(1)
  fc_siplt.plot(Th,u,'FaceColor','interp', 'EdgeColor','none','FaceAlpha',0.4)
  hold on;axis off;axis equal;colorbar
  fc_siplt.plotiso(Th,u,'LineWidth',1.5)

  figure(2)
  fc_siplt.plot(Th,u,'labels',[1,11],'FaceColor','interp', 'EdgeColor','none')
  hold on;axis off;axis equal;colorbar;
  fc_siplt.plotiso(Th,u,'labels',[1,11],'LineWidth',1.,'color','w')
  fc_siplt.plot(Th,u,'labels',[10,12],'FaceColor','interp', 'EdgeColor','none','FaceAlpha',0.4)
  fc_siplt.plotiso(Th,u,'labels',[10,12],'LineWidth',1.5)
  fc_siplt.plotmesh(Th,'d',1,'Color','k')

[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

Published with MATLAB® R2022a