Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/operators/standard_diagonal.toml @ 813:cdc2b5ebf7cb operator_storage_array_of_table

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Review: Suggested changes in standard_diagonal.toml
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 22 Sep 2021 13:09:46 +0200
parents f682e4fe3cef
children 454ba1efa644
comparison
equal deleted inserted replaced
812:8c977b336803 813:cdc2b5ebf7cb
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 = "A paper a long time ago in a galaxy far far away."
6 # Review:
7 # Suggested change:
8 # "A paper a long time ago in a galaxy far far away." -->
9 # "
10 # Ken Mattsson, Jan Nordström,
11 # Summation by parts operators for finite difference approximations of second derivatives,
12 # Journal of Computational Physics,
13 # Volume 199, Issue 2,
14 # 2004,
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.
6 20
7 [[stencil_set]] 21 [[stencil_set]]
8 22
9 order = 2 23 order = 2
10 24
28 42
29 order = 4 43 order = 4
30 H.inner = "1" 44 H.inner = "1"
31 H.closure = ["17/48", "59/48", "43/48", "49/48"] 45 H.closure = ["17/48", "59/48", "43/48", "49/48"]
32 46
47 # Review:
48 # Add missing 4th order accurate D1 operator
49 # D1.inner_stencil = ["1/12","-2/3","0","2/3","-1/12"]
50 # D1.closure_stencils = [
51 # {s = [ "-24/17", "59/34", "-4/17", "-3/34", "0", "0"], c = 1},
52 # {s = [ "-1/2", "0", "1/2", "0", "0", "0"], c = 2},
53 # {s = [ "4/43", "-59/86", "0", "59/86", "-4/43", "0"], c = 3},
54 # {s = [ "3/98", "0", "-59/98", "0", "32/49", "-4/49"], c = 4},
55 # ]
56
33 D2.inner_stencil = ["-1/12","4/3","-5/2","4/3","-1/12"] 57 D2.inner_stencil = ["-1/12","4/3","-5/2","4/3","-1/12"]
34 D2.closure_stencils = [ 58 D2.closure_stencils = [
35 {s = [ "2", "-5", "4", "-1", "0", "0"], c = 1}, 59 {s = [ "2", "-5", "4", "-1", "0", "0"], c = 1},
36 {s = [ "1", "-2", "1", "0", "0", "0"], c = 2}, 60 {s = [ "1", "-2", "1", "0", "0", "0"], c = 2},
37 {s = [ "-4/43", "59/43", "-110/43", "59/43", "-4/43", "0"], c = 3}, 61 {s = [ "-4/43", "59/43", "-110/43", "59/43", "-4/43", "0"], c = 3},