comparison src/SbpOperators/operators/standard_diagonal.toml @ 1651:707fc9761c2b feature/sbp_operators/laplace_curvilinear

Merge feature/grids/manifolds
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 26 Jun 2024 12:47:26 +0200
parents 1388149b54ad
children
comparison
equal deleted inserted replaced
1648:b7dcd3dd3181 1651:707fc9761c2b
32 D2.inner_stencil = ["1", "-2", "1"] 32 D2.inner_stencil = ["1", "-2", "1"]
33 D2.closure_stencils = [ 33 D2.closure_stencils = [
34 {s = ["1", "-2", "1"], c = 1}, 34 {s = ["1", "-2", "1"], c = 1},
35 ] 35 ]
36 36
37 D2.positivity = {theta_M = "0.3636363636", theta_R = "1.000000538455350", m_b = "2"}
38
37 D2variable.inner_stencil = [["1/2", "1/2", "0"],[ "-1/2", "-1", "-1/2"],["0", "1/2", "1/2"]] 39 D2variable.inner_stencil = [["1/2", "1/2", "0"],[ "-1/2", "-1", "-1/2"],["0", "1/2", "1/2"]]
38 D2variable.closure_stencils = [ 40 D2variable.closure_stencils = [
39 {s = [["2", "-1", "0"],["-3", "1", "0"],["1","0","0"]], c = 1}, 41 {s = [["2", "-1", "0"],["-3", "1", "0"],["1","0","0"]], c = 1},
40 ] 42 ]
43
41 44
42 [[stencil_set]] 45 [[stencil_set]]
43 46
44 order = 4 47 order = 4
45 48
62 {s = [ "2", "-5", "4", "-1", "0", "0"], c = 1}, 65 {s = [ "2", "-5", "4", "-1", "0", "0"], c = 1},
63 {s = [ "1", "-2", "1", "0", "0", "0"], c = 2}, 66 {s = [ "1", "-2", "1", "0", "0", "0"], c = 2},
64 {s = [ "-4/43", "59/43", "-110/43", "59/43", "-4/43", "0"], c = 3}, 67 {s = [ "-4/43", "59/43", "-110/43", "59/43", "-4/43", "0"], c = 3},
65 {s = [ "-1/49", "0", "59/49", "-118/49", "64/49", "-4/49"], c = 4}, 68 {s = [ "-1/49", "0", "59/49", "-118/49", "64/49", "-4/49"], c = 4},
66 ] 69 ]
70
71 D2.positivity = {theta_M = "0.2505765857", theta_R = "0.577587500088313", m_b = "4"}
67 72
68 D2variable.inner_stencil = [ 73 D2variable.inner_stencil = [
69 ["-1/8", "1/6", "-1/8", "0", "0" ], 74 ["-1/8", "1/6", "-1/8", "0", "0" ],
70 [ "1/6", "1/2", "1/2", "1/6", "0" ], 75 [ "1/6", "1/2", "1/2", "1/6", "0" ],
71 ["-1/24", "-5/6", "-3/4", "-5/6", "-1/24"], 76 ["-1/24", "-5/6", "-3/4", "-5/6", "-1/24"],
134 [ "0", "0", "0", "0", "0", "-1/8", "1/6", "-1/8" ], 139 [ "0", "0", "0", "0", "0", "-1/8", "1/6", "-1/8" ],
135 ]} 140 ]}
136 ] 141 ]
137 142
138 143
139
140 [[stencil_set]] 144 [[stencil_set]]
141 145
142 order = 6 146 order = 6
143 147
144 H.inner = "1" 148 H.inner = "1"
148 152
149 153
150 154
151 e.closure = ["1"] 155 e.closure = ["1"]
152 d1.closure = ["-25/12", "4", "-3", "4/3", "-1/4"] 156 d1.closure = ["-25/12", "4", "-3", "4/3", "-1/4"]
157
158 D2.positivity = {theta_M = "0.1878687080", theta_R = "0.3697", m_b = "7"}