Changeset 1142 for trunk/src/filter_tps.m

Timestamp:

May 5, 2024, 12:29:30 PM (5 months ago)

Author:

sommeria

Message:

filter parameter improved in filter_tps

File:

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

trunk/src/filter_tps.m

@@ -6 +6 @@
 % SubRange(NbCoord,2,NbSubdomain): range (min, max) of the coordinates x and y respectively, for each subdomain
 % NbCentre(NbSubdomain): number of source points for each subdomain
-% FF: false flags
+% FF: false flags preserved from the input, or equal to 20 for vectors excluded by the difference with the smoothed field
 % U_smooth, V_smooth: filtered velocity components at the positions of the initial data
 % Coord_tps(NbCentre,NbCoord,NbSubdomain): positions of the tps centres
-% U_tps,V_tps: weight of the tps for each subdomain
+% U_tps,V_tps: weight of the tps centers for each subdomain
 % to get the interpolated field values, use the function calc_field.m
 %
@@ -62 +62 @@
 %% smoothing parameter: CHANGED 03 May 2024 TO GET RESULTS INDEPENDENT OF SUBDOMAINSIZE
 %smoothing=Siz(1)*Siz(2)*FieldSmooth/1000%optimum smoothing increase as the area of the subdomain (division by 1000 to reach good values with the default GUI input)
-smoothing=sqrt(Siz(1)*Siz(2)/SubDomainSize)*FieldSmooth;%optimum smoothing increase as the typical mesh size =sqrt(SizX*SizY/SubDomainSize)^1/2
+NbVecSub=NbVec/NbSubDomain;% refined estimation of the nbre of vectors per subdomain
+smoothing=sqrt(Siz(1)*Siz(2)/NbVecSub)*FieldSmooth;%optimum smoothing increase as the typical mesh size =sqrt(SizX*SizY/NbVecSub)^1/2
 %% default output
 SubRange=zeros(NbCoord,2,NbSubDomain);%initialise the boundaries of subdomains
@@ -80 +81 @@
 %% 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
@@ -89 +89 @@
         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)
+        %disp([numel(ind_sel) ' vectors in subdomain #' num2str(isub)])
         % if no vector in the subdomain  #isub, skip the subdomain
         if isempty(ind_sel)
@@ -124 +124 @@
                 V_tps(1:NbCentre(isub)+3,isub)=V_tps_sub;
                 nb_select(ind_sel)=nb_select(ind_sel)+weight;
-                display(['tps done in subdomain # ' num2str(isub)  ' among ' num2str(NbSubDomain)])
+                display(['tps done with ' num2str(numel(ind_sel)) ' vectors in subdomain # ' num2str(isub)  ' among ' num2str(NbSubDomain)])
                 break
             % if too few selected vectors, increase the subrange for next iteration
@@ -146 +146 @@
                 V_tps(1:NbCentre(isub)+3,isub)=V_tps_sub;
                 nb_select(ind_sel(ind_ind_sel))=nb_select(ind_sel(ind_ind_sel))+weight;
-                display(['tps redone after elimination of erratic vectors in subdomain # ' num2str(isub) ' among ' num2str(NbSubDomain)])
+                display(['tps redone with ' num2str(numel(ind_sel)) ' vectors after elimination of ' num2str(numel(ind_ind_sel)) ' erratic vectors in subdomain # ' num2str(isub) ' among ' num2str(NbSubDomain)])
                 break
             end