Currently we use a number of different comparison mechanisms for arrays. We should check if it makes sense to provide one unified approach.

Current mechanisms:

• np.testing.assert_array_less

• np.allclose

• np.isclose

• np.linalg.norm(a - b) < tol

• np.array_equal

• np.abs(1.0 - node.coordinates[0]) < 1e-10

• ...