Changeset 3016

Timestamp:

02/10/10 17:33:48 (15 years ago)

Author:

Eric.Larour

Message:

partitioner: wrap chaco, scotch, metis
adjacency: faster computation of adjacency.

Location:

issm/trunk/src/m/partition

Files:

• adjacency.m (added)
• flagedges.m (modified) (3 diffs)
• partitioner.m (added)

issm/trunk/src/m/partition/flagedges.m

@@ -1 @@
-function flags=flagedges(partition,adjacencymatrix);
+function [flags middles]=flagedges(x,y,partition,adjacencymatrix);
 %FLAGEDGES return a vector of boolean flags that indicate whether a node is on the edge of a partition.
 %          needs partition vector and adjacencymatrix in inputs.
 %
 %
-%         Usage: flags=flagedges(partition,adjacencymatrix);
+%         Usage: [flags middlegrids]=flagedges(x,y,partition,adjacencymatrix);
 %
 
@@ -12 +12 @@
 %initialize flags
 flags=zeros(numberofgrids,1);
+
+middles=zeros(0,2);
 
 %go through grids, and figure out which grid is connected to a grid belonging to another partition, flag it.
@@ -38 +40 @@
                 grids_on_edge=grids(index);
                 flags(grids_on_edge)=1;
-                %flags(i)=1; do not flag i, otherwise, edges are doubled, results in ugly plots.
+                flags(i)=1; 
+                for j=1:length(grids_on_edge),
+                        middles(end+1,:)=[(x(grids_on_edge(j))+x(i))/2 (y(grids_on_edge(j))+y(i))/2];
+                end
         end
 end