Skip to content
GitHub Actions / Firedrake complex failed Oct 21, 2024 in 0s

8079 tests run, 6496 passed, 1582 skipped, 1 failed.

Annotations

Check failure on line 149 in tests/multigrid/test_transfer_manager.py

See this annotation in the file changed.

@github-actions github-actions / Firedrake complex

test_transfer_manager.test_cached_transfer[CG]

ValueError: Order 0 invalid for 'Lagrange' finite element.
Raw output
family = 'CG'

    @pytest.mark.parametrize("family", ("CG", "R"))
    def test_cached_transfer(family):
        # Test that we can properly reuse transfers within solve
        sp = {"mat_type": "matfree",
              "pc_type": "mg",
              "mg_coarse_pc_type": "none",
              "mg_levels_pc_type": "none"}
    
        base = UnitSquareMesh(1, 1)
        hierarchy = MeshHierarchy(base, 3)
        mesh = hierarchy[-1]
    
>       V = FunctionSpace(mesh, family, 0)

tests/multigrid/test_transfer_manager.py:149: 
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ 
petsc4py/PETSc/Log.pyx:188: in petsc4py.PETSc.Log.EventDecorator.decorator.wrapped_func
    ???
petsc4py/PETSc/Log.pyx:189: in petsc4py.PETSc.Log.EventDecorator.decorator.wrapped_func
    ???
firedrake/functionspace.py:104: in FunctionSpace
    element = make_scalar_element(mesh, family, degree, vfamily, vdegree, variant)
petsc4py/PETSc/Log.pyx:188: in petsc4py.PETSc.Log.EventDecorator.decorator.wrapped_func
    ???
petsc4py/PETSc/Log.pyx:189: in petsc4py.PETSc.Log.EventDecorator.decorator.wrapped_func
    ???
firedrake/functionspace.py:70: in make_scalar_element
    return finat.ufl.FiniteElement(family, cell=cell, degree=degree, variant=variant)
../firedrake_venv/src/FInAT/finat/ufl/finiteelement.py:145: in __init__
    ) = canonical_element_description(family, cell, degree, form_degree)
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ 

family = 'Lagrange', cell = triangle, order = 0, form_degree = None

    def canonical_element_description(family, cell, order, form_degree):
        """Given basic element information, return corresponding element information on canonical form.
    
        Input: family, cell, (polynomial) order, form_degree
        Output: family (canonical), short_name (for printing), order, value shape,
        reference value shape, sobolev_space.
    
        This is used by the FiniteElement constructor to ved input
        data against the element list and aliases defined in ufl.
        """
        # Get domain dimensions
        if cell is not None:
            tdim = cell.topological_dimension()
            gdim = cell.geometric_dimension()
            if isinstance(cell, Cell):
                cellname = cell.cellname()
            else:
                cellname = None
        else:
            tdim = None
            gdim = None
            cellname = None
    
        # Catch general FEEC notation "P" and "S"
        if form_degree is not None and family in ("P", "S"):
            family, order = feec_element(family, tdim, order, form_degree)
    
        if form_degree is not None and family in ("P L2", "S L2"):
            family, order = feec_element_l2(family, tdim, order, form_degree)
    
        # Check whether this family is an alias for something else
        while family in aliases:
            if tdim is None:
                raise ValueError("Need dimension to handle element aliases.")
            (family, order) = aliases[family](family, tdim, order, form_degree)
    
        # Check that the element family exists
        if family not in ufl_elements:
            raise ValueError(f"Unknown finite element '{family}'.")
    
        # Check that element data is valid (and also get common family
        # name)
        (family, short_name, value_rank, sobolev_space, mapping, krange, cellnames) = ufl_elements[family]
    
        # Accept CG/DG on all kind of cells, but use Q/DQ on "product" cells
        if cellname in set(cubes) - set(simplices) or isinstance(cell, TensorProductCell):
            if family == "Lagrange":
                family = "Q"
            elif family == "Discontinuous Lagrange":
                if order >= 1:
                    warnings.warn("Discontinuous Lagrange element requested on %s, creating DQ element." % cell.cellname())
                family = "DQ"
            elif family == "Discontinuous Lagrange L2":
                if order >= 1:
                    warnings.warn(f"Discontinuous Lagrange L2 element requested on {cell.cellname()}, "
                                  "creating DQ L2 element.")
                family = "DQ L2"
    
        # Validate cellname if a valid cell is specified
        if not (cellname is None or cellname in cellnames):
            raise ValueError(f"Cellname '{cellname}' invalid for '{family}' finite element.")
    
        # Validate order if specified
        if order is not None:
            if krange is None:
                raise ValueError(f"Order {order} invalid for '{family}' finite element, should be None.")
            kmin, kmax = krange
            if not (kmin is None or (asarray(order) >= kmin).all()):
>               raise ValueError(f"Order {order} invalid for '{family}' finite element.")
E               ValueError: Order 0 invalid for 'Lagrange' finite element.

../firedrake_venv/src/FInAT/finat/ufl/elementlist.py:478: ValueError