From 1f13d183eb4725db46ca4496ed0db78c7cc82bf7 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Wed, 10 Feb 2021 11:06:00 +0100
Subject: [PATCH] Make mesh more precise close to zero to have good agreement
with Stefano.
---
@Ternary_model/create_mesh.m | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/@Ternary_model/create_mesh.m b/@Ternary_model/create_mesh.m
index 8ffca12..2464527 100644
--- a/@Ternary_model/create_mesh.m
+++ b/@Ternary_model/create_mesh.m
@@ -4,7 +4,7 @@ function create_mesh(T)
if strcmp(T.mode, 'Constituent') || strcmp(T.mode, 'Const_mob') || T.v == 0
% Start with left side of mesh, at r=0/x=0, with a given step size
- x = linspace(0, -0.75*T.a, 40);
+ x = linspace(0, -0.75*T.a, 200);
x = [x, x(end)+0.001/T.precision];
while x(end)<-T.a
x_temp = x(end)-x(end-1);
--
GitLab