plot_elementstype

PURPOSE ^

PLOT_ELEMENTSTYPE - plot elements type

SYNOPSIS ^

function plot_elementstype(md,options_structure,width,i);

DESCRIPTION ^

PLOT_ELEMENTSTYPE - plot elements type

   Usage:
      plot_elementstype(md,options_structure,width,i);

   See also: PLOTMODEL

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function plot_elementstype(md,options_structure,width,i);
0002 %PLOT_ELEMENTSTYPE - plot elements type
0003 %
0004 %   Usage:
0005 %      plot_elementstype(md,options_structure,width,i);
0006 %
0007 %   See also: PLOTMODEL
0008 
0009 %process data and model
0010 [x y z elements is2d]=processmesh(md,options_structure);
0011 [data1 isongrid]=processdata(md,md.elements_type(:,1),options_structure);
0012 [data2 isongrid]=processdata(md,md.elements_type(:,2),options_structure);
0013 data=[data1 data2];
0014 
0015 %edgecolor?
0016 if ~isnan(options_structure.edgecolor),
0017     edgecolor=options_structure.edgecolor;
0018 else
0019     edgecolor='none';
0020 end
0021 
0022 %plot
0023 subplot(width,width,i);
0024 
0025 if is2d
0026     %Hutter elements
0027     posH=find(data(:,1)==hutterenum);
0028     A=elements(posH,1); B=elements(posH,2); C=elements(posH,3); 
0029     p1=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',hutterenum,'FaceColor','flat','EdgeColor',edgecolor);
0030     %MacAyeal element
0031     posH=find(data(:,1)==macayealenum);
0032     A=elements(posH,1); B=elements(posH,2); C=elements(posH,3); 
0033     p2=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',macayealenum,'FaceColor','flat','EdgeColor',edgecolor);
0034     %Pattyn element
0035     posH=find(data(:,1)==pattynenum);
0036     A=elements(posH,1); B=elements(posH,2); C=elements(posH,3); 
0037     p3=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',pattynenum,'FaceColor','flat','EdgeColor',edgecolor);
0038     %Stokes elements
0039     alpha=0.35;
0040     posS=find(data(:,2)==stokesenum);
0041     if ~isempty(posS)
0042         A=elements(posS,1); B=elements(posS,2); C=elements(posS,3);
0043     %    p4=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData',stokesenum,'FaceColor','flat','EdgeColor',edgecolor,'EdgeAlpha',alpha,'FaceAlpha',alpha);
0044     %    legend([p1 p2 p3 p4],'Hutter''s elements','MacAyeal''s elements','Pattyn''s elements','Stokes''s elements');
0045     else
0046         legend([p1 p2 p3],'Hutter''s elements','MacAyeal''s elements','Pattyn''s elements');
0047     end
0048 else
0049     %Hutter elements
0050     posH=find(data(:,1)==hutterenum);
0051     A=elements(posH,1); B=elements(posH,2); C=elements(posH,3); D=elements(posH,4); E=elements(posH,5); F=elements(posH,6);
0052     p1=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', hutterenum,'FaceColor','flat','EdgeColor',edgecolor);
0053     patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', hutterenum,'FaceColor','flat','EdgeColor',edgecolor);
0054     patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', hutterenum,'FaceColor','flat','EdgeColor',edgecolor);
0055     patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', hutterenum,'FaceColor','flat','EdgeColor',edgecolor);
0056     patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', hutterenum,'FaceColor','flat','EdgeColor',edgecolor);
0057     %MacAyeal elements
0058     posM=find(data(:,1)==macayealenum);
0059     A=elements(posM,1); B=elements(posM,2); C=elements(posM,3); D=elements(posM,4); E=elements(posM,5); F=elements(posM,6);
0060     p2=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', macayealenum,'FaceColor','flat','EdgeColor',edgecolor);
0061     patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', macayealenum,'FaceColor','flat','EdgeColor',edgecolor);
0062     patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', macayealenum,'FaceColor','flat','EdgeColor',edgecolor);
0063     patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', macayealenum,'FaceColor','flat','EdgeColor',edgecolor);
0064     patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', macayealenum,'FaceColor','flat','EdgeColor',edgecolor);
0065     %Pattyn elements
0066     posP=find(data(:,1)==pattynenum);
0067     A=elements(posP,1); B=elements(posP,2); C=elements(posP,3); D=elements(posP,4); E=elements(posP,5); F=elements(posP,6);
0068     p3=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', pattynenum,'FaceColor','flat','EdgeColor',edgecolor);
0069     patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', pattynenum,'FaceColor','flat','EdgeColor',edgecolor);
0070     patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', pattynenum,'FaceColor','flat','EdgeColor',edgecolor);
0071     patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', pattynenum,'FaceColor','flat','EdgeColor',edgecolor);
0072     patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', pattynenum,'FaceColor','flat','EdgeColor',edgecolor);
0073     %Stokes elements
0074     alpha=0.35;
0075     posS=find(data(:,2)==stokesenum);
0076     if ~isempty(posS)
0077         A=elements(posS,1); B=elements(posS,2); C=elements(posS,3); D=elements(posS,4); E=elements(posS,5); F=elements(posS,6);
0078     %    p4=patch( 'Faces', [A B C], 'Vertices', [x y z],'CData', stokesenum,'FaceColor','flat','EdgeColor',edgecolor,'EdgeAlpha',alpha,'FaceAlpha',alpha);
0079     %    patch( 'Faces', [D E F], 'Vertices', [x y z],'CData', stokesenum,'FaceColor','flat','EdgeColor',edgecolor,'EdgeAlpha',alpha,'FaceAlpha',alpha);
0080     %    patch( 'Faces', [A B E D], 'Vertices', [x y z],'CData', stokesenum,'FaceColor','flat','EdgeColor',edgecolor,'EdgeAlpha',alpha,'FaceAlpha',alpha);
0081     %    patch( 'Faces', [B E F C ], 'Vertices', [x y z],'CData', stokesenum,'FaceColor','flat','EdgeColor',edgecolor,'EdgeAlpha',alpha,'FaceAlpha',alpha);
0082     %    patch( 'Faces', [C A D F ], 'Vertices', [x y z],'CData', stokesenum,'FaceColor','flat','EdgeColor',edgecolor,'EdgeAlpha',alpha,'FaceAlpha',alpha);
0083     %    legend([p1 p2 p3 p4],'Hutter''s elements','MacAyeal''s elements','Pattyn''s elements','Stokes''s elements');
0084     else
0085         legend([p1 p2 p3],'Hutter''s elements','MacAyeal''s elements','Pattyn''s elements');
0086     end
0087 end
0088 
0089 %apply options
0090 if isnan(options_structure.title)
0091     options_structure.title='Elements type';
0092 end 
0093 if isnan(options_structure.colorbar)
0094     options_structure.colorbar=0;
0095 end
0096 applyoptions(md,[],options_structure);

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003