geofile=fc_simesh.get_geo(2,2,'condenser11');
  meshfile=fc_oogmsh.buildmesh2d(geofile,25);
  Th=fc_simesh.siMesh(meshfile);

  u=@(x,y) cos(pi*x.^2).*cos(pi*y.^2);
  U=Th.eval(u);
  w={@(x,y) y.*cos(-(x.^2+y.^2)/10),@(x,y) -x.*cos(-(x.^2+y.^2)/10)};
  W=Th.eval(w);

  figure(1)
  Th.plotmesh('d',1,'color','k', 'LineWidth',1.5)
  hold on;axis off;axis image
  Th.plotquiver(W)

  figure(2)
  Th.plotmesh('d',1,'color','k', 'LineWidth',1.5)
  hold on;axis off;axis image
  Th.plotquiver(W,'LineWidth',2,'merge',false)

  figure(3)
  Th.plotmesh('d',1,'color','k', 'LineWidth',1.5)
  hold on;axis off;axis image
  Th.plotquiver(W,'colordata',U,'labels',[2:2:8,20])
  caxis([min(U) max(U)])
  colormap('jet');colorbar

  figure(4)
  Th.plotmesh('d',1,'color','k', 'LineWidth',1.5)
  hold on;axis off;axis image
  Th.plotquiver(W,'colordata',U, 'scale',0.2)
  colormap('jet');colorbar
  fc_tools.graphics.monitors.autoGrid()

[fc-oogmsh] Input file : <fc-oogmsh>/geodir/2d/condenser11.geo
[fc-oogmsh] Mesh file <fc-oogmsh>/meshes/condenser11-25.msh [version 4.1] already exists.
  -> Use "force" flag to rebuild if needed.

Published with MATLAB® R2022a