8084 tests run, 7292 passed, 786 skipped, 6 failed.
Annotations
Check failure on line 282 in tests/vertexonly/test_interpolation_from_parent.py
github-actions / Firedrake real
test_interpolation_from_parent.test_tensor_function_interpolation[shiftedmesh-mesh-100-coords-FunctionSpace(Regge2)]
AssertionError: assert False
+ where False = <function allclose at 0x7f5cb55237f0>(array([[[ 1.43674311e-01, 1.24372901e-01],\n [ 1.24372901e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 1.49278171e-02],\n [ 1.49278171e-02, 1.51586555e-01]],\n\n [[ 3.31461468e-02, 5.06537791e-02],\n [ 5.06537791e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.76284326e-02],\n [ 1.76284326e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.93939794e-02],\n [-4.93939794e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.43758718e-02],\n [ 1.43758718e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.27013888e-01],\n [-1.27013888e-01, 4.83806987e-02]],\n\n [[ 3.67346949e-03, -1.73496851e-02],\n [-1.73496851e-02, 6.73957309e-02]],\n\n [[ 9.28985153e-03, -2.38406258e-02],\n [-2.38406258e-02, 6.89296471e-02]],\n\n [[ 4.12457457e-02, -2.77278213e-02],\n [-2.77278213e-02, 9.32473200e-03]],\n\n [[ 7.06210354e-03, -1.01052981e-02],\n [-1.01052981e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.46538578e-01],\n [-1.46538578e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.14803716e-04],\n [-3.14803716e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.37492638e-01],\n [-1.37492638e-01, 1.36037967e-01]],\n\n [[ 2.15699806e-01, -2.51198051e-01],\n [-2.51198051e-01, 1.53134491e-01]],\n\n [[ 9.80934964e-03, -2.61702863e-02],\n [-2.61702863e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, 2.76967696e-01],\n [ 2.76967696e-01, 9.74182567e-02]],\n\n [[ 7.16363728e-02, -3.19287889e-02],\n [-3.19287889e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.85390569e-02],\n [ 2.85390569e-02, 5.17205577e-02]],\n\n [[ 1.76037448e-01, -1.03851951e-02],\n [-1.03851951e-02, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.11319594e-02],\n [ 4.11319594e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.39663046e-02],\n [ 8.39663046e-02, 1.39923807e-01]],\n\n [[ 3.13298757e-01, 1.87898647e-01],\n [ 1.87898647e-01, 1.08454218e-01]]]), array([[[ 1.28112370e-01, 1.16591930e-01],\n [ 1.16591930e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 4.62988588e-02],\n [ 4.62988588e-02, 2.14328639e-01]],\n\n [[ 3.31461468e-02, 5.08404361e-02],\n [ 5.08404361e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.91865200e-02],\n [ 1.91865200e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.76152739e-02],\n [-4.76152739e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.45889613e-02],\n [ 1.45889613e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.04824895e-01],\n [-1.04824895e-01, 4.83806987e-02]],\n\n [[ 4.30518717e-03, -1.70338262e-02],\n [-1.70338262e-02, 6.73957309e-02]],\n\n [[ 8.51567131e-03, -2.42277159e-02],\n [-2.42277159e-02, 6.89296471e-02]],\n\n [[ 5.24647599e-02, -2.21183142e-02],\n [-2.21183142e-02, 9.32473200e-03]],\n\n [[ 6.21136065e-03, -1.05306696e-02],\n [-1.05306696e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.13695712e-01],\n [-1.13695712e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.93571046e-04],\n [-3.93571046e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.41610458e-01],\n [-1.41610458e-01, 1.36037967e-01]],\n\n [[ 2.93843368e-01, -2.12126270e-01],\n [-2.12126270e-01, 1.53134491e-01]],\n\n [[ 1.14137763e-02, -2.53680729e-02],\n [-2.53680729e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, -1.70571807e-01],\n [-1.70571807e-01, 9.74182567e-02]],\n\n [[ 8.78779331e-02, -2.38080087e-02],\n [-2.38080087e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.97596329e-02],\n [ 2.97596329e-02, 5.17205577e-02]],\n\n [[ 1.87161402e-01, -4.82321801e-03],\n [-4.82321801e-03, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.46684825e-02],\n [ 4.46684825e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.19259064e-02],\n [ 8.19259064e-02, 1.39923807e-01]],\n\n [[ 2.95606040e-01, 1.79052288e-01],\n [ 1.79052288e-01, 1.08454218e-01]]]))
+ where <function allclose at 0x7f5cb55237f0> = np.allclose
+ and array([[[ 1.43674311e-01, 1.24372901e-01],\n [ 1.24372901e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 1.49278171e-02],\n [ 1.49278171e-02, 1.51586555e-01]],\n\n [[ 3.31461468e-02, 5.06537791e-02],\n [ 5.06537791e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.76284326e-02],\n [ 1.76284326e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.93939794e-02],\n [-4.93939794e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.43758718e-02],\n [ 1.43758718e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.27013888e-01],\n [-1.27013888e-01, 4.83806987e-02]],\n\n [[ 3.67346949e-03, -1.73496851e-02],\n [-1.73496851e-02, 6.73957309e-02]],\n\n [[ 9.28985153e-03, -2.38406258e-02],\n [-2.38406258e-02, 6.89296471e-02]],\n\n [[ 4.12457457e-02, -2.77278213e-02],\n [-2.77278213e-02, 9.32473200e-03]],\n\n [[ 7.06210354e-03, -1.01052981e-02],\n [-1.01052981e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.46538578e-01],\n [-1.46538578e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.14803716e-04],\n [-3.14803716e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.37492638e-01],\n [-1.37492638e-01, 1.36037967e-01]],\n\n [[ 2.15699806e-01, -2.51198051e-01],\n [-2.51198051e-01, 1.53134491e-01]],\n\n [[ 9.80934964e-03, -2.61702863e-02],\n [-2.61702863e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, 2.76967696e-01],\n [ 2.76967696e-01, 9.74182567e-02]],\n\n [[ 7.16363728e-02, -3.19287889e-02],\n [-3.19287889e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.85390569e-02],\n [ 2.85390569e-02, 5.17205577e-02]],\n\n [[ 1.76037448e-01, -1.03851951e-02],\n [-1.03851951e-02, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.11319594e-02],\n [ 4.11319594e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.39663046e-02],\n [ 8.39663046e-02, 1.39923807e-01]],\n\n [[ 3.13298757e-01, 1.87898647e-01],\n [ 1.87898647e-01, 1.08454218e-01]]]) = <built-in method reshape of numpy.ndarray object at 0x7f5c7ce7d410>((23, 2, 2))
+ where <built-in method reshape of numpy.ndarray object at 0x7f5c7ce7d410> = array([[[ 1.43674311e-01, 1.24372901e-01],\n [ 1.24372901e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 1.49278171e-02],\n [ 1.49278171e-02, 1.51586555e-01]],\n\n [[ 3.31461468e-02, 5.06537791e-02],\n [ 5.06537791e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.76284326e-02],\n [ 1.76284326e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.93939794e-02],\n [-4.93939794e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.43758718e-02],\n [ 1.43758718e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.27013888e-01],\n [-1.27013888e-01, 4.83806987e-02]],\n\n [[ 3.67346949e-03, -1.73496851e-02],\n [-1.73496851e-02, 6.73957309e-02]],\n\n [[ 9.28985153e-03, -2.38406258e-02],\n [-2.38406258e-02, 6.89296471e-02]],\n\n [[ 4.12457457e-02, -2.77278213e-02],\n [-2.77278213e-02, 9.32473200e-03]],\n\n [[ 7.06210354e-03, -1.01052981e-02],\n [-1.01052981e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.46538578e-01],\n [-1.46538578e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.14803716e-04],\n [-3.14803716e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.37492638e-01],\n [-1.37492638e-01, 1.36037967e-01]],\n\n [[ 2.15699806e-01, -2.51198051e-01],\n [-2.51198051e-01, 1.53134491e-01]],\n\n [[ 9.80934964e-03, -2.61702863e-02],\n [-2.61702863e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, 2.76967696e-01],\n [ 2.76967696e-01, 9.74182567e-02]],\n\n [[ 7.16363728e-02, -3.19287889e-02],\n [-3.19287889e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.85390569e-02],\n [ 2.85390569e-02, 5.17205577e-02]],\n\n [[ 1.76037448e-01, -1.03851951e-02],\n [-1.03851951e-02, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.11319594e-02],\n [ 4.11319594e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.39663046e-02],\n [ 8.39663046e-02, 1.39923807e-01]],\n\n [[ 3.13298757e-01, 1.87898647e-01],\n [ 1.87898647e-01, 1.08454218e-01]]]).reshape
+ where array([[[ 1.43674311e-01, 1.24372901e-01],\n [ 1.24372901e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 1.49278171e-02],\n [ 1.49278171e-02, 1.51586555e-01]],\n\n [[ 3.31461468e-02, 5.06537791e-02],\n [ 5.06537791e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.76284326e-02],\n [ 1.76284326e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.93939794e-02],\n [-4.93939794e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.43758718e-02],\n [ 1.43758718e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.27013888e-01],\n [-1.27013888e-01, 4.83806987e-02]],\n\n [[ 3.67346949e-03, -1.73496851e-02],\n [-1.73496851e-02, 6.73957309e-02]],\n\n [[ 9.28985153e-03, -2.38406258e-02],\n [-2.38406258e-02, 6.89296471e-02]],\n\n [[ 4.12457457e-02, -2.77278213e-02],\n [-2.77278213e-02, 9.32473200e-03]],\n\n [[ 7.06210354e-03, -1.01052981e-02],\n [-1.01052981e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.46538578e-01],\n [-1.46538578e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.14803716e-04],\n [-3.14803716e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.37492638e-01],\n [-1.37492638e-01, 1.36037967e-01]],\n\n [[ 2.15699806e-01, -2.51198051e-01],\n [-2.51198051e-01, 1.53134491e-01]],\n\n [[ 9.80934964e-03, -2.61702863e-02],\n [-2.61702863e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, 2.76967696e-01],\n [ 2.76967696e-01, 9.74182567e-02]],\n\n [[ 7.16363728e-02, -3.19287889e-02],\n [-3.19287889e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.85390569e-02],\n [ 2.85390569e-02, 5.17205577e-02]],\n\n [[ 1.76037448e-01, -1.03851951e-02],\n [-1.03851951e-02, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.11319594e-02],\n [ 4.11319594e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.39663046e-02],\n [ 8.39663046e-02, 1.39923807e-01]],\n\n [[ 3.13298757e-01, 1.87898647e-01],\n [ 1.87898647e-01, 1.08454218e-01]]]) = Dat(DataSet(Set((np.int64(23), np.int64(23), np.int64(23)), 'set_#x7f5c76babd10'), (2, 2), 'None_nodes_dset'), None, dtype('float64'), 'function_3124').data_ro
+ where Dat(DataSet(Set((np.int64(23), np.int64(23), np.int64(23)), 'set_#x7f5c76babd10'), (2, 2), 'None_nodes_dset'), None, dtype('float64'), 'function_3124') = Coefficient(WithGeometry(FunctionSpace(<firedrake.mesh.VertexOnlyMeshTopology object at 0x7f5c76816e10>, TensorElement(FiniteElement('Discontinuous Lagrange', vertex, 0), shape=(2, 2), symmetry={}), name=None), Mesh(VectorElement(FiniteElement('Discontinuous Lagrange', Cell(vertex, 2), 0), dim=2), 3118)), 5582).dat
+ and (23, 2, 2) = array([[[ 1.28112370e-01, 1.16591930e-01],\n [ 1.16591930e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 4.62988588e-02],\n [ 4.62988588e-02, 2.14328639e-01]],\n\n [[ 3.31461468e-02, 5.08404361e-02],\n [ 5.08404361e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.91865200e-02],\n [ 1.91865200e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.76152739e-02],\n [-4.76152739e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.45889613e-02],\n [ 1.45889613e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.04824895e-01],\n [-1.04824895e-01, 4.83806987e-02]],\n\n [[ 4.30518717e-03, -1.70338262e-02],\n [-1.70338262e-02, 6.73957309e-02]],\n\n [[ 8.51567131e-03, -2.42277159e-02],\n [-2.42277159e-02, 6.89296471e-02]],\n\n [[ 5.24647599e-02, -2.21183142e-02],\n [-2.21183142e-02, 9.32473200e-03]],\n\n [[ 6.21136065e-03, -1.05306696e-02],\n [-1.05306696e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.13695712e-01],\n [-1.13695712e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.93571046e-04],\n [-3.93571046e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.41610458e-01],\n [-1.41610458e-01, 1.36037967e-01]],\n\n [[ 2.93843368e-01, -2.12126270e-01],\n [-2.12126270e-01, 1.53134491e-01]],\n\n [[ 1.14137763e-02, -2.53680729e-02],\n [-2.53680729e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, -1.70571807e-01],\n [-1.70571807e-01, 9.74182567e-02]],\n\n [[ 8.78779331e-02, -2.38080087e-02],\n [-2.38080087e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.97596329e-02],\n [ 2.97596329e-02, 5.17205577e-02]],\n\n [[ 1.87161402e-01, -4.82321801e-03],\n [-4.82321801e-03, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.46684825e-02],\n [ 4.46684825e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.19259064e-02],\n [ 8.19259064e-02, 1.39923807e-01]],\n\n [[ 2.95606040e-01, 1.79052288e-01],\n [ 1.79052288e-01, 1.08454218e-01]]]).shape
Raw output
parentmesh = Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 3113)
vertexcoords = array([[-0.35792788, -0.3257414 ],
[-0.10000695, -0.46295641],
[-0.18206083, -0.27924972],
[-0.07...-0.01114882],
[ 0.13596636, 0.32852599],
[ 0.2190158 , 0.37406391],
[ 0.54369664, 0.32932388]])
tfs = ('Regge', 2, <cyfunction FunctionSpace at 0x7f5ca7a97ac0>)
def test_tensor_function_interpolation(parentmesh, vertexcoords, tfs):
if parentmesh.name == "immersedsphere":
vertexcoords = immersed_sphere_vertexcoords(parentmesh, vertexcoords)
tfs_fam, tfs_deg, tfs_typ = tfs
vm = VertexOnlyMesh(parentmesh, vertexcoords, missing_points_behaviour=None)
vertexcoords = vm.coordinates.dat.data_ro
if (
parentmesh.coordinates.function_space().ufl_element().family()
== "Discontinuous Lagrange"
):
pytest.skip(f"Interpolating into f{tfs_fam} on a periodic mesh is not well-defined.")
V = tfs_typ(parentmesh, tfs_fam, tfs_deg)
W = TensorFunctionSpace(vm, "DG", 0)
x = SpatialCoordinate(parentmesh)
# use outer product to check Regge works
expr = outer(x, x)
assert W.shape == expr.ufl_shape
v = Function(V).interpolate(expr)
result = np.asarray([np.outer(vertexcoords[i], vertexcoords[i]) for i in range(len(vertexcoords))])
if len(result) == 0:
result = result.reshape(vertexcoords.shape + (parentmesh.geometric_dimension(),))
w_v = assemble(interpolate(v, W))
> assert np.allclose(w_v.dat.data_ro.reshape(result.shape), result)
E AssertionError: assert False
E + where False = <function allclose at 0x7f5cb55237f0>(array([[[ 1.43674311e-01, 1.24372901e-01],\n [ 1.24372901e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 1.49278171e-02],\n [ 1.49278171e-02, 1.51586555e-01]],\n\n [[ 3.31461468e-02, 5.06537791e-02],\n [ 5.06537791e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.76284326e-02],\n [ 1.76284326e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.93939794e-02],\n [-4.93939794e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.43758718e-02],\n [ 1.43758718e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.27013888e-01],\n [-1.27013888e-01, 4.83806987e-02]],\n\n [[ 3.67346949e-03, -1.73496851e-02],\n [-1.73496851e-02, 6.73957309e-02]],\n\n [[ 9.28985153e-03, -2.38406258e-02],\n [-2.38406258e-02, 6.89296471e-02]],\n\n [[ 4.12457457e-02, -2.77278213e-02],\n [-2.77278213e-02, 9.32473200e-03]],\n\n [[ 7.06210354e-03, -1.01052981e-02],\n [-1.01052981e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.46538578e-01],\n [-1.46538578e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.14803716e-04],\n [-3.14803716e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.37492638e-01],\n [-1.37492638e-01, 1.36037967e-01]],\n\n [[ 2.15699806e-01, -2.51198051e-01],\n [-2.51198051e-01, 1.53134491e-01]],\n\n [[ 9.80934964e-03, -2.61702863e-02],\n [-2.61702863e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, 2.76967696e-01],\n [ 2.76967696e-01, 9.74182567e-02]],\n\n [[ 7.16363728e-02, -3.19287889e-02],\n [-3.19287889e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.85390569e-02],\n [ 2.85390569e-02, 5.17205577e-02]],\n\n [[ 1.76037448e-01, -1.03851951e-02],\n [-1.03851951e-02, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.11319594e-02],\n [ 4.11319594e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.39663046e-02],\n [ 8.39663046e-02, 1.39923807e-01]],\n\n [[ 3.13298757e-01, 1.87898647e-01],\n [ 1.87898647e-01, 1.08454218e-01]]]), array([[[ 1.28112370e-01, 1.16591930e-01],\n [ 1.16591930e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 4.62988588e-02],\n [ 4.62988588e-02, 2.14328639e-01]],\n\n [[ 3.31461468e-02, 5.08404361e-02],\n [ 5.08404361e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.91865200e-02],\n [ 1.91865200e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.76152739e-02],\n [-4.76152739e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.45889613e-02],\n [ 1.45889613e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.04824895e-01],\n [-1.04824895e-01, 4.83806987e-02]],\n\n [[ 4.30518717e-03, -1.70338262e-02],\n [-1.70338262e-02, 6.73957309e-02]],\n\n [[ 8.51567131e-03, -2.42277159e-02],\n [-2.42277159e-02, 6.89296471e-02]],\n\n [[ 5.24647599e-02, -2.21183142e-02],\n [-2.21183142e-02, 9.32473200e-03]],\n\n [[ 6.21136065e-03, -1.05306696e-02],\n [-1.05306696e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.13695712e-01],\n [-1.13695712e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.93571046e-04],\n [-3.93571046e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.41610458e-01],\n [-1.41610458e-01, 1.36037967e-01]],\n\n [[ 2.93843368e-01, -2.12126270e-01],\n [-2.12126270e-01, 1.53134491e-01]],\n\n [[ 1.14137763e-02, -2.53680729e-02],\n [-2.53680729e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, -1.70571807e-01],\n [-1.70571807e-01, 9.74182567e-02]],\n\n [[ 8.78779331e-02, -2.38080087e-02],\n [-2.38080087e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.97596329e-02],\n [ 2.97596329e-02, 5.17205577e-02]],\n\n [[ 1.87161402e-01, -4.82321801e-03],\n [-4.82321801e-03, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.46684825e-02],\n [ 4.46684825e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.19259064e-02],\n [ 8.19259064e-02, 1.39923807e-01]],\n\n [[ 2.95606040e-01, 1.79052288e-01],\n [ 1.79052288e-01, 1.08454218e-01]]]))
E + where <function allclose at 0x7f5cb55237f0> = np.allclose
E + and array([[[ 1.43674311e-01, 1.24372901e-01],\n [ 1.24372901e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 1.49278171e-02],\n [ 1.49278171e-02, 1.51586555e-01]],\n\n [[ 3.31461468e-02, 5.06537791e-02],\n [ 5.06537791e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.76284326e-02],\n [ 1.76284326e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.93939794e-02],\n [-4.93939794e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.43758718e-02],\n [ 1.43758718e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.27013888e-01],\n [-1.27013888e-01, 4.83806987e-02]],\n\n [[ 3.67346949e-03, -1.73496851e-02],\n [-1.73496851e-02, 6.73957309e-02]],\n\n [[ 9.28985153e-03, -2.38406258e-02],\n [-2.38406258e-02, 6.89296471e-02]],\n\n [[ 4.12457457e-02, -2.77278213e-02],\n [-2.77278213e-02, 9.32473200e-03]],\n\n [[ 7.06210354e-03, -1.01052981e-02],\n [-1.01052981e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.46538578e-01],\n [-1.46538578e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.14803716e-04],\n [-3.14803716e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.37492638e-01],\n [-1.37492638e-01, 1.36037967e-01]],\n\n [[ 2.15699806e-01, -2.51198051e-01],\n [-2.51198051e-01, 1.53134491e-01]],\n\n [[ 9.80934964e-03, -2.61702863e-02],\n [-2.61702863e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, 2.76967696e-01],\n [ 2.76967696e-01, 9.74182567e-02]],\n\n [[ 7.16363728e-02, -3.19287889e-02],\n [-3.19287889e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.85390569e-02],\n [ 2.85390569e-02, 5.17205577e-02]],\n\n [[ 1.76037448e-01, -1.03851951e-02],\n [-1.03851951e-02, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.11319594e-02],\n [ 4.11319594e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.39663046e-02],\n [ 8.39663046e-02, 1.39923807e-01]],\n\n [[ 3.13298757e-01, 1.87898647e-01],\n [ 1.87898647e-01, 1.08454218e-01]]]) = <built-in method reshape of numpy.ndarray object at 0x7f5c7ce7d410>((23, 2, 2))
E + where <built-in method reshape of numpy.ndarray object at 0x7f5c7ce7d410> = array([[[ 1.43674311e-01, 1.24372901e-01],\n [ 1.24372901e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 1.49278171e-02],\n [ 1.49278171e-02, 1.51586555e-01]],\n\n [[ 3.31461468e-02, 5.06537791e-02],\n [ 5.06537791e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.76284326e-02],\n [ 1.76284326e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.93939794e-02],\n [-4.93939794e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.43758718e-02],\n [ 1.43758718e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.27013888e-01],\n [-1.27013888e-01, 4.83806987e-02]],\n\n [[ 3.67346949e-03, -1.73496851e-02],\n [-1.73496851e-02, 6.73957309e-02]],\n\n [[ 9.28985153e-03, -2.38406258e-02],\n [-2.38406258e-02, 6.89296471e-02]],\n\n [[ 4.12457457e-02, -2.77278213e-02],\n [-2.77278213e-02, 9.32473200e-03]],\n\n [[ 7.06210354e-03, -1.01052981e-02],\n [-1.01052981e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.46538578e-01],\n [-1.46538578e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.14803716e-04],\n [-3.14803716e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.37492638e-01],\n [-1.37492638e-01, 1.36037967e-01]],\n\n [[ 2.15699806e-01, -2.51198051e-01],\n [-2.51198051e-01, 1.53134491e-01]],\n\n [[ 9.80934964e-03, -2.61702863e-02],\n [-2.61702863e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, 2.76967696e-01],\n [ 2.76967696e-01, 9.74182567e-02]],\n\n [[ 7.16363728e-02, -3.19287889e-02],\n [-3.19287889e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.85390569e-02],\n [ 2.85390569e-02, 5.17205577e-02]],\n\n [[ 1.76037448e-01, -1.03851951e-02],\n [-1.03851951e-02, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.11319594e-02],\n [ 4.11319594e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.39663046e-02],\n [ 8.39663046e-02, 1.39923807e-01]],\n\n [[ 3.13298757e-01, 1.87898647e-01],\n [ 1.87898647e-01, 1.08454218e-01]]]).reshape
E + where array([[[ 1.43674311e-01, 1.24372901e-01],\n [ 1.24372901e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 1.49278171e-02],\n [ 1.49278171e-02, 1.51586555e-01]],\n\n [[ 3.31461468e-02, 5.06537791e-02],\n [ 5.06537791e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.76284326e-02],\n [ 1.76284326e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.93939794e-02],\n [-4.93939794e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.43758718e-02],\n [ 1.43758718e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.27013888e-01],\n [-1.27013888e-01, 4.83806987e-02]],\n\n [[ 3.67346949e-03, -1.73496851e-02],\n [-1.73496851e-02, 6.73957309e-02]],\n\n [[ 9.28985153e-03, -2.38406258e-02],\n [-2.38406258e-02, 6.89296471e-02]],\n\n [[ 4.12457457e-02, -2.77278213e-02],\n [-2.77278213e-02, 9.32473200e-03]],\n\n [[ 7.06210354e-03, -1.01052981e-02],\n [-1.01052981e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.46538578e-01],\n [-1.46538578e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.14803716e-04],\n [-3.14803716e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.37492638e-01],\n [-1.37492638e-01, 1.36037967e-01]],\n\n [[ 2.15699806e-01, -2.51198051e-01],\n [-2.51198051e-01, 1.53134491e-01]],\n\n [[ 9.80934964e-03, -2.61702863e-02],\n [-2.61702863e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, 2.76967696e-01],\n [ 2.76967696e-01, 9.74182567e-02]],\n\n [[ 7.16363728e-02, -3.19287889e-02],\n [-3.19287889e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.85390569e-02],\n [ 2.85390569e-02, 5.17205577e-02]],\n\n [[ 1.76037448e-01, -1.03851951e-02],\n [-1.03851951e-02, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.11319594e-02],\n [ 4.11319594e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.39663046e-02],\n [ 8.39663046e-02, 1.39923807e-01]],\n\n [[ 3.13298757e-01, 1.87898647e-01],\n [ 1.87898647e-01, 1.08454218e-01]]]) = Dat(DataSet(Set((np.int64(23), np.int64(23), np.int64(23)), 'set_#x7f5c76babd10'), (2, 2), 'None_nodes_dset'), None, dtype('float64'), 'function_3124').data_ro
E + where Dat(DataSet(Set((np.int64(23), np.int64(23), np.int64(23)), 'set_#x7f5c76babd10'), (2, 2), 'None_nodes_dset'), None, dtype('float64'), 'function_3124') = Coefficient(WithGeometry(FunctionSpace(<firedrake.mesh.VertexOnlyMeshTopology object at 0x7f5c76816e10>, TensorElement(FiniteElement('Discontinuous Lagrange', vertex, 0), shape=(2, 2), symmetry={}), name=None), Mesh(VectorElement(FiniteElement('Discontinuous Lagrange', Cell(vertex, 2), 0), dim=2), 3118)), 5582).dat
E + and (23, 2, 2) = array([[[ 1.28112370e-01, 1.16591930e-01],\n [ 1.16591930e-01, 1.06107460e-01]],\n\n [[ 1.00013901e-02, 4.62988588e-02],\n [ 4.62988588e-02, 2.14328639e-01]],\n\n [[ 3.31461468e-02, 5.08404361e-02],\n [ 5.08404361e-02, 7.79804048e-02]],\n\n [[ 6.27816587e-03, 1.91865200e-02],\n [ 1.91865200e-02, 5.86353656e-02]],\n\n [[ 2.13087119e-01, -4.76152739e-02],\n [-4.76152739e-02, 1.06398469e-02]],\n\n [[ 6.76448476e-03, 1.45889613e-02],\n [ 1.45889613e-02, 3.14640065e-02]],\n\n [[ 2.27120711e-01, -1.04824895e-01],\n [-1.04824895e-01, 4.83806987e-02]],\n\n [[ 4.30518717e-03, -1.70338262e-02],\n [-1.70338262e-02, 6.73957309e-02]],\n\n [[ 8.51567131e-03, -2.42277159e-02],\n [-2.42277159e-02, 6.89296471e-02]],\n\n [[ 5.24647599e-02, -2.21183142e-02],\n [-2.21183142e-02, 9.32473200e-03]],\n\n [[ 6.21136065e-03, -1.05306696e-02],\n [-1.05306696e-02, 1.78535764e-02]],\n\n [[ 9.88263557e-02, -1.13695712e-01],\n [-1.13695712e-01, 1.30802302e-01]],\n\n [[ 2.04039612e-03, -3.93571046e-04],\n [-3.93571046e-04, 7.59157337e-05]],\n\n [[ 1.47411213e-01, -1.41610458e-01],\n [-1.41610458e-01, 1.36037967e-01]],\n\n [[ 2.93843368e-01, -2.12126270e-01],\n [-2.12126270e-01, 1.53134491e-01]],\n\n [[ 1.14137763e-02, -2.53680729e-02],\n [-2.53680729e-02, 5.63826646e-02]],\n\n [[ 2.98657995e-01, -1.70571807e-01],\n [-1.70571807e-01, 9.74182567e-02]],\n\n [[ 8.78779331e-02, -2.38080087e-02],\n [-2.38080087e-02, 6.45009799e-03]],\n\n [[ 1.71234764e-02, 2.97596329e-02],\n [ 2.97596329e-02, 5.17205577e-02]],\n\n [[ 1.87161402e-01, -4.82321801e-03],\n [-4.82321801e-03, 1.24296098e-04]],\n\n [[ 1.84868507e-02, 4.46684825e-02],\n [ 4.46684825e-02, 1.07929325e-01]],\n\n [[ 4.79679212e-02, 8.19259064e-02],\n [ 8.19259064e-02, 1.39923807e-01]],\n\n [[ 2.95606040e-01, 1.79052288e-01],\n [ 1.79052288e-01, 1.08454218e-01]]]).shape
tests/vertexonly/test_interpolation_from_parent.py:282: AssertionError
Check failure on line 534 in tests/regression/test_interpolate_cross_mesh.py
github-actions / Firedrake real
test_interpolate_cross_mesh.test_interpolate_cross_mesh[unitsquare_Regge_source]
assert False
+ where False = <function allclose at 0x7f09c590f930>([array([[0.32470945, 0.34433209],\n [0.34433209, 0.36 ]]), array([[0.01 , 0.09404452],\n [0.09404452, 0.81 ]]), array([[0.81 , 0.14692594],\n [0.14692594, 0.01 ]]), array([[0.81, 0.81],\n [0.81, 0.81]]), array([[0.49240298, 0.44222054],\n [0.44222054, 0.43349056]]), array([[2.77555756e-17, 4.13576073e-02],\n [4.13576073e-02, 4.68375339e-17]]), ...], array([[[0.3136 , 0.336 ],\n [0.336 , 0.36 ]],\n\n [[0.01 , 0.09 ],\n [0.09 , 0.81 ]],\n\n [[0.81 , 0.09 ],\n [0.09 , 0.01 ]],\n\n [[0.81 , 0.81 ],\n [0.81 , 0.81 ]],\n\n [[0.527076 , 0.4779984 ],\n [0.4779984 , 0.43349056]],\n\n [[0. , 0. ],\n [0. , 0. ]],\n\n [[0. , 0. ],\n [0. , 0.11111111]],\n\n [[0.25 , 0. ],\n [0. , 0. ]],\n\n [[0.25 , 0.16666667],\n [0.16666667, 0.11111111]],\n\n [[1. , 0. ],\n [0. , 0. ]],\n\n [[0. , 0. ],\n [0. , 0.44444444]],\n\n [[1. , 0.33333333],\n [0.33333333, 0.11111111]],\n\n [[0.25 , 0.33333333],\n [0.33333333, 0.44444444]],\n\n [[0. , 0. ],\n [0. , 1. ]],\n\n [[1. , 0.66666667],\n [0.66666667, 0.44444444]],\n\n [[0.25 , 0.5 ],\n [0.5 , 1. ]],\n\n [[1. , 1. ],\n [1. , 1. ]],\n\n [... ]],\n\n [[0.11111111, 0.13333333],\n [0.13333333, 0.16 ]],\n\n [[1. , 0.2 ],\n [0.2 , 0.04 ]],\n\n [[1. , 0. ],\n [0. , 0. ]],\n\n [[0.44444444, 0.26666667],\n [0.26666667, 0.16 ]],\n\n [[0. , 0. ],\n [0. , 0.36 ]],\n\n [[0.11111111, 0.2 ],\n [0.2 , 0.36 ]],\n\n [[1. , 0.4 ],\n [0.4 , 0.16 ]],\n\n [[0.44444444, 0.4 ],\n [0.4 , 0.36 ]],\n\n [[0. , 0. ],\n [0. , 0.64 ]],\n\n [[0.11111111, 0.26666667],\n [0.26666667, 0.64 ]],\n\n [[1. , 0.6 ],\n [0.6 , 0.36 ]],\n\n [[0.44444444, 0.53333333],\n [0.53333333, 0.64 ]],\n\n [[0. , 0. ],\n [0. , 1. ]],\n\n [[0.11111111, 0.33333333],\n [0.33333333, 1. ]],\n\n [[1. , 0.8 ],\n [0.8 , 0.64 ]],\n\n [[0.44444444, 0.66666667],\n [0.66666667, 1. ]],\n\n [[1. , 1. ],\n [1. , 1. ]]]), atol=1e-08)
+ where <function allclose at 0x7f09c590f930> = np.allclose
Raw output
parameters = (Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 6499), Mesh(VectorElement(FiniteElement('Q', quadr...66667],
[0.66666667, 1. ]],
[[1. , 1. ],
[1. , 1. ]]]), ...)
def test_interpolate_cross_mesh(parameters):
(
m_src,
m_dest,
coords,
expr_src,
expr_dest,
expected,
V_src,
V_dest,
V_dest_2,
) = parameters
> get_expected_values(
m_src, m_dest, V_src, V_dest, coords, expected, expr_src, expr_dest
)
tests/regression/test_interpolate_cross_mesh.py:299:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
tests/regression/test_interpolate_cross_mesh.py:501: in get_expected_values
interpolate_function(
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
m_src = Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 6499)
m_dest = Mesh(VectorElement(FiniteElement('Q', quadrilateral, 1), dim=2), 6500)
V_src = WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0973e97110>, FiniteElement('Regge', triangle, 2), name=None), Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 6499))
V_dest = WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0973e967b0>, TensorElement(FiniteElement('Q', qu...al, 4), shape=(2, 2), symmetry={}), name=None), Mesh(VectorElement(FiniteElement('Q', quadrilateral, 1), dim=2), 6500))
coords = array([[0.56 , 0.6 ],
[0.1 , 0.9 ],
[0.9 , 0.1 ],
[0.9 ,...33333, 1. ],
[1. , 0.8 ],
[0.66666667, 1. ],
[1. , 1. ]])
expected = array([[[0.3136 , 0.336 ],
[0.336 , 0.36 ]],
[[0.01 , 0.09 ],
[0.09 ... 0.66666667],
[0.66666667, 1. ]],
[[1. , 1. ],
[1. , 1. ]]])
expr_src = Outer(SpatialCoordinate(Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 6499)), SpatialCoordinate(Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 6499)))
expr_dest = Outer(SpatialCoordinate(Mesh(VectorElement(FiniteElement('Q', quadrilateral, 1), dim=2), 6500)), SpatialCoordinate(Mesh(VectorElement(FiniteElement('Q', quadrilateral, 1), dim=2), 6500)))
atol = 1e-08
def interpolate_function(
m_src, m_dest, V_src, V_dest, coords, expected, expr_src, expr_dest, atol
):
f_src = Function(V_src).interpolate(expr_src)
f_dest = assemble(interpolate(f_src, V_dest))
assert extract_unique_domain(f_dest) is m_dest
got = f_dest.at(coords)
> assert np.allclose(got, expected, atol=atol)
E assert False
E + where False = <function allclose at 0x7f09c590f930>([array([[0.32470945, 0.34433209],\n [0.34433209, 0.36 ]]), array([[0.01 , 0.09404452],\n [0.09404452, 0.81 ]]), array([[0.81 , 0.14692594],\n [0.14692594, 0.01 ]]), array([[0.81, 0.81],\n [0.81, 0.81]]), array([[0.49240298, 0.44222054],\n [0.44222054, 0.43349056]]), array([[2.77555756e-17, 4.13576073e-02],\n [4.13576073e-02, 4.68375339e-17]]), ...], array([[[0.3136 , 0.336 ],\n [0.336 , 0.36 ]],\n\n [[0.01 , 0.09 ],\n [0.09 , 0.81 ]],\n\n [[0.81 , 0.09 ],\n [0.09 , 0.01 ]],\n\n [[0.81 , 0.81 ],\n [0.81 , 0.81 ]],\n\n [[0.527076 , 0.4779984 ],\n [0.4779984 , 0.43349056]],\n\n [[0. , 0. ],\n [0. , 0. ]],\n\n [[0. , 0. ],\n [0. , 0.11111111]],\n\n [[0.25 , 0. ],\n [0. , 0. ]],\n\n [[0.25 , 0.16666667],\n [0.16666667, 0.11111111]],\n\n [[1. , 0. ],\n [0. , 0. ]],\n\n [[0. , 0. ],\n [0. , 0.44444444]],\n\n [[1. , 0.33333333],\n [0.33333333, 0.11111111]],\n\n [[0.25 , 0.33333333],\n [0.33333333, 0.44444444]],\n\n [[0. , 0. ],\n [0. , 1. ]],\n\n [[1. , 0.66666667],\n [0.66666667, 0.44444444]],\n\n [[0.25 , 0.5 ],\n [0.5 , 1. ]],\n\n [[1. , 1. ],\n [1. , 1. ]],\n\n [... ]],\n\n [[0.11111111, 0.13333333],\n [0.13333333, 0.16 ]],\n\n [[1. , 0.2 ],\n [0.2 , 0.04 ]],\n\n [[1. , 0. ],\n [0. , 0. ]],\n\n [[0.44444444, 0.26666667],\n [0.26666667, 0.16 ]],\n\n [[0. , 0. ],\n [0. , 0.36 ]],\n\n [[0.11111111, 0.2 ],\n [0.2 , 0.36 ]],\n\n [[1. , 0.4 ],\n [0.4 , 0.16 ]],\n\n [[0.44444444, 0.4 ],\n [0.4 , 0.36 ]],\n\n [[0. , 0. ],\n [0. , 0.64 ]],\n\n [[0.11111111, 0.26666667],\n [0.26666667, 0.64 ]],\n\n [[1. , 0.6 ],\n [0.6 , 0.36 ]],\n\n [[0.44444444, 0.53333333],\n [0.53333333, 0.64 ]],\n\n [[0. , 0. ],\n [0. , 1. ]],\n\n [[0.11111111, 0.33333333],\n [0.33333333, 1. ]],\n\n [[1. , 0.8 ],\n [0.8 , 0.64 ]],\n\n [[0.44444444, 0.66666667],\n [0.66666667, 1. ]],\n\n [[1. , 1. ],\n [1. , 1. ]]]), atol=1e-08)
E + where <function allclose at 0x7f09c590f930> = np.allclose
tests/regression/test_interpolate_cross_mesh.py:534: AssertionError
Check failure on line 1 in tests/vertexonly/test_interpolation_from_parent.py
github-actions / Firedrake real
test_interpolation_from_parent.test_tensor_function_interpolation_parallel[shiftedmesh-mesh-100-coords-FunctionSpace(Regge2)]
subprocess.CalledProcessError: Command '['mpiexec', '-n', '1', '-genv', '_PYTEST_MPI_CHILD_PROCESS', '1', '/__w/firedrake/firedrake_venv/bin/python', '-m', 'pytest', '--runxfail', '-s', '-q', '/__w/firedrake/firedrake/tests/vertexonly/test_interpolation_from_parent.py::test_tensor_function_interpolation_parallel[shiftedmesh-mesh-100-coords-FunctionSpace(Regge2)]', ':', '-n', '2', '/__w/firedrake/firedrake_venv/bin/python', '-m', 'pytest', '--runxfail', '-s', '-q', '/__w/firedrake/firedrake/tests/vertexonly/test_interpolation_from_parent.py::test_tensor_function_interpolation_parallel[shiftedmesh-mesh-100-coords-FunctionSpace(Regge2)]', '--tb=no', '--no-summary', '--no-header', '--disable-warnings', '--show-capture=no']' returned non-zero exit status 1.
Raw output
args = ()
kwargs = {'parentmesh': Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 973), 'tfs': ('Regge', 2, <cyfunctio...1.38874478],
[ 0.97910159, 0.48091762],
[-0.04517074, 0.00871296],
[-0.38394168, 0.36883325]])}
def parallel_callback(*args, **kwargs):
> subprocess.run(cmd, check=True)
../firedrake_venv/src/pytest-mpi/pytest_mpi.py:233:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
input = None, capture_output = False, timeout = None, check = True
popenargs = (['mpiexec', '-n', '1', '-genv', '_PYTEST_MPI_CHILD_PROCESS', '1', ...],)
kwargs = {}
process = <Popen: returncode: 1 args: ['mpiexec', '-n', '1', '-genv', '_PYTEST_MPI_CHI...>
stdout = None, stderr = None, retcode = 1
def run(*popenargs,
input=None, capture_output=False, timeout=None, check=False, **kwargs):
"""Run command with arguments and return a CompletedProcess instance.
The returned instance will have attributes args, returncode, stdout and
stderr. By default, stdout and stderr are not captured, and those attributes
will be None. Pass stdout=PIPE and/or stderr=PIPE in order to capture them,
or pass capture_output=True to capture both.
If check is True and the exit code was non-zero, it raises a
CalledProcessError. The CalledProcessError object will have the return code
in the returncode attribute, and output & stderr attributes if those streams
were captured.
If timeout is given, and the process takes too long, a TimeoutExpired
exception will be raised.
There is an optional argument "input", allowing you to
pass bytes or a string to the subprocess's stdin. If you use this argument
you may not also use the Popen constructor's "stdin" argument, as
it will be used internally.
By default, all communication is in bytes, and therefore any "input" should
be bytes, and the stdout and stderr will be bytes. If in text mode, any
"input" should be a string, and stdout and stderr will be strings decoded
according to locale encoding, or by "encoding" if set. Text mode is
triggered by setting any of text, encoding, errors or universal_newlines.
The other arguments are the same as for the Popen constructor.
"""
if input is not None:
if kwargs.get('stdin') is not None:
raise ValueError('stdin and input arguments may not both be used.')
kwargs['stdin'] = PIPE
if capture_output:
if kwargs.get('stdout') is not None or kwargs.get('stderr') is not None:
raise ValueError('stdout and stderr arguments may not be used '
'with capture_output.')
kwargs['stdout'] = PIPE
kwargs['stderr'] = PIPE
with Popen(*popenargs, **kwargs) as process:
try:
stdout, stderr = process.communicate(input, timeout=timeout)
except TimeoutExpired as exc:
process.kill()
if _mswindows:
# Windows accumulates the output in a single blocking
# read() call run on child threads, with the timeout
# being done in a join() on those threads. communicate()
# _after_ kill() is required to collect that and add it
# to the exception.
exc.stdout, exc.stderr = process.communicate()
else:
# POSIX _communicate already populated the output so
# far into the TimeoutExpired exception.
process.wait()
raise
except: # Including KeyboardInterrupt, communicate handled that.
process.kill()
# We don't call process.wait() as .__exit__ does that for us.
raise
retcode = process.poll()
if check and retcode:
> raise CalledProcessError(retcode, process.args,
output=stdout, stderr=stderr)
E subprocess.CalledProcessError: Command '['mpiexec', '-n', '1', '-genv', '_PYTEST_MPI_CHILD_PROCESS', '1', '/__w/firedrake/firedrake_venv/bin/python', '-m', 'pytest', '--runxfail', '-s', '-q', '/__w/firedrake/firedrake/tests/vertexonly/test_interpolation_from_parent.py::test_tensor_function_interpolation_parallel[shiftedmesh-mesh-100-coords-FunctionSpace(Regge2)]', ':', '-n', '2', '/__w/firedrake/firedrake_venv/bin/python', '-m', 'pytest', '--runxfail', '-s', '-q', '/__w/firedrake/firedrake/tests/vertexonly/test_interpolation_from_parent.py::test_tensor_function_interpolation_parallel[shiftedmesh-mesh-100-coords-FunctionSpace(Regge2)]', '--tb=no', '--no-summary', '--no-header', '--disable-warnings', '--show-capture=no']' returned non-zero exit status 1.
/usr/lib/python3.12/subprocess.py:571: CalledProcessError
Check failure on line 1 in tests/regression/test_interpolate_cross_mesh.py
github-actions / Firedrake real
test_interpolate_cross_mesh.test_interpolate_cross_mesh_parallel[unitsquare_Regge_source]
subprocess.CalledProcessError: Command '['mpiexec', '-n', '1', '-genv', '_PYTEST_MPI_CHILD_PROCESS', '1', '/__w/firedrake/firedrake_venv/bin/python', '-m', 'pytest', '--runxfail', '-s', '-q', '/__w/firedrake/firedrake/tests/regression/test_interpolate_cross_mesh.py::test_interpolate_cross_mesh_parallel[unitsquare_Regge_source]', ':', '-n', '2', '/__w/firedrake/firedrake_venv/bin/python', '-m', 'pytest', '--runxfail', '-s', '-q', '/__w/firedrake/firedrake/tests/regression/test_interpolate_cross_mesh.py::test_interpolate_cross_mesh_parallel[unitsquare_Regge_source]', '--tb=no', '--no-summary', '--no-header', '--disable-warnings', '--show-capture=no']' returned non-zero exit status 1.
Raw output
args = ()
kwargs = {'parameters': (Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 8246), Mesh(VectorElement(FiniteEle...6667],
[0.66666667, 1. ]],
[[1. , 1. ],
[1. , 1. ]]]), ...)}
def parallel_callback(*args, **kwargs):
> subprocess.run(cmd, check=True)
../firedrake_venv/src/pytest-mpi/pytest_mpi.py:233:
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
input = None, capture_output = False, timeout = None, check = True
popenargs = (['mpiexec', '-n', '1', '-genv', '_PYTEST_MPI_CHILD_PROCESS', '1', ...],)
kwargs = {}
process = <Popen: returncode: 1 args: ['mpiexec', '-n', '1', '-genv', '_PYTEST_MPI_CHI...>
stdout = None, stderr = None, retcode = 1
def run(*popenargs,
input=None, capture_output=False, timeout=None, check=False, **kwargs):
"""Run command with arguments and return a CompletedProcess instance.
The returned instance will have attributes args, returncode, stdout and
stderr. By default, stdout and stderr are not captured, and those attributes
will be None. Pass stdout=PIPE and/or stderr=PIPE in order to capture them,
or pass capture_output=True to capture both.
If check is True and the exit code was non-zero, it raises a
CalledProcessError. The CalledProcessError object will have the return code
in the returncode attribute, and output & stderr attributes if those streams
were captured.
If timeout is given, and the process takes too long, a TimeoutExpired
exception will be raised.
There is an optional argument "input", allowing you to
pass bytes or a string to the subprocess's stdin. If you use this argument
you may not also use the Popen constructor's "stdin" argument, as
it will be used internally.
By default, all communication is in bytes, and therefore any "input" should
be bytes, and the stdout and stderr will be bytes. If in text mode, any
"input" should be a string, and stdout and stderr will be strings decoded
according to locale encoding, or by "encoding" if set. Text mode is
triggered by setting any of text, encoding, errors or universal_newlines.
The other arguments are the same as for the Popen constructor.
"""
if input is not None:
if kwargs.get('stdin') is not None:
raise ValueError('stdin and input arguments may not both be used.')
kwargs['stdin'] = PIPE
if capture_output:
if kwargs.get('stdout') is not None or kwargs.get('stderr') is not None:
raise ValueError('stdout and stderr arguments may not be used '
'with capture_output.')
kwargs['stdout'] = PIPE
kwargs['stderr'] = PIPE
with Popen(*popenargs, **kwargs) as process:
try:
stdout, stderr = process.communicate(input, timeout=timeout)
except TimeoutExpired as exc:
process.kill()
if _mswindows:
# Windows accumulates the output in a single blocking
# read() call run on child threads, with the timeout
# being done in a join() on those threads. communicate()
# _after_ kill() is required to collect that and add it
# to the exception.
exc.stdout, exc.stderr = process.communicate()
else:
# POSIX _communicate already populated the output so
# far into the TimeoutExpired exception.
process.wait()
raise
except: # Including KeyboardInterrupt, communicate handled that.
process.kill()
# We don't call process.wait() as .__exit__ does that for us.
raise
retcode = process.poll()
if check and retcode:
> raise CalledProcessError(retcode, process.args,
output=stdout, stderr=stderr)
E subprocess.CalledProcessError: Command '['mpiexec', '-n', '1', '-genv', '_PYTEST_MPI_CHILD_PROCESS', '1', '/__w/firedrake/firedrake_venv/bin/python', '-m', 'pytest', '--runxfail', '-s', '-q', '/__w/firedrake/firedrake/tests/regression/test_interpolate_cross_mesh.py::test_interpolate_cross_mesh_parallel[unitsquare_Regge_source]', ':', '-n', '2', '/__w/firedrake/firedrake_venv/bin/python', '-m', 'pytest', '--runxfail', '-s', '-q', '/__w/firedrake/firedrake/tests/regression/test_interpolate_cross_mesh.py::test_interpolate_cross_mesh_parallel[unitsquare_Regge_source]', '--tb=no', '--no-summary', '--no-header', '--disable-warnings', '--show-capture=no']' returned non-zero exit status 1.
/usr/lib/python3.12/subprocess.py:571: CalledProcessError
Check failure on line 58 in tests/regression/test_interpolate_vs_project.py
github-actions / Firedrake real
test_interpolate_vs_project.test_interpolate_vs_project[square-FunctionSpace(Regge1)]
AssertionError: assert False
+ where False = <function allclose at 0x7f09c590f930>(array([0.16666667, 0.16666667, 0. , 0.35355339, 0.20412415,\n 0. , 0. , 0.35355339, 0.20412415, 0. ,\n 0.33333333, 0.33333333, 0.35355339, 0.20412415, 0.35355339,\n 0.20412415, 0.16666667, 0.66666667, 0. , 0. ,\n 0. , 1.06066017, 0.20412415, 0.16666667, 0.66666667,\n 0. , 1.06066017, 0.20412415, 0. , 0. ,\n 0. , 0.83333333, 0.33333333, 1.06066017, 0.20412415,\n 0.5 , 0.28867513, 0. , 0.83333333, 0.33333333,\n 0.5 , 0.28867513, 1.06066017, 0.20412415, 0.66666667,\n 0.66666667, 0. , 0. , 0. , 0. ,\n 0.83333333, 0.83333333, 1.5 , 0.28867513, 1.06066017,\n 0.20412415]), array([ 1.67752016e-01, 1.67752017e-01, 4.34139936e-03, 3.55764540e-01,\n 2.08059306e-01, 2.60483878e-02, -1.44431658e-02, 3.51342246e-01,\n 2.05506089e-01, -2.01663799e-02, 3.32743962e-01, 3.15291240e-01,\n 3.52745428e-01, 2.05433671e-01, 3.98563338e-01, 2.28335034e-01,\n 1.77772838e-01, 6.65557825e-01, 7.77964870e-03, 1.11847071e-02,\n -1.96274994e-02, 1.06747413e+00, 2.05342092e-01, 1.66434174e-01,\n 6.66612260e-01, -3.95712710e-04, 1.06071258e+00, 2.04021134e-01,\n 4.98351674e-04, -8.23692783e-05, 2.56390623e-03, 8.33545034e-01,\n 3.34117401e-01, 1.06305001e+00, 2.01871348e-01, 4.93186044e-01,\n 2.89893082e-01, 1.54057892e-04, 8.33371977e-01, 3.33371805e-01,\n 4.99947592e-01, 2.88572124e-01, 1.06069796e+00, 2.04186949e-01,\n 6.67328444e-01, 6.66755900e-01, 9.29823912e-04, -1.07287349e-03,\n -1.25808343e-04, -2.52878707e-04, 8.33270113e-01, 8.33270113e-01,\n 1.49996148e+00, 2.88916374e-01, 1.06068741e+00, 2.04263275e-01]), atol=1e-06)
+ where <function allclose at 0x7f09c590f930> = np.allclose
+ and array([0.16666667, 0.16666667, 0. , 0.35355339, 0.20412415,\n 0. , 0. , 0.35355339, 0.20412415, 0. ,\n 0.33333333, 0.33333333, 0.35355339, 0.20412415, 0.35355339,\n 0.20412415, 0.16666667, 0.66666667, 0. , 0. ,\n 0. , 1.06066017, 0.20412415, 0.16666667, 0.66666667,\n 0. , 1.06066017, 0.20412415, 0. , 0. ,\n 0. , 0.83333333, 0.33333333, 1.06066017, 0.20412415,\n 0.5 , 0.28867513, 0. , 0.83333333, 0.33333333,\n 0.5 , 0.28867513, 1.06066017, 0.20412415, 0.66666667,\n 0.66666667, 0. , 0. , 0. , 0. ,\n 0.83333333, 0.83333333, 1.5 , 0.28867513, 1.06066017,\n 0.20412415]) = Dat(DataSet(Set((np.int64(56), np.int64(56), np.int64(56)), 'set_#x7f097296cb30'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8313').data
+ where Dat(DataSet(Set((np.int64(56), np.int64(56), np.int64(56)), 'set_#x7f097296cb30'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8313') = Coefficient(WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0973118bc0>, FiniteElement('Regge', triangle, 1), name=None), Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 8270)), 13974).dat
+ and array([ 1.67752016e-01, 1.67752017e-01, 4.34139936e-03, 3.55764540e-01,\n 2.08059306e-01, 2.60483878e-02, -1.44431658e-02, 3.51342246e-01,\n 2.05506089e-01, -2.01663799e-02, 3.32743962e-01, 3.15291240e-01,\n 3.52745428e-01, 2.05433671e-01, 3.98563338e-01, 2.28335034e-01,\n 1.77772838e-01, 6.65557825e-01, 7.77964870e-03, 1.11847071e-02,\n -1.96274994e-02, 1.06747413e+00, 2.05342092e-01, 1.66434174e-01,\n 6.66612260e-01, -3.95712710e-04, 1.06071258e+00, 2.04021134e-01,\n 4.98351674e-04, -8.23692783e-05, 2.56390623e-03, 8.33545034e-01,\n 3.34117401e-01, 1.06305001e+00, 2.01871348e-01, 4.93186044e-01,\n 2.89893082e-01, 1.54057892e-04, 8.33371977e-01, 3.33371805e-01,\n 4.99947592e-01, 2.88572124e-01, 1.06069796e+00, 2.04186949e-01,\n 6.67328444e-01, 6.66755900e-01, 9.29823912e-04, -1.07287349e-03,\n -1.25808343e-04, -2.52878707e-04, 8.33270113e-01, 8.33270113e-01,\n 1.49996148e+00, 2.88916374e-01, 1.06068741e+00, 2.04263275e-01]) = Dat(DataSet(Set((np.int64(56), np.int64(56), np.int64(56)), 'set_#x7f097296cb30'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8314').data
+ where Dat(DataSet(Set((np.int64(56), np.int64(56), np.int64(56)), 'set_#x7f097296cb30'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8314') = Coefficient(WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0973118bc0>, FiniteElement('Regge', triangle, 1), name=None), Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 8270)), 13976).dat
Raw output
V = WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0973118bc0>, FiniteElement('Regge', triangle, 1), name=None), Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 8270))
def test_interpolate_vs_project(V):
mesh = V.mesh()
dim = mesh.geometric_dimension()
if dim == 2:
x, y = SpatialCoordinate(mesh)
elif dim == 3:
x, y, z = SpatialCoordinate(mesh)
shape = V.ufl_element().value_shape
if dim == 2:
if len(shape) == 0:
expression = x + y
elif len(shape) == 1:
expression = as_vector([x, y])
elif len(shape) == 2:
expression = as_tensor(([x, y], [x, y]))
elif dim == 3:
if len(shape) == 0:
expression = x + y + z
elif len(shape) == 1:
expression = as_vector([x, y, z])
elif len(shape) == 2:
expression = as_tensor(([x, y, z], [x, y, z], [x, y, z]))
f = assemble(interpolate(expression, V))
expect = project(expression, V)
> assert np.allclose(f.dat.data, expect.dat.data, atol=1e-06)
E AssertionError: assert False
E + where False = <function allclose at 0x7f09c590f930>(array([0.16666667, 0.16666667, 0. , 0.35355339, 0.20412415,\n 0. , 0. , 0.35355339, 0.20412415, 0. ,\n 0.33333333, 0.33333333, 0.35355339, 0.20412415, 0.35355339,\n 0.20412415, 0.16666667, 0.66666667, 0. , 0. ,\n 0. , 1.06066017, 0.20412415, 0.16666667, 0.66666667,\n 0. , 1.06066017, 0.20412415, 0. , 0. ,\n 0. , 0.83333333, 0.33333333, 1.06066017, 0.20412415,\n 0.5 , 0.28867513, 0. , 0.83333333, 0.33333333,\n 0.5 , 0.28867513, 1.06066017, 0.20412415, 0.66666667,\n 0.66666667, 0. , 0. , 0. , 0. ,\n 0.83333333, 0.83333333, 1.5 , 0.28867513, 1.06066017,\n 0.20412415]), array([ 1.67752016e-01, 1.67752017e-01, 4.34139936e-03, 3.55764540e-01,\n 2.08059306e-01, 2.60483878e-02, -1.44431658e-02, 3.51342246e-01,\n 2.05506089e-01, -2.01663799e-02, 3.32743962e-01, 3.15291240e-01,\n 3.52745428e-01, 2.05433671e-01, 3.98563338e-01, 2.28335034e-01,\n 1.77772838e-01, 6.65557825e-01, 7.77964870e-03, 1.11847071e-02,\n -1.96274994e-02, 1.06747413e+00, 2.05342092e-01, 1.66434174e-01,\n 6.66612260e-01, -3.95712710e-04, 1.06071258e+00, 2.04021134e-01,\n 4.98351674e-04, -8.23692783e-05, 2.56390623e-03, 8.33545034e-01,\n 3.34117401e-01, 1.06305001e+00, 2.01871348e-01, 4.93186044e-01,\n 2.89893082e-01, 1.54057892e-04, 8.33371977e-01, 3.33371805e-01,\n 4.99947592e-01, 2.88572124e-01, 1.06069796e+00, 2.04186949e-01,\n 6.67328444e-01, 6.66755900e-01, 9.29823912e-04, -1.07287349e-03,\n -1.25808343e-04, -2.52878707e-04, 8.33270113e-01, 8.33270113e-01,\n 1.49996148e+00, 2.88916374e-01, 1.06068741e+00, 2.04263275e-01]), atol=1e-06)
E + where <function allclose at 0x7f09c590f930> = np.allclose
E + and array([0.16666667, 0.16666667, 0. , 0.35355339, 0.20412415,\n 0. , 0. , 0.35355339, 0.20412415, 0. ,\n 0.33333333, 0.33333333, 0.35355339, 0.20412415, 0.35355339,\n 0.20412415, 0.16666667, 0.66666667, 0. , 0. ,\n 0. , 1.06066017, 0.20412415, 0.16666667, 0.66666667,\n 0. , 1.06066017, 0.20412415, 0. , 0. ,\n 0. , 0.83333333, 0.33333333, 1.06066017, 0.20412415,\n 0.5 , 0.28867513, 0. , 0.83333333, 0.33333333,\n 0.5 , 0.28867513, 1.06066017, 0.20412415, 0.66666667,\n 0.66666667, 0. , 0. , 0. , 0. ,\n 0.83333333, 0.83333333, 1.5 , 0.28867513, 1.06066017,\n 0.20412415]) = Dat(DataSet(Set((np.int64(56), np.int64(56), np.int64(56)), 'set_#x7f097296cb30'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8313').data
E + where Dat(DataSet(Set((np.int64(56), np.int64(56), np.int64(56)), 'set_#x7f097296cb30'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8313') = Coefficient(WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0973118bc0>, FiniteElement('Regge', triangle, 1), name=None), Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 8270)), 13974).dat
E + and array([ 1.67752016e-01, 1.67752017e-01, 4.34139936e-03, 3.55764540e-01,\n 2.08059306e-01, 2.60483878e-02, -1.44431658e-02, 3.51342246e-01,\n 2.05506089e-01, -2.01663799e-02, 3.32743962e-01, 3.15291240e-01,\n 3.52745428e-01, 2.05433671e-01, 3.98563338e-01, 2.28335034e-01,\n 1.77772838e-01, 6.65557825e-01, 7.77964870e-03, 1.11847071e-02,\n -1.96274994e-02, 1.06747413e+00, 2.05342092e-01, 1.66434174e-01,\n 6.66612260e-01, -3.95712710e-04, 1.06071258e+00, 2.04021134e-01,\n 4.98351674e-04, -8.23692783e-05, 2.56390623e-03, 8.33545034e-01,\n 3.34117401e-01, 1.06305001e+00, 2.01871348e-01, 4.93186044e-01,\n 2.89893082e-01, 1.54057892e-04, 8.33371977e-01, 3.33371805e-01,\n 4.99947592e-01, 2.88572124e-01, 1.06069796e+00, 2.04186949e-01,\n 6.67328444e-01, 6.66755900e-01, 9.29823912e-04, -1.07287349e-03,\n -1.25808343e-04, -2.52878707e-04, 8.33270113e-01, 8.33270113e-01,\n 1.49996148e+00, 2.88916374e-01, 1.06068741e+00, 2.04263275e-01]) = Dat(DataSet(Set((np.int64(56), np.int64(56), np.int64(56)), 'set_#x7f097296cb30'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8314').data
E + where Dat(DataSet(Set((np.int64(56), np.int64(56), np.int64(56)), 'set_#x7f097296cb30'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8314') = Coefficient(WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0973118bc0>, FiniteElement('Regge', triangle, 1), name=None), Mesh(VectorElement(FiniteElement('Lagrange', triangle, 1), dim=2), 8270)), 13976).dat
tests/regression/test_interpolate_vs_project.py:58: AssertionError
Check failure on line 58 in tests/regression/test_interpolate_vs_project.py
github-actions / Firedrake real
test_interpolate_vs_project.test_interpolate_vs_project[cube-FunctionSpace(Regge1)]
AssertionError: assert False
+ where False = <function allclose at 0x7f09c590f930>(array([ 3. , 0.83333333, 0.66666667, 3.33333333, 5.5 ,\n 0.16666667, 0.83333333, 5. , 1.66666667, 1.44337567,\n 0.28867513, 3.46410162, 4.24264069, 0.81649658, 1.5 ,\n -0.28867513, 1.06066017, 0.20412415, 7.4246212 , 1.83711731,\n 0.5 , 0.28867513, 4. , 1.15470054, 3. ,\n 0.83333333, 0.66666667, 5. , 0.83333333, 1.66666667,\n 0.28867513, 1.44337567, 3.46410162, 1.06066017, 0.20412415,\n 1.5 , -0.28867513, 4. , -1.15470054, 5.5 ,\n 2.33333333, 0.66666667, 0.83333333, 2. , 0.16666667,\n 3.46410162, 0.57735027, 1.15470054, 2.82842712, 0.81649658,\n 1.06066017, -0.20412415, 0.5 , 0.28867513, 2. ,\n 0.16666667, 0.83333333, 0.28867513, 4.04145188, 6.92820323,\n 3.5 , 0.66666667, 0.83333333, 0.16666667, 0.33333333,\n 1. , 4. , -1.15470054, 0.35355339, -0.20412415,\n 1.5 , 0.28867513, 5.30330086, -1.83711731, 1.41421356,\n -0.81649658, 0.83333333, 2. , 0.16666667, 5.5 ,\n 2.33333333, 0.66666667, 3.46410162, ... 2.82842712, 0.81649658,\n 3. , 0.83333333, 0.66666667, 5. , 0.83333333,\n 1.66666667, 0.28867513, 3.46410162, 1.44337567, 1.06066017,\n 0.20412415, 1.5 , -0.28867513, 4. , -1.15470054,\n 0.66666667, 1.33333333, 4. , 1.66666667, 0.33333333,\n 3.5 , 5.30330086, -1.83711731, 2. , 0.83333333,\n 0.16666667, 1.15470054, 3.46410162, 0.57735027, 0.35355339,\n -0.20412415, 1.66666667, 0.33333333, 3.5 , 2. ,\n 0.83333333, 0.16666667, 1.15470054, 3.46410162, 0.57735027,\n 0.35355339, -0.20412415, 4. , 0.16666667, 2.33333333,\n 0.33333333, 0.66666667, 2. , 4. , 1.15470054,\n 0.33333333, 0.66666667, 2. , 4. , 0.16666667,\n 2.33333333, 4. , 1.15470054, 3.5 , 0.66666667,\n 0.83333333, 3.46410162, 1.44337567, 0.28867513, 0.16666667,\n 1. , 0.33333333, 1.5 , 0.28867513, 1.41421356,\n -0.81649658, 0.5 , -0.28867513, 3.46410162, 1.44337567,\n 0.28867513, 0.16666667, 1. , 0.33333333, 0.5 ,\n -0.28867513]), array([ 3.05085842, 0.86107353, 0.68690482, 3.46782652, 5.77440476,\n 0.18713935, 0.88680504, 5.31936332, 1.85744485, 1.41416709,\n 0.25946654, 3.34726724, 4.48869293, 0.50453799, 1.43858189,\n -0.24152551, 1.17622132, 0.11587108, 8.34545762, 1.4162593 ,\n 0.43928549, 0.336231 , 4.08380429, 1.15632545, 3.05085842,\n 0.86107353, 0.68690482, 5.31936333, 0.88680504, 1.85744486,\n 0.25946653, 1.41416709, 3.34726724, 1.17622132, 0.11587107,\n 1.43858189, -0.24152552, 4.08380429, -1.15632546, 5.31938635,\n 2.43400809, 0.68042695, 0.79155607, 1.89142585, 0.15524404,\n 2.96243402, 0.46879328, 1.03466801, 3.09662316, 0.4958687 ,\n 1.35357074, -0.31535984, 0.36310608, 0.22059957, 1.79674107,\n 0.12474411, 0.88337549, 0.15895804, 2.56929955, 4.49344088,\n 3.86866315, 0.67548334, 0.91686289, 0.17577899, 0.3712861 ,\n 1.13761467, 3.25491238, -1.19931831, 0.28659021, -0.22936352,\n 1.16097673, 0.22058257, 5.82865154, -2.22601608, 1.41067067,\n -0.99230848, 0.79155607, 1.89142587, 0.15524404, 5.31938637,\n 2.4340081 , 0.68042695, 2.96243402, ... 2.88245618, 0.85195414,\n 2.83284715, 0.81325083, 0.62795022, 4.69600433, 0.87761184,\n 1.55074172, 0.24006226, 3.64004362, 1.53530761, 1.00205949,\n 0.14826175, 1.41338403, -0.39455985, 3.72871958, -1.28615582,\n 0.56540966, 1.17752888, 3.99331102, 1.64832237, 0.28062424,\n 3.70645802, 5.16352959, -2.41179053, 1.80330713, 0.75497906,\n 0.13183775, 0.98857867, 3.67454032, 0.61483971, 0.41669874,\n -0.26501953, 1.64832236, 0.28062423, 3.706458 , 1.80330714,\n 0.75497905, 0.13183776, 0.98857867, 3.6745403 , 0.61483971,\n 0.41669876, -0.26501952, 4.17288651, 0.19991104, 2.34125801,\n 0.31652345, 0.62058386, 2.16684412, 2.72026236, 1.09980084,\n 0.31652346, 0.62058387, 2.16684415, 4.17288656, 0.19991104,\n 2.34125804, 2.72026246, 1.09980085, 3.34766794, 0.54552516,\n 0.89160837, 2.71838179, 1.25694571, 0.10224517, 0.17981149,\n 1.13321999, 0.38368428, 1.34894718, 0.37245404, 1.33261983,\n -0.89319085, 0.3251749 , -0.19117131, 2.71838184, 1.25694572,\n 0.10224519, 0.1798115 , 1.13322002, 0.38368428, 0.32517492,\n -0.19117132]), atol=1e-06)
+ where <function allclose at 0x7f09c590f930> = np.allclose
+ and array([ 3. , 0.83333333, 0.66666667, 3.33333333, 5.5 ,\n 0.16666667, 0.83333333, 5. , 1.66666667, 1.44337567,\n 0.28867513, 3.46410162, 4.24264069, 0.81649658, 1.5 ,\n -0.28867513, 1.06066017, 0.20412415, 7.4246212 , 1.83711731,\n 0.5 , 0.28867513, 4. , 1.15470054, 3. ,\n 0.83333333, 0.66666667, 5. , 0.83333333, 1.66666667,\n 0.28867513, 1.44337567, 3.46410162, 1.06066017, 0.20412415,\n 1.5 , -0.28867513, 4. , -1.15470054, 5.5 ,\n 2.33333333, 0.66666667, 0.83333333, 2. , 0.16666667,\n 3.46410162, 0.57735027, 1.15470054, 2.82842712, 0.81649658,\n 1.06066017, -0.20412415, 0.5 , 0.28867513, 2. ,\n 0.16666667, 0.83333333, 0.28867513, 4.04145188, 6.92820323,\n 3.5 , 0.66666667, 0.83333333, 0.16666667, 0.33333333,\n 1. , 4. , -1.15470054, 0.35355339, -0.20412415,\n 1.5 , 0.28867513, 5.30330086, -1.83711731, 1.41421356,\n -0.81649658, 0.83333333, 2. , 0.16666667, 5.5 ,\n 2.33333333, 0.66666667, 3.46410162, ... 2.82842712, 0.81649658,\n 3. , 0.83333333, 0.66666667, 5. , 0.83333333,\n 1.66666667, 0.28867513, 3.46410162, 1.44337567, 1.06066017,\n 0.20412415, 1.5 , -0.28867513, 4. , -1.15470054,\n 0.66666667, 1.33333333, 4. , 1.66666667, 0.33333333,\n 3.5 , 5.30330086, -1.83711731, 2. , 0.83333333,\n 0.16666667, 1.15470054, 3.46410162, 0.57735027, 0.35355339,\n -0.20412415, 1.66666667, 0.33333333, 3.5 , 2. ,\n 0.83333333, 0.16666667, 1.15470054, 3.46410162, 0.57735027,\n 0.35355339, -0.20412415, 4. , 0.16666667, 2.33333333,\n 0.33333333, 0.66666667, 2. , 4. , 1.15470054,\n 0.33333333, 0.66666667, 2. , 4. , 0.16666667,\n 2.33333333, 4. , 1.15470054, 3.5 , 0.66666667,\n 0.83333333, 3.46410162, 1.44337567, 0.28867513, 0.16666667,\n 1. , 0.33333333, 1.5 , 0.28867513, 1.41421356,\n -0.81649658, 0.5 , -0.28867513, 3.46410162, 1.44337567,\n 0.28867513, 0.16666667, 1. , 0.33333333, 0.5 ,\n -0.28867513]) = Dat(DataSet(Set((np.int64(556), np.int64(556), np.int64(556)), 'set_#x7f0972d98ec0'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8359').data
+ where Dat(DataSet(Set((np.int64(556), np.int64(556), np.int64(556)), 'set_#x7f0972d98ec0'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8359') = Coefficient(WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0995314560>, FiniteElement('Regge', tetrahedron, 1), name=None), Mesh(VectorElement(FiniteElement('Lagrange', tetrahedron, 1), dim=3), 8316)), 14107).dat
+ and array([ 3.05085842, 0.86107353, 0.68690482, 3.46782652, 5.77440476,\n 0.18713935, 0.88680504, 5.31936332, 1.85744485, 1.41416709,\n 0.25946654, 3.34726724, 4.48869293, 0.50453799, 1.43858189,\n -0.24152551, 1.17622132, 0.11587108, 8.34545762, 1.4162593 ,\n 0.43928549, 0.336231 , 4.08380429, 1.15632545, 3.05085842,\n 0.86107353, 0.68690482, 5.31936333, 0.88680504, 1.85744486,\n 0.25946653, 1.41416709, 3.34726724, 1.17622132, 0.11587107,\n 1.43858189, -0.24152552, 4.08380429, -1.15632546, 5.31938635,\n 2.43400809, 0.68042695, 0.79155607, 1.89142585, 0.15524404,\n 2.96243402, 0.46879328, 1.03466801, 3.09662316, 0.4958687 ,\n 1.35357074, -0.31535984, 0.36310608, 0.22059957, 1.79674107,\n 0.12474411, 0.88337549, 0.15895804, 2.56929955, 4.49344088,\n 3.86866315, 0.67548334, 0.91686289, 0.17577899, 0.3712861 ,\n 1.13761467, 3.25491238, -1.19931831, 0.28659021, -0.22936352,\n 1.16097673, 0.22058257, 5.82865154, -2.22601608, 1.41067067,\n -0.99230848, 0.79155607, 1.89142587, 0.15524404, 5.31938637,\n 2.4340081 , 0.68042695, 2.96243402, ... 2.88245618, 0.85195414,\n 2.83284715, 0.81325083, 0.62795022, 4.69600433, 0.87761184,\n 1.55074172, 0.24006226, 3.64004362, 1.53530761, 1.00205949,\n 0.14826175, 1.41338403, -0.39455985, 3.72871958, -1.28615582,\n 0.56540966, 1.17752888, 3.99331102, 1.64832237, 0.28062424,\n 3.70645802, 5.16352959, -2.41179053, 1.80330713, 0.75497906,\n 0.13183775, 0.98857867, 3.67454032, 0.61483971, 0.41669874,\n -0.26501953, 1.64832236, 0.28062423, 3.706458 , 1.80330714,\n 0.75497905, 0.13183776, 0.98857867, 3.6745403 , 0.61483971,\n 0.41669876, -0.26501952, 4.17288651, 0.19991104, 2.34125801,\n 0.31652345, 0.62058386, 2.16684412, 2.72026236, 1.09980084,\n 0.31652346, 0.62058387, 2.16684415, 4.17288656, 0.19991104,\n 2.34125804, 2.72026246, 1.09980085, 3.34766794, 0.54552516,\n 0.89160837, 2.71838179, 1.25694571, 0.10224517, 0.17981149,\n 1.13321999, 0.38368428, 1.34894718, 0.37245404, 1.33261983,\n -0.89319085, 0.3251749 , -0.19117131, 2.71838184, 1.25694572,\n 0.10224519, 0.1798115 , 1.13322002, 0.38368428, 0.32517492,\n -0.19117132]) = Dat(DataSet(Set((np.int64(556), np.int64(556), np.int64(556)), 'set_#x7f0972d98ec0'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8360').data
+ where Dat(DataSet(Set((np.int64(556), np.int64(556), np.int64(556)), 'set_#x7f0972d98ec0'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8360') = Coefficient(WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0995314560>, FiniteElement('Regge', tetrahedron, 1), name=None), Mesh(VectorElement(FiniteElement('Lagrange', tetrahedron, 1), dim=3), 8316)), 14109).dat
Raw output
V = WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0995314560>, FiniteElement('Regge', tetrahedron, 1), name=None), Mesh(VectorElement(FiniteElement('Lagrange', tetrahedron, 1), dim=3), 8316))
def test_interpolate_vs_project(V):
mesh = V.mesh()
dim = mesh.geometric_dimension()
if dim == 2:
x, y = SpatialCoordinate(mesh)
elif dim == 3:
x, y, z = SpatialCoordinate(mesh)
shape = V.ufl_element().value_shape
if dim == 2:
if len(shape) == 0:
expression = x + y
elif len(shape) == 1:
expression = as_vector([x, y])
elif len(shape) == 2:
expression = as_tensor(([x, y], [x, y]))
elif dim == 3:
if len(shape) == 0:
expression = x + y + z
elif len(shape) == 1:
expression = as_vector([x, y, z])
elif len(shape) == 2:
expression = as_tensor(([x, y, z], [x, y, z], [x, y, z]))
f = assemble(interpolate(expression, V))
expect = project(expression, V)
> assert np.allclose(f.dat.data, expect.dat.data, atol=1e-06)
E AssertionError: assert False
E + where False = <function allclose at 0x7f09c590f930>(array([ 3. , 0.83333333, 0.66666667, 3.33333333, 5.5 ,\n 0.16666667, 0.83333333, 5. , 1.66666667, 1.44337567,\n 0.28867513, 3.46410162, 4.24264069, 0.81649658, 1.5 ,\n -0.28867513, 1.06066017, 0.20412415, 7.4246212 , 1.83711731,\n 0.5 , 0.28867513, 4. , 1.15470054, 3. ,\n 0.83333333, 0.66666667, 5. , 0.83333333, 1.66666667,\n 0.28867513, 1.44337567, 3.46410162, 1.06066017, 0.20412415,\n 1.5 , -0.28867513, 4. , -1.15470054, 5.5 ,\n 2.33333333, 0.66666667, 0.83333333, 2. , 0.16666667,\n 3.46410162, 0.57735027, 1.15470054, 2.82842712, 0.81649658,\n 1.06066017, -0.20412415, 0.5 , 0.28867513, 2. ,\n 0.16666667, 0.83333333, 0.28867513, 4.04145188, 6.92820323,\n 3.5 , 0.66666667, 0.83333333, 0.16666667, 0.33333333,\n 1. , 4. , -1.15470054, 0.35355339, -0.20412415,\n 1.5 , 0.28867513, 5.30330086, -1.83711731, 1.41421356,\n -0.81649658, 0.83333333, 2. , 0.16666667, 5.5 ,\n 2.33333333, 0.66666667, 3.46410162, ... 2.82842712, 0.81649658,\n 3. , 0.83333333, 0.66666667, 5. , 0.83333333,\n 1.66666667, 0.28867513, 3.46410162, 1.44337567, 1.06066017,\n 0.20412415, 1.5 , -0.28867513, 4. , -1.15470054,\n 0.66666667, 1.33333333, 4. , 1.66666667, 0.33333333,\n 3.5 , 5.30330086, -1.83711731, 2. , 0.83333333,\n 0.16666667, 1.15470054, 3.46410162, 0.57735027, 0.35355339,\n -0.20412415, 1.66666667, 0.33333333, 3.5 , 2. ,\n 0.83333333, 0.16666667, 1.15470054, 3.46410162, 0.57735027,\n 0.35355339, -0.20412415, 4. , 0.16666667, 2.33333333,\n 0.33333333, 0.66666667, 2. , 4. , 1.15470054,\n 0.33333333, 0.66666667, 2. , 4. , 0.16666667,\n 2.33333333, 4. , 1.15470054, 3.5 , 0.66666667,\n 0.83333333, 3.46410162, 1.44337567, 0.28867513, 0.16666667,\n 1. , 0.33333333, 1.5 , 0.28867513, 1.41421356,\n -0.81649658, 0.5 , -0.28867513, 3.46410162, 1.44337567,\n 0.28867513, 0.16666667, 1. , 0.33333333, 0.5 ,\n -0.28867513]), array([ 3.05085842, 0.86107353, 0.68690482, 3.46782652, 5.77440476,\n 0.18713935, 0.88680504, 5.31936332, 1.85744485, 1.41416709,\n 0.25946654, 3.34726724, 4.48869293, 0.50453799, 1.43858189,\n -0.24152551, 1.17622132, 0.11587108, 8.34545762, 1.4162593 ,\n 0.43928549, 0.336231 , 4.08380429, 1.15632545, 3.05085842,\n 0.86107353, 0.68690482, 5.31936333, 0.88680504, 1.85744486,\n 0.25946653, 1.41416709, 3.34726724, 1.17622132, 0.11587107,\n 1.43858189, -0.24152552, 4.08380429, -1.15632546, 5.31938635,\n 2.43400809, 0.68042695, 0.79155607, 1.89142585, 0.15524404,\n 2.96243402, 0.46879328, 1.03466801, 3.09662316, 0.4958687 ,\n 1.35357074, -0.31535984, 0.36310608, 0.22059957, 1.79674107,\n 0.12474411, 0.88337549, 0.15895804, 2.56929955, 4.49344088,\n 3.86866315, 0.67548334, 0.91686289, 0.17577899, 0.3712861 ,\n 1.13761467, 3.25491238, -1.19931831, 0.28659021, -0.22936352,\n 1.16097673, 0.22058257, 5.82865154, -2.22601608, 1.41067067,\n -0.99230848, 0.79155607, 1.89142587, 0.15524404, 5.31938637,\n 2.4340081 , 0.68042695, 2.96243402, ... 2.88245618, 0.85195414,\n 2.83284715, 0.81325083, 0.62795022, 4.69600433, 0.87761184,\n 1.55074172, 0.24006226, 3.64004362, 1.53530761, 1.00205949,\n 0.14826175, 1.41338403, -0.39455985, 3.72871958, -1.28615582,\n 0.56540966, 1.17752888, 3.99331102, 1.64832237, 0.28062424,\n 3.70645802, 5.16352959, -2.41179053, 1.80330713, 0.75497906,\n 0.13183775, 0.98857867, 3.67454032, 0.61483971, 0.41669874,\n -0.26501953, 1.64832236, 0.28062423, 3.706458 , 1.80330714,\n 0.75497905, 0.13183776, 0.98857867, 3.6745403 , 0.61483971,\n 0.41669876, -0.26501952, 4.17288651, 0.19991104, 2.34125801,\n 0.31652345, 0.62058386, 2.16684412, 2.72026236, 1.09980084,\n 0.31652346, 0.62058387, 2.16684415, 4.17288656, 0.19991104,\n 2.34125804, 2.72026246, 1.09980085, 3.34766794, 0.54552516,\n 0.89160837, 2.71838179, 1.25694571, 0.10224517, 0.17981149,\n 1.13321999, 0.38368428, 1.34894718, 0.37245404, 1.33261983,\n -0.89319085, 0.3251749 , -0.19117131, 2.71838184, 1.25694572,\n 0.10224519, 0.1798115 , 1.13322002, 0.38368428, 0.32517492,\n -0.19117132]), atol=1e-06)
E + where <function allclose at 0x7f09c590f930> = np.allclose
E + and array([ 3. , 0.83333333, 0.66666667, 3.33333333, 5.5 ,\n 0.16666667, 0.83333333, 5. , 1.66666667, 1.44337567,\n 0.28867513, 3.46410162, 4.24264069, 0.81649658, 1.5 ,\n -0.28867513, 1.06066017, 0.20412415, 7.4246212 , 1.83711731,\n 0.5 , 0.28867513, 4. , 1.15470054, 3. ,\n 0.83333333, 0.66666667, 5. , 0.83333333, 1.66666667,\n 0.28867513, 1.44337567, 3.46410162, 1.06066017, 0.20412415,\n 1.5 , -0.28867513, 4. , -1.15470054, 5.5 ,\n 2.33333333, 0.66666667, 0.83333333, 2. , 0.16666667,\n 3.46410162, 0.57735027, 1.15470054, 2.82842712, 0.81649658,\n 1.06066017, -0.20412415, 0.5 , 0.28867513, 2. ,\n 0.16666667, 0.83333333, 0.28867513, 4.04145188, 6.92820323,\n 3.5 , 0.66666667, 0.83333333, 0.16666667, 0.33333333,\n 1. , 4. , -1.15470054, 0.35355339, -0.20412415,\n 1.5 , 0.28867513, 5.30330086, -1.83711731, 1.41421356,\n -0.81649658, 0.83333333, 2. , 0.16666667, 5.5 ,\n 2.33333333, 0.66666667, 3.46410162, ... 2.82842712, 0.81649658,\n 3. , 0.83333333, 0.66666667, 5. , 0.83333333,\n 1.66666667, 0.28867513, 3.46410162, 1.44337567, 1.06066017,\n 0.20412415, 1.5 , -0.28867513, 4. , -1.15470054,\n 0.66666667, 1.33333333, 4. , 1.66666667, 0.33333333,\n 3.5 , 5.30330086, -1.83711731, 2. , 0.83333333,\n 0.16666667, 1.15470054, 3.46410162, 0.57735027, 0.35355339,\n -0.20412415, 1.66666667, 0.33333333, 3.5 , 2. ,\n 0.83333333, 0.16666667, 1.15470054, 3.46410162, 0.57735027,\n 0.35355339, -0.20412415, 4. , 0.16666667, 2.33333333,\n 0.33333333, 0.66666667, 2. , 4. , 1.15470054,\n 0.33333333, 0.66666667, 2. , 4. , 0.16666667,\n 2.33333333, 4. , 1.15470054, 3.5 , 0.66666667,\n 0.83333333, 3.46410162, 1.44337567, 0.28867513, 0.16666667,\n 1. , 0.33333333, 1.5 , 0.28867513, 1.41421356,\n -0.81649658, 0.5 , -0.28867513, 3.46410162, 1.44337567,\n 0.28867513, 0.16666667, 1. , 0.33333333, 0.5 ,\n -0.28867513]) = Dat(DataSet(Set((np.int64(556), np.int64(556), np.int64(556)), 'set_#x7f0972d98ec0'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8359').data
E + where Dat(DataSet(Set((np.int64(556), np.int64(556), np.int64(556)), 'set_#x7f0972d98ec0'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8359') = Coefficient(WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0995314560>, FiniteElement('Regge', tetrahedron, 1), name=None), Mesh(VectorElement(FiniteElement('Lagrange', tetrahedron, 1), dim=3), 8316)), 14107).dat
E + and array([ 3.05085842, 0.86107353, 0.68690482, 3.46782652, 5.77440476,\n 0.18713935, 0.88680504, 5.31936332, 1.85744485, 1.41416709,\n 0.25946654, 3.34726724, 4.48869293, 0.50453799, 1.43858189,\n -0.24152551, 1.17622132, 0.11587108, 8.34545762, 1.4162593 ,\n 0.43928549, 0.336231 , 4.08380429, 1.15632545, 3.05085842,\n 0.86107353, 0.68690482, 5.31936333, 0.88680504, 1.85744486,\n 0.25946653, 1.41416709, 3.34726724, 1.17622132, 0.11587107,\n 1.43858189, -0.24152552, 4.08380429, -1.15632546, 5.31938635,\n 2.43400809, 0.68042695, 0.79155607, 1.89142585, 0.15524404,\n 2.96243402, 0.46879328, 1.03466801, 3.09662316, 0.4958687 ,\n 1.35357074, -0.31535984, 0.36310608, 0.22059957, 1.79674107,\n 0.12474411, 0.88337549, 0.15895804, 2.56929955, 4.49344088,\n 3.86866315, 0.67548334, 0.91686289, 0.17577899, 0.3712861 ,\n 1.13761467, 3.25491238, -1.19931831, 0.28659021, -0.22936352,\n 1.16097673, 0.22058257, 5.82865154, -2.22601608, 1.41067067,\n -0.99230848, 0.79155607, 1.89142587, 0.15524404, 5.31938637,\n 2.4340081 , 0.68042695, 2.96243402, ... 2.88245618, 0.85195414,\n 2.83284715, 0.81325083, 0.62795022, 4.69600433, 0.87761184,\n 1.55074172, 0.24006226, 3.64004362, 1.53530761, 1.00205949,\n 0.14826175, 1.41338403, -0.39455985, 3.72871958, -1.28615582,\n 0.56540966, 1.17752888, 3.99331102, 1.64832237, 0.28062424,\n 3.70645802, 5.16352959, -2.41179053, 1.80330713, 0.75497906,\n 0.13183775, 0.98857867, 3.67454032, 0.61483971, 0.41669874,\n -0.26501953, 1.64832236, 0.28062423, 3.706458 , 1.80330714,\n 0.75497905, 0.13183776, 0.98857867, 3.6745403 , 0.61483971,\n 0.41669876, -0.26501952, 4.17288651, 0.19991104, 2.34125801,\n 0.31652345, 0.62058386, 2.16684412, 2.72026236, 1.09980084,\n 0.31652346, 0.62058387, 2.16684415, 4.17288656, 0.19991104,\n 2.34125804, 2.72026246, 1.09980085, 3.34766794, 0.54552516,\n 0.89160837, 2.71838179, 1.25694571, 0.10224517, 0.17981149,\n 1.13321999, 0.38368428, 1.34894718, 0.37245404, 1.33261983,\n -0.89319085, 0.3251749 , -0.19117131, 2.71838184, 1.25694572,\n 0.10224519, 0.1798115 , 1.13322002, 0.38368428, 0.32517492,\n -0.19117132]) = Dat(DataSet(Set((np.int64(556), np.int64(556), np.int64(556)), 'set_#x7f0972d98ec0'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8360').data
E + where Dat(DataSet(Set((np.int64(556), np.int64(556), np.int64(556)), 'set_#x7f0972d98ec0'), (1,), 'None_nodes_dset'), None, dtype('float64'), 'function_8360') = Coefficient(WithGeometry(FunctionSpace(<firedrake.mesh.MeshTopology object at 0x7f0995314560>, FiniteElement('Regge', tetrahedron, 1), name=None), Mesh(VectorElement(FiniteElement('Lagrange', tetrahedron, 1), dim=3), 8316)), 14109).dat
tests/regression/test_interpolate_vs_project.py:58: AssertionError