1
# 1 block, 21 elems, 4 nodes-per-elem, 1 field-per-node (size 1), 0 BCs 7 Lagrange CRs
23
# numFieldsPerNode is a list of length 'numNodesPerElement'
27
# nodalFieldIDs is a table with 'numNodesPerElement' rows, and the
28
# i-th row is of length 'numFieldsPerNode[i]'
35
# The mesh looks like this:
36
# (element numbers are written 'E<n>', node numbers are written
37
# below and left of each node.)
39
# *-----*-----*-----* *-------*-------*
40
# 33| 34| 35| 36| 18| 19| 20|
42
# | E14 | E15 | E16 | | | |
43
# *-----*-----*-----* | E19 | E20 |
44
# 29| 30| 31| 32| | | |
45
# | | | | *-------*-------*
46
# | E11 | E12 | E13 | 15| 16| 17|
47
# *-----*-----*-----* | | |
48
# 25| 26| 27| 28| | | |
49
# | | | | | E17 | E18 |
50
# | E8 | E9 | E10 | | | |
51
# *-----*-----*-----* | | |
54
# *--------*------------*-------*-------*
56
# | E4 | E5 | E6 | E7 |
58
# *--------*------------*-------*-------*
60
# | E0 | E1 | E2 | E3 |
62
# *--------*------------*-------*-------*
65
# Lagrange constraints will be used to tie the following nodes:
75
# Now the element-stiffnesses and loads
76
# element connectivities (list-of-connected-nodes) are numbered counter-clockwise
77
# from the lower-left of each element.
81
# connectivity -- 4 node IDs
83
# stiffness, 4 rows, 4 entries each
93
# connectivity -- 4 node IDs
95
# stiffness, 4 rows, 4 entries each
105
# connectivity -- 4 node IDs
107
# stiffness, 4 rows, 4 entries each
117
# connectivity -- 4 node IDs
119
# stiffness, 4 rows, 4 entries each
129
# connectivity -- 4 node IDs
131
# stiffness, 4 rows, 4 entries each
141
# connectivity -- 4 node IDs
143
# stiffness, 4 rows, 4 entries each
153
# connectivity -- 4 node IDs
155
# stiffness, 4 rows, 4 entries each
165
# connectivity -- 4 node IDs
167
# stiffness, 4 rows, 4 entries each
177
# connectivity -- 4 node IDs
179
# stiffness, 4 rows, 4 entries each
189
# connectivity -- 4 node IDs
191
# stiffness, 4 rows, 4 entries each
201
# connectivity -- 4 node IDs
203
# stiffness, 4 rows, 4 entries each
213
# connectivity -- 4 node IDs
215
# stiffness, 4 rows, 4 entries each
225
# connectivity -- 4 node IDs
227
# stiffness, 4 rows, 4 entries each
237
# connectivity -- 4 node IDs
239
# stiffness, 4 rows, 4 entries each
249
# connectivity -- 4 node IDs
251
# stiffness, 4 rows, 4 entries each
261
# connectivity -- 4 node IDs
263
# stiffness, 4 rows, 4 entries each
273
# connectivity -- 4 node IDs
275
# stiffness, 4 rows, 4 entries each
285
# connectivity -- 4 node IDs
287
# stiffness, 4 rows, 4 entries each
297
# connectivity -- 4 node IDs
299
# stiffness, 4 rows, 4 entries each
309
# connectivity -- 4 node IDs
311
# stiffness, 4 rows, 4 entries each
321
# connectivity -- 4 node IDs
323
# stiffness, 4 rows, 4 entries each
339
# node IDs -- 1 row for each CR, each row has 1 entry for each node
341
# fieldID for each node
343
#weights, 1 row per node, each row has 1 entry per field-component
346
#rhs values, 1 value per CR
353
# node IDs -- 1 row for each CR, each row has 1 entry for each node
355
# fieldID for each node
357
#weights, 1 row per node, each row has 1 entry per field-component
361
#rhs values, 1 value per CR
368
# node IDs -- 1 row for each CR, each row has 1 entry for each node
370
# fieldID for each node
372
# weights, 1 row per node, each row has 1 entry per field-component
376
#rhs values, 1 value per CR
383
# node IDs -- 1 row for each CR, each row has 1 entry for each node
385
# fieldID for each node
387
#weights, 1 row per node, each row has 1 entry per field-component
390
#rhs values, 1 value per CR
397
# node IDs -- 1 row for each CR, each row has 1 entry for each node
399
# fieldID for each node
401
#weights, 1 row per node, each row has 1 entry per field-component
405
#rhs values, 1 value per CR
412
# node IDs -- 1 row for each CR, each row has 1 entry for each node
414
# fieldID for each node
416
#weights, 1 row per node, each row has 1 entry per field-component
420
#rhs values, 1 value per CR
427
# node IDs -- 1 row for each CR, each row has 1 entry for each node
429
# fieldID for each node
431
#weights, 1 row per node, each row has 1 entry per field-component
434
#rhs values, 1 value per CR