diff --git a/mlforecast/compat.py b/mlforecast/compat.py
index 4878b251..15c151b6 100644
--- a/mlforecast/compat.py
+++ b/mlforecast/compat.py
@@ -22,15 +22,6 @@ def __init__(self, *args, **kwargs):  # noqa: ARG002
             raise ImportError("Please install lightgbm to use this model.")
 
 
-try:
-    from window_ops.shift import shift_array
-except ImportError:
-    import numpy as np
-
-    def shift_array(x, offset):  # noqa: ARG002
-        return np.hstack([np.full(offset, np.nan), x[:-offset]])
-
-
 try:
     from xgboost import XGBRegressor
 except ImportError:
@@ -38,3 +29,19 @@ def shift_array(x, offset):  # noqa: ARG002
     class XGBRegressor:
         def __init__(self, *args, **kwargs):  # noqa: ARG002
             raise ImportError("Please install xgboost to use this model.")
+
+
+try:
+    from window_ops.shift import shift_array
+except ImportError:
+    import numpy as np
+    from utilsforecast.compat import njit
+
+    @njit
+    def shift_array(x, offset):
+        if offset >= x.size or offset < 0:
+            return np.full_like(x, np.nan)
+        out = np.empty_like(x)
+        out[:offset] = np.nan
+        out[offset:] = x[:-offset]
+        return out
diff --git a/nbs/compat.ipynb b/nbs/compat.ipynb
index 3f8eb958..6c1b4472 100644
--- a/nbs/compat.ipynb
+++ b/nbs/compat.ipynb
@@ -37,21 +37,28 @@
     "            )\n",
     "\n",
     "try:\n",
-    "    from window_ops.shift import shift_array\n",
-    "except ImportError:\n",
-    "    import numpy as np\n",
-    "\n",
-    "    def shift_array(x, offset):  # noqa: ARG002\n",
-    "        return np.hstack([np.full(offset, np.nan), x[:-offset]])\n",
-    "\n",
-    "try:\n",
     "    from xgboost import XGBRegressor\n",
     "except ImportError:\n",
     "    class XGBRegressor:\n",
     "        def __init__(self, *args, **kwargs):  # noqa: ARG002\n",
     "            raise ImportError(\n",
     "                \"Please install xgboost to use this model.\"\n",
-    "            )"
+    "            )\n",
+    "\n",
+    "try:\n",
+    "    from window_ops.shift import shift_array\n",
+    "except ImportError:\n",
+    "    import numpy as np\n",
+    "    from utilsforecast.compat import njit\n",
+    "\n",
+    "    @njit\n",
+    "    def shift_array(x, offset):\n",
+    "        if offset >= x.size or offset < 0:\n",
+    "            return np.full_like(x, np.nan)\n",
+    "        out = np.empty_like(x)\n",
+    "        out[:offset] = np.nan\n",
+    "        out[offset:] = x[:-offset]\n",
+    "        return out"
    ]
   }
  ],