Mercurial > repos > public > sbplib

annotate +multiblock/Contour.m @ 1227:02dfe3a56742 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add Upwind ElasticAnisotropic schemes. Seem to work really well!
author Martin Almquist <malmquist@stanford.edu>
date Sat, 16 Nov 2019 14:26:06 -0800
parents 97b9a0023d38
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
586
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
1 classdef Contour < handle
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
2 properties
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
3 grid
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
4 contours
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
5 nContours
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
7 ZData
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
8 CData
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
9
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
10 end
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
11
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
12 methods
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
13 function obj = Contour(g, gf, nContours)
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
14 obj.grid = g;
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
15 obj.nContours = nContours;
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
16
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
17 coords = obj.grid.points();
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
18 X = obj.grid.funcToPlotMatrices(coords(:,1));
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
19 Y = obj.grid.funcToPlotMatrices(coords(:,2));
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
20
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
21 V = obj.grid.funcToPlotMatrices(gf);
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
22
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
23
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
24 holdState = ishold();
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
25 hold on
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
26
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
27 contours = {1, obj.grid.nBlocks};
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
28 for i = 1:obj.grid.nBlocks
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
29 [~, contours{i}] = contour(X{i}, Y{i}, V{i},obj.nContours);
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
30 contours{i}.LevelList = contours{1}.LevelList;
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
31 end
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
32
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
33 if holdState == false
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
34 hold off
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
35 end
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
36
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
37 obj.contours = [contours{:}];
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
38
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
39 obj.ZData = gf;
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
40 obj.CData = gf;
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
41 end
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
42
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
43 function set(obj, propertyName, propertyValue)
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
44 set(obj.contours, propertyName, propertyValue);
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
45 end
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
46
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
47 function obj = set.ZData(obj, gf)
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
48 obj.ZData = gf;
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
49
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
50 V = obj.grid.funcToPlotMatrices(gf);
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
51 for i = 1:obj.grid.nBlocks
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
52 obj.contours(i).ZData = V{i};
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
53 end
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
54 end
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
55
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
56 function obj = set.CData(obj, gf)
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
57 obj.CData = gf;
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
58
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
59 V = obj.grid.funcToPlotMatrices(gf);
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
60 for i = 1:obj.grid.nBlocks
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
61 obj.contours(i).CData = V{i};
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
62 end
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
63 end
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
64 end
97b9a0023d38 Add multiblock contour plot
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
65 end