Bug#1042003: python3-dolfinx: test_vector_function.py fails with PETSc 3.19
Drew Parsons
dparsons at debian.org
Tue Jul 25 15:27:15 BST 2023
Package: python3-dolfinx
Version: 1:0.6.0-6
Severity: normal
Control: forwarded -1 https://github.com/FEniCS/dolfinx/issues/2736
The python unit test python/test/unit/fem/test_vector_function.py
fails after building dolfinx 0.6.0 against PETSc 3.19
$ pytest-3 -v python/test/unit/fem/test_vector_function.py
================================================================================================================================== test session starts ===================================================================================================================================
platform linux -- Python 3.11.4, pytest-7.4.0, pluggy-1.2.0 -- /usr/bin/python3
cachedir: .pytest_cache
hypothesis profile 'default' -> database=DirectoryBasedExampleDatabase('/projects/fenics/build/fenics-dolfinx/.hypothesis/examples')
rootdir: /projects/fenics/build/fenics-dolfinx/python
configfile: setup.cfg
plugins: xvfb-2.0.0, remotedata-0.4.0, asyncio-0.20.3, arraydiff-0.5.0, hypothesis-6.82.0, filter-subpackage-0.1.2, flaky-3.7.0, mock-3.11.1, mpi-0.6, astropy-header-0.2.2, cov-4.1.0, openfiles-0.5.0, astropy-0.10.0, doctestplus-0.13.0
asyncio: mode=Mode.STRICT
collected 10 items
python/test/unit/fem/test_vector_function.py::test_div_conforming_triangle[1-RT] FAILED [ 10%]
python/test/unit/fem/test_vector_function.py::test_div_conforming_triangle[2-RT] FAILED [ 20%]
python/test/unit/fem/test_vector_function.py::test_div_conforming_triangle[3-RT] FAILED [ 30%]
python/test/unit/fem/test_vector_function.py::test_div_conforming_triangle[4-RT] FAILED [ 40%]
python/test/unit/fem/test_vector_function.py::test_div_conforming_triangle[5-RT] FAILED [ 50%]
python/test/unit/fem/test_vector_function.py::test_div_conforming_tetrahedron[1-RT] FAILED [ 60%]
python/test/unit/fem/test_vector_function.py::test_div_conforming_tetrahedron[2-RT] FAILED [ 70%]
python/test/unit/fem/test_vector_function.py::test_div_conforming_tetrahedron[3-RT] FAILED [ 80%]
python/test/unit/fem/test_vector_function.py::test_div_conforming_tetrahedron[4-RT] FAILED [ 90%]
python/test/unit/fem/test_vector_function.py::test_div_conforming_tetrahedron[5-RT] FAILED [100%]
======================================================================================================================================== FAILURES ========================================================================================================================================
___________________________________________________________________________________________________________________________ test_div_conforming_triangle[1-RT] ___________________________________________________________________________________________________________________________
space_type = 'RT', order = 1
@pytest.mark.skip_in_parallel
@pytest.mark.parametrize('space_type', ["RT"])
@pytest.mark.parametrize('order', [1, 2, 3, 4, 5])
def test_div_conforming_triangle(space_type, order):
"""Checks that the vectors in div conforming spaces on a triangle are correctly oriented"""
# Create simple triangle mesh
def perform_test(points, cells):
domain = ufl.Mesh(ufl.VectorElement("Lagrange", "triangle", 1))
mesh = create_mesh(MPI.COMM_WORLD, cells, points, domain)
V = FunctionSpace(mesh, (space_type, order))
f = Function(V)
output = []
for dof in range(len(f.vector[:])):
f.vector[:] = np.zeros(len(f.vector[:]))
f.vector[dof] = 1
points = np.array([[.5, .5, 0], [.5, .5, 0]])
cells = np.array([0, 1])
result = f.eval(points, cells)
normal = np.array([-1., 1.])
output.append(result.dot(normal))
return output
points = np.array([[0, 0], [1, 0], [1, 1], [0, 1]])
cells = np.array([[0, 1, 2], [2, 3, 0]])
> result = perform_test(points, cells)
python/test/unit/fem/test_vector_function.py:43:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
python/test/unit/fem/test_vector_function.py:31: in perform_test
f.vector[:] = np.zeros(len(f.vector[:]))
petsc4py/PETSc/Vec.pyx:107: in petsc4py.PETSc.Vec.__getitem__
???
petsc4py/PETSc/petscvec.pxi:436: in petsc4py.PETSc.vec_getitem
???
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
> ???
E petsc4py.PETSc.Error: error code 73
E [0] VecGetValues() at ./src/vec/vec/interface/rvector.c:924
E [0] Object is in wrong state
E [0] Not for unassembled vector, did you call VecAssemblyBegin()/VecAssemblyEnd()?
petsc4py/PETSc/petscvec.pxi:388: Error
...etc
================================================================================================================================ short test summary info =================================================================================================================================
FAILED python/test/unit/fem/test_vector_function.py::test_div_conforming_triangle[1-RT] - petsc4py.PETSc.Error: error code 73
FAILED python/test/unit/fem/test_vector_function.py::test_div_conforming_triangle[2-RT] - petsc4py.PETSc.Error: error code 73
FAILED python/test/unit/fem/test_vector_function.py::test_div_conforming_triangle[3-RT] - petsc4py.PETSc.Error: error code 73
FAILED python/test/unit/fem/test_vector_function.py::test_div_conforming_triangle[4-RT] - petsc4py.PETSc.Error: error code 73
FAILED python/test/unit/fem/test_vector_function.py::test_div_conforming_triangle[5-RT] - petsc4py.PETSc.Error: error code 73
FAILED python/test/unit/fem/test_vector_function.py::test_div_conforming_tetrahedron[1-RT] - petsc4py.PETSc.Error: error code 73
FAILED python/test/unit/fem/test_vector_function.py::test_div_conforming_tetrahedron[2-RT] - petsc4py.PETSc.Error: error code 73
FAILED python/test/unit/fem/test_vector_function.py::test_div_conforming_tetrahedron[3-RT] - petsc4py.PETSc.Error: error code 73
FAILED python/test/unit/fem/test_vector_function.py::test_div_conforming_tetrahedron[4-RT] - petsc4py.PETSc.Error: error code 73
FAILED python/test/unit/fem/test_vector_function.py::test_div_conforming_tetrahedron[5-RT] - petsc4py.PETSc.Error: error code 73
=================================================================================================================================== 10 failed in 1.91s ===================================================================================================================================
More information about the debian-science-maintainers
mailing list