PLOT_MAGNETIZATION

 plot magnetization from solving Bloch-Torrey equation
 
 Input:
     1. mymesh is a structure with 10 elements:
         Nnode
         Nele
         Nface
         Pts_cmpt_reorder
         Ele_cmpt_reorder
         Pts_ind
         Pts_boundary_reorder
         Fac_boundary_reorder
         Nboundary
         Ncmpt
     2. YOUT
     3. OUT_cmpts_index
     4. ECS_cmpts_index
     5. IN_cmpts_index
     
 Output:
     figures for magnetization of inner, (outer, ECS) compartments

0001 function PLOT_MAGNETIZATION(mymesh,YOUT,OUT_cmpts_index,ECS_cmpts_index,IN_cmpts_index)
0002 
0003 % plot magnetization from solving Bloch-Torrey equation
0004 %
0005 % Input:
0006 %     1. mymesh is a structure with 10 elements:
0007 %         Nnode
0008 %         Nele
0009 %         Nface
0010 %         Pts_cmpt_reorder
0011 %         Ele_cmpt_reorder
0012 %         Pts_ind
0013 %         Pts_boundary_reorder
0014 %         Fac_boundary_reorder
0015 %         Nboundary
0016 %         Ncmpt
0017 %     2. YOUT
0018 %     3. OUT_cmpts_index
0019 %     4. ECS_cmpts_index
0020 %     5. IN_cmpts_index
0021 %
0022 % Output:
0023 %     figures for magnetization of inner, (outer, ECS) compartments
0024 
0025 xmin=0;
0026 xmax=0;
0027 ymin=0;
0028 ymax=0;
0029 zmin=0;
0030 zmax=0;
0031 for icmpt = 1:mymesh.Ncmpt
0032     xx=max(mymesh.Pts_cmpt_reorder{icmpt}(1,:));
0033     yx=max(mymesh.Pts_cmpt_reorder{icmpt}(2,:));
0034     zx=max(mymesh.Pts_cmpt_reorder{icmpt}(3,:));
0035     xn=min(mymesh.Pts_cmpt_reorder{icmpt}(1,:));
0036     yn=min(mymesh.Pts_cmpt_reorder{icmpt}(2,:));
0037     zn=min(mymesh.Pts_cmpt_reorder{icmpt}(3,:));
0038     xmin = min(xmin,xn);
0039     xmax = max(xmax,xx);
0040     ymin = min(ymin,yn);
0041     ymax = max(ymax,yx);
0042     zmin = min(zmin,zn);
0043     zmax = max(zmax,zx);
0044 end
0045 
0046 if (~isempty(OUT_cmpts_index))
0047     figure; 
0048     hold on;
0049     cmptvec = OUT_cmpts_index;
0050     for ict = 1:length(cmptvec)
0051         icmpt = cmptvec(ict);
0052         Fac = [];
0053         for iboundary = 1:mymesh.Nboundary
0054             Fac = [Fac,mymesh.Fac_boundary_reorder{icmpt}{iboundary}];
0055         end
0056         h = trisurf(Fac',mymesh.Pts_cmpt_reorder{icmpt}(1,:),mymesh.Pts_cmpt_reorder{icmpt}(2,:),...
0057             mymesh.Pts_cmpt_reorder{icmpt}(3,:),real(YOUT{end}{end}{icmpt}(:,end)));
0058         set(h,'facealpha',0.6);
0059         axis equal;
0060         axis([xmin,xmax,ymin,ymax,zmin,zmax]); colorbar('southoutside');
0061         view(3);
0062         title(['Magnetization Inner cmpts: ',num2str(OUT_cmpts_index)]);
0063     end
0064 end
0065 
0066 if (~isempty(IN_cmpts_index))
0067     figure; 
0068     hold on;
0069     cmptvec = IN_cmpts_index;
0070     for ict = 1:length(cmptvec)
0071         icmpt = cmptvec(ict);
0072         Fac = [];
0073         for iboundary = 1:mymesh.Nboundary
0074             Fac = [Fac,mymesh.Fac_boundary_reorder{icmpt}{iboundary}];
0075         end
0076         h = trisurf(Fac',mymesh.Pts_cmpt_reorder{icmpt}(1,:),mymesh.Pts_cmpt_reorder{icmpt}(2,:),...
0077             mymesh.Pts_cmpt_reorder{icmpt}(3,:),real(YOUT{end}{end}{icmpt}(:,end)));
0078         set(h,'facealpha',0.6);
0079         axis equal;
0080         axis([xmin,xmax,ymin,ymax,zmin,zmax]); colorbar('southoutside');
0081         view(3);
0082         title(['Magnetization Outer cmpts: ',num2str([IN_cmpts_index])]);
0083     end
0084 end
0085 
0086 if (~isempty(ECS_cmpts_index))
0087     figure; 
0088     hold on;
0089     cmptvec = ECS_cmpts_index;
0090     for ict = 1:length(cmptvec)
0091         icmpt = cmptvec(ict);
0092         Fac = [];
0093         for iboundary = 1:mymesh.Nboundary
0094             Fac = [Fac,mymesh.Fac_boundary_reorder{icmpt}{iboundary}];
0095         end
0096         h = trisurf(Fac',mymesh.Pts_cmpt_reorder{icmpt}(1,:),mymesh.Pts_cmpt_reorder{icmpt}(2,:),...
0097             mymesh.Pts_cmpt_reorder{icmpt}(3,:),real(YOUT{end}{end}{icmpt}(:,end)));
0098         set(h,'facealpha',0.6);
0099         axis equal;
0100         axis([xmin,xmax,ymin,ymax,zmin,zmax]); colorbar('southoutside');
0101         view(3);
0102         title(['Magnetization ECS cmpt: ',num2str([ECS_cmpts_index])]);
0103     end
0104 end