Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/operators/standard_diagonal.toml @ 833:454ba1efa644 operator_storage_array_of_table

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Review: Accept changes
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 12 Jan 2022 16:03:48 +0100
parents cdc2b5ebf7cb
children fe8fe3f01162
comparison
equal deleted inserted replaced
832:00f6bbdcd73a 833:454ba1efa644
1 [meta] 1 [meta]
2 authors = "Ken Mattson" 2 authors = "Ken Mattson"
3 description = "Standard operators for equidistant grids" 3 description = "Standard operators for equidistant grids"
4 type = "equidistant" 4 type = "equidistant"
5 cite = "A paper a long time ago in a galaxy far far away." 5 cite = """
6 # Review: 6 Ken Mattsson, Jan Nordström,
7 # Suggested change: 7 Summation by parts operators for finite difference approximations of second derivatives,
8 # "A paper a long time ago in a galaxy far far away." --> 8 Journal of Computational Physics,
9 # " 9 Volume 199, Issue 2,
10 # Ken Mattsson, Jan Nordström, 10 2004,
11 # Summation by parts operators for finite difference approximations of second derivatives, 11 Pages 503-540,
12 # Journal of Computational Physics, 12 ISSN 0021-9991,
13 # Volume 199, Issue 2, 13 https://doi.org/10.1016/j.jcp.2004.03.001.
14 # 2004, 14 """
15 # Pages 503-540,
16 # ISSN 0021-9991,
17 # https://doi.org/10.1016/j.jcp.2004.03.001.
18 # "
19 # or perhaps a shorter version.
20 15
21 [[stencil_set]] 16 [[stencil_set]]
22 17
23 order = 2 18 order = 2
24 19
42 37
43 order = 4 38 order = 4
44 H.inner = "1" 39 H.inner = "1"
45 H.closure = ["17/48", "59/48", "43/48", "49/48"] 40 H.closure = ["17/48", "59/48", "43/48", "49/48"]
46 41
47 # Review: 42 D1.inner_stencil = ["1/12","-2/3","0","2/3","-1/12"]
48 # Add missing 4th order accurate D1 operator 43 D1.closure_stencils = [
49 # D1.inner_stencil = ["1/12","-2/3","0","2/3","-1/12"] 44 {s = [ "-24/17", "59/34", "-4/17", "-3/34", "0", "0"], c = 1},
50 # D1.closure_stencils = [ 45 {s = [ "-1/2", "0", "1/2", "0", "0", "0"], c = 2},
51 # {s = [ "-24/17", "59/34", "-4/17", "-3/34", "0", "0"], c = 1}, 46 {s = [ "4/43", "-59/86", "0", "59/86", "-4/43", "0"], c = 3},
52 # {s = [ "-1/2", "0", "1/2", "0", "0", "0"], c = 2}, 47 {s = [ "3/98", "0", "-59/98", "0", "32/49", "-4/49"], c = 4},
53 # {s = [ "4/43", "-59/86", "0", "59/86", "-4/43", "0"], c = 3}, 48 ]
54 # {s = [ "3/98", "0", "-59/98", "0", "32/49", "-4/49"], c = 4},
55 # ]
56 49
57 D2.inner_stencil = ["-1/12","4/3","-5/2","4/3","-1/12"] 50 D2.inner_stencil = ["-1/12","4/3","-5/2","4/3","-1/12"]
58 D2.closure_stencils = [ 51 D2.closure_stencils = [
59 {s = [ "2", "-5", "4", "-1", "0", "0"], c = 1}, 52 {s = [ "2", "-5", "4", "-1", "0", "0"], c = 1},
60 {s = [ "1", "-2", "1", "0", "0", "0"], c = 2}, 53 {s = [ "1", "-2", "1", "0", "0", "0"], c = 2},