import math


def padding_util(order_max, acc):
    """Helper function returns number of nodes to pad at each boundary; defined by the width of FD stencils"""
    return int(math.ceil((acc + order_max)/2)+1)