git-svn-id: https://drops.fb12.tu-berlin.de/svn/bacc/trunk@89 26120e32-c555-405d...

diff --git a/src/refineQuad.m b/src/refineQuad.m
index d3bce2bedb69b6880a52e60d79ebfa16a12d9951..7ce800095debd2a96c35838976c31f64d3a6b14a 100644 (file)
--- a/src/refineQuad.m
+++ b/src/refineQuad.m
@@ -49,7 +49,7 @@ G_ref_f2s = repmat([1:c_loop]',1,4);
 clear elements coordinates neigh type
 
 % Flächeninhalt aufbauen
-updateS(1:c_loop);
+% updateS(1:c_loop);
 
 % figure(11)
 % plotShape(G_ref_C,G_ref_E,'s',G_ref_t);
@@ -151,6 +151,10 @@ while(~length(err))
   ref = find(G_ref_t>1);
   ref = reshape(ref,1,length(ref));
   if(isequal(ref,ref_old))
+        
+      
+      
+      
         figure(13)
         err_col = G_ref_t;
         if(length(err_col)<size(G_ref_E,1))
@@ -344,7 +348,7 @@ global G_ref_t;
     
         G_ref_f2s(ele,1:3)=c_ele+3:-1:c_ele+1;
 %         G_ref_f2s(ele,:)=[c_ele+3:-1:c_ele+1 ele];
-        updateS([ele,c_ele+1:c_ele+3]);
+%         updateS([ele,c_ele+1:c_ele+3]);
     elseif(G_ref_t(ele)==3)
         G_ref_C(c_coo+1,:) = (G_ref_C(el(1),:)+G_ref_C(el(4),:))/2;
         G_ref_C(c_coo+2,:) = (G_ref_C(el(2),:)+G_ref_C(el(3),:))/2;
@@ -355,7 +359,7 @@ global G_ref_t;
         G_ref_E(c_ele+1,:) = [el(1),el(2),c_coo+2,c_coo+1];
         G_ref_f2s(ele,[1 2])=c_ele+1;
 %         G_ref_f2s(ele,[3 4])=ele;
-        updateS([ele,c_ele+1]);
+%         updateS([ele,c_ele+1]);
      elseif(G_ref_t(ele)==4)
         G_ref_C(c_coo+1,:) = (G_ref_C(el(1),:)+G_ref_C(el(2),:))/2;
         G_ref_C(c_coo+2,:) = (G_ref_C(el(4),:)+G_ref_C(el(3),:))/2;
@@ -366,7 +370,7 @@ global G_ref_t;
         G_ref_E(c_ele+1,:) = [c_coo+1,el(2),el(3),c_coo+2];
         G_ref_f2s(ele,[2 3])=c_ele+1;
 %         G_ref_f2s(ele,[1 4])=ele;
-        updateS([ele,c_ele+1]);
+%         updateS([ele,c_ele+1]);
     end
     
     G_ref_t(ele) = 0;
@@ -501,21 +505,21 @@ end
 
 
 %% Berechnet die Flächeninhalte der Elemente
-function updateS(ele)
-%Globale Variabelen aufbauen
-global G_ref_E;
-global G_ref_C;
-% global G_ref_N;
-% global G_ref_f2s;
-% global G_ref_t;
-global G_ref_s;
-       for i = ele
-    % normalized Vector on every triangle
-    tri = G_ref_E(i,:);
-    a = (G_ref_C(tri(2),:)-G_ref_C(tri(1),:));
-    b = (G_ref_C(tri(4),:)-G_ref_C(tri(1),:));
-    N = cross(a',b');
-    N = N/norm(N);
-    G_ref_s(i) = sqrt(sum(N.*N));
-       end
-end
\ No newline at end of file
+% function updateS(ele)
+% %Globale Variabelen aufbauen
+% global G_ref_E;
+% global G_ref_C;
+% % global G_ref_N;
+% % global G_ref_f2s;
+% % global G_ref_t;
+% global G_ref_s;
+%      for i = ele
+%     % normalized Vector on every triangle
+%     tri = G_ref_E(i,:);
+%     a = (G_ref_C(tri(2),:)-G_ref_C(tri(1),:));
+%     b = (G_ref_C(tri(4),:)-G_ref_C(tri(1),:));
+%     N = cross(a',b');
+%     N = N/norm(N);
+%     G_ref_s(i) = sqrt(sum(N.*N));
+%      end
+% end
\ No newline at end of file