1 function fn = gen_dune_params(params)
2 %
function fn = gen_dune_params(params)
4 %
function generating
new params file from
default file
for dune-detailed
5 % simulations. The new file is written to $RBMATLABTEMP/
burgers_fem
6 % the default settings are taken from params.
8 % required fields of params:
9 % dune_params_fn: the filename of the params-file with default
10 % settings. Further the simulation parameters are required:
11 % udir_height, k, cx, cy, udir_amplitude, udir_xscale, T, nt,
12 % max_newton_steps, newton_eps, dune_grid_fn
14 % Bernard Haasdonk 18.1.2007
21 [p,f,e] = fileparts(params.dune_params_fn);
22 new_params_fn = fullfile(rbmatlabtemp,
'burgers_fem',[f,
'_new.dat']);
23 pnames = {
'mu4',
'mu1',
'mu2',
'mu3',
'mu5',
'mu6',...
24 'EndTime',
'deltaT',
'TimeStepNumber',...
25 'maxNewtonStep',
'IteratedNewtonSolverEpsilon',...
26 'Task',
'GrapeDisplay',
'VisualizeRBSpace',
'VisualizeEachNewtonStep',...
27 'VisualizeEachTimeStep',...
28 'OutputPath',
'OutputPrefix',
'Grid'};
30 % generate outputprefix
for snapshot-files, e.g.
31 % mu1_0.1,mu2_0,mu3_0,mu4_0.2,mu5_1,mu6_1
34 ['mu1_',num2str(mu(1)),...
35 ',mu2_',num2str(mu(2)),...
36 ',mu3_',num2str(mu(3)),...
37 ',mu4_',num2str(mu(4)),...
38 ',mu5_',num2str(mu(5)),...
39 ',mu6_',num2str(mu(6)) ...
42 pvalues = {[num2str(params.udir_height),' % = mu_4'],...
43 [num2str(params.k),' % = mu_1'],...
44 [num2str(params.cx),' % = mu_2'],...
45 [num2str(params.cy),' % = mu_3'],...
46 [num2str(params.udir_amplitude),' % = mu_5'],...
47 [num2str(params.udir_xscale),' % = mu_6'],...
48 num2str(params.T), ...
49 num2str(params.T/params.nt), ...
50 num2str(params.nt), ...
51 num2str(params.max_newton_steps),...
52 num2str(params.newton_eps),...
54 '0 % == GrapeDisplay',...
55 '0 % == VisualizeRBSpace',...
56 '0 % == Visualize each newton step',...
57 '0 % == Visualize each time step',...
58 [fullfile(rbmatlabtemp,'burgers_fem'),' % = outputpat for snapshots'],...
63 disp([
'generating new params file']);
65 disp([
'from old ',params.dune_params_fn]);
66 change_params_file(params.dune_params_fn, new_params_fn, ...
70 % TO BE ADJUSTED TO NEW SYNTAX