# Changeset 1059

Ignore:
Timestamp:
Nov 29, 2018, 5:44:44 PM (2 years ago)
Message:

various bugs repaired

Location:
trunk/src
Files:
11 edited

Unmodified
Removed

• ## trunk/src/filter_tps.m

 r1027 %% 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 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)

• ## trunk/src/phys_XYZ.m

 r1027 testangle=0; test_refraction=0; if exist('Zindex','var')&& isequal(Zindex,round(Zindex))&& Zindex>0 && isfield(Calib,'SliceCoord')&&length(Calib.SliceCoord)>=Zindex if exist('Zindex','var')&& isequal(Zindex,round(Zindex))&& Zindex>0 && isfield(Calib,'SliceCoord')&&size(Calib.SliceCoord,1)>=Zindex if isfield(Calib, 'SliceAngle') && ~isequal(Calib.SliceAngle,[0 0 0]) && ~isequal(Calib.SliceAngle(Zindex,:),[0 0 0]) testangle=1; %         Z0=norm_plane*Calib.SliceCoord(Zindex,:)'/norm_plane(3); end Z0=Calib.SliceCoord(Zindex,3);%horizontal plane z=cte Z0=Calib.SliceCoord(Zindex,3);%horizontal plane z=cte %     end Z0virt=Z0;