| 1 |
|
|---|
| 2 | ###########################################################
|
|---|
| 3 |
|
|---|
| 4 | # GridCreate: ndim nparts
|
|---|
| 5 | GridCreate: 3 8
|
|---|
| 6 |
|
|---|
| 7 | # GridSetExtents: part ilower(ndim) iupper(ndim)
|
|---|
| 8 | GridSetExtents: 0 (1- 1- 1-) (4+ 3+ 4+)
|
|---|
| 9 | GridSetExtents: 1 (1- 1- 1-) (4+ 3+ 4+)
|
|---|
| 10 | GridSetExtents: 2 (1- 1- 1-) (4+ 3+ 4+)
|
|---|
| 11 | GridSetExtents: 3 (1- 1- 1-) (4+ 3+ 4+)
|
|---|
| 12 | GridSetExtents: 4 (1- 1- 1-) (2+ 6+ 4+)
|
|---|
| 13 | GridSetExtents: 5 (1- 1- 1-) (6+ 2+ 4+)
|
|---|
| 14 | GridSetExtents: 6 (1- 1- 1-) (6+ 2+ 4+)
|
|---|
| 15 | GridSetExtents: 7 (1- 1- 1-) (4+ 3+ 4+)
|
|---|
| 16 |
|
|---|
| 17 | # GridSetVariables: part nvars vartypes[nvars]
|
|---|
| 18 | # CELL = 0
|
|---|
| 19 | # NODE = 1
|
|---|
| 20 | # XFACE = 2
|
|---|
| 21 | # YFACE = 3
|
|---|
| 22 | # ZFACE = 4
|
|---|
| 23 | # XEDGE = 5
|
|---|
| 24 | # YEDGE = 6
|
|---|
| 25 | # ZEDGE = 7
|
|---|
| 26 | GridSetVariables: 0 1 [0]
|
|---|
| 27 | GridSetVariables: 1 1 [0]
|
|---|
| 28 | GridSetVariables: 2 1 [0]
|
|---|
| 29 | GridSetVariables: 3 1 [0]
|
|---|
| 30 | GridSetVariables: 4 1 [0]
|
|---|
| 31 | GridSetVariables: 5 1 [0]
|
|---|
| 32 | GridSetVariables: 6 1 [0]
|
|---|
| 33 | GridSetVariables: 7 1 [0]
|
|---|
| 34 |
|
|---|
| 35 | # GridAddVariables: part index(ndim) vartype
|
|---|
| 36 |
|
|---|
| 37 | # GridSetNeighborBox: part ilower(ndim) iupper(ndim) \
|
|---|
| 38 | # nbor_part nbor_ilower(ndim) nbor_iupper(ndim) \
|
|---|
| 39 | # index_map[ndim]
|
|---|
| 40 | GridSetNeighborBox: 0 (5- 1- 1-) (5- 3+ 4+) 1 (1- 1- 1-) (1- 3+ 4+) [0 1 2]
|
|---|
| 41 | GridSetNeighborBox: 1 (0+ 1- 1-) (0+ 3+ 4+) 0 (4+ 1- 1-) (4+ 3+ 4+) [0 1 2]
|
|---|
| 42 |
|
|---|
| 43 | GridSetNeighborBox: 1 (5- 1- 1-) (5- 3+ 4+) 2 (1- 1- 1-) (1- 3+ 4+) [0 1 2]
|
|---|
| 44 | GridSetNeighborBox: 2 (0+ 1- 1-) (0+ 3+ 4+) 1 (4+ 1- 1-) (4+ 3+ 4+) [0 1 2]
|
|---|
| 45 |
|
|---|
| 46 | GridSetNeighborBox: 2 (5- 1- 1-) (5- 3+ 4+) 3 (1- 1- 1-) (1- 3+ 4+) [0 1 2]
|
|---|
| 47 | GridSetNeighborBox: 3 (0+ 1- 1-) (0+ 3+ 4+) 2 (4+ 1- 1-) (4+ 3+ 4+) [0 1 2]
|
|---|
| 48 |
|
|---|
| 49 | GridSetNeighborBox: 3 (5- 1- 1-) (5- 4+ 4+) 4 (1- 1- 1-) (1- 4+ 4+) [0 1 2]
|
|---|
| 50 | GridSetNeighborBox: 4 (0+ 1- 1-) (0+ 3+ 4+) 3 (4+ 1- 1-) (4+ 3+ 4+) [0 1 2]
|
|---|
| 51 |
|
|---|
| 52 | GridSetNeighborBox: 4 (3- 5- 1-) (3- 6+ 4+) 5 (1- 1- 1-) (1- 2+ 4+) [0 1 2]
|
|---|
| 53 | GridSetNeighborBox: 5 (0+ 1- 1-) (0+ 2+ 4+) 4 (2+ 5- 1-) (2+ 6+ 4+) [0 1 2]
|
|---|
| 54 |
|
|---|
| 55 | GridSetNeighborBox: 4 (3- 1- 1-) (3- 2+ 4+) 6 (1- 1- 1-) (1- 2+ 4+) [0 1 2]
|
|---|
| 56 | GridSetNeighborBox: 6 (0+ 1- 1-) (0+ 2+ 4+) 4 (2+ 1- 1-) (2+ 2+ 4+) [0 1 2]
|
|---|
| 57 |
|
|---|
| 58 | GridSetNeighborBox: 6 (1- 4- 1-) (2+ 4- 4+) 7 (3- 1- 1-) (4+ 1- 4+) [0 1 2]
|
|---|
| 59 | GridSetNeighborBox: 7 (3- 0+ 1-) (4+ 0+ 4+) 6 (1- 3+ 1-) (2+ 3+ 4+) [0 1 2]
|
|---|
| 60 |
|
|---|
| 61 | # GridSetPeriodic: part periodic[ndim]
|
|---|
| 62 |
|
|---|
| 63 | ###########################################################
|
|---|
| 64 |
|
|---|
| 65 | # StencilCreate: nstencils sizes[nstencils]
|
|---|
| 66 | StencilCreate: 7 [9 9 9 9 8 7 7]
|
|---|
| 67 |
|
|---|
| 68 | # StencilSetEntry: stencil_num entry offset[ndim] var value
|
|---|
| 69 | StencilSetEntry: 0 0 [ 0 0 0] 0 6.0
|
|---|
| 70 | StencilSetEntry: 0 1 [-1 0 0] 0 0.0
|
|---|
| 71 | StencilSetEntry: 0 2 [ 1 0 0] 0 0.0
|
|---|
| 72 | StencilSetEntry: 0 3 [ 0 -1 0] 0 -1.0
|
|---|
| 73 | StencilSetEntry: 0 4 [ 0 1 0] 0 -1.0
|
|---|
| 74 | StencilSetEntry: 0 5 [ 0 0 -1] 0 -1.0
|
|---|
| 75 | StencilSetEntry: 0 6 [ 0 0 1] 0 -1.0
|
|---|
| 76 | StencilSetEntry: 0 7 [-1 1 0] 0 -1.0
|
|---|
| 77 | StencilSetEntry: 0 8 [ 1 -1 0] 0 -1.0
|
|---|
| 78 |
|
|---|
| 79 | # StencilSetEntry: stencil_num entry offset[ndim] var value
|
|---|
| 80 | StencilSetEntry: 1 0 [ 0 0 0] 0 600.0
|
|---|
| 81 | StencilSetEntry: 1 1 [-1 0 0] 0 0.0
|
|---|
| 82 | StencilSetEntry: 1 2 [ 1 0 0] 0 0.0
|
|---|
| 83 | StencilSetEntry: 1 3 [ 0 -1 0] 0 -100.0
|
|---|
| 84 | StencilSetEntry: 1 4 [ 0 1 0] 0 -100.0
|
|---|
| 85 | StencilSetEntry: 1 5 [ 0 0 -1] 0 -100.0
|
|---|
| 86 | StencilSetEntry: 1 6 [ 0 0 1] 0 -100.0
|
|---|
| 87 | StencilSetEntry: 1 7 [-1 -1 0] 0 -100.0
|
|---|
| 88 | StencilSetEntry: 1 8 [ 1 1 0] 0 -100.0
|
|---|
| 89 |
|
|---|
| 90 | # StencilSetEntry: stencil_num entry offset[ndim] var value
|
|---|
| 91 | StencilSetEntry: 2 0 [ 0 0 0] 0 0.6
|
|---|
| 92 | StencilSetEntry: 2 1 [-1 0 0] 0 0.0
|
|---|
| 93 | StencilSetEntry: 2 2 [ 1 0 0] 0 0.0
|
|---|
| 94 | StencilSetEntry: 2 3 [ 0 -1 0] 0 -0.1
|
|---|
| 95 | StencilSetEntry: 2 4 [ 0 1 0] 0 -0.1
|
|---|
| 96 | StencilSetEntry: 2 5 [ 0 0 -1] 0 -0.1
|
|---|
| 97 | StencilSetEntry: 2 6 [ 0 0 1] 0 -0.1
|
|---|
| 98 | StencilSetEntry: 2 7 [-1 1 0] 0 -0.1
|
|---|
| 99 | StencilSetEntry: 2 8 [ 1 -1 0] 0 -0.1
|
|---|
| 100 |
|
|---|
| 101 | # StencilSetEntry: stencil_num entry offset[ndim] var value
|
|---|
| 102 | StencilSetEntry: 3 0 [ 0 0 0] 0 6.0
|
|---|
| 103 | StencilSetEntry: 3 1 [-1 0 0] 0 0.0
|
|---|
| 104 | StencilSetEntry: 3 2 [ 1 0 0] 0 0.0
|
|---|
| 105 | StencilSetEntry: 3 3 [ 0 -1 0] 0 -1.0
|
|---|
| 106 | StencilSetEntry: 3 4 [ 0 1 0] 0 -1.0
|
|---|
| 107 | StencilSetEntry: 3 5 [ 0 0 -1] 0 -1.0
|
|---|
| 108 | StencilSetEntry: 3 6 [ 0 0 1] 0 -1.0
|
|---|
| 109 | StencilSetEntry: 3 7 [-1 -1 0] 0 -1.0
|
|---|
| 110 | StencilSetEntry: 3 8 [ 1 1 0] 0 -1.0
|
|---|
| 111 |
|
|---|
| 112 | # StencilSetEntry: stencil_num entry offset[ndim] var value
|
|---|
| 113 | StencilSetEntry: 4 0 [ 0 0 0] 0 6.0
|
|---|
| 114 | StencilSetEntry: 4 1 [-1 0 0] 0 -1.0
|
|---|
| 115 | StencilSetEntry: 4 2 [ 1 0 0] 0 -1.0
|
|---|
| 116 | StencilSetEntry: 4 3 [ 0 -1 0] 0 -1.0
|
|---|
| 117 | StencilSetEntry: 4 4 [ 0 1 0] 0 -1.0
|
|---|
| 118 | StencilSetEntry: 4 5 [ 0 0 -1] 0 -1.0
|
|---|
| 119 | StencilSetEntry: 4 6 [ 0 0 1] 0 -1.0
|
|---|
| 120 | StencilSetEntry: 4 7 [-1 -1 0] 0 0.0
|
|---|
| 121 |
|
|---|
| 122 | # StencilSetEntry: stencil_num entry offset[ndim] var value
|
|---|
| 123 | StencilSetEntry: 5 0 [ 0 0 0] 0 6.0
|
|---|
| 124 | StencilSetEntry: 5 1 [-1 0 0] 0 -1.0
|
|---|
| 125 | StencilSetEntry: 5 2 [ 1 0 0] 0 -1.0
|
|---|
| 126 | StencilSetEntry: 5 3 [ 0 -1 0] 0 -1.0
|
|---|
| 127 | StencilSetEntry: 5 4 [ 0 1 0] 0 -1.0
|
|---|
| 128 | StencilSetEntry: 5 5 [ 0 0 -1] 0 -1.0
|
|---|
| 129 | StencilSetEntry: 5 6 [ 0 0 1] 0 -1.0
|
|---|
| 130 |
|
|---|
| 131 | # StencilSetEntry: stencil_num entry offset[ndim] var value
|
|---|
| 132 | StencilSetEntry: 6 0 [ 0 0 0] 0 20.4
|
|---|
| 133 | StencilSetEntry: 6 1 [-1 0 0] 0 -0.1
|
|---|
| 134 | StencilSetEntry: 6 2 [ 1 0 0] 0 -0.1
|
|---|
| 135 | StencilSetEntry: 6 3 [ 0 -1 0] 0 -10.0
|
|---|
| 136 | StencilSetEntry: 6 4 [ 0 1 0] 0 -10.0
|
|---|
| 137 | StencilSetEntry: 6 5 [ 0 0 -1] 0 -0.1
|
|---|
| 138 | StencilSetEntry: 6 6 [ 0 0 1] 0 -0.1
|
|---|
| 139 |
|
|---|
| 140 | ###########################################################
|
|---|
| 141 |
|
|---|
| 142 | # GraphSetStencil: part var stencil_num
|
|---|
| 143 | GraphSetStencil: 0 0 0
|
|---|
| 144 | GraphSetStencil: 1 0 1
|
|---|
| 145 | GraphSetStencil: 2 0 2
|
|---|
| 146 | GraphSetStencil: 3 0 3
|
|---|
| 147 | GraphSetStencil: 4 0 4
|
|---|
| 148 | GraphSetStencil: 5 0 5
|
|---|
| 149 | GraphSetStencil: 6 0 5
|
|---|
| 150 | GraphSetStencil: 7 0 6
|
|---|
| 151 |
|
|---|
| 152 | # GraphAddEntries: \
|
|---|
| 153 | # part ilower(ndim) iupper(ndim) stride[ndim] var \
|
|---|
| 154 | # to_part to_ilower(ndim) to_iupper(ndim) to_stride[ndim] to_var \
|
|---|
| 155 | # index_map[ndim] entry value
|
|---|
| 156 |
|
|---|
| 157 | ###########################################################
|
|---|
| 158 |
|
|---|
| 159 | # MatrixCreate:
|
|---|
| 160 |
|
|---|
| 161 | # MatrixSetValues: \
|
|---|
| 162 | # part ilower(ndim) iupper(ndim) stride[ndim] var entry value
|
|---|
| 163 |
|
|---|
| 164 | ##### Dirichlet boundary conditions are set automatically
|
|---|
| 165 |
|
|---|
| 166 | ##### Fix up stencils along part boundaries
|
|---|
| 167 |
|
|---|
| 168 | MatrixSetValues: 0 (4+ 1- 1-) (4+ 3+ 4+) [1 1 1] 0 8 0.0
|
|---|
| 169 | MatrixSetValues: 0 (4+ 1- 1-) (4+ 3+ 4+) [1 1 1] 0 2 -100.0
|
|---|
| 170 | MatrixSetValues: 0 (4+ 1- 1-) (4+ 3+ 4+) [1 1 1] 0 0 105.0
|
|---|
| 171 |
|
|---|
| 172 | MatrixSetValues: 1 (1- 1- 1-) (1- 3+ 4+) [1 1 1] 0 7 0.0
|
|---|
| 173 | MatrixSetValues: 1 (1- 1- 1-) (1- 3+ 4+) [1 1 1] 0 1 -1.0
|
|---|
| 174 | MatrixSetValues: 1 (1- 1- 1-) (1- 3+ 4+) [1 1 1] 0 0 501.0
|
|---|
| 175 | MatrixSetValues: 1 (4+ 1- 1-) (4+ 3+ 4+) [1 1 1] 0 8 0.0
|
|---|
| 176 | MatrixSetValues: 1 (4+ 1- 1-) (4+ 3+ 4+) [1 1 1] 0 2 -0.1
|
|---|
| 177 | MatrixSetValues: 1 (4+ 1- 1-) (4+ 3+ 4+) [1 1 1] 0 0 500.1
|
|---|
| 178 |
|
|---|
| 179 | MatrixSetValues: 2 (1- 1- 1-) (1- 3+ 4+) [1 1 1] 0 7 0.0
|
|---|
| 180 | MatrixSetValues: 2 (1- 1- 1-) (1- 3+ 4+) [1 1 1] 0 1 -100.0
|
|---|
| 181 | MatrixSetValues: 2 (1- 1- 1-) (1- 3+ 4+) [1 1 1] 0 0 100.5
|
|---|
| 182 | MatrixSetValues: 2 (4+ 1- 1-) (4+ 3+ 4+) [1 1 1] 0 8 0.0
|
|---|
| 183 | MatrixSetValues: 2 (4+ 1- 1-) (4+ 3+ 4+) [1 1 1] 0 2 -1.0
|
|---|
| 184 | MatrixSetValues: 2 (4+ 1- 1-) (4+ 3+ 4+) [1 1 1] 0 0 1.5
|
|---|
| 185 |
|
|---|
| 186 | MatrixSetValues: 3 (1- 1- 1-) (1- 3+ 4+) [1 1 1] 0 7 0.0
|
|---|
| 187 | MatrixSetValues: 3 (1- 1- 1-) (1- 3+ 4+) [1 1 1] 0 1 -0.1
|
|---|
| 188 | MatrixSetValues: 3 (1- 1- 1-) (1- 3+ 4+) [1 1 1] 0 1 5.1
|
|---|
| 189 | MatrixSetValues: 3 (1- 1- 1-) (1- 3+ 4+) [1 1 1] 0 7 0.0
|
|---|
| 190 | MatrixSetValues: 3 (1- 1- 1-) (1- 3+ 4+) [1 1 1] 0 1 -1.0
|
|---|
| 191 |
|
|---|
| 192 | MatrixSetValues: 4 (1- 1- 1-) (1- 4+ 4+) [1 1 1] 0 1 0.0
|
|---|
| 193 | MatrixSetValues: 4 (1- 1- 1-) (1- 4+ 4+) [1 1 1] 0 7 -1.0
|
|---|
| 194 |
|
|---|
| 195 | # May need to get the right discretization between parts 9 and 10
|
|---|
| 196 | # MatrixSetValues: 9 (1- 3+ 1-) (2+ 3+ 4+) [1 1 1] 0 0 4.0
|
|---|
| 197 | # MatrixSetValues: 9 (1- 3+ 1-) (2+ 3+ 4+) [1 1 1] 0 4 -1.0
|
|---|
| 198 |
|
|---|
| 199 | MatrixSetValues: 7 (3- 1- 1-) (4+ 1- 4+) [1 1 1] 0 0 11.4
|
|---|
| 200 | MatrixSetValues: 7 (3- 1- 1-) (4+ 1- 4+) [1 1 1] 0 3 -1.0
|
|---|
| 201 |
|
|---|
| 202 | # MatrixAddToValues: \
|
|---|
| 203 | # part ilower(ndim) iupper(ndim) var nentries entries[nentries] values[nentries]
|
|---|
| 204 | #MatrixAddToValues: 0 (1- 1- 1-) (4+ 4+ 4+) 0 4 [0 1 2 3] [6.0 -1.0 -1.0 -1.0]
|
|---|
| 205 |
|
|---|
| 206 | ###########################################################
|
|---|
| 207 |
|
|---|
| 208 | # VectorCreate:
|
|---|
| 209 |
|
|---|
| 210 | ###########################################################
|
|---|
| 211 |
|
|---|
| 212 | #
|
|---|
| 213 | # part: 0 1 2 3 4 5 6 7
|
|---|
| 214 | # nvar: 12 12 12 12 12 12 12 12
|
|---|
| 215 | # total: 96
|
|---|
| 216 | #
|
|---|
| 217 | # ProcessPoolCreate: num_pools
|
|---|
| 218 | ProcessPoolCreate: 1
|
|---|
| 219 |
|
|---|
| 220 | # ProcessPoolSetPart: pool part
|
|---|
| 221 | ProcessPoolSetPart: 0 0
|
|---|
| 222 | ProcessPoolSetPart: 0 1
|
|---|
| 223 | ProcessPoolSetPart: 0 2
|
|---|
| 224 | ProcessPoolSetPart: 0 3
|
|---|
| 225 | ProcessPoolSetPart: 0 4
|
|---|
| 226 | ProcessPoolSetPart: 0 5
|
|---|
| 227 | ProcessPoolSetPart: 0 6
|
|---|
| 228 | ProcessPoolSetPart: 0 7
|
|---|
| 229 |
|
|---|
| 230 | # ProcessPoolCreate: num_pools
|
|---|
| 231 | ProcessPoolCreate: 8
|
|---|
| 232 |
|
|---|
| 233 | # ProcessPoolSetPart: pool part
|
|---|
| 234 | ProcessPoolSetPart: 0 0
|
|---|
| 235 | ProcessPoolSetPart: 1 1
|
|---|
| 236 | ProcessPoolSetPart: 2 2
|
|---|
| 237 | ProcessPoolSetPart: 3 3
|
|---|
| 238 | ProcessPoolSetPart: 4 4
|
|---|
| 239 | ProcessPoolSetPart: 5 5
|
|---|
| 240 | ProcessPoolSetPart: 6 6
|
|---|
| 241 | ProcessPoolSetPart: 7 7
|
|---|
| 242 |
|
|---|
| 243 | ###########################################################
|
|---|
| 244 |
|
|---|