Changeset 905 for trunk/src/calc_field_tps.m

Timestamp:

May 30, 2015, 8:44:20 PM (9 years ago)

Author:

sommeria

Message:

projection with tps corrected + minor changes

File:

• trunk/src/calc_field_tps.m (modified) (1 diff)

trunk/src/calc_field_tps.m

@@ -113 +113 @@
     nbvec_sub=NbCentre(isub);
     check_range=(Coord_interp >=ones(nb_sites,1)*SubRange(:,1,isub)' & Coord_interp<=ones(nb_sites,1)*SubRange(:,2,isub)');
-    ind_sel=find(sum(check_range,2)==nb_coord);
+    ind_sel=find(sum(check_range,2)==nb_coord);% select points whose all coordinates are in the prescribed range
     nbval(ind_sel)=nbval(ind_sel)+1;% records the number of values for eacn interpolation point (in case of subdomain overlap)
     if check_grid