Changeset 2987

Timestamp:

02/08/10 17:16:45 (15 years ago)

Author:

Mathieu Morlighem

Message:

now use FlagElements to compute flag list

File:

• issm/trunk/src/m/classes/public/modelextract.m (modified) (1 diff)

issm/trunk/src/m/classes/public/modelextract.m

@@ -34 +34 @@
 
 %first
-if ischar(area),
-        if isempty(area),
-                flag_elem=zeros(md1.numberofelements,1);
-                invert=0;
-        elseif strcmpi(area,'all')
-                flag_elem=ones(md1.numberofelements,1);
-                invert=0;
-        else
-                %make sure that we actually don't want the elements outside the domain outline!
-                if strcmpi(area(1),'~'),
-                        area=area(2:length(area));
-                        invert=1;
-                else
-                        invert=0;
-                end
-
-                %ok, flag_elem elements
-                A=expread(area,1);
-                flag_elem=ContourToMesh(md1.elements(:,1:3),md1.x,md1.y,A,'element',1);
-        end
-        if invert, flag_elem=~flag_elem; end
-
-elseif isfloat(area),
-        if size(area,1)~=md1.numberofelements,
-                setelementstypeusage();
-                error('Flags must be of same size as number of elements in model');
-        end
-        flag_elem=area;
-else
-        error('Invalide option');
-end
+flag_elem=FlagElements(md,area);
 
 %kick out all elements with 3 dirichlets