Homec4science
Diffusion yearII_reports (master) Browse

rSTICAZZI/physnum/rap7/Convergence.m
0cf006019435master

Last Change
Raw File

Convergence.m
View Options

% Ce script Matlab automatise la production de resultats
% lorsqu'on doit faire une serie de simulations en
% variant un des parametres d'entree.
%
% Il utilise les arguments du programme (voir ConfigFile.h)
% pour remplacer la valeur d'un parametre du fichier d'input
% par la valeur scannee.
%
%% Parametres %%
%%%%%%%%%%%%%%%%
repertoire = './'; % Chemin d'acces au code compile (NB: enlever le ./ sous Windows)
executable = 'exercice7'; % Nom de l'executable (NB: ajouter .exe sous Windows)
input = 'configuration.in'; % Nom du fichier d'entree de base
out_mode = 3; % target mode, do not change
L = 20;
tfin = 1.5;
% number of samples
n = 1500;
% point
target = 5;
beta = 0.9;
u = 6;
dt = zeros(n, 1);
sampling = zeros(n, 1);
for i=1:n
dt(i) = tfin / (i + 50);
sampling(i) = i + 50;
end
dx = dt * u / beta;
omega = 5;
f_ref = - sin(omega * (target / u - tfin));
b_left_type = 2;
% border A * sin(omega * t)
b_left = sprintf('{%.15g, %.15g, %.15g}', 1, 5, 0);
b_right_type = 3;
% border follow
b_right = '';
% border fixed
%b_right = sprintf('%.15g', 0);
%% Simulations %%
%%%%%%%%%%%%%%%%%
output = cell(n);
for i=1:n
output{i} = sprintf('m_outs/conv_dt%.15g', dt(i));
cmd = sprintf('%s%s %s N=%d L=%.15g u=%.15g t=%.15g dt=%.15g "b_left=%s" b_left_type=%d "b_right=%s" b_right_type=%d out_mode=%d target=%.15g sampling=%d output=%s', ...
repertoire, executable, input,...
floor(L/dx(i)), L, u, ...
tfin, dt(i), ...
b_left, b_left_type, ...
b_right, b_right_type, ...
out_mode, target, sampling(i), output{i});
disp(cmd);
system(cmd);
end
%% Loading data %%
%%%%%%%%%%%%%
% pre-build f
f = zeros(n, 1);
dtnew = zeros(n, 1);
j= 1;
for i=1:n
data = load(sprintf('%s.out', output{i}));
if abs(data(end, 1) - tfin) < 1e-14
args = [target-2*dx(i) target-dx(i) target target+dx(i) target+2*dx(i)];
fit = polyfit(args, data(end, 2:6), 1);
val = abs(fit(1) * target + fit(2) - f_ref); % center value
if (val < 1e-3)
f(j) = val;
dtnew(j) = dt(i);
fprintf('%d -> f = %.15g\n', j, f(j));
else
j = j - 1;
end
else
j = j - 1;
end
j = j + 1;
end
% redefine n
m = j-1;
f = f(1:m);
dtnew = dtnew(1:m);
%% PLOT %%
figure
grid on
xlabel('t')
ylabel(sprintf('f(x = %.15g', target))
%plot(dt(1:restr), abs(f(1:restr) - last), 'k+');
loglog(dtnew, f, 'k+');
%% LINEAR FIT %%
fit = polyfit(log(dtnew), log(f), 1);
fprintf('fit: %g * x + %g\n', fit(1), fit(2));
%% WRITE ON FILE %%
outfile = 'outs/conv5m_alt.dat';
% print data on file
out = fopen(outfile, 'w');
fprintf(out, 'dt f\n');
for i = 1:m
fprintf(out, '%.15g %.15g\n', dtnew(i), f(i));
end
fclose(out);
c4science · Help