From 7fa776c751858aa14afe65729f7cb0830af15e67 Mon Sep 17 00:00:00 2001
From: Olivier Sauter <olivier.sauter@epfl.ch>
Date: Tue, 16 Apr 2002 20:59:28 +0000
Subject: [PATCH] to go faster
git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@1872 d63d8f72-b253-0410-a779-e742ad2e26cf
---
AUG/sxrbgeometry.m | 12 ++++--------
1 file changed, 4 insertions(+), 8 deletions(-)
diff --git a/AUG/sxrbgeometry.m b/AUG/sxrbgeometry.m
index b8145297..9d631c2b 100644
--- a/AUG/sxrbgeometry.m
+++ b/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;
--
GitLab