simonbb 0.0.4__tar.gz → 0.0.6__tar.gz

{simonbb-0.0.4 → simonbb-0.0.6}/PKG-INFO

@@ -1,7 +1,7 @@
 Metadata-Version: 2.4
 Name: simonbb
-Version: 0.0.4
-Summary: A small package to facilitate teaching AI courses
+Version: 0.0.6
+Summary: A small package to facilitate teaching AI courses at business schools
 Project-URL: Homepage, https://simonbb.com
 Author-email: Huaxia Rui <huaxia.rui@simon.rochester.edu>
 License-Expression: MIT
@@ -20,7 +20,9 @@ Requires-Dist: matplotlib~=3.10.8
 Requires-Dist: numpy~=2.4.0
 Requires-Dist: pytorch-ignite~=0.5.4
 Requires-Dist: pyyaml~=6.0.3
+Requires-Dist: ray~=2.56.0
 Requires-Dist: sentence-transformers~=5.3.0
+Requires-Dist: stable-baselines3~=2.9.0
 Requires-Dist: torchvision~=0.24.1
 Requires-Dist: torch~=2.9.1
 Requires-Dist: transformers~=4.57.3

{simonbb-0.0.4 → simonbb-0.0.6}/pyproject.toml

@@ -4,11 +4,11 @@ build-backend = "hatchling.build"
 
 [project]
 name = "simonbb"
-version = "0.0.4"
+version = "0.0.6"
 authors = [
 	{ name="Huaxia Rui", email="huaxia.rui@simon.rochester.edu" },
 ]
-description = "A small package to facilitate teaching AI courses"
+description = "A small package to facilitate teaching AI courses at business schools"
 readme = "README.md"
 requires-python = ">=3.9"
 classifiers = [
@@ -32,6 +32,8 @@ dependencies = [
 	"sentence-transformers~=5.3.0",
 	"gymnasium~=1.3.0",
 	"pytorch-ignite~=0.5.4",
+	"stable_baselines3~=2.9.0",
+	"ray~=2.56.0",
 	"pyyaml~=6.0.3"
 ]
 [project.urls]

{simonbb-0.0.4 → simonbb-0.0.6}/src/simonbb/cis433/actions.py

@@ -61,10 +61,7 @@ class EpsilonTracker:
     """
     Updates epsilon according to linear schedule
     """
-    def __init__(self, selector: EpsilonGreedyActionSelector,
-                 eps_start: tt.Union[int, float],
-                 eps_final: tt.Union[int, float],
-                 eps_frames: int):
+    def __init__(self, selector: EpsilonGreedyActionSelector, eps_start: tt.Union[int, float], eps_final: tt.Union[int, float], eps_frames: int):
         self.selector = selector
         self.eps_start = eps_start
         self.eps_final = eps_final

simonbb-0.0.6/src/simonbb/cis433/dqn.py

@@ -0,0 +1,312 @@
+import random
+import argparse
+import numpy as np
+import torch
+import torch.nn as nn
+import dataclasses
+from datetime import timedelta
+import typing as tt
+import os
+import glob
+
+from ignite.engine import Engine
+from ignite.metrics import RunningAverage
+from ignite.contrib.handlers import tensorboard_logger as tb_logger
+from ray import tune
+
+from .experience import ExperienceFirstLast, ExperienceSourceFirstLast, ExperienceReplayBuffer
+from .ignite import EndOfEpisodeHandler, EpisodeEvents, EpisodeFPSHandler, PeriodicEvents, PeriodEvents
+
+#%%
+
+SEED = 123
+
+@dataclasses.dataclass
+class Hyperparams:
+    env_name: str
+    stop_reward: float
+    run_name: str
+    replay_size: int
+    replay_initial: int
+    target_net_sync: int
+    epsilon_frames: int
+
+    learning_rate: float = 0.0001
+    batch_size: int = 32
+    gamma: float = 0.99
+    epsilon_start: float = 1.0
+    epsilon_final: float = 0.1
+
+    tuner_mode: bool = False
+    episodes_to_solve: int = 500
+
+
+GAME_PARAMS = {
+    'pong': Hyperparams(
+        env_name="PongNoFrameskip-v4",
+        stop_reward=18.0,
+        run_name="pong",
+        replay_size=100_000,
+        replay_initial=10_000,
+        target_net_sync=1000,
+        epsilon_frames=100_000,
+        epsilon_final=0.02,
+    ),
+    'breakout-small': Hyperparams(
+        env_name="BreakoutNoFrameskip-v4",
+        stop_reward=500.0,
+        run_name="breakout-small",
+        replay_size=300_000,
+        replay_initial=20_000,
+        target_net_sync=1000,
+        epsilon_frames=1_000_000,
+        batch_size=64,
+    ),
+    'breakout': Hyperparams(
+        env_name="BreakoutNoFrameskip-v4",
+        stop_reward=500.0,
+        run_name='breakout',
+        replay_size=1_000_000,
+        replay_initial=50_000,
+        target_net_sync=10_000,
+        epsilon_frames=10_000_000,
+        learning_rate=0.00025,
+    ),
+    'invaders': Hyperparams(
+        env_name="SpaceInvadersNoFrameskip-v4",
+        stop_reward=500.0,
+        run_name='invaders',
+        replay_size=10_000_000,
+        replay_initial=50_000,
+        target_net_sync=10_000,
+        epsilon_frames=10_000_000,
+        learning_rate=0.00025,
+    ),
+}
+
+
+def unpack_batch( batch: tt.List[ExperienceFirstLast] ):
+    """
+    Convert a batch of transitions to Numpy arrays suitable for training
+    """
+    states, actions, rewards, dones, last_states = [],[],[],[],[]
+    for exp in batch:
+        states.append(exp.state)
+        actions.append(exp.action)
+        rewards.append(exp.reward)
+        dones.append(exp.last_state is None)
+        if exp.last_state is None: # end of episode
+            lstate = exp.state  # the result will be masked anyway, but this avoid special handling of None
+        else:
+            lstate = exp.last_state
+        last_states.append(lstate)
+    return np.asarray(states), np.array(actions), np.array(rewards, dtype=np.float32), np.array(dones, dtype=bool), np.asarray(last_states)
+
+
+def calc_loss_dqn( batch: tt.List[ExperienceFirstLast], net: nn.Module, tgt_net: nn.Module, gamma: float, device: torch.device) -> torch.Tensor:
+    states, actions, rewards, dones, next_states = unpack_batch(batch)
+    states_v = torch.as_tensor(states).to(device)
+    next_states_v = torch.as_tensor(next_states).to(device)
+    actions_v = torch.LongTensor(actions).to(device)
+    rewards_v = torch.FloatTensor(rewards).to(device)
+    done_mask = torch.BoolTensor(dones).to(device)
+
+    actions_v = actions_v.unsqueeze(-1)
+    state_action_vals = net(states_v).gather(1, actions_v).squeeze(-1)
+    with torch.no_grad():
+        next_state_vals = tgt_net(next_states_v).max(1)[0]
+        next_state_vals[done_mask] = 0.0
+
+    bellman_vals = next_state_vals.detach() * gamma + rewards_v
+    return nn.MSELoss()(state_action_vals, bellman_vals)
+
+
+def calc_loss_double_dqn( batch: tt.List[ExperienceFirstLast], net: nn.Module, tgt_net: nn.Module, gamma: float, device: torch.device):
+    states, actions, rewards, dones, next_states = unpack_batch(batch)
+    states_v = torch.as_tensor(states).to(device)
+    actions_v = torch.tensor(actions).to(device)
+    rewards_v = torch.tensor(rewards).to(device)
+    done_mask = torch.BoolTensor(dones).to(device)
+
+    actions_v = actions_v.unsqueeze(-1)
+    state_action_vals = net(states_v).gather(1, actions_v).squeeze(-1)
+    with torch.no_grad():
+        next_states_v = torch.as_tensor(next_states).to(device)
+        next_state_acts = net(next_states_v).max(1)[1] # use net() to get the best action for the next state
+        next_state_acts = next_state_acts.unsqueeze(-1)
+        next_state_vals = tgt_net(next_states_v).gather(1, next_state_acts).squeeze(-1) # use tgt_net to get the best action value corresponding to the selected action
+        next_state_vals[done_mask] = 0.0
+        exp_sa_vals = next_state_vals.detach() * gamma + rewards_v
+    return nn.MSELoss()(state_action_vals, exp_sa_vals)
+
+
+def batch_generator(buffer: ExperienceReplayBuffer, initial: int, batch_size: int) -> tt.Generator[tt.List[ExperienceFirstLast], None, None]:
+    buffer.populate(initial) # get initial fresh samples
+    while True:
+        buffer.populate(1)
+        yield buffer.sample(batch_size)
+
+
+@torch.no_grad()
+def calc_values_of_states(states: np.ndarray, net: nn.Module, device: torch.device):
+    mean_vals = []
+    for batch in np.array_split(states, 64):
+        states_v = torch.tensor(batch).to(device)
+        action_values_v = net(states_v)
+        best_action_values_v = action_values_v.max(1)[0]
+        mean_vals.append(best_action_values_v.mean().item())
+    return np.mean(mean_vals)
+
+
+def setup_ignite(
+        engine: Engine, params: Hyperparams, exp_source: ExperienceSourceFirstLast,
+        run_name: str, log_dir: str, extra_metrics: tt.Iterable[str] = (),
+        tuner_reward_episode: int = 100, tuner_reward_min: float = -19, 
+):
+    handler = EndOfEpisodeHandler( exp_source, bound_avg_reward=params.stop_reward )
+    handler.attach(engine)
+    EpisodeFPSHandler().attach(engine)
+
+    @engine.on(EpisodeEvents.EPISODE_COMPLETED)
+    def episode_completed(trainer: Engine):
+        passed = trainer.state.metrics.get('time_passed', 0)
+        print("Episode %d: reward=%.0f, steps=%s, speed=%.1f f/s, elapsed=%s" % (
+            trainer.state.episode, trainer.state.episode_reward,
+            trainer.state.episode_steps, trainer.state.metrics.get('avg_fps', 0),
+            timedelta(seconds=int(passed))))
+
+    @engine.on(EpisodeEvents.BOUND_REWARD_REACHED)
+    def game_solved(trainer: Engine):
+        passed = trainer.state.metrics['time_passed']
+        print("Game solved in %s, after %d episodes and %d iterations!" % (
+            timedelta(seconds=int(passed)), trainer.state.episode,
+            trainer.state.iteration))
+        trainer.should_terminate = True
+        trainer.state.solved = True
+
+    tb = tb_logger.TensorboardLogger(log_dir=log_dir)
+    run_avg = RunningAverage(output_transform=lambda v: v['loss'])
+    run_avg.attach(engine, "avg_loss")
+
+    metrics = ['reward', 'steps', 'avg_reward']
+    handler = tb_logger.OutputHandler(tag="episodes", metric_names=metrics)
+    event = EpisodeEvents.EPISODE_COMPLETED
+    tb.attach(engine, log_handler=handler, event_name=event)
+
+    # write to tensorboard every 100 iterations
+    PeriodicEvents().attach(engine)
+    metrics = ['avg_loss', 'avg_fps']
+    metrics.extend(extra_metrics)
+    handler = tb_logger.OutputHandler(tag="train", metric_names=metrics, output_transform=lambda a: a)
+    event = PeriodEvents.ITERS_100_COMPLETED
+    tb.attach(engine, log_handler=handler, event_name=event)
+
+    if params.tuner_mode:
+        @engine.on(EpisodeEvents.EPISODE_COMPLETED)
+        def episode_completed(trainer: Engine):
+            avg_reward = trainer.state.metrics.get('avg_reward')
+            max_episodes = params.episodes_to_solve * 1.1
+            if trainer.state.episode > tuner_reward_episode and \
+                    avg_reward < tuner_reward_min:
+                trainer.should_terminate = True
+                trainer.state.solved = False
+            elif trainer.state.episode > max_episodes:
+                trainer.should_terminate = True
+                trainer.state.solved = False
+            if trainer.should_terminate:
+                print(f"Episode {trainer.state.episode}, "
+                      f"avg_reward {avg_reward:.2f}, terminating")
+
+
+# hyperparams tuner
+TrainFunc = tt.Callable[  [Hyperparams, torch.device, dict],  tt.Optional[int] ]
+
+BASE_SPACE = {
+    "learning_rate": tune.loguniform(1e-5, 1e-4),
+    "gamma": tune.choice([0.9, 0.92, 0.95, 0.98, 0.99, 0.995]),
+}
+
+def tune_params(
+        base_params: Hyperparams, train_func: TrainFunc, device: torch.device,
+        samples: int = 10, extra_space: tt.Optional[tt.Dict[str, tt.Any]] = None,
+):
+    """
+    Perform hyperparameters tune.
+    :param train_func: Train function, has to return "episodes" key with metric
+    :param device: torch device
+    :param samples: count of samples to perform
+    :param extra_space: additional search space
+    """
+    search_space = dict(BASE_SPACE)
+    if extra_space is not None:
+        search_space.update(extra_space)
+    config = tune.TuneConfig(num_samples=samples)
+
+    def objective(config: dict, device: torch.device) -> dict:
+        keys = dataclasses.asdict(base_params).keys()
+        upd = {"tuner_mode": True}
+        for k, v in config.items():
+            if k in keys:
+                upd[k] = v
+        params = dataclasses.replace(base_params, **upd)
+        res = train_func(params, device, config)
+        return {"episodes": res if res is not None else 10**6}
+
+    obj = tune.with_parameters(objective, device=device)
+    if device.type == "cuda":
+        obj = tune.with_resources(obj, {"gpu": 1})
+    tuner = tune.Tuner(obj, param_space=search_space, tune_config=config)
+    results = tuner.fit()
+    best = results.get_best_result(metric="episodes", mode="min")
+    print(best.config)
+    print(best.metrics)
+
+
+def train_or_tune(
+        args: argparse.Namespace,
+        train_func: TrainFunc,
+	device: torch.device,
+        best_params: Hyperparams,
+        extra_params: tt.Optional[dict] = None,
+        extra_space: tt.Optional[dict] = None,
+):
+    random.seed(SEED)
+    torch.manual_seed(SEED)
+    if args.params == "common":
+        params = GAME_PARAMS['pong']
+    else:
+        params = best_params
+
+    if extra_params is None:
+        extra_params = {}
+    if args.tune is None:
+        train_func(params, device, extra_params)
+    else:
+        tune_params(params, train_func, device, samples=args.tune, extra_space=extra_space)
+
+
+def save_model(model, path, global_step, start_iteration):
+    torch.save({
+        'model_state_dict': model.state_dict(),
+        'global_step': global_step,
+        'start_iteration': start_iteration
+    }, path)
+
+
+def load_model(model, path):    
+    if not os.path.exists(path):
+        raise FileNotFoundError(f"No checkpoint found at {path}")   
+    checkpoint = torch.load(path) # Load the checkpoint    
+    if 'model_state_dict' in checkpoint:
+        model.load_state_dict(checkpoint['model_state_dict'])
+    else:
+        raise KeyError(f"'model_state_dict' not found in the checkpoint at {path}")
+    model.eval()  # Set the model to evaluation mode
+    return checkpoint  # Return the checkpoint to access global_step
+
+def get_latest_checkpoint(dir_path):    
+    checkpoint_files = glob.glob(os.path.join(dir_path, '*.pth'))
+    if not checkpoint_files:
+        return None    
+    latest_checkpoint = max(checkpoint_files, key=os.path.getmtime)
+    return latest_checkpoint

simonbb-0.0.6/src/simonbb/cis433/wrappers.py

@@ -0,0 +1,100 @@
+import typing as tt
+import gymnasium as gym
+from gymnasium import spaces
+import collections
+import numpy as np
+from stable_baselines3.common import atari_wrappers
+
+class ImageToPyTorch(gym.ObservationWrapper):
+    def __init__(self, env):
+        super(ImageToPyTorch, self).__init__(env)
+        obs = self.observation_space
+        assert isinstance(obs, gym.spaces.Box)
+        assert len(obs.shape) == 3
+        new_shape = (obs.shape[-1], obs.shape[0], obs.shape[1])
+        self.observation_space = gym.spaces.Box(
+            low=obs.low.min(), high=obs.high.max(),
+            shape=new_shape, dtype=obs.dtype)
+
+    def observation(self, observation):
+        return np.moveaxis(observation, 2, 0)
+
+
+class FrameStack1D(gym.Wrapper):
+    """
+    Stacks observations into 1D array
+    """
+    def __init__(self, env, k):
+        assert isinstance(env, gym.Env)
+        assert isinstance(env.observation_space, gym.spaces.Box)
+        assert len(env.observation_space.shape) == 1
+        super(FrameStack1D, self).__init__(env)
+        self.k = k
+        self.frames = collections.deque(maxlen=k)
+        shp = env.observation_space.shape
+        self.observation_space = gym.spaces.Box(
+            low=np.min(env.observation_space.low),
+            high=np.max(env.observation_space.high),
+            shape=(env.observation_space.shape[0]*k,),
+            dtype=env.observation_space.dtype)
+
+    def reset(self):
+        ob, extra = self.env.reset()
+        for _ in range(self.k):
+            self.frames.append(ob)
+        return self._get_ob(), extra
+
+    def step(self, action):
+        ob, reward, done, trunc, info = self.env.step(action)
+        self.frames.append(ob)
+        return self._get_ob(), reward, done, trunc, info
+
+    def _get_ob(self):
+        assert len(self.frames) == self.k
+        return np.concatenate(self.frames)
+
+
+
+class BufferWrapper(gym.ObservationWrapper):
+    def __init__(self, env, n_steps):
+        super(BufferWrapper, self).__init__(env)
+        obs = env.observation_space
+        assert isinstance(obs, spaces.Box)
+        new_obs = gym.spaces.Box(
+            obs.low.repeat(n_steps, axis=0),
+            obs.high.repeat(n_steps, axis=0), dtype=obs.dtype)
+        self.observation_space = new_obs
+        self.buffer = collections.deque(maxlen=n_steps)
+
+    def reset(self, *, seed: tt.Optional[int] = None,
+              options: tt.Optional[dict[str, tt.Any]] = None):
+        for _ in range(self.buffer.maxlen-1):
+            self.buffer.append(self.env.observation_space.low)
+        obs, extra = self.env.reset()
+        return self.observation(obs), extra
+
+    def observation(self, observation: np.ndarray) -> np.ndarray:
+        self.buffer.append(observation)
+        return np.concatenate(self.buffer)
+
+
+def wrap_dqn(env: gym.Env, stack_frames: int = 4,
+             episodic_life: bool = True, clip_reward: bool = True,
+             noop_max: int = 0) -> gym.Env:
+    """
+    Apply a common set of wrappers for Atari games.
+    :param env: Environment to wrap
+    :param stack_frames: count of frames to stack, default=4
+    :param episodic_life: convert life to end of episode
+    :param clip_reward: reward clipping
+    :param noop_max: how many NOOP actions to execute
+    :return: wrapped environment
+    """
+    assert 'NoFrameskip' in env.spec.id
+    env = atari_wrappers.AtariWrapper(
+        env, clip_reward=clip_reward, noop_max=noop_max,
+        terminal_on_life_loss=episodic_life)
+    env = ImageToPyTorch(env)
+    env = BufferWrapper(env, stack_frames)
+    return env
+

{simonbb-0.0.4 → simonbb-0.0.6}/src/simonbb/utils.py

@@ -1,11 +1,32 @@
+import time
 import re
+import collections
+import sys
 import pickle
 import json
 from collections import Counter
 import urllib.request
 import numpy as np
+import torch
+import torch.nn as nn
 
+class WeightedMSELoss(nn.Module):
+    def __init__(self, size_average=True):
+        super(WeightedMSELoss, self).__init__()
+        self.size_average = size_average
 
+    def forward(self, input, target, weights=None):
+        if weights is None:
+            return nn.MSELoss(self.size_average)(input, target)
+
+        loss_rows = (input - target) ** 2
+        if len(loss_rows.size()) != 1:
+            loss_rows = torch.sum(loss_rows, dim=1)
+        res = (weights * loss_rows).sum()
+        if self.size_average:
+            res /= len(weights)
+        return res
+    
 class TextVectorizer:
     def __init__(self, max_tokens=20000, output_mode="int", output_sequence_length=None, ngrams=None, standardize=None):
         self.max_tokens = max_tokens
@@ -123,6 +144,97 @@ class TextVectorizer:
         return vectorizer
 
 
+
+class TBMeanTracker:
+    """
+    TensorBoard value tracker: allows to batch fixed amount of historical values and write their mean into TB
+
+    Designed and tested with pytorch-tensorboard in mind
+    """
+    def __init__(self, writer, batch_size):
+        """
+        :param writer: writer with close() and add_scalar() methods
+        :param batch_size: integer size of batch to track
+        """
+        assert isinstance(batch_size, int)
+        assert writer is not None
+        self.writer = writer
+        self.batch_size = batch_size
+
+    def __enter__(self):
+        self._batches = collections.defaultdict(list)
+        return self
+
+    def __exit__(self, exc_type, exc_val, exc_tb):
+        self.writer.close()
+
+    @staticmethod
+    def _as_float(value):
+        assert isinstance(value, (float, int, np.ndarray, np.generic, torch.autograd.Variable)) or torch.is_tensor(value)
+        tensor_val = None
+        if isinstance(value, torch.autograd.Variable):
+            tensor_val = value.data
+        elif torch.is_tensor(value):
+            tensor_val = value
+
+        if tensor_val is not None:
+            return tensor_val.float().mean().item()
+        elif isinstance(value, np.ndarray):
+            return float(np.mean(value))
+        else:
+            return float(value)
+
+    def track(self, param_name, value, iter_index):
+        assert isinstance(param_name, str)
+        assert isinstance(iter_index, int)
+
+        data = self._batches[param_name]
+        data.append(self._as_float(value))
+
+        if len(data) >= self.batch_size:
+            self.writer.add_scalar(param_name, np.mean(data), iter_index)
+            data.clear()
+
+
+class RewardTracker:
+    def __init__(self, writer, min_ts_diff=1.0):
+        """
+        Constructs RewardTracker
+        :param writer: writer to use for writing stats
+        :param min_ts_diff: minimal time difference to track speed
+        """
+        self.writer = writer
+        self.min_ts_diff = min_ts_diff
+
+    def __enter__(self):
+        self.ts = time.time()
+        self.ts_frame = 0
+        self.total_rewards = []
+        return self
+
+    def __exit__(self, *args):
+        self.writer.close()
+
+    def reward(self, reward, frame, epsilon=None):
+        self.total_rewards.append(reward)
+        mean_reward = np.mean(self.total_rewards[-100:])
+        ts_diff = time.time() - self.ts
+        if ts_diff > self.min_ts_diff:
+            speed = (frame - self.ts_frame) / ts_diff
+            self.ts_frame = frame
+            self.ts = time.time()
+            epsilon_str = "" if epsilon is None else ", eps %.2f" % epsilon
+            print("%d: done %d episodes, mean reward %.3f, speed %.2f f/s%s" % (
+                frame, len(self.total_rewards), mean_reward, speed, epsilon_str
+            ))
+            sys.stdout.flush()
+            self.writer.add_scalar("speed", speed, frame)
+        if epsilon is not None:
+            self.writer.add_scalar("epsilon", epsilon, frame)
+        self.writer.add_scalar("reward_100", mean_reward, frame)
+        self.writer.add_scalar("reward", reward, frame)
+        return mean_reward if len(self.total_rewards) > 30 else None
+    
 def query_ollama(prompt, model="llama3", url="http://localhost:11434/api/chat"):
     # Create the data payload as a dictionary
     data = {