Changeset 5861

Timestamp:

09/17/10 08:28:07 (15 years ago)

Author:

Eric.Larour

Message:

added prescribedVertices options

File:

• issm/trunk/src/m/classes/public/bamg.m (modified) (2 diffs)

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

@@ -10 +10 @@
 %   - hminVertices: minimum edge length for each vertex (mesh)
 %   - hmaxVertices: maximum edge length for each vertex (mesh)
+%   - prescribedVertices:  vertices that should be kept by bamg. [x,y,marker]; marker is optional
 %
 %   - anisomax    : maximum ration between the smallest and largest edges (default is 10^30)
@@ -223 +224 @@
                         count=count+nods;
                 end
+        end
+
+        %Deal with vertices that need to be kept by mesher
+        if exist(options,'prescribedVertices'),
+
+                %recover xvertices yvertices
+                prescribedVertices=getfieldvalue(options,'prescribedVertices');
+                xvertices=prescribedVertices(:,1);
+                yvertices=prescribedVertices(:,2);
+                if(size(prescribedVertices,2)==3),
+                        markers=prescribedVertices(:,3);
+                else 
+                        markers=[];
+                end
+        
+                %only keep those inside
+                flags=ContourToNodes(xvertices,yvertices,domain(1),0);
+                pos=find(flags);
+                xvertices=xvertices(pos);
+                yvertices=yvertices(pos);
+                if ~isempty(markers),
+                        markers=markers(pos);
+                end
+
+                %Add all points to bamg_geometry
+                nods=length(xvertices);
+                if isempty(markers),
+                        markers=4*ones(nods,1);
+                end
+
+                bamg_geometry.Vertices=[bamg_geometry.Vertices; [xvertices yvertices markers]]; 
+
         end