Changeset 3039

Timestamp:

02/16/10 11:02:14 (15 years ago)

Author:

Eric.Larour

Message:

Added partition wrapper to John's partitioners.
Added flagedges routine.

Location:

issm/trunk/src/m/partition

Files:

• flagedges.m (modified) (2 diffs)
• partitioner.m (modified) (3 diffs)

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

@@ -20 +20 @@
         %some feedback
         if i==1,
-                fprintf('    progress report: %5.2f %s',i/numberofgrids*100,'%');
+                fprintf('    flagedges progress report: %5.2f %s',i/numberofgrids*100,'%');
         elseif mod(i,1000),
                 fprintf('\b\b\b\b\b\b\b');
@@ -31 +31 @@
         end
         
-        grids=adjacencymatrix(i,:); grids=find(grids);
+        grids=adjacencymatrix(:,i); grids=find(grids);
         
         part_i=partition(i);

issm/trunk/src/m/partition/partitioner.m

@@ -7 +7 @@
 %   npart: number of partitions.
 %   weighting: 'on' or 'off': default off
+%   section:  1 by defaults(1=bisection, 2=quadrisection, 3=octasection)
 %   
 %   Usage:
@@ -19 +20 @@
 options=addfielddefault(options,'package','chaco');
 options=addfielddefault(options,'npart',10);
-options=addfielddefault(options,'weighting','off');
+options=addfielddefault(options,'weighting','on');
+options=addfielddefault(options,'section',1);
 
 %get package: 
@@ -30 +32 @@
         method=[1 1 0 0 1 1 50 0 .001 7654321];
         method(1)=3;    %  global method (3=inertial (geometric))
-        method(6)=1;    %  ndims (1=bisection, 2=quadrisection, 3=octasection)
         method(3)=1;    %  vertex weights (0=off, 1=on)
+        
+        %specify bisection
+        method(6)=getfieldvalue(options,'section');%  ndims (1=bisection, 2=quadrisection, 3=octasection)
 
         %are we using weights?