From 6932c1f7e13c03862d050c14ef72efd5db64da7b Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Mon, 16 Nov 2020 14:58:44 +0100
Subject: [PATCH] Finer mesh around x0 also for client.
---
@Ternary_model/create_mesh.m | 7 +++++++
1 file changed, 7 insertions(+)
diff --git a/@Ternary_model/create_mesh.m b/@Ternary_model/create_mesh.m
index 35884d8..6646f2f 100644
--- a/@Ternary_model/create_mesh.m
+++ b/@Ternary_model/create_mesh.m
@@ -33,5 +33,12 @@ elseif strcmp(T.mode, 'Client') && T.v ~= 0
% T.x = linspace(0, T.system_size, 24000);
T.x = [linspace(0, -T.a-1, 1000), linspace(-T.a-1+0.0001, -T.a+1, 4800),...
linspace(-T.a+1+0.0001, T.system_size, 1000)];
+ % Fine mesh close to x0 if ic = 'Gauss'
+ [~, ind] = min(abs(T.x-T.x0));
+ [~, ind_plus] = min(abs(T.x-(T.x0+0.6)));
+ [~, ind_minus] = min(abs(T.x-(T.x0-0.6)));
+ T.x = [T.x(1:ind_minus),...
+ linspace(T.x(ind_minus)+0.0001, T.x(ind_plus)-0.0001, 1000),...
+ T.x(ind_plus:end)];
end
end
--
GitLab