def test_order(problem, stabilization):
    """Assert the correct discretization order.
    """
    mesh_sizes = [16, 32, 64]
    errors, hmax = _compute_errors(problem, mesh_sizes, stabilization)

    # Compute the numerical order of convergence.
    order = helpers.compute_numerical_order_of_convergence(hmax, errors)

    # The test is considered passed if the numerical order of convergence
    # matches the expected order in at least the first step in the coarsest
    # spatial discretization, and is not getting worse as the spatial
    # discretizations are refining.
    tol = 0.1
    expected_order = 2.0
    assert (order > expected_order - tol).all()
    return