u=@(x,y,z) 2*cos(x).*sin(y+z); [q3,me3,toGlobal3]=fc_meshtools.simplicial.getMeshArrays('3D',3); U3=fc_meshtools.eval(u,q3); [q1,me1,toGlobal1]=fc_meshtools.simplicial.getMeshArrays('3D',1); figure(1) P=fc_tools.graphics.PlaneCoefs([0 0 1],[1 -1 1]); fc_graphics4mesh.slicemesh(q3,me3,P) hold on fc_graphics4mesh.plotmesh(q1,me1,'color','k') fc_graphics4mesh.sliceiso(q3,me3,U3,P) axis equal;axis image view(-11,15) figure(2) P=fc_tools.graphics.PlaneCoefs([0 0 1],[1 -1 1]); fc_graphics4mesh.slicemesh(q3,me3,P) hold on fc_graphics4mesh.plotmesh(q1,me1,'color','k') fc_graphics4mesh.sliceiso(q3,me3,U3,P,'Linewidth',2, 'isocolorbar',true,'LineWidth',2,'format','%.3f'); P=fc_tools.graphics.PlaneCoefs([0 0 1],[1 -1 0]); fc_graphics4mesh.slice(q3,me3,U3,P) fc_graphics4mesh.sliceiso(q3,me3,U3,P,'color','w','Linewidth',2); axis equal;axis image colorbar('Location','westoutside') caxis([min(U3),max(U3)]) view(-11,15) figure(3) P=[fc_tools.graphics.PlaneCoefs([0 0 1],[1 0 0]); ... fc_tools.graphics.PlaneCoefs([0 0 1],[0 1 0]); ... fc_tools.graphics.PlaneCoefs([0 0 1],[0 0 1])]; fc_graphics4mesh.slicemesh(q3,me3,P) hold on fc_graphics4mesh.plotmesh(q1,me1,'color','k') fc_graphics4mesh.sliceiso(q3,me3,U3,P,'isocolorbar',true,'LineWidth',2, 'format','%.3f'); view(3);axis equal;axis image figure(4) P=[fc_tools.graphics.PlaneCoefs([0 0 1],[0 1 0]); ... fc_tools.graphics.PlaneCoefs([0 0 1],[0 -1 1])]; fc_graphics4mesh.slice(q3,me3,U3,P) hold on fc_graphics4mesh.plotmesh(q1,me1,'color','k') fc_graphics4mesh.sliceiso(q3,me3,U3,P,'Color','w','LineWidth',2); caxis([min(U3),max(U3)]) view(3);axis equal;axis image colorbar
