Mercurial > repos > public > sbplib

comparison convergencePlot.m @ 755:14f0058356f2 feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with feature/grids
author Martin Almquist <malmquist@stanford.edu>
date Fri, 15 Jun 2018 18:10:26 -0700
parents 38446922c32a
children
comparison
equal deleted inserted replaced
754:5264ce57b573 755:14f0058356f2
1 function hand = convergencePlot(orders, h, e)
2 N = length(orders);
3
4 fh = figure();
5 ah = axes();
6 ah.XScale = 'log';
7 ah.YScale = 'log';
8 hold on
9 ph = {};
10 phc = {};
11 legends = {};
12 for i = 1:N
13 ph{i} = loglog(h{i}, e{i});
14 phc{i} = plotConvergenceFit(orders{i}, h{i}, e{i});
15
16 ph{i}.LineStyle = 'none';
17 ph{i}.Marker = Color.solidMarkers{i};
18 ph{i}.MarkerSize = 12;
19 ph{i}.Color = Color.colors{i};
20 ph{i}.MarkerFaceColor = Color.colors{i};
21
22 legends{i} = sprintf('$o = %d$', orders{i});
23 end
24 hold off
25
26 lh = legend([ph{:}], legends);
27 lh.Interpreter = 'latex';
28 lh.Location = 'SouthEast';
29
30 for i = 1:N
31 uistack(phc{i}, 'bottom');
32 end
33
34 xlabel('$h$', 'interpreter', 'latex')
35 ylabel('Error', 'interpreter', 'latex')
36
37 % xlim([0.7e-2, 1e-1])
38 % ylim([3e-5, 4])
39
40 grid on
41
42 ah = gca();
43 ah.TickLabelInterpreter = 'latex';
44 setFontSize(fh);
45
46 % if savePngs
47 % savepng(fh, 'fig/conv/conv',600)
48 % end
49
50 hand = struct();
51 hand.fig = fh;
52 hand.data = ph;
53 hand.fits = phc;
54 hand.legend = lh;
55 end