From b037bc9438b57e5bc8c6402411ed14e67dfacdbb Mon Sep 17 00:00:00 2001
From: Johannes Pfeifer <jpfeifer@gmx.de>
Date: Tue, 17 May 2022 18:29:42 +0200
Subject: [PATCH] GSA sampling from prior range: replace infinity by huge
 number to avoid NaN

---
 matlab/gsa/stab_map_.m | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/matlab/gsa/stab_map_.m b/matlab/gsa/stab_map_.m
index 86f3aef98..94a49a20f 100644
--- a/matlab/gsa/stab_map_.m
+++ b/matlab/gsa/stab_map_.m
@@ -175,7 +175,9 @@ if fload==0
         end
         if opt_gsa.prior_range
             for j=1:np
-                lpmat(:,j)=lpmat(:,j).*(bounds.ub(j+nshock)-bounds.lb(j+nshock))+bounds.lb(j+nshock);
+                lower_bound=max(-options_.huge_number,bounds.lb(j+nshock));
+                upper_bound=min(options_.huge_number,bounds.ub(j+nshock));
+                lpmat(:,j)=lpmat(:,j).*(upper_bound-lower_bound)+lower_bound;
             end
         else
             xx=prior_draw_gsa(0,[lpmat0 lpmat]);