 Timestamp:
 Dec 11, 2011, 11:08:21 PM (12 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/src/civ_matlab.m
r326 r331 181 181 isy2=iby2+2; 182 182 183 %get the previous guess for displacement 184 185 if ~check_patch1 186 [Data,VelType]=read_civdata(ObjectName,InputField,ParamIn.VelType);%TO DEVELOP 187 end 188 189 Guess =interp_tps(Data,X,Y) 183 190 184 % Data.Civ1_U_Diff=zeros(size(Data.Civ1_X)); 191 185 % Data.Civ1_V_Diff=zeros(size(Data.Civ1_X)); … … 209 203 PointCoord(:,1)=reshape(GridX,[],1); 210 204 PointCoord(:,2)=reshape(GridY,[],1); 205 206 Guess = tps_eval(PointCoord,[Data.Civ1_X_tps Data.Civ1_Y_tps]) 207 Shiftx=Guess*Data.Civ1_U_tps; 208 Shifty=Guess*Data.Civ1_V_tps; 209 211 210 if ~isempty(par_civ2.maskname)&& ~strcmp(maskname,par_civ2.maskname)% mask exist, not already read in civ1 212 211 mask=imread(par_civ2.maskname); … … 673 672 FF=zeros(length(X),1); 674 673 check_empty=zeros(1,NbSubDomain); 674 SubRangx=zeros(NbSubDomain,2);%initialise the positions of subdomains 675 SubRangy=zeros(NbSubDomain,2); 675 676 for isub=1:NbSubDomain 676 677 SubRangx(isub,:)=[CentreX(isub)0.55*SizX CentreX(isub)+0.55*SizX];
Note: See TracChangeset
for help on using the changeset viewer.