comparison +noname/printSolutions.m @ 142:484b48e95c83

Removed ylim from setup1dPlot added some comments and fixed timestepper paramters.
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 29 Feb 2016 15:00:52 +0100
parents c7efff913935
children
comparison
equal deleted inserted replaced
141:cb2b12246b7e 142:484b48e95c83
1 function printSolutions(filename) 1 function printSolutions(filename)
2 sf = SolutionFile(filename); 2 sf = SolutionFile(filename);
3 3
4 method = {}; 4 name = {};
5 order = []; 5 order = [];
6 m = []; 6 m = [];
7 T = []; 7 T = [];
8 t = []; 8 t = [];
9 runtime = []; 9 runtime = [];
11 11
12 for i = 1:length(sf.keys) 12 for i = 1:length(sf.keys)
13 key = sf.keys{i}; 13 key = sf.keys{i};
14 entry = sf.get(key); 14 entry = sf.get(key);
15 15
16 16 name = [name key.name];
17 method = [method key.method];
18 order = [order key.order]; 17 order = [order key.order];
19 m = [m key.m]; 18 m = [m key.m];
20 T = [T key.T]; 19 T = [T key.T];
21 t = [t entry.repr.t]; 20 t = [t entry.repr.t];
22 runtime = [runtime entry.runtime]; 21 runtime = [runtime entry.runtime];
23 k = [k entry.k]; 22 k = [k entry.k];
24 end 23 end
25 24
26 methodW = findFieldWidth('%s',method); 25 nameW = findFieldWidth('%s', name);
27 orderW = findFieldWidth('%d',order); 26 orderW = findFieldWidth('%d', order);
28 mW = findFieldWidth('%d',m); 27 mW = findFieldWidth('%d', m);
29 TW = findFieldWidth('%d',T); 28 TW = findFieldWidth('%d', T);
30 tW = findFieldWidth('%.3e',t); 29 tW = findFieldWidth('%.3e', t);
31 runtimeW = findFieldWidth('%.3f',runtime); 30 runtimeW = findFieldWidth('%.3f', runtime);
32 kW = findFieldWidth('%.4f',k); 31 kW = findFieldWidth('%.4f', k);
33 32
34 for i = 1:length(sf.keys) 33 for i = 1:length(sf.keys)
35 fprintf('[%*s: o=%-*d, m=%-*d, T=%-*d]: t=%-*.3e, runtime=%*.3f, k=%*.4f\n',methodW, method{i}, orderW,order(i),mW,m(i),TW,T(i), tW, t(i), runtimeW,runtime(i), kW, k(i)); 34 fprintf('[%*s: o=%-*d, m=%-*d, T=%-*d]: t=%-*.3e, runtime=%*.3f, k=%*.4f\n', nameW, name{i}, orderW,order(i),mW,m(i),TW,T(i), tW, t(i), runtimeW,runtime(i), kW, k(i));
36 end 35 end
37 36
38 end 37 end