if np.product(hfg.shape) == 0:
        hfg = np.ones((1,), int)