code: def pegelsfcast_C(h, obj, npaths=None, level=None, bootstrap=None):
      forecast = np.full(h, fill_value=np.nan)
      m = obj["m"]
      n = obj["n"]
      states = obj["states"]
      cesforecast(
          states=states,
          n=n,
          m=m,
          season=switch_ces(obj["seasontype"]),
          h=h,
          f=forecast,
          **obj["par"],
      )
      return forecast
  
imports: import math
  from typing import Tuple
  import numpy
  from numba import njit
  from statsmodels.tsa.seasonal import seasonal_decompose
  from .ets import restrict_to_bounds, results
  from .utils import CACHE, NOGIL
label_desc: Use image node_id 13 for calling a global function with example usage: pegelsfcast_C(h, obj, npaths, level, bootstrap) and returns: forecast