Unisolvence of Symmetric Node Patterns for Polynomial Spaces on the Simplex

Finite elements with polynomial basis functions on the simplex with a symmetric distribution of nodes should have a unique polynomial representation. Unisolvence not only requires that the number of nodes equals the number of independent polynomials spanning a polynomial space of a given degree, but also that the Vandermonde matrix controlling their mapping to the Lagrange interpolating polynomials can be inverted. Here, a necessary condition for unisolvence is presented for polynomial spaces that have non-decreasing degrees when going from the edges and the various faces to the interior of the simplex. It leads to a proof of a conjecture on a necessary condition for unisolvence, requiring the node pattern to be the same as that of the regular simplex.


Introduction
Mass-lumped finite elements with a symmetric distribution of nodes on the simplex enable explicit time stepping when solving the wave equation by avoiding the inversion of a large sparse mass matrix. The elements of the lumped mass matrix are proportional to quadrature weights. To avoid zero or negative weights that make time stepping unstable, the elements are enriched with higher-degree polynomials in the interior. Elements of degree 2 [1][2][3][4], 3 [5,6], 4 [7], 5 [8], 6 [9], 7 and 8 [10,11] have been found. The last two papers also contain degree-9 elements, but the one in [11] has degree 10 instead of 9 on the edges. In 3D, elements of degree 2 [7], 3 [8], and 4 [12] were found. In the last paper, some of the elements only involve a subset of the higher-degree polynomials on the faces or in the interior of the tetrahedron. Elements of that type will not be considered in this paper.
The construction of these elements involves several steps. Nodes in a symmetric node pattern can be classified according to their equivalence class. For a given degree on the edges, a node pattern has to be chosen, as well as degrees for the faces and interior. All elements found so far include the vertices. Then, a set of equations has to be formed that require quadrature to be exact for a certain set of polynomials. The resulting system may be inconsistent, zero-dimensional with a finite number of solutions, or infinitely many solutions may exist. Among these solutions, only the real-valued ones with positive weights and nodes that do not lie outside the element, or move outside their symmetry class or coincide with another node are of interest, if they exist and can be found. Finally, the basis polynomials have to be mapped to Lagrange interpolating polynomials. It may happen that the Vandermonde matrix of the linear system that governs the mapping may be rank deficient and the element has to be discarded because it is not unisolvent.
To avoid that problem, it helps to know in advance which node patterns for given degrees on edges and interior, and faces in 3D, can provide unisolvence. As degenerate cases such as coinciding nodes may still occur, a necessary condition is all one can ask for. Conjecture 3.1 in [11] states that a node pattern or orbits of a symmetric set of nodes can be unisolvent only if it is the same as the pattern of the simplex with a regular symmetric distribution of nodes, placed on a uniform grid or lattice. A proof of this conjecture as a necessary condition for unisolvence on triangles was presented in [13], for arbitrary polynomial degrees on edges and in the interior, the latter not smaller then the degree on the edges. Triangles and tetrahedra are considered up to degree 10 and 7, respectively, in [14], which only came to my attention when the current paper was almost finished-all part of the game.
The basic principle is that the Lagrange interpolating polynomial of a node of an equivalence class or orbit with a higher degree of symmetry should be invariant under the same coordinate permutations that leave that node invariant. A Lagrange interpolating polynomial can be said to belong to the node where it is one, while zero on all other nodes. This principle was applied in [13] for 2D and allows a straightforward generalization to higher dimensions-the subject of the current paper. Section 2 describes the equivalence classes or orbits that define the symmetric node patterns on the simplex. Section 3 defines the polynomial spaces in terms of monomials, possibly with different degrees on the m-faces. It is shown that the invertibility of the Vandermonde matrix that defines the Lagrange interpolating polynomials can be decomposed into subproblems that only involve the interior nodes for simplices in m dimensions, with m up to the d space dimensions of the d-simplex.
A necessary condition for unisolvence is that the number of unknown coefficients in the Lagrange interpolating polynomials should equal the number of equations obtained by requiring the interpolants to be one on one node and zero on the others. Section 4 contains an example that illustrates how the number of independent coefficients for a given polynomial degree depends on the symmetry class of the node where the interpolant is one and how many equations are available for each class. This leads to a set of linear constraints on the node pattern.
Section 5 provides a formal definition of the linear system, with a matrix depending on the dimension and acting on a vector that defines the node pattern, the latter describing the number of generating nodes per equivalence class or orbit. The right-hand side depends on the maximum degree of the polynomial space. It is subsequently proven that this system has a unique solution. Since the regular symmetric distribution of nodes is known to be unisolvent [15], this proves Conjecture 3.1 of [11]. Section 6 summarizes the main result.
x k and 0 ≤ x k ≤ 1. Symmetric node sets are represented by the symmetry group S d+1 containing the (d + 1)! permutations of the coordinates of a generating or reference node x g . They can be divided onto equivalence classes by selecting subsets for which two or more of the x k are equal.
Keast [16] denotes the equivalence classes by n = [n 1 , n 2 , . . . , n k ], with n 1 ≥ n 2 ≥ . . . ≥ n k > 0, where n k is the number of repetitions of a node parameter a k for a generating interior node with barycentric coordinate x g (n) = (a 1 , a 1 , . . . , a 1 , a 2 , . . . , a 2 , . . . , a k , . . . , a k ), where all a , = 1, . . . , k, are distinct and positive and k =1 n = d + 1. The number of nodes per class characterised by n for an interior node is The equivalence classes on the m-face can be obtained from the lower-dimensional msimplex by zero padding of the barycentric coordinates of the interior nodes to a length d + 1 and by replacing The other facets then follow from the coordinate permutations. For the m-faces, the number of generated nodes for a class n is By using m +1 = k =1 n , the right-hand side does not explicitly refer to the face dimension m.
In what follows, it is assumed that the equivalence classes n j for the interior nodes are sorted in descending order, starting in the first column and moving to the next in case of equality. Below, the equivalence class or orbit defined by n j will sometimes be loosely referred to as class j.
Plane through 2 vertices and centroid The generating nodes are expressed in barycentric coordinates. The number of generated nodes is ν j = ν(3, n j ) Figure 1 shows an example of the interior nodes for d = 2 and Table 1 lists the equivalence classes and generating nodes for d = 3. For each class, the symmetric set of nodes consists of the ν(d, n j ) permutations of the generating or reference node x g (n j ).
A symmetric node pattern defines the number of generating nodes for each equivalence class. For the interior, the pattern is denoted by a vector K d of length n c (d), where K d j , j = 1, . . . , n c (d), is the number of generating nodes x g (n j ) in class n j (d). Note that the notation K d j rather than k d j is used to avoid confusion with an index k. Concatenation of the interior node patterns K m for the m-faces, with m increasing from 0 to d, provides the overall node pattern (K 0 , K 1 , . . . , K d ) for the d-simplex. For the mass-lumped finite elements mentioned in the introduction, the vertices were always included, implying that K 0 = (1).

Polynomials
Polynomials basis functions are chosen to have a degree p d in the interior and degree p m when restricted to the m-face. The mass-lumped finite elements mentioned in the introduction generally have p k ≥ p k−1 for k = 1, . . . , d. On a vertex, p 0 = 1. The overall degree of the element is taken as p 1 on the edges, with the vertices at the endpoints included.
Define an operator S[ f ] that, when acting on a function f (x 0 , x 1 , . . . , x d ), produces a set of functions obtained by permuting its arguments. The same notation will be used if f is not a single function but a set of functions. In that case, the operator is applied to each element of the set and the results are combined into a single set with duplicates removed. Alternatively, let S k be a permutation matrix, with k = 1, . . . , (d + 1)!, corresponding to an element S k of the set of permutations operators S, in the sense that S k [ f (x)] = f (S k x) when the arguments are expressed as a vector x. Another approach is a set of index arrays s k that lets S k, x = x s k, for each row of the matrix S k . Then, S k, ,s k, = 1 and the other elements are zero.
The bubble function on one of the m-faces, the one with x k = 0 for m < k ≤ d, is η m = Consider the polynomials containing products of the bubble function with the monomials in the set Note that (x 1 , . . . , x m ) defines a point on the m-face, but that The polynomial space is given by resulting in a total of n total The element characterized by a set of nodes and polynomials is unisolvent if A is not singular, i.e., det(A) = 0 or rank(A) = n total x (d). Proof If the nodes are sorted according the m-faces that contain them, with increasing m, and the polynomials as well, the Vandermonde matrix is block diagonal:

Corollary 3.1 In Theorem 1, the polynomial space V m on the m-simplex can be replaced by
Proof On the right-hand side of equation (7), the polynomials φ m j for the m-face have a bubble functions η m as factor, which is shared by the polynomials ψ m j on the m -faces with m > m. Because of the block upper triangular structure, the bubble function η m is inherited by the ψ m i on the m-face, that is, ψ m i = η m ψ m i , consistent with Lemma 3.1.10 in [17]. Factoring out η m provides the result.
Given a generating node x j = x g (n j ) for equivalence class n j , its symmetric counterparts are obtained by S k x j for all permutation matrices S k . Depending on the amount of symmetry, a subset of these permutation matrices may already suffice.
The Lagrange interpolating polynomial that is one at x j and zero at all other nodes on the simplex can be expressed as It can be loosely referred to as the Lagrange polynomial for node x j . After applying a permutation matrix S k , ψ j (S k x) will be one for S k x = S k x j and zero at the other nodes, defined by If x j is one of the generating nodes for equivalence class j and x j is one of its symmetric counterparts, then equation (9b) provides the Lagrange interpolating polynomial ψ j (x) for that node if ψ j (x) is known. This relation is convenient when coding up a finite element, because it reduces the number of explicit functions and their gradients, needed in for instance [18], to functions for the generating nodes only, which then can be called for the other nodes by rearranging the function arguments.

Counting Coefficients
Assuming that the Vandermonde matrix A d,d for the interior nodes is invertible, the Lagrange polynomial for x j with ψ j (x k ) = δ j,k can be expressed as where the coefficients are taken from row j of B = A d,d −1 . If a permutation S k is applied to a node that stays invariant, equation (9b) with j = j states that ψ j (x) should also be invariant under S k . The monomials x , however, may or may not stay invariant. Equating ψ j (x) = ψ j (S k x) for all permutations that leave x j invariant provides additional equations for the b j,m and decreases the number of independent coefficients.
The number of coefficients per class can be summarized by a vector r = (3, 6, 10) T . The terms of the Lagrange interpolants that contain only one coefficient per monomial are c 1,1 x 0 x 1 x 2 for class 1, c 2,1 x 0 x 1 x 2 and c 2,6 x 3 2 for class 2, and all terms for class 3. These numbers can be summarized by a vector f = (1, 2, 10) T .
The coefficients for each Lagrange polynomial should follow from the requirement that it is one on one node and zero on the others. For class 2, substitution of a class-2 node of the form x = (a, a, b) and its permutations into ψ 2 produces 3 equations, ψ 2 (a, a, b) = 1, ψ 2 (a, b, a) = 0 and ψ 2 (b, a, a) = 0, of which the last 2 are the same. Likewise a class-3 node of the form x = (a, b, c) provides 3 equations for 6 nodes. Repeating this for ψ 1 produces 1 equation for a node of class 1, 2 or 3. For ψ 5 , the number of equations are 1, 3, and 6, respectively. Substitution of the centroid provides one equation for each ψ k .
Requiring the number of equations to be the same as the number of coefficients as a necessary condition for unisolvence provides WK = r, with where r counts the number of coefficients for a given degree and W counts the number of available equations for the nodes of each class. The current example has K = (1, 1, 1) T as solution, one generating node for each class.
A subset is formed by the coefficients of the Lagrange interpolants that involve a single monomial. These are counted by f. The corresponding subset of W is denoted by V and leads to the set of conditions VK = f. The above term c 1,1 x 0 x 1 x 2 has class-1 symmetry and results in v 1,1 = 1, since there is only one such term. The term c 2,1 x 0 x 1 x 2 also has class-1 symmetry and is counted by v 2,1 , and c 2,6 x 3 2 has class-2 symmetry is accounted for by v 2,2 . The third row of V count how the monomials for ψ 5 are distributed over the equivalence classes and is the same as the last row of W. All together, A necessary condition for unisolvence is that the total number of nodes equals the number of independent polynomials spanning a polynomial space, in this case the monomials of equation (4) with m = d. The example reveals additional requirements for each equivalence class, as a number of linear constraints on K of the form or WK = r( p).
Formal definitions will be given in the next section. Simple cases are d = 0 and 1. For d = 0, K 1 = 1 if the vertices are included and K 1 = 0 if not.
In the 1-D case, there is a unique polynomial of degree p that interpolates p + 1 distinct points. The number of interior nodes n x = K 1 + 2K 2 , where subscript 1 corresponds to the midpoint of an edge and subscript 2 to a different point, excluding the endpoints. The factor 2 is due to symmetry. The number of nodes should equal the number of coefficients in the polynomial: (11) and (12) become

The expressions in equations
and

Necessary Conditions for Unisolvence
According to Corollary 3.1, it suffices to consider unisolvence for polynomials of degree p = p m − m − 1 and the node pattern K m for the interior nodes of the simplex in m dimensions, for m = 0, . . . , d. From here onwards, a fixed spatial dimension m = d will be considered and the related sub-or superscript or argument will often be dropped. One necessary condition for unisolvence is that the number of nodes n x (d) in the interior of the d-simplex equals the number of monomials n p (d, p) in the set defined by equation (4) for m = d. This provides one linear constraint, detailed as Property A.1 in Appendix A. Other necessary conditions in the form of constraints are the related to the symmetry of the equivalence classes, one condition for each class, leading to matrices V and W of size n c ×n c and vectors f and r of length n c . The conditions (11) are obtained by counting terms in the Lagrange interpolant with a single coefficient for a monomial. The conditions (12) involve a count of these coefficients as well as the number of terms with more than one monomial per coefficient. Before going to formal definitions of the matrices and right-hand sides, two more examples may help to understand what they represent.

Example 2
Consider the Lagrange interpolant for the centroid. According to equation (9), it should be symmetric with respect to the centroid and, therefore, can only depend on polynomials that are symmetric in x 0 to x d . The number of coefficients required to represent these polynomials in terms of the monomials of equation (4) with m = d is counted by r 1 ( p), for class 1 with n 1 = [d + 1] that describes the symmetry of the centroid.
The first row of W describes how these polynomials are distributed over the nodes of the various symmetry classes. Symmetric polynomials for each class can be obtained by either taking monomials with the same symmetry or by adding the least number of monomials that provide the desired symmetry. For class 1, this results in only one coefficient per class, and hence, the first row of W has w 1,k = 1 for k = 1, . . . , n c . Property B.3 lists additional details.
A subset of these polynomials share the symmetry of class 1 and do not have to be added. Their number is counted by f 1 ( p) and the way they are distributed over the classes by the first row of V, which then will have 1 on the diagonal and zeros elsewhere, spelled out as Property A.2.

Example 3
Another example is the hyperplane defined by x 0 = x 1 = · · · = x d−1 . A node on that plane either can be taken as a generating node for class n c − 1, or it belongs to a class j < n c − 1 with a higher degree of symmetry. The centroid is among them. The right-hand side f n c −1 ( p) counts how many monomials have the same invariance as the nodes on the chosen hyperplane and row n c − 1 of V describes how these are distributed over the classes these nodes belong to.
Nodes outside the plane have mirror symmetry w.r.t. the chosen hyperplane. Therefore, only half of the number of nodes outside the hyperplane can contribute to define the coefficients of the Lagrange interpolants that are one on a node of the hyperplane and have the same mirror symmetry, as stated by equation (9). This determines row number n c − 1 of W, whereas r n c −1 ( p) counts the number of contributing monomials. Details are given in Properties A.5, B.4 and B.5 and their proofs.
For unisolvence, the necessary conditions on the node pattern K are summarized by VK = f( p), where V is a square matrix with n c (d) rows and columns and f( p) is a vector that depends on the degree p = p d − d − 1. The matrix is defined by the following.
Recall that K j is the number of generating nodes x g (n j ) in equivalence class n j for the interior of the d-simplex. The other nodes are obtained as the permutations S i x g (n j ) of the generating node.
A subset L j of the set S of all the permutations on the simplex is defined by those that leave the nodes of equivalence class n j invariant. The entry v j,k is defined as the number of nodes of the same or another equivalence class n k that are also invariant under the permutations of L j . Stated in another way: if is the set of nodes obtained from node x by applying the permutations of L j , then the set contains the nodes of class n k that remain invariant under L j and v j,k = X j,k is defined as its cardinality. (a, a, b, b) Table 1 for 3D belongs to class 3 with n 3 = [2, 2]. Among its 6 permutations, (a, a, b, b), (a, b, a, b), (a, b, b, a), (b, a, a, b),  (b, a, b, a) and (b, b, a, a), there are two that remain invariant under L 3 , namely the node itself and (b, b, a, a), implying v 3,3 = 2. From the other classes, only the centroid is invariant under L 3 . Therefore, row 3 of V is (1, 0, 2, 0, 0).

Example 4 The node
The components f j ( p) on the right-hand side of equation (11) are defined as the number of monomials in the set P p (x 0 , . . . , x d ) of equation (4) that are invariant under the coordinate permutation L j . Before going to a formal definition, the following is worth noting.

Theorem 2 The powers of the monomials in the set P p of equation (4) in d dimensions share their equivalence classes with a regular, unisolvent set of nodes on the d-simplex.
Proof The distribution of powers in equation (4) for m = d follows that of the regular dsimplex with sides of length p in natural coordinates, after multiplying the length of its edges by p. In this way, k / p is identified with the x k -axis for k = 0, . . . , d. This regular simplex with nodes on a uniform grid is unisolvent [15] and has the same number of nodes as the number of monomials, specified by equation (A1). An index permutation s j,k corresponding to S k ∈ L j for class j, applied to both the node coordinates and the powers, leaves the monomial the same: Because interchanging coordinates simultaneously with their powers does not alter the monomial, the generating nodes for the coordinates and for the powers belong to the same equivalence class. If some of the coordinates of a node in class j are the same, the corresponding powers appear as a sum in the monomials, which is invariant under the permutations from the set L j . If some of the powers are the same in agreement with class j, the product of the coordinates is invariant under the permutations from the set L j .
Note that this match between the equivalence classes of the powers of monomials and those of the interpolating nodes has a wider applicability to other symmetric nodes sets in combination with symmetric sets of polynomials or other types of symmetric basis functions.
With Theorem 2, f j (d) can be defined analogously to equations (15) and (16). Let characterize the monomials by their powers. The set of sets characterizes monomials obtained from the powers ∈ O d p by applying the permutations of L j and The matrix W and the related right-hand side vector r( p) are not needed in what follows. Appendix B contains their definitions as well as various properties.

Theorem 3 The matrix V can be inverted.
Proof Entry v j,k of V represents the number of nodes of equivalence class n k that are invariant under the permutations L j that leave nodes of class n j invariant.
A node of class n j is by definition invariant under L j , implying that v j, j ≥ 1. Classes with a higher degree of symmetry have k < j if the equivalence classes are ordered according to Sect. 2, implying v j,k ≥ 0, whereas nodes with a lower degree of symmetry are excluded, implying v j,k = 0 for k > j. Therefore, the matrix is lower triangular with non-negative integers on the diagonal. Because its determinant is the product of the diagonal values, the linear system VK = f has a unique solution.

Theorem 4 A necessary condition for unisolvence of a node pattern for degree p on the d-simplex is that it is the same as the node pattern of the regular degree-p element on the d-simplex.
Proof The regular simplex is unisolvent [15] and satisfies equation (11). Because this system has a unique solution according to Theorem 3, this proves Conjecture 3.1 of [11], stating that a necessary condition for unisolvence in the interior of the d-simplex is a node pattern or orbit that is the same as that of a regular symmetric distribution of nodes.
From the monomial powers for degree p in d dimensions, not only f can be determined, but also K. The Mathematica [19] code fragment in Listing 1 computes these node patterns as Kreg on a regular element from degree 0 to pmax on the simplex in ndim dimensions, without referring to f. More properties of V and f can be found in Appendix A. Appendix C lists examples of V and f for dimensions 1 to 6 up to degree 23. The Mathematica kernel script sympat.m that generated these results with the option to include W and r, for limited values of the degree p and a small number of dimensions d, is included as supplemental file.

Conclusions
For a polynomial space on the d-simplex with non-decreasing degrees on its m-faces as m increases, a necessary condition on the symmetric node pattern for unisolvence has been presented. It leads to a proof of the conjecture that the node pattern should be the same as that of the regular d-simplex for a given degree, with nodes on a uniform grid. The line of reasoning, in particular of Theorem 2 stating that the symmetry of the equivalence classes or orbits for the nodes should match that of the monomial powers in the polynomial basis, may be generalized to other symmetric distributions of nodes and polynomial spaces or other types of basis functions with analogous symmetry properties.

Conflict of interest
The author declares that he have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Appendix A Properties of V and f
Some additional properties of the system (11) are presented.
Identifying this with the last row, numbered n c (d), of the linear system (11) provides v n c (d),k = ν(d, n k ) for k = 1, . . . , n c and f n c (d) ( p) = p+d d . Note that the only permutation that leaves a node of class n n c invariant is the identity, and therefore, all the nodes are involved.  max(2, d)).
Proof The coordinates of the centroid are invariant under all permutations and there is only one centroid, so v j,1 = 1 for j = 1, . . . , n c (d). For degree p = 0, there is only one degree of freedom. The only node that has ν(d, n j ) = 1 in equation (1) is the centroid in class j = 1 with n 1 = [d + 1], implying K j = δ 1, j and r j (0) = 1 for all j = 1, . . . , n c (d).
Property A. 5 The number of monomials f n c (d)−1 ( p) obeys the recurrence relation Proof Class n c (d) − 1 with n n c (d)−1 = [2, 1, . . . , 1] has a generating node with x 0 = x 1 . Monomials with this symmetry are elements of

Appendix B Properties of W and r
In [13], a different criterion was used to prove the 2-D case on triangles. Although not explicitly stated in that paper, the linear system is that of equation (12) instead of equation (11).
The matrix W has n c (d) rows and columns and the vector r( p) depends on the degree p = p d −d −1. The matrix is defined as follows. Recall the definition of N j (x) = S i x S i ∈ L j in equation (15). Then, class n k , i = 1, . . . , ν(d, n k ) .
The right-hand side r is given by with O d p from equation (18) and M j ( ) = S i T S i ∈ L j from (19). The main difference with the system (11) is that not only those nodes are counted that are invariant under the permutations that leave the generating node invariant, but also polynomials with the same invariance, constructed by adding the least number of monomials to obtain that invariance.
Some properties of W are listed next. Proof For degree zero ( p = 0), there is only one degree of freedom. The only node that has ν(d, n j ) = 1 in equation (1) is the centroid of class j = 1 with n 1 = [d + 1], implying K j = δ 1, j and r j (d, 0) = 1 for all j = 1, . . . , n c . A degree one element with p = 1 requires d + 1 nodes. This only can be obtained for class 2 characterized by n 2 = [d − 1, 1] and having ν(d, n 2 ) = d + 1. Then, K j = δ 2, j leading to r j (d, 1) = w j,2 .

Property B.2 The last row of W is the same as that of V and last row of r obeys r n
Proof Because the only permutation for n n c that leave a node invariant is the identity, N n c (x i ) = x i for any node i and, therefore, X j,k in equation (16) and Z j,k (B5) coincide for j = n c , proving v n c ,k = w n c ,k for k = 1, . . . , n c . The same argument provides r n c (d) (d, p) = f n c (d) (d, p).

Property B.3
The first row of r has entries r 1 (d, p) = p ≤ p (d + 1) for p = 0, 1, . . ., where p ≤k (n) is the number of partitions of n into at most k parts.
Proof The Lagrange interpolating polynomial related to the centroid in class n 1 should be invariant under all permutations of the nodes and, therefore, only depend on polynomials that are symmetric in x 0 to x d . Let these be given by λ k , with λ 0 = 1, Then, the Lagrange interpolant can be expressed as where m k ≥ 0 for k = 1, . . . d and with starting values from property B.1 and equation (14). The number of partitions of n into at most k parts has a recursion with special cases p ≤2 (n) = ceil((n + 1)/2) = 1 + floor(n/2) and p ≤3 (n) = round (n + 3) 2 /12 [20, e.g.]. Clearly, Because each reference node contributes only once and its symmetric counterparts are included by the use of symmetric polynomials, the number of coefficients should obey Lagrange interpolant for nodes on the hyperplane should be symmetric inx 1 and only contain monomials with even powers ofx 1 . Therefore, effectively only half the nodes of the simplex outside the hyperplanex 1 = 0 contribute, either forx 1 < 0 orx 1 > 0. The total number of nodes is k w n c ,k K k = k v n c −1,k K k + 2u n c −1,k K k for all K and w n c −1,k = u n c −1,k + v n c −1,k by definition. Eliminating u n c −1,k provides the result.
Property B.5 The number of monomials r n c −1 (d, p) that can contribute to the Lagrange interpolant of a node of class n c − 1 obeys the recurrence relation The distribution of powers in equation (4) for m = d follows that of the regular d)-simplex with sides of length p in natural coordinates, after multiplying the length of its edges by p. In this way, k / p is identified with the x k -axis for k = 0, . . . , d. According to Property B.4, only even powers of 1 occur. If the nodes with odd 1 are removed, there should be r n c −1 (d, p) nodes left. The nodes that have been removed cover a simplex of similar shape, but lack one node for the longest edge in the 1 direction, leaving r n c −1 (d, p − 1) of them. With a total of r n c (d, p), this geometrical approach proves the property.
Alternatively, the relevant monomials are elements of the set Its number of elements is the same as that of P n c −1 p−1 , with r n c −1 (d, p − 1) elements, since the only difference is a factorx 1 in the monomials.
The generating function for the total number of nodes in equation (A1) is The following is based on observations, but remains to be proven:

Conjecture 1
The matrix W is symmetric and invertible.

Conjecture 2
The node patterns obtained for the system (12) are the same as those resulting from (11).
This alternative approach was useful for the 2-D case and its generalization to 3D, but did not help to treat the problem in an arbitrary number of dimensions.