Changeset 15990

Timestamp:

08/28/13 09:06:10 (12 years ago)

Author:

Mathieu Morlighem

Message:

NEW: added vertex list support to FlagElements.m

Location:

issm/trunk-jpl/src/m/geometry

Files:

• FlagElements.m (modified) (2 diffs)
• FlagElements.py (modified) (2 diffs)

issm/trunk-jpl/src/m/geometry/FlagElements.m

@@ -2 +2 @@
 %FLAGELEMENTS - flag the elements in an region
 %
-%   The region can be given with an exp file, a list of elements.
+%   The region can be given with an exp file, a list of elements or vertices
 %
 %   Usage: 
@@ -46 +46 @@
                 end
         elseif isfloat(region) | islogical(region),
-                if size(region,1)~=md.mesh.numberofelements,
+                if size(region,1)==md.mesh.numberofelements,
+                        flag=region;
+                elseif size(region,1)==md.mesh.numberofvertices,
+                        flag=logical(sum(region(md.mesh.elements)>0,2)==size(md.mesh.elements,2));
+                else
                         help FlagElements
                         error('Flaglist for region must be of same size as number of elements in model');
                 end
-                flag=region;
         else
                 error('Invalid region option');

issm/trunk-jpl/src/m/geometry/FlagElements.py

@@ -10 +10 @@
         FLAGELEMENTS - flag the elements in an region
 
-           The region can be given with an exp file, a list of elements.
+           The region can be given with an exp file, a list of elements or vertices
 
            Usage: 
@@ -54 +54 @@
 
         elif isinstance(region,numpy.ndarray) or isinstance(region,bool):
-                if not numpy.size(region,0)==md.mesh.numberofelements:
+                if numpy.size(region,0)==md.mesh.numberofelements:
+                        flag=region
+                else if numpy.size(region,0)==md.mesh.numberofvertices:
+                        raise TypeError("look at FlagElements.m and implement the same...")
+                else:
                         raise TypeError("Flaglist for region must be of same size as number of elements in model.")
-                flag=region
 
         else: