From e1fd51fc530779beae2749a3684ec3038bd636c2 Mon Sep 17 00:00:00 2001
From: Olivier Sauter <olivier.sauter@epfl.ch>
Date: Tue, 8 Aug 2017 16:12:09 +0000
Subject: [PATCH] use isfinite instead of ~isnan since more general and avoids
NaNs and Infs, ready to move function up for all machines
git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@8147 d63d8f72-b253-0410-a779-e742ad2e26cf
---
crpptbx/JET/get_grids_1d.m | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/crpptbx/JET/get_grids_1d.m b/crpptbx/JET/get_grids_1d.m
index 497eb6c5..57d2be4f 100644
--- a/crpptbx/JET/get_grids_1d.m
+++ b/crpptbx/JET/get_grids_1d.m
@@ -52,9 +52,9 @@ for it=1:length(gdat_data.t)
it_rt_eff = it_rt(it);
it_vol_eff = it_vol(it);
if (nbdim_x == 1)
- ii=find(~isnan(gdat_data.grids_1d.rhopolnorm));
+ ii=find(isfinite(gdat_data.grids_1d.rhopolnorm));
else
- ii=find(~isnan(gdat_data.grids_1d.rhopolnorm(:,it)));
+ ii=find(isfinite(gdat_data.grids_1d.rhopolnorm(:,it)));
end
if (nbdim_x == 1)
if length(ii)==length(gdat_data.grids_1d.rhopolnorm)
--
GitLab