diff --git a/FIAT/fdm_element.py b/FIAT/fdm_element.py
index 875a9843..ab82c2dc 100644
--- a/FIAT/fdm_element.py
+++ b/FIAT/fdm_element.py
@@ -35,7 +35,6 @@ def __init__(self, ref_el, degree, bc_order=1, formdegree=0):
         # Define the generalized eigenproblem on a GLL element
         gll_degree = degree + formdegree
         gll = GaussLobattoLegendre(ref_el, gll_degree)
-        xhat = numpy.array([list(x.get_point_dict().keys())[0][0] for x in gll.dual_basis()])
         E = numpy.eye(gll.space_dimension())
 
         bdof = []
@@ -83,7 +82,7 @@ def __init__(self, ref_el, degree, bc_order=1, formdegree=0):
         if formdegree == 0:
             basis = numpy.dot(S.T, E0)
             # Eigenfunctions in the Lagrange basis
-            self._points = xhat
+            self._points = numpy.array([list(x.get_point_dict().keys())[0][0] for x in gll.dual_basis()])
             self._tabulation = numpy.dot(S.T, E.T)
         else:
             # Take the derivative of the eigenbasis and normalize