Skip to content
Snippets Groups Projects
Commit 6932c1f7 authored by Lars Hubatsch's avatar Lars Hubatsch
Browse files

Finer mesh around x0 also for client.

parent 17d69652
No related branches found
No related tags found
No related merge requests found
......@@ -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
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment