LINE Solver
MATLAB API documentation
Loading...
Searching...
No Matches
test_qrf_bas_1.m
1%% QRF BAS test from code release
2f=1;
3N=10;
4M=5;
5ZM=2;
6MR=5;
7F=[5,10,10,10,10];
8K=[2,2,2,2,2];
9%% r
10r=[];
11r(1,1)=0;
12r(1,2)=5.000000e-01;
13r(1,3)=0;
14r(1,4)=0;
15r(1,5)=5.000000e-01;
16r(2,1)=5.000000e-01;
17r(2,2)=0;
18r(2,3)=5.000000e-01;
19r(2,4)=0;
20r(2,5)=0;
21r(3,1)=0;
22r(3,2)=5.000000e-01;
23r(3,3)=0;
24r(3,4)=5.000000e-01;
25r(3,5)=0;
26r(4,1)=0;
27r(4,2)=0;
28r(4,3)=5.000000e-01;
29r(4,4)=0;
30r(4,5)=5.000000e-01;
31r(5,1)=5.000000e-01;
32r(5,2)=0;
33r(5,3)=0;
34r(5,4)=5.000000e-01;
35r(5,5)=0;
36%% mu
37mu =[];
38mu(1,1,1)=1.016186e+00;
39mu(1,1,2)=2.585708e-05;
40mu(1,2,1)=1.569888e-03;
41mu(1,2,2)=1.413298e-02;
42mu(2,1,1)=1.016186e+00;
43mu(2,1,2)=2.585708e-05;
44mu(2,2,1)=1.569888e-03;
45mu(2,2,2)=1.413298e-02;
46mu(3,1,1)=1.016186e+00;
47mu(3,1,2)=2.585708e-05;
48mu(3,2,1)=1.569888e-03;
49mu(3,2,2)=1.413298e-02;
50mu(4,1,1)=1.016186e+00;
51mu(4,1,2)=2.585708e-05;
52mu(4,2,1)=1.569888e-03;
53mu(4,2,2)=1.413298e-02;
54mu(5,1,1)=1.016186e+00;
55mu(5,1,2)=2.585708e-05;
56mu(5,2,1)=1.569888e-03;
57mu(5,2,2)=1.413298e-02;
58%% v
59v =[];
60v(1,1,1)=0;
61v(1,1,2)=0;
62v(1,2,1)=0;
63v(1,2,2)=0;
64v(2,1,1)=0;
65v(2,1,2)=0;
66v(2,2,1)=0;
67v(2,2,2)=0;
68v(3,1,1)=0;
69v(3,1,2)=0;
70v(3,2,1)=0;
71v(3,2,2)=0;
72v(4,1,1)=0;
73v(4,1,2)=0;
74v(4,2,1)=0;
75v(4,2,2)=0;
76v(5,1,1)=0;
77v(5,1,2)=0;
78v(5,2,1)=0;
79v(5,2,2)=0;
80%% BB
81BB =[];
82BB(1,1)=0;
83BB(1,2)=0;
84BB(1,3)=0;
85BB(1,4)=0;
86BB(1,5)=0;
87BB(2,1)=0;
88BB(2,2)=1;
89BB(2,3)=0;
90BB(2,4)=0;
91BB(2,5)=0;
92BB(3,1)=0;
93BB(3,2)=1;
94BB(3,3)=0;
95BB(3,4)=0;
96BB(3,5)=1;
97BB(4,1)=0;
98BB(4,2)=0;
99BB(4,3)=0;
100BB(4,4)=0;
101BB(4,5)=1;
102BB(5,1)=0;
103BB(5,2)=1;
104BB(5,3)=0;
105BB(5,4)=0;
106BB(5,5)=1;
107%% MM
108MM=[];
109MM(1,1)=0;
110MM(1,2)=0;
111MM(2,1)=2;
112MM(2,2)=0;
113MM(3,1)=2;
114MM(3,2)=5;
115MM(4,1)=5;
116MM(4,2)=0;
117MM(5,1)=5;
118MM(5,2)=2;
119%% ZZ
120ZZ = [];
121ZZ(1)=0;
122ZZ(2)=1;
123ZZ(3)=2;
124ZZ(4)=1;
125ZZ(5)=2;
126%% MM1
127MM1 =[];
128MM1(1,1)=0;
129MM1(1,2)=0;
130MM1(1,3)=0;
131MM1(1,4)=0;
132MM1(1,5)=0;
133MM1(2,1)=-1;
134MM1(2,2)=-1;
135MM1(2,3)=-1;
136MM1(2,4)=-1;
137MM1(2,5)=3;
138MM1(3,1)=0;
139MM1(3,2)=0;
140MM1(3,3)=0;
141MM1(3,4)=0;
142MM1(3,5)=0;
143MM1(4,1)=-1;
144MM1(4,2)=5;
145MM1(4,3)=-1;
146MM1(4,4)=-1;
147MM1(4,5)=-1;
148%%
149p2est=qrf_bas_mmi(f,M,MR,MM,MM1,ZZ,ZM,BB,K,F,N,mu,v,r)