QuadTreeCode/QuadTreeTest.m at master · tjssmy/QuadTreeCode · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51

QTGlobalsInit

CreateQTkFct(300,200,200.0,300.0,'r',[10,0.3,1],0)
CreateInitQTBlocks(2,3,0);

RefineReg(1) = QTCirRefineReg(30,30,20,1);
RefineReg(2) = QTCirRefineReg(170,260,10,1);
RefineReg(3) = QTCirRefineReg(30,275,10,1);

CreateQTSmap(RefineReg);
RefineQTBlocks(5,0.1)

figure; DrawQTBlocks(1,0,1,1,0); title('Refined blocks');

% GetQTFct(KVAL,'b',1,'block k no interp')
GetQTSurfaceAndFct(KVAL,1,1,'block k with interp',[]);


% IntBC(1) = QTRectIntBC(10,40,10,40,1);
% IntBC(2) = QTRectIntBC(160,180,250,270,0);

IntBC(1) = QTCirIntBC(20,20,10,[0 1e3],'g');
IntBC(end+1) = QTCirIntBC(170,260,10,[0 8e3 -1000],'m');
IntBC(end+1) = QTCirIntBC(30,275,10,[0 8e3 -0000],'m');
% IntBC(end+1) = QTCirIntBC(52,148,5,[0 0e3 -000],'m');


IntBC = SolveQTDCSol(['f','f','f','f'],IntBC);

% GetQTFct(VVAL,'b',1,'Solution no interp')
[~,h] = GetQTSurfaceAndFct(VVAL,1,1,'Solution with interp',[]);


GetQTFctGradAndFlow(1,'flow');


PrintQTIntBCVals(IntBC)

IntBC(3).trans = 1;
IntBC(3).TransMode(1) = 'c';
IntBC(3).valTrans{1} = [0];
IntBC(3).TransMode(2) = 'c';
IntBC(3).valTrans{2} = [8e3]';
IntBC(3).TransMode(3) = 'l';
IntBC(3).valTrans{3} = [0.0,-3,1000]';

figure

SolveQTTransSol(1500,1000,10,IntBC)