Changeset 1059 for trunk/src/filter_tps.m

Timestamp:

Nov 29, 2018, 5:44:44 PM (6 years ago)

Author:

sommeria

Message:

various bugs repaired

File:

• trunk/src/filter_tps.m (modified) (2 diffs)

trunk/src/filter_tps.m

@@ -79 +79 @@
 %% calculate tps coeff in each subdomain
 for isub=1:NbSubDomain
+    isub
     SubRange(1,:,isub)=[CentreX(isub)-0.55*Siz(1) CentreX(isub)+0.55*Siz(1)];%bounds of subdomain #isub in x coordinate
     SubRange(2,:,isub)=[CentreY(isub)-0.55*Siz(2) CentreY(isub)+0.55*Siz(2)];%bounds of subdomain #isub in y coordinate
@@ -87 +88 @@
         ind_sel_previous=ind_sel;% record the set of selected vector indices for next iteration
         ind_sel=find(Coord(:,1)>=SubRange(1,1,isub) & Coord(:,1)<=SubRange(1,2,isub) & Coord(:,2)>=SubRange(2,1,isub) & Coord(:,2)<=SubRange(2,2,isub));
+        numel(ind_sel)
         % if no vector in the subdomain  #isub, skip the subdomain
         if isempty(ind_sel)