Homec4science
Diffusion yearII_reports (master) Browse

rSTICAZZI/physnum/rap7/Convnot
ba5e3b30d249master

Last Change
Raw File

Convnot
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 = 1; % target mode, do not change
method = 1; % order 1 or 2
L = 20;
tfin = 1.5;
% number of samples
n = 1000;
% point
target = 5;
beta = 0.9;
u = 6;
dt = zeros(n, 1);
N = zeros(n, 1);
dx = zeros(n, 1);
omega = 5;
f_ref = - sin(omega * (target / u - tfin));
for i=1:n
if (i == n)
dt(i) = tfin / 10000;
else
dt(i) = tfin / (i + 50);
end
dx(i) = dt(i) * u / beta;
N(i) = floor(L/dx(i));
end
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 = {};
dtnew = {};
M = 0;
for i=1:n
j = floor(target * N(i) / L); % approximation
if ((j*dx(i)) == target) || (i == n)
output{M+1} = sprintf('m_outs/conv_dt%.15g', dt(i));
dtnew{M+1} = 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 output=%s', ...
repertoire, executable, input,...
N(i), L, u, ...
tfin, dt(i), ...
b_left, b_left_type, ...
b_right, b_right_type, ...
out_mode, target, output{M+1});
disp(cmd);
system(cmd);
M = M + 1;
end
end
%% Analyse %%
%%%%%%%%%%%%%
% pre-build f
f = zeros(M, 1);
for i=1:M
data = load(sprintf('%s.out', output{i}));
f(i) = data(end, 2);
%args = [dx(i)*j, dx(i)*(j+1)];
%fit = polyfit(args, data(M, 2:end), 1);
%f(i) = fit(1) * target + fit(2);
end
%% PLOT %%
dtnew_ = cell2mat(dtnew)';
f_ = [];
dtnew__ = [];
if method == 0
last = f(M);
dtnew__ = dtnew_(1:end-1);
f_ = abs(f(1:end-1) - last);
elseif method == 1
dtnew__ = dtnew_;
f_ = abs(f - f_ref);
end
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(dtnew__, f_, 1);
fprintf('fit: %g * x + %g\n', fit(1), fit(2));
%% WRITE ON FILE %%
outfile = sprintf('outs/conv5m%d.dat', method);
% print data on file
out = fopen(outfile, 'w');
fprintf(out, 'dt f\n');
for i = 1:length(dtnew__)
fprintf(out, '%.15g %.15g\n', dtnew__(i), f_(i));
end
fclose(out);
c4science · Help