to go faster

git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@1872 d63d8f72-b253-0410-a779-e742ad2e26cf

AUG/sxrbgeometry.m

@@ -59,12 +59,8 @@ end
 % compute R
 zmag=varargin{1};
 varargout{1}=NaN*ones(length(ang_B),length(zmag));
-for i=1:length(zmag)
-  if zmag(i)>Z_B
-    ii=find(ang_B<((pi-atan((zmag(i)-Z_B)/(R_B-0.)))*180/pi));
-    varargout{1}(ii,i) = R_B - (zmag(i)-Z_B)./tan(pi-ang_B(ii).*pi/180);
-  elseif zmag(i)<Z_B
-    ii=find(ang_B>((pi-atan((zmag(i)-Z_B)/(R_B-0.)))*180/pi));
-    varargout{1}(ii,i) = R_B - (zmag(i)-Z_B)./tan(pi-ang_B(ii).*pi/180);
-  end
+for i=1:length(ang_B)
+  varargout{1}(i,:) = R_B - (zmag'-Z_B)./tan(pi-ang_B(i).*pi/180);
 end
+iiout=find(varargout{1}<=0 | varargout{1}>=2.5);
+varargout{1}(iiout)=NaN;