Constructors for first two columns of FE periodic table

fuse/__init__.py

@@ -1,5 +1,5 @@
 from fuse.cells import Point, Edge, polygon, make_tetrahedron, constructCellComplex
-from fuse.groups import S1, S2, S3, D4, Z3, Z4, C3, C4, S4, A4, diff_C3, tet_edges, tet_faces, sq_edges, GroupRepresentation, PermutationSetRepresentation, get_cyc_group, get_sym_group
+from fuse.groups import S1, S2, S3, D4, Z3, Z4, C3, C4, S4, A4, diff_C3, new_C3, tet_edges, tet_faces, sq_edges, GroupRepresentation, PermutationSetRepresentation, get_cyc_group, get_sym_group
 from fuse.dof import DeltaPairing, DOF, L2Pairing, FuseFunction, PointKernel, VectorKernel, BarycentricPolynomialKernel, PolynomialKernel, ComponentKernel
 from fuse.triples import ElementTriple, DOFGenerator, immerse
 from fuse.traces import TrH1, TrGrad, TrHess, TrHCurl, TrHDiv

fuse/dof.py

@@ -210,14 +210,14 @@ def _from_dict(obj_dict):
 class BarycentricPolynomialKernel(BaseKernel):
 
     def __init__(self, fn, g=None, symbols=[]):
-        if hasattr(fn, "__iter__"):
+        if hasattr(fn, "__iter__") or isinstance(fn, sp.Matrix):
             # if len(symbols) != 0 and any(not sp.sympify(fn[i]).as_poly() for i in range(len(fn))):
             #     raise ValueError("Function components must be able to be interpreted as a sympy polynomial")
             self.fn = [sp.Poly(fn[i], symbols) for i in range(len(fn))]
             self.shape = len(fn)
         else:
-            if len(symbols) != 0 and not sp.sympify(fn).as_poly():
-                raise ValueError("Function must be able to be interpreted as a sympy polynomial")
+            # if len(symbols) != 0 and not sp.sympify(fn).as_poly():
+            #     raise ValueError("Function must be able to be interpreted as a sympy polynomial")
             # self.fn = sp.sympify(fn)
             self.fn = sp.Poly(fn, symbols)
             self.shape = 0