Deployed 6a66208 with MkDocs version: 1.6.1

scripts/plots/analyze_comparisons.py

@@ -6,18 +6,17 @@ Rate, Distance Remaining).
 """
 
 import os
-import pandas as pd
-import numpy as np
-import matplotlib.pyplot as plt
 
+import matplotlib.pyplot as plt
+import numpy as np
+import pandas as pd
 from plot_config import (
-    COLORS,
-    apply_style,
-    BEST_PERFORMER_MARKER,
-    BEST_PERFORMER_TEXT,
     BEST_PERFORMER_COLOR,
-    create_common_parser,
+    BEST_PERFORMER_TEXT,
+    COLORS,
     LEGEND_KWARGS,
+    apply_style,
+    create_common_parser,
 )
 
 
@@ -84,7 +83,8 @@ def plot_grouped_bar(
 
     fig, ax = plt.subplots(figsize=figsize)
     bar_width = 0.35
-    x_indices = np.arange(len(morphologies))
+    group_spacing = 1.3
+    x_indices = np.arange(len(morphologies)) * group_spacing
     all_bars = {}
     all_means = []
 
@@ -118,7 +118,7 @@ def plot_grouped_bar(
         )
         all_bars[arch] = (x_pos, means, stds, bars)
 
-    for m_idx, m in enumerate(morphologies):
+    for m_idx, _ in enumerate(morphologies):
         m_means = {arch: all_bars[arch][1][m_idx] for arch in architectures}
         best_arch = (
             max(m_means, key=m_means.get) if higher_is_better else min(m_means, key=m_means.get)
@@ -144,12 +144,13 @@ def plot_grouped_bar(
 
     x_ticks_pos = (
         x_indices
+        + bar_width  # center the label in the 3 bars
         + (bar_width / 2 if len(architectures) % 2 == 0 else 0)
         - (bar_width / 2 if len(architectures) == 2 else 0)
     )
     ax.set_xticks(x_ticks_pos)
     ax.set_xticklabels([f"{m} Arms" for m in morphologies])
-    ax.tick_params(axis="x", pad=25)  # More padding for the squares
+    ax.tick_params(axis="x")  # More padding for the squares
 
     # X-axis at zero
     ax.axhline(0, color="black", linewidth=1.5)
@@ -172,20 +173,7 @@ def plot_grouped_bar(
             plt.FuncFormatter(lambda x, _: f"{x:.2f}" if abs(x) < 10 else f"{x:.0f}")
         )
 
-    _add_square_placeholders(ax, x_ticks_pos, [f"{m} Arms" for m in morphologies])
-
-    # Add custom legend entry for best performer
-    ax.plot(
-        [],
-        [],
-        marker=BEST_PERFORMER_MARKER,
-        color="w",
-        markerfacecolor=BEST_PERFORMER_COLOR,
-        markersize=15,
-        label="Best Performance",
-        ls="",
-    )
-    ax.legend(**LEGEND_KWARGS, ncol=len(architectures) + 1)
+    ax.legend(**LEGEND_KWARGS, ncol=len(architectures))
     ax.set_facecolor("white")
     fig.patch.set_facecolor("white")
 
@@ -305,22 +293,7 @@ def plot_grouped_bar_alt(
             plt.FuncFormatter(lambda x, _: f"{x:.2f}" if abs(x) < 10 else f"{x:.0f}")
         )
 
-    # In this alt plot, placeholders might be per architecture
-    _add_square_placeholders(
-        ax, x_indices, [arch.replace("_", "\n").title() for arch in architectures]
-    )
-
-    ax.plot(
-        [],
-        [],
-        marker=BEST_PERFORMER_MARKER,
-        color="w",
-        markerfacecolor=BEST_PERFORMER_COLOR,
-        markersize=15,
-        label="Best Performance",
-        ls="",
-    )
-    ax.legend(**LEGEND_KWARGS, ncol=len(morphologies) + 1)
+    ax.legend(**LEGEND_KWARGS, ncol=len(morphologies))
     ax.set_facecolor("white")
     fig.patch.set_facecolor("white")
 
@@ -358,8 +331,8 @@ if __name__ == "__main__":
         plot_grouped_bar(
             df=df,
             metric_col="approx_max_velocity",
-            ylabel="Max Forward Velocity (cm/s)",
-            title="Graceful Degradation: Velocity Across Morphologies",
+            ylabel="",
+            title="Maximal forward velocity (in cm/s)",
             output_filename="poster_plot_velocity.png",
             output_dir=OUTPUT_DIR,
             higher_is_better=True,

scripts/plots/analyze_convergence.py

@@ -44,20 +44,24 @@ class Columns(str, Enum):
     # ... (rest of the file remains same, just need to update plotting functions and obtain_data)
     """Column names expected in every evaluation CSV."""
 
+    CHECKPOINT = "checkpoint"
     ARCH = "architecture"
-    TIMESTEPS = "total_trained_timesteps"
-    REWARD = "accumulated_reward"
+    TIMESTEPS = "trained_timesteps"
+    REWARD = "eval_return"
     VELOCITY = "velocity"
+    EVAL_STEPS = "eval_steps"
+    FINAL_XY_DIST = "final_xy_dist"
+    INITIAL_XY_DIST = "initial_xy_dist"
+    REACHED_TARGET = "reached_target"
 
 
 # Maps architecture display names to the path of their evaluation CSV.
 # Update these paths once real evaluation data is available.
 FILE_MAPPING: dict[str, str] = {
-    "centralized 2 arms": "runs/dummy/dummy_centralized_2_arms.csv",
-    "centralized 5 arms": "runs/dummy/dummy_centralized_5_arms.csv",
-    "decentralized fully connected": "runs/dummy/dummy_decentralized_fully_connected.csv",
-    "decentralized ring-level": "runs/dummy/dummy_decentralized_ring-level.csv",
-    "decentralized segment-level": "runs/dummy/dummy_decentralized_segment-level.csv",
+    # "centralized 2 arms": "runs/dummy/dummy_centralized_2_arms.csv",
+    "centralized 5 arms": "runs/final-v2-centralized/checkpoint_evaluation.csv",
+    "decentralized fully connected": "runs/final-v2-fully-conn/checkpoint_evaluation.csv",
+    "decentralized ring-level": "runs/final-v2-ring/checkpoint_evaluation.csv",
 }
 
 # Architecture profiles for dummy data generation: (max_reward, max_velocity, sigmoid_speed)
@@ -108,7 +112,14 @@ def load_metrics(file_mapping: dict[str, str]) -> pd.DataFrame:
     Loads one CSV per architecture, injects the architecture name as a column,
     and returns the combined DataFrame with only the required columns.
     """
-    required = [Columns.TIMESTEPS, Columns.REWARD, Columns.VELOCITY]
+    required = [
+        Columns.CHECKPOINT,
+        Columns.TIMESTEPS,
+        Columns.REWARD,
+        Columns.INITIAL_XY_DIST,
+        Columns.FINAL_XY_DIST,
+        Columns.EVAL_STEPS,
+    ]
     dfs = []
 
     for arch_name, filepath in file_mapping.items():
@@ -124,17 +135,30 @@ def load_metrics(file_mapping: dict[str, str]) -> pd.DataFrame:
             continue
 
         df = df[required].copy()
+        df[Columns.VELOCITY] = (df[Columns.INITIAL_XY_DIST] - df[Columns.FINAL_XY_DIST]) / df[
+            Columns.EVAL_STEPS
+        ]
         df[Columns.ARCH] = arch_name
+        df[Columns.VELOCITY] = (df[Columns.INITIAL_XY_DIST] - df[Columns.FINAL_XY_DIST]) / df[
+            Columns.EVAL_STEPS
+        ]
+
         dfs.append(df)
 
     return pd.concat(dfs, ignore_index=True) if dfs else pd.DataFrame()
 
 
-def _convergence_timestep(series: pd.Series, timesteps: pd.Series) -> float:
+def _convergence_timestep(
+    series: pd.Series, timesteps: pd.Series, checkpoints: pd.Series
+) -> tuple[float, int, int]:
     """Returns the first timestep where the smoothed series reaches 95% of its peak."""
     smoothed = series.rolling(window=SMOOTHING_WINDOW, min_periods=1).mean()
     threshold = smoothed.max() * CONVERGENCE_THRESHOLD
-    return timesteps[smoothed >= threshold].iloc[0]
+
+    mask = smoothed >= threshold
+    first_idx = mask.idxmax()
+
+    return timesteps.loc[first_idx], first_idx, checkpoints.loc[first_idx]
 
 
 def analyze_convergence(df: pd.DataFrame) -> pd.DataFrame:
@@ -144,21 +168,40 @@ def analyze_convergence(df: pd.DataFrame) -> pd.DataFrame:
     """
     results = []
 
+    centralized_base = 0
+
     for arch in df[Columns.ARCH].unique():
         arch_data = df[df[Columns.ARCH] == arch].sort_values(Columns.TIMESTEPS)
 
+        reward_timestep, reward_checkpoint_idx, reward_checkpoint = _convergence_timestep(
+            arch_data[Columns.REWARD],
+            arch_data[Columns.TIMESTEPS],
+            arch_data[Columns.CHECKPOINT],
+        )
+
+        velocity_timestep, velocity_checkpoint_idx, velocity_checkpoint = _convergence_timestep(
+            arch_data[Columns.VELOCITY],
+            arch_data[Columns.TIMESTEPS],
+            arch_data[Columns.CHECKPOINT],
+        )
+
         results.append(
             {
                 "Architecture": arch,
-                "Reward_Convergence_Timestep": _convergence_timestep(
-                    arch_data[Columns.REWARD], arch_data[Columns.TIMESTEPS]
-                ),
-                "Velocity_Convergence_Timestep": _convergence_timestep(
-                    arch_data[Columns.VELOCITY], arch_data[Columns.TIMESTEPS]
-                ),
+                "Reward_Convergence_Timestep": reward_timestep,
+                "Reward_Convergence_Checkpoint_Idx": reward_checkpoint_idx,
+                "Reward_Convergence_Checkpoint": reward_checkpoint,
+                "Velocity_Convergence_Timestep": velocity_timestep,
+                "Velocity_Convergence_Checkpoint_Idx": velocity_checkpoint_idx,
+                "Velocity_Convergence_Checkpoint": velocity_checkpoint,
             }
         )
 
+        if arch == "centralized 5 arms":
+            centralized_base = reward_checkpoint
+        else:
+            print(arch, "speedup:", 1 - reward_checkpoint / centralized_base)
+
     return pd.DataFrame(results)
 
 
@@ -303,6 +346,7 @@ def plot_results(df: pd.DataFrame, results: pd.DataFrame, output_dir: str, **kwa
 def obtain_data() -> pd.DataFrame:
     """Resolves the file mapping, falling back to generated dummy CSVs if needed."""
     global USING_DUMMY_DATA
+
     if not any(os.path.exists(p) for p in FILE_MAPPING.values()):
         logger.info("No real evaluation files found. Generating dummy CSVs at expected locations.")
         generate_dummy_csvs(FILE_MAPPING)
@@ -319,6 +363,12 @@ def run_analysis(output_dir: str, **kwargs):
         return
 
     results = analyze_convergence(df)
+    print(
+        results[
+            ["Architecture", "Reward_Convergence_Checkpoint_Idx", "Reward_Convergence_Checkpoint"]
+        ]
+    )
+
     plot_results(df, results, output_dir, **kwargs)
     logger.info("Analysis complete. Plots saved to disk.")
 

scripts/plots/plot_config.py

@@ -4,26 +4,24 @@ import matplotlib.pyplot as plt
 # Shared Color Palette (Colorblind friendly, high contrast)
 # Matches poster design
 COLORS = {
-    "CENTRALIZED": "#2B4162",  # Deep Slate Blue
-    "FULLY_CONNECTED": "#FA9F42",  # Vibrant Orange
-    "RING_LEVEL": "#4E937A",  # Muted Teal
-    "SEGMENT_LEVEL": "#B4436C",  # Soft Red
-    "DECENTRALIZED": "#4E937A",  # Default decentralized fallback
+    "CENTRALIZED": "#0D567C",  # Blue
+    "FULLY_CONNECTED": "#8C0E0F",  # Reddish
+    "RING": "#FCB305",  # Pale Yellow
 }
 
 
-def apply_style(font_size=28):
+def apply_style(font_size=36):
     """
     Applies the shared typography and aesthetic settings to Matplotlib.
     """
     plt.rcParams.update(
         {
             "font.size": font_size,
-            "axes.labelsize": font_size + 4,
-            "axes.titlesize": font_size + 8,
-            "xtick.labelsize": font_size - 4,
-            "ytick.labelsize": font_size - 4,
-            "legend.fontsize": font_size - 6,
+            "axes.labelsize": font_size,
+            "axes.titlesize": font_size,
+            "xtick.labelsize": font_size,
+            "ytick.labelsize": font_size,
+            "legend.fontsize": font_size,
             "axes.linewidth": 2,
             "axes.spines.top": False,
             "axes.spines.right": False,
@@ -44,7 +42,7 @@ BEST_PERFORMER_COLOR = "#D4AF37"  # Gold
 # Centralized Legend Configuration
 LEGEND_KWARGS = {
     "loc": "upper center",
-    "bbox_to_anchor": (0.5, -0.5),
+    "bbox_to_anchor": (0.5, -0.12),
     "frameon": False,
 }
 

scripts/poster_visualisations/render_poster_videos.py

@@ -0,0 +1,150 @@
+from __future__ import annotations
+
+import argparse
+from pathlib import Path
+
+from brittle_star_project.evaluation.checkpoint import load_metadata, metadata_to_configs
+from brittle_star_project.evaluation.eval_env_builder import build_eval_env
+from brittle_star_project.evaluation.video import record_episode_multi_camera
+
+_ARCH_DIR_MAP = {
+    "CENTRALIZED": "centralized",
+    "FULLY_CONNECTED": "fully-connected",
+    "RING": "ring",
+    "SEGMENT": "segment",
+}
+
+ROBOT_COLOR_MAP = {
+    "CENTRALIZED": "#0D567C",  # Blue
+    "FULLY_CONNECTED": "#8C0E0F",  # Reddish
+    "RING": "#FCB304",  # Pale Yellow
+}
+
+
+def _arch_dir(name: str) -> str:
+    return _ARCH_DIR_MAP.get(name, name.lower())
+
+
+def _resolve_overrides(overrides: list[str], count: int) -> list[str | None]:
+    if not overrides:
+        return [None] * count
+    if len(overrides) == 1 and count > 1:
+        return overrides * count
+    if len(overrides) != count:
+        raise ValueError("morphology overrides must match the number of models")
+    return overrides
+
+
+def main() -> None:
+    parser = argparse.ArgumentParser(description="Render top-down and follow videos for poster.")
+    parser.add_argument("models", nargs="+", help="Paths to .flax checkpoints")
+    parser.add_argument(
+        "--morphology-override",
+        action="append",
+        default=[],
+        help="Override morphology YAML path (repeat to match models)",
+    )
+    parser.add_argument("--output-root", default="vids/poster")
+    parser.add_argument("--seed", type=int, default=0)
+    parser.add_argument("--max-steps", type=int, default=5000)
+    parser.add_argument("--topdown-camera", type=int, default=0)
+    parser.add_argument("--follow-camera", type=int, default=1)
+    parser.add_argument("--topdown-camera-x", type=float, default=-3.0)
+    parser.add_argument("--topdown-camera-y", type=float, default=0.0)
+    parser.add_argument("--topdown-camera-z", type=float, default=4.5)
+    parser.add_argument("--topdown-camera-fovy", type=float, default=None)
+    parser.add_argument("--target-x", type=float, default=-6.0)
+    parser.add_argument("--target-y", type=float, default=0.0)
+    parser.add_argument("--width", type=int, default=2160)
+    parser.add_argument("--height", type=int, default=960)
+    parser.add_argument("--fps", type=int, default=60)
+    parser.add_argument(
+        "--robot-color",
+        default="#2B4162",
+        help="Hex color for the brittle star robot",
+    )
+    args = parser.parse_args()
+
+    if (args.target_x is None) != (args.target_y is None):
+        raise ValueError("target-x and target-y must be provided together")
+
+    target_xy = None
+    if args.target_x is not None:
+        target_xy = (float(args.target_x), float(args.target_y))
+
+    camera_fovy = None
+    if args.topdown_camera_fovy is not None:
+        camera_fovy = {args.topdown_camera: float(args.topdown_camera_fovy)}
+
+    camera_x = None
+    if args.topdown_camera_x is not None:
+        camera_x = {args.topdown_camera: float(args.topdown_camera_x)}
+
+    camera_y = None
+    if args.topdown_camera_y is not None:
+        camera_y = {args.topdown_camera: float(args.topdown_camera_y)}
+
+    camera_z = None
+    if args.topdown_camera_z is not None:
+        camera_z = {args.topdown_camera: float(args.topdown_camera_z)}
+
+    camera_xyz = (camera_x, camera_y, camera_z)
+
+    overrides = _resolve_overrides(args.morphology_override, len(args.models))
+    output_root = Path(args.output_root)
+
+    for model_path_str, override in zip(args.models, overrides):
+        model_path = Path(model_path_str)
+        metadata = load_metadata(model_path, None)
+        training = metadata_to_configs(metadata)
+
+        bundle = build_eval_env(
+            model_path=model_path,
+            training=training,
+            metadata=metadata,
+            morphology_override_path=override,
+        )
+
+        arch_dir = _arch_dir(bundle.architecture)
+        arms_dir = f"{bundle.num_active_arms}arms"
+        out_dir = output_root / arms_dir / arch_dir
+        out_dir.mkdir(parents=True, exist_ok=True)
+
+        output_paths = {
+            args.topdown_camera: out_dir / "topdown.mp4",
+            args.follow_camera: out_dir / "follow.mp4",
+        }
+
+        print(bundle.architecture)
+        color = ROBOT_COLOR_MAP.get(bundle.architecture, args.robot_color)
+
+        result = record_episode_multi_camera(
+            env=bundle.env,
+            policy=bundle.policy,
+            seed=args.seed,
+            max_steps=args.max_steps,
+            action_low=bundle.action_low,
+            action_high=bundle.action_high,
+            action_mask=bundle.action_mask,
+            output_paths=output_paths,
+            camera_ids=[args.topdown_camera, args.follow_camera],
+            camera_fovy=camera_fovy,
+            camera_xyz=camera_xyz,
+            target_xy=target_xy,
+            robot_color=color,
+            width=args.width,
+            height=args.height,
+            fps=args.fps,
+        )
+
+        final_dist = "n/a" if result.final_xy_dist is None else f"{result.final_xy_dist:.3f}"
+        print(
+            f"{arms_dir}/{arch_dir}: return={result.return_:.6f}, len={result.length}, "
+            f"target_reached={result.reached_target}, final_xy_dist={final_dist}"
+        )
+
+        bundle.env.close()
+
+
+if __name__ == "__main__":
+    main()

scripts/poster_visualisations/render_poster_videos.sh

@@ -0,0 +1,10 @@
+#!/usr/bin/env bash
+
+# Multi-camera renders per model -> vids/poster/{arms}arms/{arch}/topdown.mp4 + follow.mp4
+
+path=$1
+
+uv run scripts/poster_visualisations/render_poster_videos.py \
+  "$path"/centralized.flax \
+  "$path"/fully-connected.flax \
+  "$path"/ring.flax \

scripts/poster_visualisations/render_static_path_image.py

@@ -0,0 +1,241 @@
+from __future__ import annotations
+
+import argparse
+from pathlib import Path
+
+import numpy as np
+
+from brittle_star_project.evaluation.checkpoint import load_metadata, metadata_to_configs
+from brittle_star_project.evaluation.eval_env_builder import build_eval_env
+from brittle_star_project.evaluation.rollout import (
+    _get_observations,
+    _maybe_clip_action,
+    _target_reached,
+)
+from brittle_star_project.evaluation.video import (
+    _apply_camera_overrides,
+    _ensure_offscreen_size,
+    hex_to_rgba,
+)
+
+ROBOT_COLOR_MAP = {
+    "CENTRALIZED": "#0D567C",  # Blue
+    "FULLY_CONNECTED": "#8C0E0F",  # Reddish
+    "RING": "#FCB304",  # Pale Yellow
+}
+
+
+def _enum_value(enum_obj, *names: str) -> int:
+    for name in names:
+        if hasattr(enum_obj, name):
+            return int(getattr(enum_obj, name))
+    raise AttributeError(f"Could not find any of {names!r} on {enum_obj!r}")
+
+
+def _append_sphere(scene, mujoco, center: np.ndarray, radius: float, rgba: np.ndarray) -> None:
+    geom = scene.geoms[scene.ngeom]
+    mujoco.mjv_initGeom(
+        geom,
+        mujoco.mjtGeom.mjGEOM_SPHERE,
+        np.asarray([radius, 0.0, 0.0], dtype=np.float32),
+        center,
+        np.eye(3, dtype=np.float32).reshape(-1),
+        rgba,
+    )
+    scene.ngeom += 1
+
+
+def main() -> None:
+    parser = argparse.ArgumentParser(description="Render a static path image from a rollout.")
+    parser.add_argument("model", help="Path to .flax checkpoint")
+    parser.add_argument("--morphology-override", default=None)
+    parser.add_argument("--output-path", required=True)
+    parser.add_argument("--seed", type=int, default=0)
+    parser.add_argument("--max-steps", type=int, default=5000)
+    parser.add_argument("--body-name", default="BrittleStarMorphology/central_disk")
+    parser.add_argument("--camera-id", type=int, default=0)
+    parser.add_argument("--camera-x", type=float, default=-3.0)
+    parser.add_argument("--camera-y", type=float, default=0.0)
+    parser.add_argument("--camera-z", type=float, default=4.5)
+    parser.add_argument("--camera-fovy", type=float, default=None)
+    parser.add_argument("--target-x", type=float, default=-6.0)
+    parser.add_argument("--target-y", type=float, default=0.0)
+    parser.add_argument("--width", type=int, default=2160)
+    parser.add_argument("--height", type=int, default=960)
+    parser.add_argument("--frame-stride", type=int, default=15)
+    parser.add_argument(
+        "--path-color", default="#FA9F42"
+    )  # ring = #888888, centralized = #2B4162, fully connected = FA9F42
+    parser.add_argument(
+        "--robot-color",
+        default="#FA9F42",
+        help="Hex color for brittle star robot (e.g. #ff0000)",
+    )
+    args = parser.parse_args()
+
+    model_path = Path(args.model)
+    metadata = load_metadata(model_path, None)
+    training = metadata_to_configs(metadata)
+
+    bundle = build_eval_env(
+        model_path=model_path,
+        training=training,
+        metadata=metadata,
+        morphology_override_path=args.morphology_override,
+    )
+
+    if (args.target_x is None) != (args.target_y is None):
+        raise ValueError("target-x and target-y must be provided together")
+
+    target_xy = None
+    if args.target_x is not None:
+        target_xy = (float(args.target_x), float(args.target_y))
+
+    try:
+        import imageio
+        import mujoco
+    except ImportError as e:
+        raise ImportError(
+            "Static image rendering requires 'mujoco' and 'imageio'. "
+            "Please install the evaluation dependencies: `uv pip install .[evaluation]`"
+        ) from e
+
+    reset_kwargs = {}
+    if target_xy is not None:
+        reset_kwargs["target_position"] = (target_xy[0], target_xy[1], 0.0)
+
+    state = bundle.env.reset(seed=args.seed, **reset_kwargs)
+    model = state.mj_model
+    data = state.mj_data
+
+    _apply_camera_overrides(
+        model,
+        camera_fovy={args.camera_id: float(args.camera_fovy)}
+        if args.camera_fovy is not None
+        else None,
+        camera_xyz=(
+            {args.camera_id: float(args.camera_x)} if args.camera_x is not None else None,
+            {args.camera_id: float(args.camera_y)} if args.camera_y is not None else None,
+            {args.camera_id: float(args.camera_z)} if args.camera_z is not None else None,
+        ),
+    )
+    _ensure_offscreen_size(model, args.width, args.height)
+
+    body_id = mujoco.mj_name2id(
+        model, mujoco.mjtObj.mjOBJ_BODY, "BrittleStarMorphology/central_disk"
+    )
+
+    robot_rgba = hex_to_rgba(ROBOT_COLOR_MAP.get(bundle.architecture, args.robot_color), 1.0)
+
+    # Optionally override robot color by recoloring geoms belonging to the robot's body subtree.
+    if args.robot_color is not None:
+        # Collect body IDs in the subtree rooted at `body_id` by walking parent links.
+        nbody = int(model.nbody)
+        body_parent = model.body_parentid
+        robot_body_ids = set([int(body_id)])
+        for i in range(1, nbody):
+            cur = int(i)
+            # walk up until root (0) or until we hit the robot root
+            while cur not in (-1, 0, int(body_id)):
+                cur = int(body_parent[cur])
+            if cur == int(body_id):
+                robot_body_ids.add(i)
+
+        # Recolor geoms whose body id is in the robot subtree
+        for g in range(int(model.ngeom)):
+            if int(model.geom_bodyid[g]) in robot_body_ids:
+                model.geom_rgba[g][:] = robot_rgba
+
+    # Optionally override robot color by recoloring geoms belonging to the robot's body subtree.
+    if args.robot_color is not None:
+        # Collect body IDs in the subtree rooted at `body_id` by walking parent links.
+        nbody = int(model.nbody)
+        body_parent = model.body_parentid
+        robot_body_ids = set([int(body_id)])
+        for i in range(1, nbody):
+            cur = int(i)
+            # walk up until root (0) or until we hit the robot root
+            while cur not in (-1, 0, int(body_id)):
+                cur = int(body_parent[cur])
+            if cur == int(body_id):
+                robot_body_ids.add(i)
+
+        # Recolor geoms whose body id is in the robot subtree
+        for g in range(int(model.ngeom)):
+            if int(model.geom_bodyid[g]) in robot_body_ids:
+                model.geom_rgba[g][:] = robot_rgba
+
+    positions = []
+    observations = _get_observations(state)
+
+    for _ in range(int(args.max_steps)):
+        positions.append(np.asarray(data.xpos[body_id], dtype=np.float32))
+
+        obs_dict = observations or {}
+        action = bundle.policy.act(observations=obs_dict)
+        if bundle.action_mask is not None:
+            action = action[bundle.action_mask]
+        action = _maybe_clip_action(action, bundle.action_low, bundle.action_high)
+
+        state = bundle.env.step(state=state, action=action)
+        data = state.mj_data
+        observations = _get_observations(state)
+
+        if _target_reached(state=state):
+            break
+
+    positions_arr = np.vstack(positions)
+    if len(positions_arr) < 2:
+        raise ValueError("Need at least two rollout positions to render a path")
+
+    path_points = positions_arr.copy()
+    path_points[:, 2] -= 0.02
+
+    path_step = max(1, int(args.frame_stride))
+    path_points_visible = path_points[::path_step]
+    path_rgba = hex_to_rgba(ROBOT_COLOR_MAP.get(bundle.architecture, args.path_color), 0.92)
+
+    ctx = mujoco.GLContext(args.width, args.height)
+    ctx.make_current()
+    try:
+        catmask = _enum_value(mujoco.mjtCatBit, "mjCAT_ALL")
+        camera_type = _enum_value(mujoco.mjtCamera, "mjCAMERA_FIXED")
+        font_scale = _enum_value(mujoco.mjtFontScale, "mjFONTSCALE_100")
+
+        maxgeom = int(model.ngeom + len(path_points_visible) + 8)
+        scene = mujoco.MjvScene(model, maxgeom=maxgeom)
+        option = mujoco.MjvOption()
+        perturb = mujoco.MjvPerturb()
+        camera = mujoco.MjvCamera()
+        mujoco.mjv_defaultOption(option)
+        mujoco.mjv_defaultPerturb(perturb)
+        mujoco.mjv_defaultCamera(camera)
+        camera.type = camera_type
+        camera.fixedcamid = int(args.camera_id)
+        if hasattr(camera, "trackbodyid"):
+            camera.trackbodyid = -1
+
+        context = mujoco.MjrContext(model, font_scale)
+        viewport = mujoco.MjrRect(0, 0, args.width, args.height)
+
+        mujoco.mjv_updateScene(model, data, option, perturb, camera, catmask, scene)
+
+        for idx, path_point in enumerate(path_points_visible):
+            path_rgba[3] = 0.10 + 0.70 * (idx / max(len(path_points_visible) - 1, 1))
+            _append_sphere(scene, mujoco, path_point, 0.03, path_rgba)
+
+        rgb = np.empty((args.height, args.width, 3), dtype=np.uint8)
+        depth = np.empty((args.height, args.width), dtype=np.float32)
+        mujoco.mjr_render(viewport, scene, context)
+        mujoco.mjr_readPixels(rgb, depth, viewport, context)
+        imageio.imwrite(args.output_path, np.flipud(rgb))
+
+        context.free()
+    finally:
+        ctx.free()
+
+    bundle.env.close()
+
+
+if __name__ == "__main__":
+    main()

scripts/poster_visualisations/render_static_path_image.sh

@@ -0,0 +1,10 @@
+#!/usr/bin/env bash
+
+# Static path image + optional ghost render
+
+path=$1 # path to .flax model with metadata.yaml alongside it
+
+uv run scripts/poster_visualisations/render_static_path_image.py \
+  "$path" \
+  --output-path vids/poster/5arms/centralized/path.png \
+  --ghost-overlay

scripts/simulate.py

@@ -107,6 +107,9 @@ def main(dict_cfg: DictConfig) -> None:
             action_mask=action_mask,
             output_path=output_path,
             camera_id=sim_cfg.camera_id,
+            width=sim_cfg.video_width,
+            height=sim_cfg.video_height,
+            fps=sim_cfg.video_fps,
         )
 
         save_evaluation_metadata(

scripts/tools/download_wandb_project.py

@@ -0,0 +1,142 @@
+from pathlib import Path
+from concurrent.futures import ThreadPoolExecutor, as_completed
+import threading
+import wandb
+import argparse
+
+# tune these depending on network / W&B limits
+MAX_RUN_WORKERS = 8
+MAX_FILE_WORKERS = 16
+MAX_ARTIFACT_WORKERS = 8
+
+api = wandb.Api()
+
+print_lock = threading.Lock()
+
+
+def safe_print(*args, **kwargs):
+    with print_lock:
+        print(*args, **kwargs)
+
+
+def download_file(file, run_dir):
+    target = run_dir / file.name
+
+    try:
+        # skip existing files
+        if target.exists():
+            return f"SKIP FILE {target}"
+
+        target.parent.mkdir(parents=True, exist_ok=True)
+
+        file.download(root=run_dir, replace=False)
+
+        return f"DONE FILE {target}"
+
+    except Exception as e:
+        return f"FAIL FILE {target}: {e}"
+
+
+def sanitize_artifact_name(name: str):
+    return name.replace(":", "_")
+
+
+def download_artifact(artifact, artifact_root):
+    try:
+        artifact_name = sanitize_artifact_name(artifact.name)
+        artifact_dir = artifact_root / artifact_name
+
+        if artifact_dir.exists() and any(artifact_dir.iterdir()):
+            return f"SKIP ARTIFACT {artifact.name}"
+
+        artifact_dir.mkdir(parents=True, exist_ok=True)
+
+        artifact.download(root=artifact_dir)
+
+        return f"DONE ARTIFACT {artifact.name}"
+
+    except Exception as e:
+        return f"FAIL ARTIFACT {artifact.name}: {e}"
+
+
+def download_run(run, root):
+    run_dir = root / f"{run.name}"
+    run_dir.mkdir(parents=True, exist_ok=True)
+
+    safe_print(f"\n=== {run.name} ({run.id}) ===")
+
+    # -------------------------
+    # Download regular run files
+    # -------------------------
+    files = list(run.files())
+
+    with ThreadPoolExecutor(max_workers=MAX_FILE_WORKERS) as executor:
+        futures = [executor.submit(download_file, file, run_dir) for file in files]
+
+        for future in as_completed(futures):
+            safe_print(future.result())
+
+    # -------------------------
+    # Download logged artifacts
+    # -------------------------
+    artifact_root = run_dir / "artifacts"
+
+    try:
+        artifacts = list(run.logged_artifacts())
+        safe_print(f"Found {len(artifacts)} artifacts for {run.name}")
+
+        with ThreadPoolExecutor(max_workers=MAX_ARTIFACT_WORKERS) as executor:
+            futures = [
+                executor.submit(download_artifact, artifact, artifact_root)
+                for artifact in artifacts
+            ]
+
+            for future in as_completed(futures):
+                safe_print(future.result())
+
+    except Exception as e:
+        safe_print(f"Artifact download failed for {run.name}: {e}")
+
+    # -------------------------
+    # OPTIONAL: download used/input artifacts
+    # -------------------------
+    # try:
+    #     used_artifacts = list(run.used_artifacts())
+    #     used_root = run_dir / "used_artifacts"
+    #
+    #     for artifact in used_artifacts:
+    #         download_artifact(artifact, used_root)
+    # except Exception as e:
+    #     safe_print(f"Used artifact download failed: {e}")
+
+    safe_print(f"Finished {run.name}")
+
+
+def main(entity: str, project: str, root: Path):
+    root.mkdir(exist_ok=True)
+
+    runs = list(api.runs(f"{entity}/{project}"))
+
+    safe_print(f"Found {len(runs)} runs")
+
+    with ThreadPoolExecutor(max_workers=MAX_RUN_WORKERS) as executor:
+        futures = [executor.submit(download_run, run, root) for run in runs]
+
+        for future in as_completed(futures):
+            try:
+                future.result()
+            except Exception as e:
+                safe_print("RUN FAILED:", e)
+
+    safe_print("\nAll downloads complete.")
+
+
+if __name__ == "__main__":
+    parser = argparse.ArgumentParser()
+    parser.add_argument("--entity", type=str, default="SEL3-2026-Groep-4")
+    parser.add_argument("--project", type=str, required=True)
+    parser.add_argument("--root", type=str, default="runs")
+    args = parser.parse_args()
+
+    root = Path(args.root)
+    main(entity=args.entity, project=args.project, root=root)