Changeset 1131

Timestamp:

06/26/09 13:48:40 (16 years ago)

Author:

Mathieu Morlighem

Message:

Added refining

File:

• issm/trunk/src/m/utils/Exp/expcoarsen.m (modified) (1 diff)

issm/trunk/src/m/utils/Exp/expcoarsen.m

@@ -37 +37 @@
 
         %stop if we have reached end of profile (always keep the last point)
-        while j<numpoints-1,,
+        while j<numpoints,
 
                 %See whether we keep this point or not
-                if sqrt((A(i).x(j)-A(i).x(j+1))^2+(A(i).y(j)-A(i).y(j+1))^2)<resolution
+                distance=sqrt((A(i).x(j)-A(i).x(j+1))^2+(A(i).y(j)-A(i).y(j+1))^2);
+                if distance<resolution & j<numpoints-1  %do not remove last point
                         A(i).x(j+1)=[];
                         A(i).y(j+1)=[];
                         numpoints=numpoints-1;
                 else
-                        %update current point
-                        j=j+1;
+                        division=floor(distance/resolution)+2;
+                        if division>2,
+                                x=linspace(A(i).x(j),A(i).x(j+1),division)';
+                                y=linspace(A(i).y(j),A(i).y(j+1),division)';
+                                A(i).x=[A(i).x(1:j);x(2:end-1); A(i).x(j+1:end)];
+                                A(i).y=[A(i).y(1:j);y(2:end-1); A(i).y(j+1:end)];
+
+                                %update current point
+                                j=j+1+division-2;
+                                numpoints=numpoints+division-2;
+                        else
+                                %update current point
+                                j=j+1;
+                        end
                 end
         end