meshfile=fc_oogmsh.gmsh.buildmesh3d('cylinder3holes',10);
  Th=fc_simesh.siMesh(meshfile);
  u=@(x,y,z) cos(2*x-y-z).*sin(x-2*y+z);
  U=Th.eval(u);

  fc_tools.graphics.monitors.onGrid(2,2,'figures',1:3);

  figure(1)
  fc_siplt.plot(Th,U);
  axis image

  figure(2)
  fc_siplt.plot(Th,U,'d',2);
  axis image

  figure(3)
  fc_siplt.plot(Th,U,'d',1,'Linewidth',3);
  axis image
  colormap('jet')
  colorbar

[fc-oogmsh] Input file : <fc-oogmsh>/geodir/3d/cylinder3holes.geo
[fc-oogmsh] Mesh file <fc-oogmsh>/meshes/cylinder3holes-10.msh [version 4.1] already exists.
  -> Use "force" flag to rebuild if needed.
  Mesh cylinder3holes-10.msh is a 3-dimensional mesh
    Force dimension to 3

Published with MATLAB® R2022a