From 475f7d5bfc76d82da4372610fd817c4a055e4596 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Fri, 18 Sep 2020 10:05:43 +0200
Subject: [PATCH] Adapt FRAP flux to new function signature of int_prob.
---
prob_laplace.m | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/prob_laplace.m b/prob_laplace.m
index 321695f..910b2da 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -177,8 +177,8 @@ s_dot(i) = sum(diff(T_mov.x).*f.^2./(g0.*(1-pt).*u_interp));
% plot(x_interp, f)
% pause()
end
-csvwrite([s, 'FRAP_Flux.csv'], [x_interp; f])
-csvwrite([s, 'FRAP_entropy.csv'], [T_mov.t; s_dot]);
+% csvwrite([s, 'FRAP_Flux.csv'], [x_interp; f])
+% csvwrite([s, 'FRAP_entropy.csv'], [T_mov.t; s_dot]);
%% Check integral of solution, should be mass conserving and sum to 1.
% integrate to 50, to avoid right boundary
[~, ind] = min(abs(T_prec(4).x-50));
@@ -262,7 +262,7 @@ p = nan(length(ls), n_T);
for j = 1:n_T
tic
parfor i = 1:length(ls)
- p(i, j) = int_prob(ls(i), F, x0, direction, T_mov, j, 5);
+ p(i, j) = int_prob(ls(i), F, x0, direction, j, T_mov, 5);
end
toc
end
--
GitLab