From 07a05989f77f486a9737c4e741d40ebe75aa5b94 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Sun, 7 Feb 2021 17:21:12 +0100
Subject: [PATCH] nansum for probs in prob_laplace.
---
prob_laplace.m | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/prob_laplace.m b/prob_laplace.m
index 7d35a56..8658102 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -219,7 +219,7 @@ end
toc
%% Normalization factor
% N = normalization(T, x0, 0, t_ind, direc, -params{1})
-N=sum(p)/length(p)*l_max;
+N=nansum(p)/length(p)*l_max;
% N = sum(p)*0.001;
m = sum(ls.*p/N)/length(p)*l_max;
% should sum to one
@@ -263,6 +263,9 @@ parfor i = 1:length(ls)
end
toc
save(['prob_7_7_bp', num2str(bp(j))]);
+pa = '/Users/hubatsch/Nextcloud/Langevin_vs_MeanField/data/constD/';
+csvwrite([pa, 'jump_length_7_7_bp', num2str(bp(j)), '.csv'], ls)
+csvwrite([pa, 'prob_7_7_bp', num2str(bp(j)), '.csv'], p);
end
%% Mean cross jump length chi 7/3
run_jump_lengths(-5, 1, 7/3, 1, 4, '73_1_01', 10^-6, 3, 0.5);
--
GitLab