-
Notifications
You must be signed in to change notification settings - Fork 1
/
SBP_quad_data.jl
114 lines (100 loc) · 20 KB
/
SBP_quad_data.jl
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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
# SBP Node data #
#N=1
Q_GLegendre_N1 = 6;
Q_GLegendre_rq_N1 = [-0.577350269189626,0.577350269189626,-1.000000000000000,-1.000000000000000,-0.577350269189626,0.577350269189626];
Q_GLegendre_sq_N1 = [-1.000000000000000,-1.000000000000000,-0.577350269189626,0.577350269189626,0.577350269189626,-0.577350269189626];
Q_GLegendre_wq_N1 = [0.333333333333333,0.333333333333333,0.333333333333333,0.333333333333333,0.333333333333333,0.333333333333333];
GQ_N1 = (Q_GLegendre_rq_N1, Q_GLegendre_sq_N1, Q_GLegendre_wq_N1)
#N=2
Q_GLegendre_N2 = 10;
Q_GLegendre_rq_N2 = [-0.774596669241483,0.000000000000000,0.774596669241483,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.774596669241483,0.000000000000000,0.774596669241483,-0.333333333333333];
Q_GLegendre_sq_N2 = [-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.774596669241483,0.000000000000000,0.774596669241483,0.774596669241483,0.000000000000000,-0.774596669241483,-0.333333333333333];
Q_GLegendre_wq_N2 = [0.083333333333333,0.200000000000000,0.083333333333333,0.083333333333333,0.200000000000000,0.083333333333333,0.083333333333333,0.200000000000000,0.083333333333333,0.900000000000000];
GQ_N2 = (Q_GLegendre_rq_N2, Q_GLegendre_sq_N2, Q_GLegendre_wq_N2)
#N=3
Q_GLegendre_N3 = 18;
Q_GLegendre_rq_N3 = [-0.861136311594054,-0.339981043584857,0.339981043584856,0.861136311594053,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.861136311594054,-0.339981043584857,0.339981043584856,0.861136311594053,0.168314227951316,-0.542461986333870,0.168314227951316,-0.625852241617446,-0.542461986333870,-0.625852241617446];
Q_GLegendre_sq_N3 = [-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.861136311594054,-0.339981043584857,0.339981043584856,0.861136311594053,0.861136311594054,0.339981043584857,-0.339981043584856,-0.861136311594053,-0.542461986333870,0.168314227951316,-0.625852241617446,0.168314227951316,-0.625852241617446,-0.542461986333870];
Q_GLegendre_wq_N3 = [0.030198029745131,0.080913081365980,0.080913081365980,0.030198029745131,0.030198029745131,0.080913081365980,0.080913081365980,0.030198029745131,0.030198029745131,0.080913081365980,0.080913081365980,0.030198029745131,0.222222222222222,0.222222222222222,0.222222222222222,0.222222222222222,0.222222222222222,0.222222222222222];
GQ_N3 = (Q_GLegendre_rq_N3, Q_GLegendre_sq_N3, Q_GLegendre_wq_N3)
#N=4
Q_GLegendre_N4 = 22;
Q_GLegendre_rq_N4 = [-0.906179845938664,-0.538469310105683,0.000000000000000,0.538469310105683,0.906179845938664,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.906179845938664,-0.538469310105683,0.000000000000000,0.538469310105683,0.906179845938664,-0.333333333333333,-0.721132537169100,-0.721132537169100,0.442265074338200,-0.123152095118386,-0.123152095118386,-0.753695809763228];
Q_GLegendre_sq_N4 = [-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.906179845938664,-0.538469310105683,0.000000000000000,0.538469310105683,0.906179845938664,0.906179845938664,0.538469310105683,0.000000000000000,-0.538469310105683,-0.906179845938664,-0.333333333333333,-0.721132537169100,0.442265074338200,-0.721132537169100,-0.123152095118386,-0.753695809763228,-0.123152095118386];
Q_GLegendre_wq_N4 = [0.013202630162003,0.041060919360860,0.037074169667907,0.041060919360860,0.013202630162003,0.013202630162003,0.041060919360860,0.037074169667907,0.041060919360860,0.013202630162003,0.013202630162003,0.041060919360860,0.037074169667907,0.041060919360860,0.013202630162003,0.182199822395302,0.210858659241677,0.210858659241677,0.210858659241677,0.249473464579591,0.249473464579591,0.249473464579591];
GQ_N4 = (Q_GLegendre_rq_N4, Q_GLegendre_sq_N4, Q_GLegendre_wq_N4)
#N=5
Q_GLegendre_N5 = 34;
Q_GLegendre_rq_N5 = [-0.932469514203153,-0.661209386466264,-0.238619186083197,0.238619186083197,0.661209386466263,0.932469514203152,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.932469514203153,-0.661209386466264,-0.238619186083197,0.238619186083197,0.661209386466263,0.932469514203152,-0.333333333333333,-0.550932095205921,-0.550932095205921,0.101864190411842,0.611382635695146,-0.744884752682278,0.611382635695146,-0.866497883012868,-0.744884752682278,-0.866497883012868,0.154088161163644,-0.813213547833420,0.154088161163644,-0.340874613330224,-0.813213547833420,-0.340874613330224];
Q_GLegendre_sq_N5 = [-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.932469514203153,-0.661209386466264,-0.238619186083197,0.238619186083197,0.661209386466263,0.932469514203152,0.932469514203153,0.661209386466264,0.238619186083197,-0.238619186083197,-0.661209386466263,-0.932469514203152,-0.333333333333333,-0.550932095205921,0.101864190411842,-0.550932095205921,-0.744884752682278,0.611382635695146,-0.866497883012868,0.611382635695146,-0.866497883012868,-0.744884752682278,-0.813213547833420,0.154088161163644,-0.340874613330224,0.154088161163644,-0.340874613330224,-0.813213547833420];
Q_GLegendre_wq_N5 = [0.006658503841347,0.016313739546177,0.025165807590991,0.025165807590991,0.016313739546177,0.006658503841347,0.006658503841347,0.016313739546177,0.025165807590991,0.025165807590991,0.016313739546177,0.006658503841347,0.006658503841347,0.016313739546177,0.025165807590991,0.025165807590991,0.016313739546177,0.006658503841347,0.209054393645172,0.095225751929022,0.095225751929022,0.095225751929022,0.060381677758004,0.060381677758004,0.060381677758004,0.060381677758004,0.060381677758004,0.060381677758004,0.142358329691442,0.142358329691442,0.142358329691442,0.142358329691442,0.142358329691442,0.142358329691442];
GQ_N5 = (Q_GLegendre_rq_N5, Q_GLegendre_sq_N5, Q_GLegendre_wq_N5)
#N=6
Q_GLegendre_N6 = 42;
Q_GLegendre_rq_N6 = [-0.949107912342759,-0.741531185599395,-0.405845151377397,0.000000000000000,0.405845151377397,0.741531185599396,0.949107912342757,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.949107912342759,-0.741531185599395,-0.405845151377397,0.000000000000000,0.405845151377397,0.741531185599396,0.949107912342757,-0.896694477209711,-0.896694477209711,0.793388954419422,-0.210196994443282,-0.210196994443282,-0.579606011113436,-0.593408612745712,-0.593408612745712,0.186817225491424,-0.273069434548136,-0.867319561224946,-0.273069434548136,0.140388995773082,-0.867319561224946,0.140388995773082,-0.874214255458302,0.514759945734612,-0.874214255458302,-0.640545690276310,0.514759945734612,-0.640545690276310];
Q_GLegendre_sq_N6 = [-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.949107912342759,-0.741531185599395,-0.405845151377397,0.000000000000000,0.405845151377397,0.741531185599396,0.949107912342757,0.949107912342759,0.741531185599395,0.405845151377397,0.000000000000000,-0.405845151377397,-0.741531185599396,-0.949107912342757,-0.896694477209711,0.793388954419422,-0.896694477209711,-0.210196994443282,-0.579606011113436,-0.210196994443282,-0.593408612745712,0.186817225491424,-0.593408612745712,-0.867319561224946,-0.273069434548136,0.140388995773082,-0.273069434548136,0.140388995773082,-0.867319561224946,0.514759945734612,-0.874214255458302,-0.640545690276310,-0.874214255458302,-0.640545690276310,0.514759945734612];
Q_GLegendre_wq_N6 = [0.002716104484117,0.009406550665867,0.013917356931396,0.015633693986814,0.013917356931396,0.009406550665867,0.002716104484117,0.002716104484117,0.009406550665867,0.013917356931396,0.015633693986814,0.013917356931396,0.009406550665867,0.002716104484117,0.002716104484117,0.009406550665867,0.013917356931396,0.015633693986814,0.013917356931396,0.009406550665867,0.002716104484117,0.032730621299045,0.032730621299045,0.032730621299045,0.140317170156143,0.140317170156143,0.140317170156143,0.117404101260644,0.117404101260644,0.117404101260644,0.087979462909756,0.087979462909756,0.087979462909756,0.087979462909756,0.087979462909756,0.087979462909756,0.066271064990874,0.066271064990874,0.066271064990874,0.066271064990874,0.066271064990874,0.066271064990874];
GQ_N6 = (Q_GLegendre_rq_N6, Q_GLegendre_sq_N6, Q_GLegendre_wq_N6)
#N=1
Q_GLobatto_N1 = 6;
Q_GLobatto_rq_N1 = [-1.000000000000000,0.000000000000000,1.000000000000000,0.000000000000000,-1.000000000000000,-1.000000000000000];
Q_GLobatto_sq_N1 = [-1.000000000000000,-1.000000000000000,-1.000000000000000,0.000000000000000,1.000000000000000,0.000000000000000];
Q_GLobatto_wq_N1 = [0.333333333333333,0.333333333333333,0.333333333333333,0.333333333333333,0.333333333333333,0.333333333333333];
GL_N1 = (Q_GLobatto_rq_N1, Q_GLobatto_sq_N1, Q_GLobatto_wq_N1)
#N=2
Q_GLobatto_N2 = 10;
Q_GLobatto_rq_N2 = [-1.000000000000000,-0.447213595499958,0.447213595499958,1.000000000000000,0.447213595499958,-0.447213595499958,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.333333333333333];
Q_GLobatto_sq_N2 = [-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.447213595499958,0.447213595499958,1.000000000000000,0.447213595499958,-0.447213595499958,-0.333333333333333];
Q_GLobatto_wq_N2 = [0.033333333333333,0.166666666666667,0.166666666666667,0.033333333333333,0.166666666666667,0.166666666666667,0.033333333333333,0.166666666666667,0.166666666666667,0.900000000000000];
GL_N2 = (Q_GLobatto_rq_N2, Q_GLobatto_sq_N2, Q_GLobatto_wq_N2)
#N=3
Q_GLobatto_N3 = 15;
Q_GLobatto_rq_N3 = [-1.000000000000000,-0.654653670707977,0.000000000000000,0.654653670707977,1.000000000000000,0.654653670707977,0.000000000000000,-0.654653670707977,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.585309648672816,-0.585309648672816,0.170619297345632];
Q_GLobatto_sq_N3 = [-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.654653670707977,0.000000000000000,0.654653670707977,1.000000000000000,0.654653670707977,0.000000000000000,-0.654653670707977,-0.585309648672816,0.170619297345632,-0.585309648672816];
Q_GLobatto_wq_N3 = [0.014698618394802,0.056443964486202,0.097526003618641,0.056443964486202,0.014698618394802,0.056443964486202,0.097526003618641,0.056443964486202,0.014698618394802,0.056443964486202,0.097526003618641,0.056443964486202,0.441554115680819,0.441554115680819,0.441554115680819];
GL_N3 = (Q_GLobatto_rq_N3, Q_GLobatto_sq_N3, Q_GLobatto_wq_N3)
#N=4
Q_GLobatto_N4 = 24;
Q_GLobatto_rq_N4 = [-1.000000000000000,-0.765055323929465,-0.285231516480645,0.285231516480645,0.765055323929466,1.000000000000000,0.765055323929465,0.285231516480645,-0.285231516480645,-0.765055323929466,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.153862961347394,-0.153862961347394,-0.692274077305212,0.462899656119068,-0.801694541056790,0.462899656119068,-0.661205115062278,-0.801694541056790,-0.661205115062278];
Q_GLobatto_sq_N4 = [-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.765055323929465,-0.285231516480645,0.285231516480645,0.765055323929466,1.000000000000000,0.765055323929465,0.285231516480645,-0.285231516480645,-0.765055323929466,-0.153862961347394,-0.692274077305212,-0.153862961347394,-0.801694541056790,0.462899656119068,-0.661205115062278,0.462899656119068,-0.661205115062278,-0.801694541056790];
Q_GLobatto_wq_N4 = [0.007531613457659,0.021086271203965,0.045504525339300,0.045504525339300,0.021086271203965,0.007531613457659,0.021086271203965,0.045504525339300,0.045504525339300,0.021086271203965,0.007531613457659,0.021086271203965,0.045504525339300,0.045504525339300,0.021086271203965,0.304805308340987,0.304805308340987,0.304805308340987,0.110574075890745,0.110574075890745,0.110574075890745,0.110574075890745,0.110574075890745,0.110574075890745];
GL_N4 = (Q_GLobatto_rq_N4, Q_GLobatto_sq_N4, Q_GLobatto_wq_N4)
#N=5
Q_GLobatto_N5 = 33;
Q_GLobatto_rq_N5 = [-1.000000000000000,-0.830223896278567,-0.468848793470714,0.000000000000000,0.468848793470714,0.830223896278567,1.000000000000000,0.830223896278567,0.468848793470714,0.000000000000000,-0.468848793470714,-0.830223896278567,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.481290025456663,-0.481290025456663,-0.037419949086674,-0.085394918585138,-0.085394918585138,-0.829210162829724,-0.826342421973033,-0.826342421973033,0.652684843946066,-0.834845895799092,0.320885646218972,-0.834845895799092,-0.486039750419880,0.320885646218972,-0.486039750419880];
Q_GLobatto_sq_N5 = [-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.830223896278567,-0.468848793470714,0.000000000000000,0.468848793470714,0.830223896278567,1.000000000000000,0.830223896278567,0.468848793470714,0.000000000000000,-0.468848793470714,-0.830223896278567,-0.481290025456663,-0.037419949086674,-0.481290025456663,-0.085394918585138,-0.829210162829724,-0.085394918585138,-0.826342421973033,0.652684843946066,-0.826342421973033,0.320885646218972,-0.834845895799092,-0.486039750419880,-0.834845895799092,-0.486039750419880,0.320885646218972];
Q_GLobatto_wq_N5 = [0.002314778534840,0.013700973144398,0.020117879303463,0.023426958130206,0.020117879303463,0.013700973144398,0.002314778534840,0.013700973144398,0.020117879303463,0.023426958130206,0.020117879303463,0.013700973144398,0.002314778534840,0.013700973144398,0.020117879303463,0.023426958130206,0.020117879303463,0.013700973144398,0.176274129326032,0.176274129326032,0.176274129326032,0.113088210164246,0.113088210164246,0.113088210164246,0.076346374759077,0.076346374759077,0.076346374759077,0.103789255428272,0.103789255428272,0.103789255428272,0.103789255428272,0.103789255428272,0.103789255428272];
GL_N5 = (Q_GLobatto_rq_N5, Q_GLobatto_sq_N5, Q_GLobatto_wq_N5)
#N=6
Q_GLobatto_N6 = 40;
Q_GLobatto_rq_N6 = [-1.000000000000000,-0.871740148509606,-0.591700181433142,-0.209299217902479,0.209299217902479,0.591700181433143,0.871740148509605,1.000000000000000,0.871740148509606,0.591700181433142,0.209299217902479,-0.209299217902479,-0.591700181433143,-0.871740148509605,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.333333333333333,-0.832518933286299,-0.832518933286299,0.665037866572598,-0.018986264273055,-0.018986264273055,-0.962027471453890,-0.130967625384817,-0.130967625384817,-0.738064749230366,-0.591781036586590,-0.591781036586590,0.183562073173180,-0.465390307260470,0.362915057330052,-0.465390307260470,-0.897524750069582,0.362915057330052,-0.897524750069582];
Q_GLobatto_sq_N6 = [-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-1.000000000000000,-0.871740148509606,-0.591700181433142,-0.209299217902479,0.209299217902479,0.591700181433143,0.871740148509605,1.000000000000000,0.871740148509606,0.591700181433142,0.209299217902479,-0.209299217902479,-0.591700181433143,-0.871740148509605,-0.333333333333333,-0.832518933286299,0.665037866572598,-0.832518933286299,-0.018986264273055,-0.962027471453890,-0.018986264273055,-0.130967625384817,-0.738064749230366,-0.130967625384817,-0.591781036586590,0.183562073173180,-0.591781036586590,0.362915057330052,-0.465390307260470,-0.897524750069582,-0.465390307260470,-0.897524750069582,0.362915057330052];
Q_GLobatto_wq_N6 = [0.000886028699507,0.010480045161785,0.011623548297153,0.001211835242253,0.001211835242253,0.011623548297153,0.010480045161785,0.000886028699507,0.010480045161785,0.011623548297153,0.001211835242253,0.001211835242253,0.011623548297153,0.010480045161785,0.000886028699507,0.010480045161785,0.011623548297153,0.001211835242253,0.001211835242253,0.011623548297153,0.010480045161785,0.171525496248007,0.074048529130209,0.074048529130209,0.074048529130209,0.049562169655280,0.049562169655280,0.049562169655280,0.136384241605594,0.136384241605594,0.136384241605594,0.139937548771021,0.139937548771021,0.139937548771021,0.081021062993335,0.081021062993335,0.081021062993335,0.081021062993335,0.081021062993335,0.081021062993335];
GL_N6 = (Q_GLobatto_rq_N6, Q_GLobatto_sq_N6, Q_GLobatto_wq_N6)
#N=1
Tri_diage_N1 = 7;
Tri_diage_rq_N1 = [0.0, 1.0, 0.0, 0.5, 0.5, 0.0, 0.3333333333333333];
Tri_diage_sq_N1 = [0.0, 0.0, 1.0, 0.0, 0.5, 0.5, 0.3333333333333333];
Tri_diage_wq_N1 = [0.024999999999999998, 0.024999999999999998, 0.024999999999999998, 0.06666666666666667, 0.06666666666666667, 0.06666666666666667, 0.22500000000000006];
Tri_diage_N1 = (2*Tri_diage_rq_N1 .-1, 2*Tri_diage_sq_N1 .-1, 4*Tri_diage_wq_N1)
#N=2
Tri_diage_N2 = 12;
Tri_diage_rq_N2 = [0.0, 1.0, 0.0, 0.21285435711180825, 0.21285435711180825, 0.5742912857763836, 0.27639320225002106, 0.7236067977499789, 0.7236067977499789, 0.27639320225002106, 0.0, 0.0];
Tri_diage_sq_N2 = [0.0, 0.0, 1.0, 0.5742912857763836, 0.21285435711180825, 0.21285435711180825, 0.0, 0.0, 0.27639320225002106, 0.7236067977499789, 0.7236067977499789, 0.27639320225002106];
Tri_diage_wq_N2 = [0.006261126504899741, 0.006261126504899741, 0.006261126504899741, 0.10675793966098839, 0.10675793966098839, 0.10675793966098839, 0.026823800250389242, 0.026823800250389242, 0.026823800250389242, 0.026823800250389242, 0.026823800250389242, 0.026823800250389242];
Tri_diage_N2 = (2*Tri_diage_rq_N2 .-1, 2*Tri_diage_sq_N2 .-1, 4*Tri_diage_wq_N2)
#N=3
Tri_diage_N3 = 18;
Tri_diage_rq_N3 = [0.0, 1.0, 0.0, 0.5, 0.5, 0.0, 0.4243860251718814, 0.4243860251718814, 0.1512279496562372, 0.14200508409677795, 0.14200508409677795, 0.7159898318064442, 0.17267316464601146, 0.8273268353539885, 0.8273268353539885, 0.17267316464601146, 0.0, 0.0];
Tri_diage_sq_N3 = [0.0, 0.0, 1.0, 0.0, 0.5, 0.5, 0.1512279496562372, 0.4243860251718814, 0.4243860251718814, 0.7159898318064442, 0.14200508409677795, 0.14200508409677795, 0.0, 0.0, 0.17267316464601146, 0.8273268353539885, 0.8273268353539885, 0.17267316464601146];
Tri_diage_wq_N3 = [0.0022825661430496253, 0.0022825661430496253, 0.0022825661430496253, 0.015504052643022513, 0.015504052643022513, 0.015504052643022513, 0.07467669469983994, 0.07467669469983994, 0.07467669469983994, 0.051518167995569394, 0.051518167995569394, 0.051518167995569394, 0.011342592592592586, 0.011342592592592586, 0.011342592592592586, 0.011342592592592586, 0.011342592592592586, 0.011342592592592586];
Tri_diage_N3 = (2*Tri_diage_rq_N3 .-1, 2*Tri_diage_sq_N3 .-1, 4*Tri_diage_wq_N3)
#N=4
Tri_diage_N4 = 27;
Tri_diage_rq_N4 = [0.0, 1.0, 0.0, 0.10367750814280517, 0.10367750814280517, 0.7926449837143896, 0.2653313804842097, 0.2653313804842097, 0.46933723903158064, 0.35738424175967753, 0.6426157582403225, 0.6426157582403225, 0.35738424175967753, 0.0, 0.0, 0.11747233803526758, 0.8825276619647324, 0.8825276619647324, 0.11747233803526758, 0.0, 0.0, 0.5870855671333673, 0.3246404722650515, 0.3246404722650515, 0.5870855671333673, 0.0882739606015811, 0.0882739606015811];
Tri_diage_sq_N4 = [0.0, 0.0, 1.0, 0.7926449837143896, 0.10367750814280517, 0.10367750814280517, 0.46933723903158064, 0.2653313804842097, 0.2653313804842097, 0.0, 0.0, 0.35738424175967753, 0.6426157582403225, 0.6426157582403225, 0.35738424175967753, 0.0, 0.0, 0.11747233803526758, 0.8825276619647324, 0.8825276619647324, 0.11747233803526758, 0.0882739606015811, 0.0882739606015811, 0.5870855671333673, 0.3246404722650515, 0.3246404722650515, 0.5870855671333673];
Tri_diage_wq_N4 = [0.001090393904993471, 0.001090393904993471, 0.001090393904993471, 0.02839719066391149, 0.02839719066391149, 0.02839719066391149, 0.03996004802785181, 0.03996004802785181, 0.03996004802785181, 0.0069669428714637, 0.0069669428714637, 0.0069669428714637, 0.0069669428714637, 0.0069669428714637, 0.0069669428714637, 0.005519747637357106, 0.005519747637357106, 0.005519747637357106, 0.005519747637357106, 0.005519747637357106, 0.005519747637357106, 0.03612282652613417, 0.03612282652613417, 0.03612282652613417, 0.03612282652613417, 0.03612282652613417, 0.03612282652613417];
Tri_diage_N4 = (2*Tri_diage_rq_N4 .-1, 2*Tri_diage_sq_N4 .-1, 4*Tri_diage_wq_N4)
GQ_SBP = (GQ_N1, GQ_N2, GQ_N3, GQ_N4, GQ_N5, GQ_N6);
GL_SBP = (GL_N1, GL_N2, GL_N3, GL_N4, GL_N5, GL_N6);
Tri_diage = (Tri_diage_N1, Tri_diage_N2, Tri_diage_N3, Tri_diage_N4);