From f77666c5b0c7b0d098fb7bb8f9ae85ac8db0086d Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Tue, 17 Nov 2020 14:32:42 +0100
Subject: [PATCH] Reduced difference in normalization to 2.5%.
---
prob_laplace.m | 14 +++++++-------
1 file changed, 7 insertions(+), 7 deletions(-)
diff --git a/prob_laplace.m b/prob_laplace.m
index 0ae9bb8..6ae42f3 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -304,22 +304,22 @@ end
%% MOVING BOUNDARY JUMP LENGTH DISTRIBUTION
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% load entropy_mov_bound_in_to_out.mat
-direc = -1; % 1: jump from left to right, -1: jump from right to left.
-x0 = sort(10-direc*(0:0.02:4.01));
+direc = 1; % 1: jump from left to right, -1: jump from right to left.
+x0 = sort(10-direc*(0:0.005:4.01));
% t_snap = repmat((0.25:0.5:9.75)', 1, length(x0)); % midpoint
t_snap = repmat((0.:0.5:9.5)', 1, length(x0)); % starting point
si_t = size(t_snap);
-x0 = sort(10-T_mov.v*t_snap-direc*repmat(0:0.02:4.01, [si_t(1), 1]), 2);
+x0 = sort(10-T_mov.v*t_snap-direc*repmat(0:0.005:4.01, [si_t(1), 1]), 2);
%% Run simulations for 'delta' IC across outside
F = {};
t = linspace(0, 1, 21);
-params = params_mov_bound;
for j = 1:si_t(1)
tic
parfor i = 1:si_t(2)
+ params = params_mov_bound;
+ params{8} = x0(j, i);
+ params{1} = -10+params{9}*t_snap(j);
F{i, j} = Ternary_model(0, 'Gauss', params, t, 0.1);
- F{i, j}.x0 = x0(j, i);
- F{i, j}.a = -10+F{i, j}.v*t_snap(j);
F{i, j}.solve_tern_frap();
end
toc
@@ -336,7 +336,7 @@ for j = 1:n_T
tic
F1 = F(:, j);
bp = bound(j);
- for i = 1:length(ls)
+ parfor i = 1:length(ls)
p(i, j) = int_prob(ls(i), F1, x0(j, :), direc, j, bp, 10, T_mov);
end
toc
--
GitLab