@@ -19,13 +19,13 @@ class LagrangeDualSet(dual_set.DualSet):
19
19
20
20
:arg ref_el: The simplicial complex.
21
21
:arg degree: The polynomial degree.
22
- :point_variant: The point distribution variant passed on to recursivenodes.
23
- :lexsort_entities : A flag to sort entities by lexicographically by coordinate,
24
- if false then entities are sorted first by dimension and then by
25
- entity id. The DOFs are always sorted by the entity ordering
26
- and then lexicographically by coordinate .
22
+ :arg point_variant: The point distribution variant passed on to recursivenodes.
23
+ :arg sort_entities : A flag to sort entities by support vertex ids.
24
+ If false then entities are sorted first by dimension and then by
25
+ entity id. The DOFs are always sorted by the entity ordering
26
+ and then lexicographically by lattice multiindex .
27
27
"""
28
- def __init__ (self , ref_el , degree , point_variant = "equispaced" , lexsort_entities = False ):
28
+ def __init__ (self , ref_el , degree , point_variant = "equispaced" , sort_entities = False ):
29
29
nodes = []
30
30
entity_ids = {}
31
31
entity_permutations = {}
@@ -38,8 +38,8 @@ def __init__(self, ref_el, degree, point_variant="equispaced", lexsort_entities=
38
38
entity_permutations [dim ][entity ] = perms
39
39
40
40
entities = [(dim , entity ) for dim in sorted (top ) for entity in sorted (top [dim ])]
41
- if lexsort_entities :
42
- # sort the entities by their support vertex ids
41
+ if sort_entities :
42
+ # sort the entities by support vertex ids
43
43
support = [top [dim ][entity ] for dim , entity in entities ]
44
44
entities = [entity for verts , entity in sorted (zip (support , entities ))]
45
45
@@ -67,16 +67,16 @@ class Lagrange(finite_element.CiarletElement):
67
67
variant='equispaced,Iso(2)' with degree=1 gives the P2:P1 iso element.
68
68
variant='Alfeld' can be used to obtain a barycentrically refined
69
69
macroelement for Scott-Vogelius.
70
- :arg lexsort_entities : A flag to sort entities lexicographically by coordinate,
71
- if false then entities are sorted first by dimension and then by
72
- entity id. The DOFs are always sorted by the entity ordering
73
- and then lexicographically by coordinate .
70
+ :arg sort_entities : A flag to sort entities by support vertex ids.
71
+ If false then entities are sorted first by dimension and then by
72
+ entity id. The DOFs are always sorted by the entity ordering
73
+ and then lexicographically by lattice multiindex .
74
74
"""
75
- def __init__ (self , ref_el , degree , variant = "equispaced" , lexsort_entities = False ):
75
+ def __init__ (self , ref_el , degree , variant = "equispaced" , sort_entities = False ):
76
76
splitting , point_variant = parse_lagrange_variant (variant )
77
77
if splitting is not None :
78
78
ref_el = splitting (ref_el )
79
- dual = LagrangeDualSet (ref_el , degree , point_variant = point_variant , lexsort_entities = lexsort_entities )
79
+ dual = LagrangeDualSet (ref_el , degree , point_variant = point_variant , sort_entities = sort_entities )
80
80
if ref_el .shape == LINE :
81
81
# In 1D we can use the primal basis as the expansion set,
82
82
# avoiding any round-off coming from a basis transformation
0 commit comments