next fix

backend/services/planning_engine.py

@@ -41,6 +41,9 @@ GE_MIN_EXPORT_W      = 1.0
 # tento počet 15min slotů, držíme plánovací spodek na rezervě (arb_base_wh) místo planner floor —
 # priorita: beze „ztráty na prodeji“ (sell >= 0) držet buffer, hluboký vývoz až těsně před záporným prodejem.
 DEFAULT_PLANNER_DISCHARGE_RELAX_PREWINDOW_SLOTS = 8
+# Měkká kotva: chceme být u planner floor už v posledním slotu před prvním sell < 0.
+# Penalizace je v Kč/Wh (např. 0.02 = 20 Kč/kWh).
+PRENEG_SELL_SOC_ANCHOR_SLACK_PENALTY_CZK_PER_WH = 0.02
 CORRECTION_WINDOW_H  = 1        # hodina zpět pro výpočet korekčního faktoru
 CORRECTION_MIN_CLAMP = 0.5      # spodní limit korekčního faktoru
 CORRECTION_MAX_CLAMP = 1.5      # horní limit korekčního faktoru
@@ -519,6 +522,8 @@ def solve_dispatch(
     ca  = [pulp.LpVariable(f"ca_{t}",  0, slots[t].pv_a_forecast_w)        for t in range(T)]
     hp  = [pulp.LpVariable(f"hp_{t}",  0, heat_pump.rated_heating_power_w) for t in range(T)]
     soc_deficit_24h = pulp.LpVariable("soc_deficit_24h", 0, battery.usable_capacity_wh)
+    soc_anchor_slack = None
+    t_anchor = None
 
     # GEN port cut-off (BA81): binární proměnná pouze pokud je feature povolená v konfiguraci site/invertoru.
     gen_cutoff_enabled = bool(getattr(grid, "deye_gen_microinverter_cutoff_enabled", False))
@@ -547,6 +552,13 @@ def solve_dispatch(
         avg_buy_terminal * TERMINAL_SOC_VALUE_FACTOR / 1000.0
     )
 
+    # Kotva: poslední slot před prvním sell<0 by měl končit u planner floor (pokud relaxace existuje).
+    # Slack penalizujeme v objective; samotné omezení přidáme až po definici soc.
+    first_neg_sell_idx = next((i for i, s in enumerate(slots) if float(s.sell_price) < 0), None)
+    if first_neg_sell_idx is not None and first_neg_sell_idx > 0 and floor_pct is not None:
+        t_anchor = first_neg_sell_idx - 1
+        soc_anchor_slack = pulp.LpVariable("soc_anchor_slack_wh", 0, float(battery.usable_capacity_wh))
+
     # --- Účelová funkce (jen OTE sloty; terminal SoC shadow price na konci horizontu) ---
     prob += (
         pulp.lpSum(
@@ -564,9 +576,18 @@ def solve_dispatch(
         )
         + soc_deficit_24h * soc_deficit_penalty_czk_kwh / 1000
         - terminal_soc_kcz_per_wh * soc[T - 1]
+        + (
+            soc_anchor_slack * PRENEG_SELL_SOC_ANCHOR_SLACK_PENALTY_CZK_PER_WH
+            if soc_anchor_slack is not None
+            else 0
+        )
     )
 
     # --- Omezení ---
+    if t_anchor is not None and soc_anchor_slack is not None:
+        target_floor_wh = float(soc_min_series[t_anchor])
+        prob += soc[t_anchor] <= target_floor_wh + soc_anchor_slack
+
     for t in range(T):
         s = slots[t]
         pv_a_net = s.pv_a_forecast_w - ca[t]