Harness context codex item completed
{
"type": "codex_item_completed",
"item": {
"type": "CommandExecution",
"id": "exec-f8c50588-162c-4f0b-803a-7985bed4d1dd",
"process_id": "98962",
"command": [
"/bin/sh",
"-lc",
"sed -n '1080,1515p' /app/solution/train_current_model.py; sed -n '280,346p' /app/solution/two_tower_model.py; sed -n '320,540p' /app/solution/train_current_model.py; sed -n '880,1085p' /app/solution/train_current_model.py"
],
"cwd": "file:///app",
"parsed_cmd": [
{
"type": "read",
"cmd": "sed -n '1080,1515p' /app/solution/train_current_model.py",
"name": "train_current_model.py",
"path": "/app/solution/train_current_model.py"
},
{
"type": "read",
"cmd": "sed -n '280,346p' /app/solution/two_tower_model.py",
"name": "two_tower_model.py",
"path": "/app/solution/two_tower_model.py"
},
{
"type": "read",
"cmd": "sed -n '320,540p' /app/solution/train_current_model.py",
"name": "train_current_model.py",
"path": "/app/solution/train_current_model.py"
},
{
"type": "read",
"cmd": "sed -n '880,1085p' /app/solution/train_current_model.py",
"name": "train_current_model.py",
"path": "/app/solution/train_current_model.py"
}
],
"source": "unified_exec_startup",
"status": "completed",
"stdout": " event_features[row, :count] = prepared.event_features[positions]\n activity_features[row, :count] = prepared.activity_features[positions]\n state_ids[row, :count] = prepared.state_ids[positions]\n targets[row, :count] = prepared.sequence_items[target_positions]\n weights[row, :count] = target_weights[target_positions]\n score_from = min(count, max(int(score_starts[row]) - start, 0))\n scored[row, score_from:count] = True\n processed[row] = count\n user_sequence = prepared.sequence_items[\n offset : offset + int(lengths[user_row])\n ]\n histories[row, :count] = history_windows(\n user_sequence,\n absolute_start,\n count,\n config.bptt_steps,\n padding_item,\n )\n return (\n input_items,\n event_features,\n activity_features,\n state_ids,\n targets,\n scored,\n histories,\n processed,\n weights,\n )\n\n\ndef score_target_microbatches(\n model: base.SequentialTwoTower,\n prepared: SequenceData,\n config: Config,\n sequence_vectors: torch.Tensor,\n flat_states: np.ndarray,\n flat_activity: np.ndarray,\n positive_items: torch.Tensor,\n negatives: torch.Tensor,\n negative_weights: torch.Tensor,\n sample_weights: np.ndarray,\n device: torch.device,\n backward: bool,\n) -> tuple[dict[str, float], torch.Tensor | None]:\n \"\"\"Compute one logical block's exact BCE in bounded-memory slices.\"\"\"\n target_count = len(positive_items)\n total_weight = float(sample_weights.sum())\n if total_weight <= 0:\n raise RuntimeError(\"target block has no positive sample weight\")\n sequence_gradient = torch.zeros_like(sequence_vectors) if backward else None\n statistics = {\n \"loss_sum\": 0.0,\n \"positive_sum\": 0.0,\n \"negative_sum\": 0.0,\n \"pair_correct\": 0.0,\n \"weight_sum\": 0.0,\n }\n for begin in range(0, target_count, config.loss_microbatch_size):\n end = min(begin + config.loss_microbatch_size, target_count)\n selected_states = flat_states[begin:end]\n if backward:\n sequence_input = (\n sequence_vectors[begin:end].detach().requires_grad_(True)\n )\n else:\n sequence_input = sequence_vectors[begin:end]\n liked_vectors, disliked_vectors, _ = pool_explicit_states(\n model, prepared, selected_states, device\n )\n counts = prepared.state_counts[selected_states]\n numeric = np.zeros((end - begin, 6), dtype=np.float32)\n numeric[:, :4] = flat_activity[begin:end]\n numeric[:, 4] = np.log1p(counts[:, 0]) / 6.0\n numeric[:, 5] = np.log1p(counts[:, 1]) / 4.0\n query_vectors = user_queries(\n model,\n sequence_input,\n liked_vectors,\n disliked_vectors,\n torch.from_numpy(numeric).to(device),\n )\n positive_vectors, negative_vectors = candidate_vectors(\n model,\n positive_items[begin:end],\n negatives[begin:end],\n )\n positive_logits = (query_vectors * positive_vectors).sum(dim=1)\n negative_logits = torch.einsum(\n \"bd,bkd->bk\", query_vectors, negative_vectors\n )\n positive_loss = F.softplus(-positive_logits)\n negative_loss = corrected_negative_loss(\n negative_logits, negative_weights[begin:end], config\n )\n selected_weights = torch.from_numpy(\n sample_weights[begin:end].astype(np.float32, copy=False)\n ).to(device)\n loss_sum = (\n selected_weights * (positive_loss + negative_loss) / 2.0\n ).sum()\n if not bool(torch.isfinite(loss_sum)):\n raise RuntimeError(\"dense next-item training produced non-finite loss\")\n if backward:\n (loss_sum / total_weight).backward()\n if sequence_input.grad is None:\n raise RuntimeError(\"sequence microbatch did not produce a gradient\")\n sequence_gradient[begin:end] = sequence_input.grad\n statistics[\"loss_sum\"] += float(loss_sum.detach().cpu())\n statistics[\"positive_sum\"] += float(positive_logits.detach().sum().cpu())\n statistics[\"negative_sum\"] += float(negative_logits.detach().sum().cpu())\n statistics[\"pair_correct\"] += float(\n (positive_logits.detach()[:, None] > negative_logits.detach()).sum().cpu()\n )\n statistics[\"weight_sum\"] += float(selected_weights.sum().detach().cpu())\n return statistics, sequence_gradient\n\n\ndef release_accelerator_cache(device: torch.device) -> None:\n if device.type == \"mps\":\n torch.mps.empty_cache()\n\n\ndef run_training_epoch(\n model: base.SequentialTwoTower,\n optimizer: torch.optim.Optimizer,\n prepared: SequenceData,\n config: Config,\n target_weights: np.ndarray,\n window_starts: np.ndarray,\n target_counts: np.ndarray,\n item_correction: torch.Tensor,\n known_positive_packed: np.ndarray | None,\n generator: torch.Generator,\n device: torch.device,\n epoch: int,\n) -> dict[str, float]:\n model.train()\n train_users = np.flatnonzero(target_counts > 0)\n sorted_users = train_users[\n np.argsort(target_counts[train_users], kind=\"stable\")\n ]\n groups = [\n sorted_users[start : start + config.user_batch_size]\n for start in range(0, len(sorted_users), config.user_batch_size)\n ]\n group_order = torch.randperm(len(groups), generator=generator).numpy()\n negative_count = config.uniform_negatives + config.inbatch_negatives\n total_loss = 0.0\n total_targets = 0\n total_weight = 0.0\n total_positive = 0.0\n total_negative = 0.0\n total_pair_correct = 0.0\n batches = 0\n started = time.monotonic()\n\n for group_number in group_order:\n users = groups[int(group_number)]\n group_counts = target_counts[users]\n hidden = torch.zeros(\n 1, len(users), config.embedding_dim, device=device\n )\n max_targets = int(group_counts.max())\n for start in range(0, max_targets, config.bptt_steps):\n steps = min(config.bptt_steps, max_targets - start)\n block = build_dense_block(\n prepared,\n target_weights,\n users,\n window_starts[users],\n group_counts,\n np.zeros(len(users), dtype=np.int64),\n start,\n steps,\n config,\n )\n sequence_outputs, _ = event_sequence(\n model,\n torch.from_numpy(block[0].astype(np.int64)).to(device),\n torch.from_numpy(block[1].astype(np.float32)).to(device),\n hidden,\n )\n hidden = torch.stack(\n [\n sequence_outputs[row, count - 1]\n if count\n else hidden[0, row]\n for row, count in enumerate(block[7])\n ]\n ).unsqueeze(0).detach()\n flat_scored = block[5].reshape(-1)\n if not np.any(flat_scored):\n continue\n mask = torch.from_numpy(flat_scored).to(device)\n sequence_vectors = sequence_outputs.reshape(\n -1, config.embedding_dim\n )[mask]\n flat_states = block[3].reshape(-1)[flat_scored]\n positive_items = torch.from_numpy(\n block[4].reshape(-1)[flat_scored].astype(np.int64)\n ).to(device)\n known_history = torch.from_numpy(\n block[6].reshape(-1, config.bptt_steps)[flat_scored].astype(\n np.int64\n )\n ).to(device)\n liked_items = torch.from_numpy(\n prepared.liked_states[flat_states].astype(np.int64)\n ).to(device)\n flat_user_rows = np.broadcast_to(\n users[:, None], block[5].shape\n ).reshape(-1)[flat_scored]\n negatives, negative_weights = draw_training_negatives(\n positive_items,\n known_history,\n liked_items,\n flat_user_rows,\n known_positive_packed,\n prepared,\n config,\n item_correction,\n generator,\n device,\n )\n optimizer.zero_grad(set_to_none=True)\n statistics, sequence_gradient = score_target_microbatches(\n model,\n prepared,\n config,\n sequence_vectors,\n flat_states,\n block[2].reshape(-1, 4)[flat_scored].astype(np.float32),\n positive_items,\n negatives,\n negative_weights,\n block[8].reshape(-1)[flat_scored],\n device,\n backward=True,\n )\n if sequence_gradient is None:\n raise RuntimeError(\"training block omitted its sequence gradient\")\n sequence_vectors.backward(sequence_gradient)\n nn.utils.clip_grad_norm_(model.parameters(), config.gradient_clip)\n optimizer.step()\n\n count = len(positive_items)\n total_loss += statistics[\"loss_sum\"]\n total_targets += count\n total_weight += statistics[\"weight_sum\"]\n total_positive += statistics[\"positive_sum\"]\n total_negative += statistics[\"negative_sum\"]\n total_pair_correct += statistics[\"pair_correct\"]\n batches += 1\n if batches % 25 == 0:\n release_accelerator_cache(device)\n if batches % 250 == 0:\n batch_bce = statistics[\"loss_sum\"] / max(\n statistics[\"weight_sum\"], 1.0e-12\n )\n running_bce = total_loss / max(total_weight, 1.0e-12)\n log(\n f\"epoch={epoch}/{config.epochs} batches={batches:,} \"\n f\"targets={total_targets:,} \"\n f\"batch_bce={batch_bce:.6f} \"\n f\"running_bce={running_bce:.6f} \"\n f\"elapsed={time.monotonic() - started:.1f}s\"\n )\n\n return {\n \"epoch\": float(epoch),\n \"loss\": total_loss / max(total_weight, 1.0e-12),\n \"positive_logit\": total_positive / max(total_targets, 1),\n \"negative_logit\": total_negative / max(total_targets * negative_count, 1),\n \"pair_accuracy\": total_pair_correct\n / max(total_targets * negative_count, 1),\n \"targets\": float(total_targets),\n \"target_weight_sum\": total_weight,\n \"batches\": float(batches),\n \"seconds\": time.monotonic() - started,\n }\n\n\[email protected]_grad()\ndef evaluate_validation_bce(\n model: base.SequentialTwoTower,\n prepared: SequenceData,\n config: Config,\n target_weights: np.ndarray,\n window_starts: np.ndarray,\n training_counts: np.ndarray,\n validation_counts: np.ndarray,\n item_correction: torch.Tensor,\n known_positive_packed: np.ndarray | None,\n device: torch.device,\n) -> dict[str, float]:\n \"\"\"Measure deterministic BCE on each user's chronological target suffix.\"\"\"\n model.eval()\n totals = training_counts + validation_counts\n users_with_validation = np.flatnonzero(validation_counts > 0)\n ordered = users_with_validation[\n np.argsort(totals[users_with_validation], kind=\"stable\")\n ]\n groups = [\n ordered[start : start + config.user_batch_size]\n for start in range(0, len(ordered), config.user_batch_size)\n ]\n generator = torch.Generator(device=\"cpu\").manual_seed(config.seed + 20_000)\n total_loss = 0.0\n total_targets = 0\n total_weight = 0.0\n started = time.monotonic()\n\n for group_number, users in enumerate(groups, start=1):\n group_totals = totals[users]\n group_starts = training_counts[users]\n hidden = torch.zeros(\n 1, len(users), config.embedding_dim, device=device\n )\n max_targets = int(group_totals.max())\n for start in range(0, max_targets, config.bptt_steps):\n steps = min(config.bptt_steps, max_targets - start)\n block = build_dense_block(\n prepared,\n target_weights,\n users,\n window_starts[users],\n group_totals,\n group_starts,\n start,\n steps,\n config,\n )\n sequence_outputs, _ = event_sequence(\n model,\n torch.from_numpy(block[0].astype(np.int64)).to(device),\n torch.from_numpy(block[1].astype(np.float32)).to(device),\n hidden,\n )\n hidden = torch.stack(\n [\n sequence_outputs[row, count - 1]\n if count\n else hidden[0, row]\n for row, count in enumerate(block[7])\n ]\n ).unsqueeze(0)\n flat_scored = block[5].reshape(-1)\n if not np.any(flat_scored):\n continue\n mask = torch.from_numpy(flat_scored).to(device)\n sequence_vectors = sequence_outputs.reshape(\n -1, config.embedding_dim\n )[mask]\n flat_states = block[3].reshape(-1)[flat_scored]\n positive_items = torch.from_numpy(\n block[4].reshape(-1)[flat_scored].astype(np.int64)\n ).to(device)\n known_history = torch.from_numpy(\n block[6].reshape(-1, config.bptt_steps)[flat_scored].astype(\n np.int64\n )\n ).to(device)\n liked_items = torch.from_numpy(\n prepared.liked_states[flat_states].astype(np.int64)\n ).to(device)\n flat_user_rows = np.broadcast_to(\n users[:, None], block[5].shape\n ).reshape(-1)[flat_scored]\n negatives, negative_weights = draw_training_negatives(\n positive_items,\n known_history,\n liked_items,\n flat_user_rows,\n known_positive_packed,\n prepared,\n config,\n item_correction,\n generator,\n device,\n )\n statistics, _ = score_target_microbatches(\n model,\n prepared,\n config,\n sequence_vectors,\n flat_states,\n block[2].reshape(-1, 4)[flat_scored].astype(np.float32),\n positive_items,\n negatives,\n negative_weights,\n block[8].reshape(-1)[flat_scored],\n device,\n backward=False,\n )\n total_loss += statistics[\"loss_sum\"]\n total_targets += len(positive_items)\n total_weight += statistics[\"weight_sum\"]\n if group_number % 25 == 0:\n release_accelerator_cache(device)\n return {\n \"validation_loss\": total_loss / max(total_weight, 1.0e-12),\n \"validation_targets\": float(total_targets),\n \"validation_target_weight_sum\": total_weight,\n \"validation_seconds\": time.monotonic() - started,\n }\n\n\ndef snapshot_model(model: nn.Module) -> dict[str, torch.Tensor]:\n return {\n name: value.detach().cpu().clone()\n for name, value in model.state_dict().items()\n }\n\n\ndef save_epoch_checkpoint(\n checkpoint_dir: Path,\n epoch: int,\n model_state: dict[str, torch.Tensor],\n config: Config,\n diagnostics: dict[str, float],\n optimizer: torch.optim.Optimizer,\n generator: torch.Generator,\n) -> Path:\n checkpoint_dir.mkdir(parents=True, exist_ok=True)\n path = checkpoint_dir / f\"epoch-{epoch:02d}.pt\"\n torch.save(\n {\n \"epoch\": epoch,\n \"config\": asdict(config),\n \"diagnostics\": diagnostics,\n \"model_state\": model_state,\n \"optimizer_state\": optimizer.state_dict(),\n \"generator_state\": generator.get_state(),\n },\n path,\n destination[row, : len(retained)] = retained.astype(np.uint16)\n return liked, disliked, like_counts, dislike_counts\n\n\nclass SequentialTwoTower(nn.Module):\n \"\"\"Order-sensitive user tower and independent catalog-aware item tower.\"\"\"\n\n def __init__(self, item_features: np.ndarray, config: TowerConfig) -> None:\n super().__init__()\n n_items, n_item_features = item_features.shape\n dim = config.embedding_dim\n self.n_items = n_items\n self.context_embedding = nn.Embedding(\n n_items + 1, dim, padding_idx=n_items\n )\n self.event_projection = nn.Sequential(\n nn.Linear(dim + 3, config.hidden_dim),\n nn.GELU(),\n )\n self.sequence_encoder = nn.GRU(\n input_size=config.hidden_dim,\n hidden_size=dim,\n num_layers=1,\n batch_first=True,\n )\n self.user_tower = nn.Sequential(\n nn.Linear(dim * 3 + 6, config.user_hidden_1),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.user_hidden_1, config.user_hidden_2),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.user_hidden_2, dim),\n )\n self.candidate_embedding = nn.Embedding(n_items, dim)\n self.register_buffer(\n \"item_features\", torch.from_numpy(item_features), persistent=True\n )\n self.item_tower = nn.Sequential(\n nn.Linear(dim + n_item_features, config.item_hidden_1),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.item_hidden_1, config.item_hidden_2),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.item_hidden_2, dim),\n )\n nn.init.normal_(self.context_embedding.weight, std=0.02)\n nn.init.normal_(self.candidate_embedding.weight, std=0.02)\n with torch.no_grad():\n self.context_embedding.weight[n_items].zero_()\n\n def _mean_context_channel(self, items: torch.Tensor) -> torch.Tensor:\n mask = items.ne(self.n_items)\n vectors = self.context_embedding(items)\n return (vectors * mask.unsqueeze(-1)).sum(dim=1) / mask.sum(\n dim=1, keepdim=True\n ).clamp_min(1)\n\n def encode_items(self, items: torch.Tensor) -> torch.Tensor:\n inputs = torch.cat(\n [self.candidate_embedding(items), self.item_features[items]], dim=1\n )\n return self.item_tower(inputs)\n\n\ndef invalid_negative_mask(\n dislike_count: int,\n) -> np.ndarray:\n start = int(\n np.searchsorted(\n timestamps,\n history_end - base.PROFILE_LOOKBACK_DAYS * base.DAY,\n side=\"left\",\n )\n )\n result = np.zeros(6, dtype=np.float32)\n recent_items = items[start:]\n if len(recent_items):\n result[0] = np.log1p(len(recent_items)) / 8.0\n result[1] = np.log1p(len(np.unique(recent_items))) / 6.0\n result[2] = float(np.mean(event_features[start:, 0]))\n result[3] = float(np.mean(event_features[start:, 1]))\n result[4] = np.log1p(like_count) / 6.0\n result[5] = np.log1p(dislike_count) / 4.0\n return result\n\n\ndef prepare_data(\n public: Path, stage: str, config: Config\n) -> tuple[SequenceData, CandidateSet]:\n data = public / \"data\"\n catalog = pd.read_parquet(data / \"candidate_catalog.parquet\").sort_values(\n \"popularity_rank\"\n )\n item_ids = catalog.item_id.to_numpy(dtype=np.int64, copy=True)\n item_index = pd.Index(item_ids)\n n_items = len(item_ids)\n if n_items >= np.iinfo(np.uint16).max:\n raise RuntimeError(\"uint16 item storage is too small for this catalog\")\n\n if stage == \"validation\":\n prefixes = (\"train\",)\n history_end = base.TRAIN_END\n candidate_name = \"validation_eligible_candidates.npz\"\n state_name = \"validation_preference_state.parquet\"\n else:\n prefixes = (\"train\", \"val\")\n history_end = base.VALIDATION_END\n candidate_name = \"test_eligible_candidates.npz\"\n state_name = \"test_preference_state.parquet\"\n eligibility = CandidateSet.read(data / candidate_name)\n if not np.array_equal(eligibility.item_ids, item_ids):\n raise RuntimeError(\"catalog and candidate artifact item order differ\")\n\n explicit_by_user = base.load_explicit_events(data, prefixes, item_index)\n inference_likes, inference_dislikes, like_counts, dislike_counts = (\n base.load_inference_state(\n data / state_name,\n item_index,\n eligibility.uids,\n n_items,\n )\n )\n inference_rows = {int(uid): row for row, uid in enumerate(eligibility.uids)}\n inference_user_features = np.zeros((len(eligibility.uids), 6), dtype=np.float32)\n inference_sequence_rows = np.full(len(eligibility.uids), -1, dtype=np.int32)\n\n # Snapshot zero is the shared empty explicit state.\n liked_states = [np.full(base.LIKE_ITEMS, n_items, dtype=np.uint16)]\n disliked_states = [\n np.full(base.DISLIKE_ITEMS, n_items, dtype=np.uint16)\n ]\n state_counts: list[tuple[int, int]] = [(0, 0)]\n user_ids: list[int] = []\n offsets = [0]\n timestamp_parts: list[np.ndarray] = []\n item_parts: list[np.ndarray] = []\n event_parts: list[np.ndarray] = []\n activity_parts: list[np.ndarray] = []\n state_parts: list[np.ndarray] = []\n\n listen_paths = [data / f\"{prefix}_listens.parquet\" for prefix in prefixes]\n started = time.monotonic()\n for seen_users, (uid, listens) in enumerate(\n base.merge_listen_streams(listen_paths, item_index), start=1\n ):\n parts = build_sequence_parts(\n listens,\n explicit_by_user.get(uid, base.empty_explicit_events()),\n n_items,\n liked_states,\n disliked_states,\n state_counts,\n )\n if parts is None:\n continue\n sequence_row = len(user_ids)\n user_ids.append(uid)\n timestamp_parts.append(parts[0].astype(np.uint32, copy=False))\n item_parts.append(parts[1])\n event_parts.append(parts[2])\n activity_parts.append(parts[3])\n state_parts.append(parts[4])\n offsets.append(offsets[-1] + len(parts[1]))\n inference_row = inference_rows.get(uid)\n if inference_row is not None:\n inference_sequence_rows[inference_row] = sequence_row\n inference_user_features[inference_row] = final_user_features(\n parts[0],\n parts[1],\n parts[2],\n history_end,\n int(like_counts[inference_row]),\n int(dislike_counts[inference_row]),\n )\n if seen_users % 1_000 == 0:\n log(\n f\"prepared users={seen_users:,} tokens={offsets[-1]:,} \"\n f\"states={len(liked_states):,} \"\n f\"elapsed={time.monotonic() - started:.1f}s\"\n )\n\n missing = np.flatnonzero(inference_sequence_rows < 0)\n if len(missing):\n inference_user_features[missing, 4] = np.log1p(like_counts[missing]) / 6.0\n inference_user_features[missing, 5] = (\n np.log1p(dislike_counts[missing]) / 4.0\n )\n log(f\"target users without qualifying catalog sequence={len(missing):,}\")\n\n prepared = SequenceData(\n item_ids=item_ids,\n item_features=base.build_item_features(catalog),\n user_ids=np.asarray(user_ids, dtype=np.int64),\n offsets=np.asarray(offsets, dtype=np.int64),\n sequence_timestamps=np.concatenate(timestamp_parts),\n sequence_items=np.concatenate(item_parts),\n event_features=np.concatenate(event_parts),\n activity_features=np.concatenate(activity_parts),\n state_ids=np.concatenate(state_parts),\n liked_states=np.stack(liked_states),\n disliked_states=np.stack(disliked_states),\n state_counts=np.asarray(state_counts, dtype=np.int32),\n inference_user_rows=inference_sequence_rows,\n inference_liked_items=inference_likes,\n inference_disliked_items=inference_dislikes,\n inference_user_features=inference_user_features,\n )\n log(\n f\"prepared stage={stage} sequence_users={len(prepared.user_ids):,} \"\n f\"tokens={len(prepared.sequence_items):,} \"\n f\"next_item_targets={prepared.target_count:,} \"\n f\"states={len(prepared.liked_states):,}\"\n )\n return prepared, eligibility\n\n\ndef cache_metadata(stage: str, config: Config) -> dict[str, object]:\n return {\n \"cache_version\": CACHE_VERSION,\n \"stage\": stage,\n \"positive_listen_threshold\": base.LISTEN_POSITIVE_THRESHOLD,\n \"profile_lookback_days\": base.PROFILE_LOOKBACK_DAYS,\n \"event_time_feature\": \"exp-negative-gap-over-60-days\",\n }\n\n\ndef save_cache(\n path: Path, prepared: SequenceData, stage: str, config: Config\n) -> None:\n path.parent.mkdir(parents=True, exist_ok=True)\n np.savez(\n path,\n metadata=np.asarray(json.dumps(cache_metadata(stage, config), sort_keys=True)),\n **prepared.__dict__,\n )\n log(f\"saved public-derived full-sequence cache: {path}\")\n\n\ndef load_cache(path: Path, stage: str, config: Config) -> SequenceData:\n with np.load(path, allow_pickle=False) as artifact:\n actual = json.loads(str(artifact[\"metadata\"].item()))\n expected = cache_metadata(stage, config)\n if actual != expected:\n raise RuntimeError(f\"sequence cache differs: {actual} != {expected}\")\n values = {\n field: artifact[field].copy()\n for field in SequenceData.__dataclass_fields__\n }\n prepared = SequenceData(**values)\n log(\n f\"loaded full-sequence cache: tokens={len(prepared.sequence_items):,} \"\n f\"targets={prepared.target_count:,}\"\n )\n return prepared\n\n\ndef load_or_prepare(\n public: Path,\n stage: str,\n config: Config,\n cache: Path | None,\n) -> tuple[SequenceData, CandidateSet]:\n candidate_name = (\n \"validation_eligible_candidates.npz\"\n if stage == \"validation\"\n else \"test_eligible_candidates.npz\"\n )\n eligibility = CandidateSet.read(public / \"data\" / candidate_name)\n if cache is not None and cache.is_file():\n prepared = load_cache(cache, stage, config)\n if not np.array_equal(prepared.item_ids, eligibility.item_ids):\n raise RuntimeError(\"cached catalog differs from candidate artifact\")\n return prepared, eligibility\n prepared, built_eligibility = prepare_data(public, stage, config)\n if cache is not None:\n save_cache(cache, prepared, stage, config)\n return prepared, built_eligibility\n\n\ndef build_item_correction(\n prepared: SequenceData,\n config: Config,\n target_starts: np.ndarray,\n target_counts: np.ndarray,\n) -> torch.Tensor:\n \"\"\"Estimate in-batch propensities from the selected training targets.\"\"\"\n raise RuntimeError(\"chronological target partition produced negative counts\")\n if not np.array_equal(\n window_starts + training_counts + validation_counts, total_targets\n ):\n raise RuntimeError(\"chronological target window is not a final contiguous suffix\")\n return window_starts\n\n\ndef chronological_target_partition(\n prepared: SequenceData,\n validation_fraction: float,\n max_train_targets_per_user: int,\n max_validation_targets_per_user: int,\n) -> TargetPartition:\n \"\"\"Select a recent contiguous window with a final out-of-time BCE suffix.\"\"\"\n totals = np.maximum(prepared.lengths - 1, 0).astype(np.int64)\n validation = np.floor(totals * validation_fraction).astype(np.int64)\n validation[(totals >= 2) & (validation == 0)] = 1\n if max_validation_targets_per_user:\n validation = np.minimum(validation, max_validation_targets_per_user)\n\n history_training = totals - validation\n training = history_training.copy()\n if max_train_targets_per_user:\n training = np.minimum(training, max_train_targets_per_user)\n window_starts = retained_target_window_starts(totals, validation, training)\n return TargetPartition(\n window_starts=window_starts,\n training_counts=training,\n validation_counts=validation,\n history_training_counts=history_training,\n )\n\n\ndef load_organic_like_outcomes(\n public: Path, stage: str, item_ids: np.ndarray\n) -> dict[int, dict[int, np.ndarray]]:\n \"\"\"Load mapped organic-like timestamps from public fitting history.\"\"\"\n prefixes = (\"train\",) if stage == \"validation\" else (\"train\", \"val\")\n frames = [\n pd.read_parquet(\n public / \"data\" / f\"{prefix}_likes.parquet\",\n columns=[\"uid\", \"timestamp\", \"item_id\", \"is_organic\"],\n )\n for prefix in prefixes\n ]\n likes = pd.concat(frames, ignore_index=True)\n likes = likes[likes.is_organic == 1].copy()\n mapped = pd.Index(item_ids).get_indexer(\n likes.item_id.to_numpy(dtype=np.int64)\n )\n likes = likes.loc[mapped >= 0, [\"uid\", \"timestamp\"]].copy()\n likes[\"item\"] = mapped[mapped >= 0].astype(np.int32, copy=False)\n likes = likes.sort_values([\"uid\", \"item\", \"timestamp\"], kind=\"stable\")\n nested: dict[int, dict[int, list[int]]] = {}\n for raw_uid, raw_timestamp, raw_item in likes.itertuples(\n index=False, name=None\n ):\n nested.setdefault(int(raw_uid), {}).setdefault(int(raw_item), []).append(\n int(raw_timestamp)\n )\n result = {\n uid: {\n item: np.asarray(timestamps, dtype=np.int64)\n for item, timestamps in items.items()\n }\n for uid, items in nested.items()\n }\n log(\n f\"loaded organic-like outcomes: rows={len(likes):,} \"\n f\"users={len(result):,}\"\n )\n return result\n\n\ndef build_target_weights(\n public: Path,\n stage: str,\n prepared: SequenceData,\n training_counts: np.ndarray,\n config: Config,\n) -> np.ndarray:\n \"\"\"Build nonzero per-token weights without dropping adjacent targets.\"\"\"\n if config.target_weighting == \"uniform\":\n weights = np.ones(len(prepared.sequence_items), dtype=np.float32)\n else:\n weights = np.full(\n len(prepared.sequence_items),\n config.ordinary_target_weight,\n dtype=np.float32,\n )\n organic_strong = (\n (prepared.event_features[:, 0] > 0.5)\n & (prepared.event_features[:, 1] >= 0.8)\n )\n weights[organic_strong] = config.organic_strong_target_weight\n if config.target_weighting == \"discovery-graded\":\n like_outcomes = load_organic_like_outcomes(\n public, stage, prepared.item_ids\n )\n history_end = (\n base.TRAIN_END if stage == \"validation\" else base.VALIDATION_END\n )\n horizon = config.discovery_horizon_days * base.DAY\n grade_three = 0\n for user_row, uid in enumerate(prepared.user_ids):\n start = int(prepared.offsets[user_row])\n end = int(prepared.offsets[user_row + 1])\n if end - start < 2:\n continue\n last_training_target = start + int(training_counts[user_row])\n first_validation_target = last_training_target + 1\n split_timestamp = (\n int(prepared.sequence_timestamps[first_validation_target])\n if first_validation_target < end\n else history_end\n )\n user_likes = like_outcomes.get(int(uid), {})\n seen = {int(prepared.sequence_items[start])}\n for position in range(start + 1, end):\n item = int(prepared.sequence_items[position])\n novel = item not in seen\n seen.add(item)\n if not novel:\n # Repeats remain training examples, but they cannot be\n # eligible discovery outcomes at the corresponding\n # cutoff and therefore keep only the ordinary weight.\n weights[position] = config.ordinary_target_weight\n continue\n timestamps = user_likes.get(item)\n if timestamps is None:\n continue\n timestamp = int(prepared.sequence_timestamps[position])\n upper = min(\n timestamp + horizon,\n split_timestamp\n if position <= last_training_target\n else history_end,\n )\n outcome = int(np.searchsorted(timestamps, timestamp, side=\"left\"))\n if outcome < len(timestamps) and int(timestamps[outcome]) < upper:\n weights[position] = config.organic_like_target_weight\n grade_three += 1\n log(\n f\"grade-3 first-encounter targets={grade_three:,} \"\n f\"horizon_days={config.discovery_horizon_days}\"\n )\n log(\n f\"target weighting={config.target_weighting} \"\n f\"min={weights.min():.3f} mean={weights.mean():.3f} \"\n f\"max={weights.max():.3f}\"\n )\n return weights\n\n\ndef build_dense_block(\n prepared: SequenceData,\n target_weights: np.ndarray,\n users: np.ndarray,\n window_starts: np.ndarray,\n process_counts: np.ndarray,\n score_starts: np.ndarray,\n start: int,\n steps: int,\n config: Config,\n) -> tuple[np.ndarray, ...]:\n \"\"\"Build one causal sequence block and mark the target rows to score.\"\"\"\n batch_size = len(users)\n padding_item = len(prepared.item_ids)\n input_items = np.full((batch_size, steps), padding_item, dtype=np.uint16)\n event_features = np.zeros((batch_size, steps, 3), dtype=np.float16)\n activity_features = np.zeros((batch_size, steps, 4), dtype=np.float16)\n state_ids = np.zeros((batch_size, steps), dtype=np.uint32)\n targets = np.zeros((batch_size, steps), dtype=np.uint16)\n weights = np.ones((batch_size, steps), dtype=np.float32)\n scored = np.zeros((batch_size, steps), dtype=bool)\n histories = np.full(\n (batch_size, steps, config.bptt_steps),\n padding_item,\n dtype=np.uint16,\n )\n processed = np.zeros(batch_size, dtype=np.int32)\n\n lengths = prepared.lengths\n for row, user_row in enumerate(users):\n count = min(steps, max(int(process_counts[row]) - start, 0))\n if not count:\n continue\n offset = int(prepared.offsets[user_row])\n window_start = int(window_starts[row])\n absolute_start = window_start + start\n positions = slice(\n offset + absolute_start,\n offset + absolute_start + count,\n )\n target_positions = slice(\n offset + absolute_start + 1,\n offset + absolute_start + count + 1,\n )\n input_items[row, :count] = prepared.sequence_items[positions]\n event_features[row, :count] = prepared.event_features[positions]\n activity_features[row, :count] = prepared.activity_features[positions]\n state_ids[row, :count] = prepared.state_ids[positions]\n targets[row, :count] = prepared.sequence_items[target_positions]\n weights[row, :count] = target_weights[target_positions]\n score_from = min(count, max(int(score_starts[row]) - start, 0))\n",
"stderr": "",
"aggregated_output": " event_features[row, :count] = prepared.event_features[positions]\n activity_features[row, :count] = prepared.activity_features[positions]\n state_ids[row, :count] = prepared.state_ids[positions]\n targets[row, :count] = prepared.sequence_items[target_positions]\n weights[row, :count] = target_weights[target_positions]\n score_from = min(count, max(int(score_starts[row]) - start, 0))\n scored[row, score_from:count] = True\n processed[row] = count\n user_sequence = prepared.sequence_items[\n offset : offset + int(lengths[user_row])\n ]\n histories[row, :count] = history_windows(\n user_sequence,\n absolute_start,\n count,\n config.bptt_steps,\n padding_item,\n )\n return (\n input_items,\n event_features,\n activity_features,\n state_ids,\n targets,\n scored,\n histories,\n processed,\n weights,\n )\n\n\ndef score_target_microbatches(\n model: base.SequentialTwoTower,\n prepared: SequenceData,\n config: Config,\n sequence_vectors: torch.Tensor,\n flat_states: np.ndarray,\n flat_activity: np.ndarray,\n positive_items: torch.Tensor,\n negatives: torch.Tensor,\n negative_weights: torch.Tensor,\n sample_weights: np.ndarray,\n device: torch.device,\n backward: bool,\n) -> tuple[dict[str, float], torch.Tensor | None]:\n \"\"\"Compute one logical block's exact BCE in bounded-memory slices.\"\"\"\n target_count = len(positive_items)\n total_weight = float(sample_weights.sum())\n if total_weight <= 0:\n raise RuntimeError(\"target block has no positive sample weight\")\n sequence_gradient = torch.zeros_like(sequence_vectors) if backward else None\n statistics = {\n \"loss_sum\": 0.0,\n \"positive_sum\": 0.0,\n \"negative_sum\": 0.0,\n \"pair_correct\": 0.0,\n \"weight_sum\": 0.0,\n }\n for begin in range(0, target_count, config.loss_microbatch_size):\n end = min(begin + config.loss_microbatch_size, target_count)\n selected_states = flat_states[begin:end]\n if backward:\n sequence_input = (\n sequence_vectors[begin:end].detach().requires_grad_(True)\n )\n else:\n sequence_input = sequence_vectors[begin:end]\n liked_vectors, disliked_vectors, _ = pool_explicit_states(\n model, prepared, selected_states, device\n )\n counts = prepared.state_counts[selected_states]\n numeric = np.zeros((end - begin, 6), dtype=np.float32)\n numeric[:, :4] = flat_activity[begin:end]\n numeric[:, 4] = np.log1p(counts[:, 0]) / 6.0\n numeric[:, 5] = np.log1p(counts[:, 1]) / 4.0\n query_vectors = user_queries(\n model,\n sequence_input,\n liked_vectors,\n disliked_vectors,\n torch.from_numpy(numeric).to(device),\n )\n positive_vectors, negative_vectors = candidate_vectors(\n model,\n positive_items[begin:end],\n negatives[begin:end],\n )\n positive_logits = (query_vectors * positive_vectors).sum(dim=1)\n negative_logits = torch.einsum(\n \"bd,bkd->bk\", query_vectors, negative_vectors\n )\n positive_loss = F.softplus(-positive_logits)\n negative_loss = corrected_negative_loss(\n negative_logits, negative_weights[begin:end], config\n )\n selected_weights = torch.from_numpy(\n sample_weights[begin:end].astype(np.float32, copy=False)\n ).to(device)\n loss_sum = (\n selected_weights * (positive_loss + negative_loss) / 2.0\n ).sum()\n if not bool(torch.isfinite(loss_sum)):\n raise RuntimeError(\"dense next-item training produced non-finite loss\")\n if backward:\n (loss_sum / total_weight).backward()\n if sequence_input.grad is None:\n raise RuntimeError(\"sequence microbatch did not produce a gradient\")\n sequence_gradient[begin:end] = sequence_input.grad\n statistics[\"loss_sum\"] += float(loss_sum.detach().cpu())\n statistics[\"positive_sum\"] += float(positive_logits.detach().sum().cpu())\n statistics[\"negative_sum\"] += float(negative_logits.detach().sum().cpu())\n statistics[\"pair_correct\"] += float(\n (positive_logits.detach()[:, None] > negative_logits.detach()).sum().cpu()\n )\n statistics[\"weight_sum\"] += float(selected_weights.sum().detach().cpu())\n return statistics, sequence_gradient\n\n\ndef release_accelerator_cache(device: torch.device) -> None:\n if device.type == \"mps\":\n torch.mps.empty_cache()\n\n\ndef run_training_epoch(\n model: base.SequentialTwoTower,\n optimizer: torch.optim.Optimizer,\n prepared: SequenceData,\n config: Config,\n target_weights: np.ndarray,\n window_starts: np.ndarray,\n target_counts: np.ndarray,\n item_correction: torch.Tensor,\n known_positive_packed: np.ndarray | None,\n generator: torch.Generator,\n device: torch.device,\n epoch: int,\n) -> dict[str, float]:\n model.train()\n train_users = np.flatnonzero(target_counts > 0)\n sorted_users = train_users[\n np.argsort(target_counts[train_users], kind=\"stable\")\n ]\n groups = [\n sorted_users[start : start + config.user_batch_size]\n for start in range(0, len(sorted_users), config.user_batch_size)\n ]\n group_order = torch.randperm(len(groups), generator=generator).numpy()\n negative_count = config.uniform_negatives + config.inbatch_negatives\n total_loss = 0.0\n total_targets = 0\n total_weight = 0.0\n total_positive = 0.0\n total_negative = 0.0\n total_pair_correct = 0.0\n batches = 0\n started = time.monotonic()\n\n for group_number in group_order:\n users = groups[int(group_number)]\n group_counts = target_counts[users]\n hidden = torch.zeros(\n 1, len(users), config.embedding_dim, device=device\n )\n max_targets = int(group_counts.max())\n for start in range(0, max_targets, config.bptt_steps):\n steps = min(config.bptt_steps, max_targets - start)\n block = build_dense_block(\n prepared,\n target_weights,\n users,\n window_starts[users],\n group_counts,\n np.zeros(len(users), dtype=np.int64),\n start,\n steps,\n config,\n )\n sequence_outputs, _ = event_sequence(\n model,\n torch.from_numpy(block[0].astype(np.int64)).to(device),\n torch.from_numpy(block[1].astype(np.float32)).to(device),\n hidden,\n )\n hidden = torch.stack(\n [\n sequence_outputs[row, count - 1]\n if count\n else hidden[0, row]\n for row, count in enumerate(block[7])\n ]\n ).unsqueeze(0).detach()\n flat_scored = block[5].reshape(-1)\n if not np.any(flat_scored):\n continue\n mask = torch.from_numpy(flat_scored).to(device)\n sequence_vectors = sequence_outputs.reshape(\n -1, config.embedding_dim\n )[mask]\n flat_states = block[3].reshape(-1)[flat_scored]\n positive_items = torch.from_numpy(\n block[4].reshape(-1)[flat_scored].astype(np.int64)\n ).to(device)\n known_history = torch.from_numpy(\n block[6].reshape(-1, config.bptt_steps)[flat_scored].astype(\n np.int64\n )\n ).to(device)\n liked_items = torch.from_numpy(\n prepared.liked_states[flat_states].astype(np.int64)\n ).to(device)\n flat_user_rows = np.broadcast_to(\n users[:, None], block[5].shape\n ).reshape(-1)[flat_scored]\n negatives, negative_weights = draw_training_negatives(\n positive_items,\n known_history,\n liked_items,\n flat_user_rows,\n known_positive_packed,\n prepared,\n config,\n item_correction,\n generator,\n device,\n )\n optimizer.zero_grad(set_to_none=True)\n statistics, sequence_gradient = score_target_microbatches(\n model,\n prepared,\n config,\n sequence_vectors,\n flat_states,\n block[2].reshape(-1, 4)[flat_scored].astype(np.float32),\n positive_items,\n negatives,\n negative_weights,\n block[8].reshape(-1)[flat_scored],\n device,\n backward=True,\n )\n if sequence_gradient is None:\n raise RuntimeError(\"training block omitted its sequence gradient\")\n sequence_vectors.backward(sequence_gradient)\n nn.utils.clip_grad_norm_(model.parameters(), config.gradient_clip)\n optimizer.step()\n\n count = len(positive_items)\n total_loss += statistics[\"loss_sum\"]\n total_targets += count\n total_weight += statistics[\"weight_sum\"]\n total_positive += statistics[\"positive_sum\"]\n total_negative += statistics[\"negative_sum\"]\n total_pair_correct += statistics[\"pair_correct\"]\n batches += 1\n if batches % 25 == 0:\n release_accelerator_cache(device)\n if batches % 250 == 0:\n batch_bce = statistics[\"loss_sum\"] / max(\n statistics[\"weight_sum\"], 1.0e-12\n )\n running_bce = total_loss / max(total_weight, 1.0e-12)\n log(\n f\"epoch={epoch}/{config.epochs} batches={batches:,} \"\n f\"targets={total_targets:,} \"\n f\"batch_bce={batch_bce:.6f} \"\n f\"running_bce={running_bce:.6f} \"\n f\"elapsed={time.monotonic() - started:.1f}s\"\n )\n\n return {\n \"epoch\": float(epoch),\n \"loss\": total_loss / max(total_weight, 1.0e-12),\n \"positive_logit\": total_positive / max(total_targets, 1),\n \"negative_logit\": total_negative / max(total_targets * negative_count, 1),\n \"pair_accuracy\": total_pair_correct\n / max(total_targets * negative_count, 1),\n \"targets\": float(total_targets),\n \"target_weight_sum\": total_weight,\n \"batches\": float(batches),\n \"seconds\": time.monotonic() - started,\n }\n\n\[email protected]_grad()\ndef evaluate_validation_bce(\n model: base.SequentialTwoTower,\n prepared: SequenceData,\n config: Config,\n target_weights: np.ndarray,\n window_starts: np.ndarray,\n training_counts: np.ndarray,\n validation_counts: np.ndarray,\n item_correction: torch.Tensor,\n known_positive_packed: np.ndarray | None,\n device: torch.device,\n) -> dict[str, float]:\n \"\"\"Measure deterministic BCE on each user's chronological target suffix.\"\"\"\n model.eval()\n totals = training_counts + validation_counts\n users_with_validation = np.flatnonzero(validation_counts > 0)\n ordered = users_with_validation[\n np.argsort(totals[users_with_validation], kind=\"stable\")\n ]\n groups = [\n ordered[start : start + config.user_batch_size]\n for start in range(0, len(ordered), config.user_batch_size)\n ]\n generator = torch.Generator(device=\"cpu\").manual_seed(config.seed + 20_000)\n total_loss = 0.0\n total_targets = 0\n total_weight = 0.0\n started = time.monotonic()\n\n for group_number, users in enumerate(groups, start=1):\n group_totals = totals[users]\n group_starts = training_counts[users]\n hidden = torch.zeros(\n 1, len(users), config.embedding_dim, device=device\n )\n max_targets = int(group_totals.max())\n for start in range(0, max_targets, config.bptt_steps):\n steps = min(config.bptt_steps, max_targets - start)\n block = build_dense_block(\n prepared,\n target_weights,\n users,\n window_starts[users],\n group_totals,\n group_starts,\n start,\n steps,\n config,\n )\n sequence_outputs, _ = event_sequence(\n model,\n torch.from_numpy(block[0].astype(np.int64)).to(device),\n torch.from_numpy(block[1].astype(np.float32)).to(device),\n hidden,\n )\n hidden = torch.stack(\n [\n sequence_outputs[row, count - 1]\n if count\n else hidden[0, row]\n for row, count in enumerate(block[7])\n ]\n ).unsqueeze(0)\n flat_scored = block[5].reshape(-1)\n if not np.any(flat_scored):\n continue\n mask = torch.from_numpy(flat_scored).to(device)\n sequence_vectors = sequence_outputs.reshape(\n -1, config.embedding_dim\n )[mask]\n flat_states = block[3].reshape(-1)[flat_scored]\n positive_items = torch.from_numpy(\n block[4].reshape(-1)[flat_scored].astype(np.int64)\n ).to(device)\n known_history = torch.from_numpy(\n block[6].reshape(-1, config.bptt_steps)[flat_scored].astype(\n np.int64\n )\n ).to(device)\n liked_items = torch.from_numpy(\n prepared.liked_states[flat_states].astype(np.int64)\n ).to(device)\n flat_user_rows = np.broadcast_to(\n users[:, None], block[5].shape\n ).reshape(-1)[flat_scored]\n negatives, negative_weights = draw_training_negatives(\n positive_items,\n known_history,\n liked_items,\n flat_user_rows,\n known_positive_packed,\n prepared,\n config,\n item_correction,\n generator,\n device,\n )\n statistics, _ = score_target_microbatches(\n model,\n prepared,\n config,\n sequence_vectors,\n flat_states,\n block[2].reshape(-1, 4)[flat_scored].astype(np.float32),\n positive_items,\n negatives,\n negative_weights,\n block[8].reshape(-1)[flat_scored],\n device,\n backward=False,\n )\n total_loss += statistics[\"loss_sum\"]\n total_targets += len(positive_items)\n total_weight += statistics[\"weight_sum\"]\n if group_number % 25 == 0:\n release_accelerator_cache(device)\n return {\n \"validation_loss\": total_loss / max(total_weight, 1.0e-12),\n \"validation_targets\": float(total_targets),\n \"validation_target_weight_sum\": total_weight,\n \"validation_seconds\": time.monotonic() - started,\n }\n\n\ndef snapshot_model(model: nn.Module) -> dict[str, torch.Tensor]:\n return {\n name: value.detach().cpu().clone()\n for name, value in model.state_dict().items()\n }\n\n\ndef save_epoch_checkpoint(\n checkpoint_dir: Path,\n epoch: int,\n model_state: dict[str, torch.Tensor],\n config: Config,\n diagnostics: dict[str, float],\n optimizer: torch.optim.Optimizer,\n generator: torch.Generator,\n) -> Path:\n checkpoint_dir.mkdir(parents=True, exist_ok=True)\n path = checkpoint_dir / f\"epoch-{epoch:02d}.pt\"\n torch.save(\n {\n \"epoch\": epoch,\n \"config\": asdict(config),\n \"diagnostics\": diagnostics,\n \"model_state\": model_state,\n \"optimizer_state\": optimizer.state_dict(),\n \"generator_state\": generator.get_state(),\n },\n path,\n destination[row, : len(retained)] = retained.astype(np.uint16)\n return liked, disliked, like_counts, dislike_counts\n\n\nclass SequentialTwoTower(nn.Module):\n \"\"\"Order-sensitive user tower and independent catalog-aware item tower.\"\"\"\n\n def __init__(self, item_features: np.ndarray, config: TowerConfig) -> None:\n super().__init__()\n n_items, n_item_features = item_features.shape\n dim = config.embedding_dim\n self.n_items = n_items\n self.context_embedding = nn.Embedding(\n n_items + 1, dim, padding_idx=n_items\n )\n self.event_projection = nn.Sequential(\n nn.Linear(dim + 3, config.hidden_dim),\n nn.GELU(),\n )\n self.sequence_encoder = nn.GRU(\n input_size=config.hidden_dim,\n hidden_size=dim,\n num_layers=1,\n batch_first=True,\n )\n self.user_tower = nn.Sequential(\n nn.Linear(dim * 3 + 6, config.user_hidden_1),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.user_hidden_1, config.user_hidden_2),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.user_hidden_2, dim),\n )\n self.candidate_embedding = nn.Embedding(n_items, dim)\n self.register_buffer(\n \"item_features\", torch.from_numpy(item_features), persistent=True\n )\n self.item_tower = nn.Sequential(\n nn.Linear(dim + n_item_features, config.item_hidden_1),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.item_hidden_1, config.item_hidden_2),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.item_hidden_2, dim),\n )\n nn.init.normal_(self.context_embedding.weight, std=0.02)\n nn.init.normal_(self.candidate_embedding.weight, std=0.02)\n with torch.no_grad():\n self.context_embedding.weight[n_items].zero_()\n\n def _mean_context_channel(self, items: torch.Tensor) -> torch.Tensor:\n mask = items.ne(self.n_items)\n vectors = self.context_embedding(items)\n return (vectors * mask.unsqueeze(-1)).sum(dim=1) / mask.sum(\n dim=1, keepdim=True\n ).clamp_min(1)\n\n def encode_items(self, items: torch.Tensor) -> torch.Tensor:\n inputs = torch.cat(\n [self.candidate_embedding(items), self.item_features[items]], dim=1\n )\n return self.item_tower(inputs)\n\n\ndef invalid_negative_mask(\n dislike_count: int,\n) -> np.ndarray:\n start = int(\n np.searchsorted(\n timestamps,\n history_end - base.PROFILE_LOOKBACK_DAYS * base.DAY,\n side=\"left\",\n )\n )\n result = np.zeros(6, dtype=np.float32)\n recent_items = items[start:]\n if len(recent_items):\n result[0] = np.log1p(len(recent_items)) / 8.0\n result[1] = np.log1p(len(np.unique(recent_items))) / 6.0\n result[2] = float(np.mean(event_features[start:, 0]))\n result[3] = float(np.mean(event_features[start:, 1]))\n result[4] = np.log1p(like_count) / 6.0\n result[5] = np.log1p(dislike_count) / 4.0\n return result\n\n\ndef prepare_data(\n public: Path, stage: str, config: Config\n) -> tuple[SequenceData, CandidateSet]:\n data = public / \"data\"\n catalog = pd.read_parquet(data / \"candidate_catalog.parquet\").sort_values(\n \"popularity_rank\"\n )\n item_ids = catalog.item_id.to_numpy(dtype=np.int64, copy=True)\n item_index = pd.Index(item_ids)\n n_items = len(item_ids)\n if n_items >= np.iinfo(np.uint16).max:\n raise RuntimeError(\"uint16 item storage is too small for this catalog\")\n\n if stage == \"validation\":\n prefixes = (\"train\",)\n history_end = base.TRAIN_END\n candidate_name = \"validation_eligible_candidates.npz\"\n state_name = \"validation_preference_state.parquet\"\n else:\n prefixes = (\"train\", \"val\")\n history_end = base.VALIDATION_END\n candidate_name = \"test_eligible_candidates.npz\"\n state_name = \"test_preference_state.parquet\"\n eligibility = CandidateSet.read(data / candidate_name)\n if not np.array_equal(eligibility.item_ids, item_ids):\n raise RuntimeError(\"catalog and candidate artifact item order differ\")\n\n explicit_by_user = base.load_explicit_events(data, prefixes, item_index)\n inference_likes, inference_dislikes, like_counts, dislike_counts = (\n base.load_inference_state(\n data / state_name,\n item_index,\n eligibility.uids,\n n_items,\n )\n )\n inference_rows = {int(uid): row for row, uid in enumerate(eligibility.uids)}\n inference_user_features = np.zeros((len(eligibility.uids), 6), dtype=np.float32)\n inference_sequence_rows = np.full(len(eligibility.uids), -1, dtype=np.int32)\n\n # Snapshot zero is the shared empty explicit state.\n liked_states = [np.full(base.LIKE_ITEMS, n_items, dtype=np.uint16)]\n disliked_states = [\n np.full(base.DISLIKE_ITEMS, n_items, dtype=np.uint16)\n ]\n state_counts: list[tuple[int, int]] = [(0, 0)]\n user_ids: list[int] = []\n offsets = [0]\n timestamp_parts: list[np.ndarray] = []\n item_parts: list[np.ndarray] = []\n event_parts: list[np.ndarray] = []\n activity_parts: list[np.ndarray] = []\n state_parts: list[np.ndarray] = []\n\n listen_paths = [data / f\"{prefix}_listens.parquet\" for prefix in prefixes]\n started = time.monotonic()\n for seen_users, (uid, listens) in enumerate(\n base.merge_listen_streams(listen_paths, item_index), start=1\n ):\n parts = build_sequence_parts(\n listens,\n explicit_by_user.get(uid, base.empty_explicit_events()),\n n_items,\n liked_states,\n disliked_states,\n state_counts,\n )\n if parts is None:\n continue\n sequence_row = len(user_ids)\n user_ids.append(uid)\n timestamp_parts.append(parts[0].astype(np.uint32, copy=False))\n item_parts.append(parts[1])\n event_parts.append(parts[2])\n activity_parts.append(parts[3])\n state_parts.append(parts[4])\n offsets.append(offsets[-1] + len(parts[1]))\n inference_row = inference_rows.get(uid)\n if inference_row is not None:\n inference_sequence_rows[inference_row] = sequence_row\n inference_user_features[inference_row] = final_user_features(\n parts[0],\n parts[1],\n parts[2],\n history_end,\n int(like_counts[inference_row]),\n int(dislike_counts[inference_row]),\n )\n if seen_users % 1_000 == 0:\n log(\n f\"prepared users={seen_users:,} tokens={offsets[-1]:,} \"\n f\"states={len(liked_states):,} \"\n f\"elapsed={time.monotonic() - started:.1f}s\"\n )\n\n missing = np.flatnonzero(inference_sequence_rows < 0)\n if len(missing):\n inference_user_features[missing, 4] = np.log1p(like_counts[missing]) / 6.0\n inference_user_features[missing, 5] = (\n np.log1p(dislike_counts[missing]) / 4.0\n )\n log(f\"target users without qualifying catalog sequence={len(missing):,}\")\n\n prepared = SequenceData(\n item_ids=item_ids,\n item_features=base.build_item_features(catalog),\n user_ids=np.asarray(user_ids, dtype=np.int64),\n offsets=np.asarray(offsets, dtype=np.int64),\n sequence_timestamps=np.concatenate(timestamp_parts),\n sequence_items=np.concatenate(item_parts),\n event_features=np.concatenate(event_parts),\n activity_features=np.concatenate(activity_parts),\n state_ids=np.concatenate(state_parts),\n liked_states=np.stack(liked_states),\n disliked_states=np.stack(disliked_states),\n state_counts=np.asarray(state_counts, dtype=np.int32),\n inference_user_rows=inference_sequence_rows,\n inference_liked_items=inference_likes,\n inference_disliked_items=inference_dislikes,\n inference_user_features=inference_user_features,\n )\n log(\n f\"prepared stage={stage} sequence_users={len(prepared.user_ids):,} \"\n f\"tokens={len(prepared.sequence_items):,} \"\n f\"next_item_targets={prepared.target_count:,} \"\n f\"states={len(prepared.liked_states):,}\"\n )\n return prepared, eligibility\n\n\ndef cache_metadata(stage: str, config: Config) -> dict[str, object]:\n return {\n \"cache_version\": CACHE_VERSION,\n \"stage\": stage,\n \"positive_listen_threshold\": base.LISTEN_POSITIVE_THRESHOLD,\n \"profile_lookback_days\": base.PROFILE_LOOKBACK_DAYS,\n \"event_time_feature\": \"exp-negative-gap-over-60-days\",\n }\n\n\ndef save_cache(\n path: Path, prepared: SequenceData, stage: str, config: Config\n) -> None:\n path.parent.mkdir(parents=True, exist_ok=True)\n np.savez(\n path,\n metadata=np.asarray(json.dumps(cache_metadata(stage, config), sort_keys=True)),\n **prepared.__dict__,\n )\n log(f\"saved public-derived full-sequence cache: {path}\")\n\n\ndef load_cache(path: Path, stage: str, config: Config) -> SequenceData:\n with np.load(path, allow_pickle=False) as artifact:\n actual = json.loads(str(artifact[\"metadata\"].item()))\n expected = cache_metadata(stage, config)\n if actual != expected:\n raise RuntimeError(f\"sequence cache differs: {actual} != {expected}\")\n values = {\n field: artifact[field].copy()\n for field in SequenceData.__dataclass_fields__\n }\n prepared = SequenceData(**values)\n log(\n f\"loaded full-sequence cache: tokens={len(prepared.sequence_items):,} \"\n f\"targets={prepared.target_count:,}\"\n )\n return prepared\n\n\ndef load_or_prepare(\n public: Path,\n stage: str,\n config: Config,\n cache: Path | None,\n) -> tuple[SequenceData, CandidateSet]:\n candidate_name = (\n \"validation_eligible_candidates.npz\"\n if stage == \"validation\"\n else \"test_eligible_candidates.npz\"\n )\n eligibility = CandidateSet.read(public / \"data\" / candidate_name)\n if cache is not None and cache.is_file():\n prepared = load_cache(cache, stage, config)\n if not np.array_equal(prepared.item_ids, eligibility.item_ids):\n raise RuntimeError(\"cached catalog differs from candidate artifact\")\n return prepared, eligibility\n prepared, built_eligibility = prepare_data(public, stage, config)\n if cache is not None:\n save_cache(cache, prepared, stage, config)\n return prepared, built_eligibility\n\n\ndef build_item_correction(\n prepared: SequenceData,\n config: Config,\n target_starts: np.ndarray,\n target_counts: np.ndarray,\n) -> torch.Tensor:\n \"\"\"Estimate in-batch propensities from the selected training targets.\"\"\"\n raise RuntimeError(\"chronological target partition produced negative counts\")\n if not np.array_equal(\n window_starts + training_counts + validation_counts, total_targets\n ):\n raise RuntimeError(\"chronological target window is not a final contiguous suffix\")\n return window_starts\n\n\ndef chronological_target_partition(\n prepared: SequenceData,\n validation_fraction: float,\n max_train_targets_per_user: int,\n max_validation_targets_per_user: int,\n) -> TargetPartition:\n \"\"\"Select a recent contiguous window with a final out-of-time BCE suffix.\"\"\"\n totals = np.maximum(prepared.lengths - 1, 0).astype(np.int64)\n validation = np.floor(totals * validation_fraction).astype(np.int64)\n validation[(totals >= 2) & (validation == 0)] = 1\n if max_validation_targets_per_user:\n validation = np.minimum(validation, max_validation_targets_per_user)\n\n history_training = totals - validation\n training = history_training.copy()\n if max_train_targets_per_user:\n training = np.minimum(training, max_train_targets_per_user)\n window_starts = retained_target_window_starts(totals, validation, training)\n return TargetPartition(\n window_starts=window_starts,\n training_counts=training,\n validation_counts=validation,\n history_training_counts=history_training,\n )\n\n\ndef load_organic_like_outcomes(\n public: Path, stage: str, item_ids: np.ndarray\n) -> dict[int, dict[int, np.ndarray]]:\n \"\"\"Load mapped organic-like timestamps from public fitting history.\"\"\"\n prefixes = (\"train\",) if stage == \"validation\" else (\"train\", \"val\")\n frames = [\n pd.read_parquet(\n public / \"data\" / f\"{prefix}_likes.parquet\",\n columns=[\"uid\", \"timestamp\", \"item_id\", \"is_organic\"],\n )\n for prefix in prefixes\n ]\n likes = pd.concat(frames, ignore_index=True)\n likes = likes[likes.is_organic == 1].copy()\n mapped = pd.Index(item_ids).get_indexer(\n likes.item_id.to_numpy(dtype=np.int64)\n )\n likes = likes.loc[mapped >= 0, [\"uid\", \"timestamp\"]].copy()\n likes[\"item\"] = mapped[mapped >= 0].astype(np.int32, copy=False)\n likes = likes.sort_values([\"uid\", \"item\", \"timestamp\"], kind=\"stable\")\n nested: dict[int, dict[int, list[int]]] = {}\n for raw_uid, raw_timestamp, raw_item in likes.itertuples(\n index=False, name=None\n ):\n nested.setdefault(int(raw_uid), {}).setdefault(int(raw_item), []).append(\n int(raw_timestamp)\n )\n result = {\n uid: {\n item: np.asarray(timestamps, dtype=np.int64)\n for item, timestamps in items.items()\n }\n for uid, items in nested.items()\n }\n log(\n f\"loaded organic-like outcomes: rows={len(likes):,} \"\n f\"users={len(result):,}\"\n )\n return result\n\n\ndef build_target_weights(\n public: Path,\n stage: str,\n prepared: SequenceData,\n training_counts: np.ndarray,\n config: Config,\n) -> np.ndarray:\n \"\"\"Build nonzero per-token weights without dropping adjacent targets.\"\"\"\n if config.target_weighting == \"uniform\":\n weights = np.ones(len(prepared.sequence_items), dtype=np.float32)\n else:\n weights = np.full(\n len(prepared.sequence_items),\n config.ordinary_target_weight,\n dtype=np.float32,\n )\n organic_strong = (\n (prepared.event_features[:, 0] > 0.5)\n & (prepared.event_features[:, 1] >= 0.8)\n )\n weights[organic_strong] = config.organic_strong_target_weight\n if config.target_weighting == \"discovery-graded\":\n like_outcomes = load_organic_like_outcomes(\n public, stage, prepared.item_ids\n )\n history_end = (\n base.TRAIN_END if stage == \"validation\" else base.VALIDATION_END\n )\n horizon = config.discovery_horizon_days * base.DAY\n grade_three = 0\n for user_row, uid in enumerate(prepared.user_ids):\n start = int(prepared.offsets[user_row])\n end = int(prepared.offsets[user_row + 1])\n if end - start < 2:\n continue\n last_training_target = start + int(training_counts[user_row])\n first_validation_target = last_training_target + 1\n split_timestamp = (\n int(prepared.sequence_timestamps[first_validation_target])\n if first_validation_target < end\n else history_end\n )\n user_likes = like_outcomes.get(int(uid), {})\n seen = {int(prepared.sequence_items[start])}\n for position in range(start + 1, end):\n item = int(prepared.sequence_items[position])\n novel = item not in seen\n seen.add(item)\n if not novel:\n # Repeats remain training examples, but they cannot be\n # eligible discovery outcomes at the corresponding\n # cutoff and therefore keep only the ordinary weight.\n weights[position] = config.ordinary_target_weight\n continue\n timestamps = user_likes.get(item)\n if timestamps is None:\n continue\n timestamp = int(prepared.sequence_timestamps[position])\n upper = min(\n timestamp + horizon,\n split_timestamp\n if position <= last_training_target\n else history_end,\n )\n outcome = int(np.searchsorted(timestamps, timestamp, side=\"left\"))\n if outcome < len(timestamps) and int(timestamps[outcome]) < upper:\n weights[position] = config.organic_like_target_weight\n grade_three += 1\n log(\n f\"grade-3 first-encounter targets={grade_three:,} \"\n f\"horizon_days={config.discovery_horizon_days}\"\n )\n log(\n f\"target weighting={config.target_weighting} \"\n f\"min={weights.min():.3f} mean={weights.mean():.3f} \"\n f\"max={weights.max():.3f}\"\n )\n return weights\n\n\ndef build_dense_block(\n prepared: SequenceData,\n target_weights: np.ndarray,\n users: np.ndarray,\n window_starts: np.ndarray,\n process_counts: np.ndarray,\n score_starts: np.ndarray,\n start: int,\n steps: int,\n config: Config,\n) -> tuple[np.ndarray, ...]:\n \"\"\"Build one causal sequence block and mark the target rows to score.\"\"\"\n batch_size = len(users)\n padding_item = len(prepared.item_ids)\n input_items = np.full((batch_size, steps), padding_item, dtype=np.uint16)\n event_features = np.zeros((batch_size, steps, 3), dtype=np.float16)\n activity_features = np.zeros((batch_size, steps, 4), dtype=np.float16)\n state_ids = np.zeros((batch_size, steps), dtype=np.uint32)\n targets = np.zeros((batch_size, steps), dtype=np.uint16)\n weights = np.ones((batch_size, steps), dtype=np.float32)\n scored = np.zeros((batch_size, steps), dtype=bool)\n histories = np.full(\n (batch_size, steps, config.bptt_steps),\n padding_item,\n dtype=np.uint16,\n )\n processed = np.zeros(batch_size, dtype=np.int32)\n\n lengths = prepared.lengths\n for row, user_row in enumerate(users):\n count = min(steps, max(int(process_counts[row]) - start, 0))\n if not count:\n continue\n offset = int(prepared.offsets[user_row])\n window_start = int(window_starts[row])\n absolute_start = window_start + start\n positions = slice(\n offset + absolute_start,\n offset + absolute_start + count,\n )\n target_positions = slice(\n offset + absolute_start + 1,\n offset + absolute_start + count + 1,\n )\n input_items[row, :count] = prepared.sequence_items[positions]\n event_features[row, :count] = prepared.event_features[positions]\n activity_features[row, :count] = prepared.activity_features[positions]\n state_ids[row, :count] = prepared.state_ids[positions]\n targets[row, :count] = prepared.sequence_items[target_positions]\n weights[row, :count] = target_weights[target_positions]\n score_from = min(count, max(int(score_starts[row]) - start, 0))\n",
"exit_code": 0,
"duration": {
"secs": 0,
"nanos": 4460
},
"formatted_output": " event_features[row, :count] = prepared.event_features[positions]\n activity_features[row, :count] = prepared.activity_features[positions]\n state_ids[row, :count] = prepared.state_ids[positions]\n targets[row, :count] = prepared.sequence_items[target_positions]\n weights[row, :count] = target_weights[target_positions]\n score_from = min(count, max(int(score_starts[row]) - start, 0))\n scored[row, score_from:count] = True\n processed[row] = count\n user_sequence = prepared.sequence_items[\n offset : offset + int(lengths[user_row])\n ]\n histories[row, :count] = history_windows(\n user_sequence,\n absolute_start,\n count,\n config.bptt_steps,\n padding_item,\n )\n return (\n input_items,\n event_features,\n activity_features,\n state_ids,\n targets,\n scored,\n histories,\n processed,\n weights,\n )\n\n\ndef score_target_microbatches(\n model: base.SequentialTwoTower,\n prepared: SequenceData,\n config: Config,\n sequence_vectors: torch.Tensor,\n flat_states: np.ndarray,\n flat_activity: np.ndarray,\n positive_items: torch.Tensor,\n negatives: torch.Tensor,\n negative_weights: torch.Tensor,\n sample_weights: np.ndarray,\n device: torch.device,\n backward: bool,\n) -> tuple[dict[str, float], torch.Tensor | None]:\n \"\"\"Compute one logical block's exact BCE in bounded-memory slices.\"\"\"\n target_count = len(positive_items)\n total_weight = float(sample_weights.sum())\n if total_weight <= 0:\n raise RuntimeError(\"target block has no positive sample weight\")\n sequence_gradient = torch.zeros_like(sequence_vectors) if backward else None\n statistics = {\n \"loss_sum\": 0.0,\n \"positive_sum\": 0.0,\n \"negative_sum\": 0.0,\n \"pair_correct\": 0.0,\n \"weight_sum\": 0.0,\n }\n for begin in range(0, target_count, config.loss_microbatch_size):\n end = min(begin + config.loss_microbatch_size, target_count)\n selected_states = flat_states[begin:end]\n if backward:\n sequence_input = (\n sequence_vectors[begin:end].detach().requires_grad_(True)\n )\n else:\n sequence_input = sequence_vectors[begin:end]\n liked_vectors, disliked_vectors, _ = pool_explicit_states(\n model, prepared, selected_states, device\n )\n counts = prepared.state_counts[selected_states]\n numeric = np.zeros((end - begin, 6), dtype=np.float32)\n numeric[:, :4] = flat_activity[begin:end]\n numeric[:, 4] = np.log1p(counts[:, 0]) / 6.0\n numeric[:, 5] = np.log1p(counts[:, 1]) / 4.0\n query_vectors = user_queries(\n model,\n sequence_input,\n liked_vectors,\n disliked_vectors,\n torch.from_numpy(numeric).to(device),\n )\n positive_vectors, negative_vectors = candidate_vectors(\n model,\n positive_items[begin:end],\n negatives[begin:end],\n )\n positive_logits = (query_vectors * positive_vectors).sum(dim=1)\n negative_logits = torch.einsum(\n \"bd,bkd->bk\", query_vectors, negative_vectors\n )\n positive_loss = F.softplus(-positive_logits)\n negative_loss = corrected_negative_loss(\n negative_logits, negative_weights[begin:end], config\n )\n selected_weights = torch.from_numpy(\n sample_weights[begin:end].astype(np.float32, copy=False)\n ).to(device)\n loss_sum = (\n selected_weights * (positive_loss + negative_loss) / 2.0\n ).sum()\n if not bool(torch.isfinite(loss_sum)):\n raise RuntimeError(\"dense next-item training produced non-finite loss\")\n if backward:\n (loss_sum / total_weight).backward()\n if sequence_input.grad is None:\n raise RuntimeError(\"sequence microbatch did not produce a gradient\")\n sequence_gradient[begin:end] = sequence_input.grad\n statistics[\"loss_sum\"] += float(loss_sum.detach().cpu())\n statistics[\"positive_sum\"] += float(positive_logits.detach().sum().cpu())\n statistics[\"negative_sum\"] += float(negative_logits.detach().sum().cpu())\n statistics[\"pair_correct\"] += float(\n (positive_logits.detach()[:, None] > negative_logits.detach()).sum().cpu()\n )\n statistics[\"weight_sum\"] += float(selected_weights.sum().detach().cpu())\n return statistics, sequence_gradient\n\n\ndef release_accelerator_cache(device: torch.device) -> None:\n if device.type == \"mps\":\n torch.mps.empty_cache()\n\n\ndef run_training_epoch(\n model: base.SequentialTwoTower,\n optimizer: torch.optim.Optimizer,\n prepared: SequenceData,\n config: Config,\n target_weights: np.ndarray,\n window_starts: np.ndarray,\n target_counts: np.ndarray,\n item_correction: torch.Tensor,\n known_positive_packed: np.ndarray | None,\n generator: torch.Generator,\n device: torch.device,\n epoch: int,\n) -> dict[str, float]:\n model.train()\n train_users = np.flatnonzero(target_counts > 0)\n sorted_users = train_users[\n np.argsort(target_counts[train_users], kind=\"stable\")\n ]\n groups = [\n sorted_users[start : start + config.user_batch_size]\n for start in range(0, len(sorted_users), config.user_batch_size)\n ]\n group_order = torch.randperm(len(groups), generator=generator).numpy()\n negative_count = config.uniform_negatives + config.inbatch_negatives\n total_loss = 0.0\n total_targets = 0\n total_weight = 0.0\n total_positive = 0.0\n total_negative = 0.0\n total_pair_correct = 0.0\n batches = 0\n started = time.monotonic()\n\n for group_number in group_order:\n users = groups[int(group_number)]\n group_counts = target_counts[users]\n hidden = torch.zeros(\n 1, len(users), config.embedding_dim, device=device\n )\n max_targets = int(group_counts.max())\n for start in range(0, max_targets, config.bptt_steps):\n steps = min(config.bptt_steps, max_targets - start)\n block = build_dense_block(\n prepared,\n target_weights,\n users,\n window_starts[users],\n group_counts,\n np.zeros(len(users), dtype=np.int64),\n start,\n steps,\n config,\n )\n sequence_outputs, _ = event_sequence(\n model,\n torch.from_numpy(block[0].astype(np.int64)).to(device),\n torch.from_numpy(block[1].astype(np.float32)).to(device),\n hidden,\n )\n hidden = torch.stack(\n [\n sequence_outputs[row, count - 1]\n if count\n else hidden[0, row]\n for row, count in enumerate(block[7])\n ]\n ).unsqueeze(0).detach()\n flat_scored = block[5].reshape(-1)\n if not np.any(flat_scored):\n continue\n mask = torch.from_numpy(flat_scored).to(device)\n sequence_vectors = sequence_outputs.reshape(\n -1, config.embedding_dim\n )[mask]\n flat_states = block[3].reshape(-1)[flat_scored]\n positive_items = torch.from_numpy(\n block[4].reshape(-1)[flat_scored].astype(np.int64)\n ).to(device)\n known_history = torch.from_numpy(\n block[6].reshape(-1, config.bptt_steps)[flat_scored].astype(\n np.int64\n )\n ).to(device)\n liked_items = torch.from_numpy(\n prepared.liked_states[flat_states].astype(np.int64)\n ).to(device)\n flat_user_rows = np.broadcast_to(\n users[:, None], block[5].shape\n ).reshape(-1)[flat_scored]\n negatives, negative_weights = draw_training_negatives(\n positive_items,\n known_history,\n liked_items,\n flat_user_rows,\n known_positive_packed,\n prepared,\n config,\n item_correction,\n generator,\n device,\n )\n optimizer.zero_grad(set_to_none=True)\n statistics, sequence_gradient = score_target_microbatches(\n model,\n prepared,\n config,\n sequence_vectors,\n flat_states,\n block[2].reshape(-1, 4)[flat_scored].astype(np.float32),\n positive_items,\n negatives,\n negative_weights,\n block[8].reshape(-1)[flat_scored],\n device,\n backward=True,\n )\n if sequence_gradient is None:\n raise RuntimeError(\"training block omitted its sequence gradient\")\n sequence_vectors.backward(sequence_gradient)\n nn.utils.clip_grad_norm_(model.parameters(), config.gradient_clip)\n optimizer.step()\n\n count = len(positive_items)\n total_loss += statistics[\"loss_sum\"]\n total_targets += count\n total_weight += statistics[\"weight_sum\"]\n total_positive += statistics[\"positive_sum\"]\n total_negative += statistics[\"negative_sum\"]\n total_pair_correct += statistics[\"pair_correct\"]\n batches += 1\n if batches % 25 == 0:\n release_accelerator_cache(device)\n if batches % 250 == 0:\n batch_bce = statistics[\"loss_sum\"] / max(\n statistics[\"weight_sum\"], 1.0e-12\n )\n running_bce = total_loss / max(total_weight, 1.0e-12)\n log(\n f\"epoch={epoch}/{config.epochs} batches={batches:,} \"\n f\"targets={total_targets:,} \"\n f\"batch_bce={batch_bce:.6f} \"\n f\"running_bce={running_bce:.6f} \"\n f\"elapsed={time.monotonic() - started:.1f}s\"\n )\n\n return {\n \"epoch\": float(epoch),\n \"loss\": total_loss / max(total_weight, 1.0e-12),\n \"positive_logit\": total_positive / max(total_targets, 1),\n \"negative_logit\": total_negative / max(total_targets * negative_count, 1),\n \"pair_accuracy\": total_pair_correct\n / max(total_targets * negative_count, 1),\n \"targets\": float(total_targets),\n \"target_weight_sum\": total_weight,\n \"batches\": float(batches),\n \"seconds\": time.monotonic() - started,\n }\n\n\[email protected]_grad()\ndef evaluate_validation_bce(\n model: base.SequentialTwoTower,\n prepared: SequenceData,\n config: Config,\n target_weights: np.ndarray,\n window_starts: np.ndarray,\n training_counts: np.ndarray,\n validation_counts: np.ndarray,\n item_correction: torch.Tensor,\n known_positive_packed: np.ndarray | None,\n device: torch.device,\n) -> dict[str, float]:\n \"\"\"Measure deterministic BCE on each user's chronological target suffix.\"\"\"\n model.eval()\n totals = training_counts + validation_counts\n users_with_validation = np.flatnonzero(validation_counts > 0)\n ordered = users_with_validation[\n np.argsort(totals[users_with_validation], kind=\"stable\")\n ]\n groups = [\n ordered[start : start + config.user_batch_size]\n for start in range(0, len(ordered), config.user_batch_size)\n ]\n generator = torch.Generator(device=\"cpu\").manual_seed(config.seed + 20_000)\n total_loss = 0.0\n total_targets = 0\n total_weight = 0.0\n started = time.monotonic()\n\n for group_number, users in enumerate(groups, start=1):\n group_totals = totals[users]\n group_starts = training_counts[users]\n hidden = torch.zeros(\n 1, len(users), config.embedding_dim, device=device\n )\n max_targets = int(group_totals.max())\n for start in range(0, max_targets, config.bptt_steps):\n steps = min(config.bptt_steps, max_targets - start)\n block = build_dense_block(\n prepared,\n target_weights,\n users,\n window_starts[users],\n group_totals,\n group_starts,\n start,\n steps,\n config,\n )\n sequence_outputs, _ = event_sequence(\n model,\n torch.from_numpy(block[0].astype(np.int64)).to(device),\n torch.from_numpy(block[1].astype(np.float32)).to(device),\n hidden,\n )\n hidden = torch.stack(\n [\n sequence_outputs[row, count - 1]\n if count\n else hidden[0, row]\n for row, count in enumerate(block[7])\n ]\n ).unsqueeze(0)\n flat_scored = block[5].reshape(-1)\n if not np.any(flat_scored):\n continue\n mask = torch.from_numpy(flat_scored).to(device)\n sequence_vectors = sequence_outputs.reshape(\n -1, config.embedding_dim\n )[mask]\n flat_states = block[3].reshape(-1)[flat_scored]\n positive_items = torch.from_numpy(\n block[4].reshape(-1)[flat_scored].astype(np.int64)\n ).to(device)\n known_history = torch.from_numpy(\n block[6].reshape(-1, config.bptt_steps)[flat_scored].astype(\n np.int64\n )\n ).to(device)\n liked_items = torch.from_numpy(\n prepared.liked_states[flat_states].astype(np.int64)\n ).to(device)\n flat_user_rows = np.broadcast_to(\n users[:, None], block[5].shape\n ).reshape(-1)[flat_scored]\n negatives, negative_weights = draw_training_negatives(\n positive_items,\n known_history,\n liked_items,\n flat_user_rows,\n known_positive_packed,\n prepared,\n config,\n item_correction,\n generator,\n device,\n )\n statistics, _ = score_target_microbatches(\n model,\n prepared,\n config,\n sequence_vectors,\n flat_states,\n block[2].reshape(-1, 4)[flat_scored].astype(np.float32),\n positive_items,\n negatives,\n negative_weights,\n block[8].reshape(-1)[flat_scored],\n device,\n backward=False,\n )\n total_loss += statistics[\"loss_sum\"]\n total_targets += len(positive_items)\n total_weight += statistics[\"weight_sum\"]\n if group_number % 25 == 0:\n release_accelerator_cache(device)\n return {\n \"validation_loss\": total_loss / max(total_weight, 1.0e-12),\n \"validation_targets\": float(total_targets),\n \"validation_target_weight_sum\": total_weight,\n \"validation_seconds\": time.monotonic() - started,\n }\n\n\ndef snapshot_model(model: nn.Module) -> dict[str, torch.Tensor]:\n return {\n name: value.detach().cpu().clone()\n for name, value in model.state_dict().items()\n }\n\n\ndef save_epoch_checkpoint(\n checkpoint_dir: Path,\n epoch: int,\n model_state: dict[str, torch.Tensor],\n config: Config,\n diagnostics: dict[str, float],\n optimizer: torch.optim.Optimizer,\n generator: torch.Generator,\n) -> Path:\n checkpoint_dir.mkdir(parents=True, exist_ok=True)\n path = checkpoint_dir / f\"epoch-{epoch:02d}.pt\"\n torch.save(\n {\n \"epoch\": epoch,\n \"config\": asdict(config),\n \"diagnostics\": diagnostics,\n \"model_state\": model_state,\n \"optimizer_state\": optimizer.state_dict(),\n \"generator_state\": generator.get_state(),\n },\n path,\n destination[row, : len(retained)] = retained.astype(np.uint16)\n return liked, disliked, like_counts, dislike_counts\n\n\nclass SequentialTwoTower(nn.Module):\n \"\"\"Order-sensitive user tower and independent catalog-aware item tower.\"\"\"\n\n def __init__(self, item_features: np.ndarray, config: TowerConfig) -> None:\n super().__init__()\n n_items, n_item_features = item_features.shape\n dim = config.embedding_dim\n self.n_items = n_items\n self.context_embedding = nn.Embedding(\n n_items + 1, dim, padding_idx=n_items\n )\n self.event_projection = nn.Sequential(\n nn.Linear(dim + 3, config.hidden_dim),\n nn.GELU(),\n )\n self.sequence_encoder = nn.GRU(\n input_size=config.hidden_dim,\n hidden_size=dim,\n num_layers=1,\n batch_first=True,\n )\n self.user_tower = nn.Sequential(\n nn.Linear(dim * 3 + 6, config.user_hidden_1),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.user_hidden_1, config.user_hidden_2),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.user_hidden_2, dim),\n )\n self.candidate_embedding = nn.Embedding(n_items, dim)\n self.register_buffer(\n \"item_features\", torch.from_numpy(item_features), persistent=True\n )\n self.item_tower = nn.Sequential(\n nn.Linear(dim + n_item_features, config.item_hidden_1),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.item_hidden_1, config.item_hidden_2),\n nn.GELU(),\n nn.Dropout(config.dropout),\n nn.Linear(config.item_hidden_2, dim),\n )\n nn.init.normal_(self.context_embedding.weight, std=0.02)\n nn.init.normal_(self.candidate_embedding.weight, std=0.02)\n with torch.no_grad():\n self.context_embedding.weight[n_items].zero_()\n\n def _mean_context_channel(self, items: torch.Tensor) -> torch.Tensor:\n mask = items.ne(self.n_items)\n vectors = self.context_embedding(items)\n return (vectors * mask.unsqueeze(-1)).sum(dim=1) / mask.sum(\n dim=1, keepdim=True\n ).clamp_min(1)\n\n def encode_items(self, items: torch.Tensor) -> torch.Tensor:\n inputs = torch.cat(\n [self.candidate_embedding(items), self.item_features[items]], dim=1\n )\n return self.item_tower(inputs)\n\n\ndef invalid_negative_mask(\n dislike_count: int,\n) -> np.ndarray:\n start = int(\n np.searchsorted(\n timestamps,\n history_end - base.PROFILE_LOOKBACK_DAYS * base.DAY,\n side=\"left\",\n )\n )\n result = np.zeros(6, dtype=np.float32)\n recent_items = items[start:]\n if len(recent_items):\n result[0] = np.log1p(len(recent_items)) / 8.0\n result[1] = np.log1p(len(np.unique(recent_items))) / 6.0\n result[2] = float(np.mean(event_features[start:, 0]))\n result[3] = float(np.mean(event_features[start:, 1]))\n result[4] = np.log1p(like_count) / 6.0\n result[5] = np.log1p(dislike_count) / 4.0\n return result\n\n\ndef prepare_data(\n public: Path, stage: str, config: Config\n) -> tuple[SequenceData, CandidateSet]:\n data = public / \"data\"\n catalog = pd.read_parquet(data / \"candidate_catalog.parquet\").sort_values(\n \"popularity_rank\"\n )\n item_ids = catalog.item_id.to_numpy(dtype=np.int64, copy=True)\n item_index = pd.Index(item_ids)\n n_items = len(item_ids)\n if n_items >= np.iinfo(np.uint16).max:\n raise RuntimeError(\"uint16 item storage is too small for this catalog\")\n\n if stage == \"validation\":\n prefixes = (\"train\",)\n history_end = base.TRAIN_END\n candidate_name = \"validation_eligible_candidates.npz\"\n state_name = \"validation_preference_state.parquet\"\n else:\n prefixes = (\"train\", \"val\")\n history_end = base.VALIDATION_END\n candidate_name = \"test_eligible_candidates.npz\"\n state_name = \"test_preference_state.parquet\"\n eligibility = CandidateSet.read(data / candidate_name)\n if not np.array_equal(eligibility.item_ids, item_ids):\n raise RuntimeError(\"catalog and candidate artifact item order differ\")\n\n explicit_by_user = base.load_explicit_events(data, prefixes, item_index)\n inference_likes, inference_dislikes, like_counts, dislike_counts = (\n base.load_inference_state(\n data / state_name,\n item_index,\n eligibility.uids,\n n_items,\n )\n )\n inference_rows = {int(uid): row for row, uid in enumerate(eligibility.uids)}\n inference_user_features = np.zeros((len(eligibility.uids), 6), dtype=np.float32)\n inference_sequence_rows = np.full(len(eligibility.uids), -1, dtype=np.int32)\n\n # Snapshot zero is the shared empty explicit state.\n liked_states = [np.full(base.LIKE_ITEMS, n_items, dtype=np.uint16)]\n disliked_states = [\n np.full(base.DISLIKE_ITEMS, n_items, dtype=np.uint16)\n ]\n state_counts: list[tuple[int, int]] = [(0, 0)]\n user_ids: list[int] = []\n offsets = [0]\n timestamp_parts: list[np.ndarray] = []\n item_parts: list[np.ndarray] = []\n event_parts: list[np.ndarray] = []\n activity_parts: list[np.ndarray] = []\n state_parts: list[np.ndarray] = []\n\n listen_paths = [data / f\"{prefix}_listens.parquet\" for prefix in prefixes]\n started = time.monotonic()\n for seen_users, (uid, listens) in enumerate(\n base.merge_listen_streams(listen_paths, item_index), start=1\n ):\n parts = build_sequence_parts(\n listens,\n explicit_by_user.get(uid, base.empty_explicit_events()),\n n_items,\n liked_states,\n disliked_states,\n state_counts,\n )\n if parts is None:\n continue\n sequence_row = len(user_ids)\n user_ids.append(uid)\n timestamp_parts.append(parts[0].astype(np.uint32, copy=False))\n item_parts.append(parts[1])\n event_parts.append(parts[2])\n activity_parts.append(parts[3])\n state_parts.append(parts[4])\n offsets.append(offsets[-1] + len(parts[1]))\n inference_row = inference_rows.get(uid)\n if inference_row is not None:\n inference_sequence_rows[inference_row] = sequence_row\n inference_user_features[inference_row] = final_user_features(\n parts[0],\n parts[1],\n parts[2],\n history_end,\n int(like_counts[inference_row]),\n int(dislike_counts[inference_row]),\n )\n if seen_users % 1_000 == 0:\n log(\n f\"prepared users={seen_users:,} tokens={offsets[-1]:,} \"\n f\"states={len(liked_states):,} \"\n f\"elapsed={time.monotonic() - started:.1f}s\"\n )\n\n missing = np.flatnonzero(inference_sequence_rows < 0)\n if len(missing):\n inference_user_features[missing, 4] = np.log1p(like_counts[missing]) / 6.0\n inference_user_features[missing, 5] = (\n np.log1p(dislike_counts[missing]) / 4.0\n )\n log(f\"target users without qualifying catalog sequence={len(missing):,}\")\n\n prepared = SequenceData(\n item_ids=item_ids,\n item_features=base.build_item_features(catalog),\n user_ids=np.asarray(user_ids, dtype=np.int64),\n offsets=np.asarray(offsets, dtype=np.int64),\n sequence_timestamps=np.concatenate(timestamp_parts),\n sequence_items=np.concatenate(item_parts),\n event_features=np.concatenate(event_parts),\n activity_features=np.concatenate(activity_parts),\n state_ids=np.concatenate(state_parts),\n liked_states=np.stack(liked_states),\n disliked_states=np.stack(disliked_states),\n state_counts=np.asarray(state_counts, dtype=np.int32),\n inference_user_rows=inference_sequence_rows,\n inference_liked_items=inference_likes,\n inference_disliked_items=inference_dislikes,\n inference_user_features=inference_user_features,\n )\n log(\n f\"prepared stage={stage} sequence_users={len(prepared.user_ids):,} \"\n f\"tokens={len(prepared.sequence_items):,} \"\n f\"next_item_targets={prepared.target_count:,} \"\n f\"states={len(prepared.liked_states):,}\"\n )\n return prepared, eligibility\n\n\ndef cache_metadata(stage: str, config: Config) -> dict[str, object]:\n return {\n \"cache_version\": CACHE_VERSION,\n \"stage\": stage,\n \"positive_listen_threshold\": base.LISTEN_POSITIVE_THRESHOLD,\n \"profile_lookback_days\": base.PROFILE_LOOKBACK_DAYS,\n \"event_time_feature\": \"exp-negative-gap-over-60-days\",\n }\n\n\ndef save_cache(\n path: Path, prepared: SequenceData, stage: str, config: Config\n) -> None:\n path.parent.mkdir(parents=True, exist_ok=True)\n np.savez(\n path,\n metadata=np.asarray(json.dumps(cache_metadata(stage, config), sort_keys=True)),\n **prepared.__dict__,\n )\n log(f\"saved public-derived full-sequence cache: {path}\")\n\n\ndef load_cache(path: Path, stage: str, config: Config) -> SequenceData:\n with np.load(path, allow_pickle=False) as artifact:\n actual = json.loads(str(artifact[\"metadata\"].item()))\n expected = cache_metadata(stage, config)\n if actual != expected:\n raise RuntimeError(f\"sequence cache differs: {actual} != {expected}\")\n values = {\n field: artifact[field].copy()\n for field in SequenceData.__dataclass_fields__\n }\n prepared = SequenceData(**values)\n log(\n f\"loaded full-sequence cache: tokens={len(prepared.sequence_items):,} \"\n f\"targets={prepared.target_count:,}\"\n )\n return prepared\n\n\ndef load_or_prepare(\n public: Path,\n stage: str,\n config: Config,\n cache: Path | None,\n) -> tuple[SequenceData, CandidateSet]:\n candidate_name = (\n \"validation_eligible_candidates.npz\"\n if stage == \"validation\"\n else \"test_eligible_candidates.npz\"\n )\n eligibility = CandidateSet.read(public / \"data\" / candidate_name)\n if cache is not None and cache.is_file():\n prepared = load_cache(cache, stage, config)\n if not np.array_equal(prepared.item_ids, eligibility.item_ids):\n raise RuntimeError(\"cached catalog differs from candidate artifact\")\n return prepared, eligibility\n prepared, built_eligibility = prepare_data(public, stage, config)\n if cache is not None:\n save_cache(cache, prepared, stage, config)\n return prepared, built_eligibility\n\n\ndef build_item_correction(\n prepared: SequenceData,\n config: Config,\n target_starts: np.ndarray,\n target_counts: np.ndarray,\n) -> torch.Tensor:\n \"\"\"Estimate in-batch propensities from the selected training targets.\"\"\"\n raise RuntimeError(\"chronological target partition produced negative counts\")\n if not np.array_equal(\n window_starts + training_counts + validation_counts, total_targets\n ):\n raise RuntimeError(\"chronological target window is not a final contiguous suffix\")\n return window_starts\n\n\ndef chronological_target_partition(\n prepared: SequenceData,\n validation_fraction: float,\n max_train_targets_per_user: int,\n max_validation_targets_per_user: int,\n) -> TargetPartition:\n \"\"\"Select a recent contiguous window with a final out-of-time BCE suffix.\"\"\"\n totals = np.maximum(prepared.lengths - 1, 0).astype(np.int64)\n validation = np.floor(totals * validation_fraction).astype(np.int64)\n validation[(totals >= 2) & (validation == 0)] = 1\n if max_validation_targets_per_user:\n validation = np.minimum(validation, max_validation_targets_per_user)\n\n history_training = totals - validation\n training = history_training.copy()\n if max_train_targets_per_user:\n training = np.minimum(training, max_train_targets_per_user)\n window_starts = retained_target_window_starts(totals, validation, training)\n return TargetPartition(\n window_starts=window_starts,\n training_counts=training,\n validation_counts=validation,\n history_training_counts=history_training,\n )\n\n\ndef load_organic_like_outcomes(\n public: Path, stage: str, item_ids: np.ndarray\n) -> dict[int, dict[int, np.ndarray]]:\n \"\"\"Load mapped organic-like timestamps from public fitting history.\"\"\"\n prefixes = (\"train\",) if stage == \"validation\" else (\"train\", \"val\")\n frames = [\n pd.read_parquet(\n public / \"data\" / f\"{prefix}_likes.parquet\",\n columns=[\"uid\", \"timestamp\", \"item_id\", \"is_organic\"],\n )\n for prefix in prefixes\n ]\n likes = pd.concat(frames, ignore_index=True)\n likes = likes[likes.is_organic == 1].copy()\n mapped = pd.Index(item_ids).get_indexer(\n likes.item_id.to_numpy(dtype=np.int64)\n )\n likes = likes.loc[mapped >= 0, [\"uid\", \"timestamp\"]].copy()\n likes[\"item\"] = mapped[mapped >= 0].astype(np.int32, copy=False)\n likes = likes.sort_values([\"uid\", \"item\", \"timestamp\"], kind=\"stable\")\n nested: dict[int, dict[int, list[int]]] = {}\n for raw_uid, raw_timestamp, raw_item in likes.itertuples(\n index=False, name=None\n ):\n nested.setdefault(int(raw_uid), {}).setdefault(int(raw_item), []).append(\n int(raw_timestamp)\n )\n result = {\n uid: {\n item: np.asarray(timestamps, dtype=np.int64)\n for item, timestamps in items.items()\n }\n for uid, items in nested.items()\n }\n log(\n f\"loaded organic-like outcomes: rows={len(likes):,} \"\n f\"users={len(result):,}\"\n )\n return result\n\n\ndef build_target_weights(\n public: Path,\n stage: str,\n prepared: SequenceData,\n training_counts: np.ndarray,\n config: Config,\n) -> np.ndarray:\n \"\"\"Build nonzero per-token weights without dropping adjacent targets.\"\"\"\n if config.target_weighting == \"uniform\":\n weights = np.ones(len(prepared.sequence_items), dtype=np.float32)\n else:\n weights = np.full(\n len(prepared.sequence_items),\n config.ordinary_target_weight,\n dtype=np.float32,\n )\n organic_strong = (\n (prepared.event_features[:, 0] > 0.5)\n & (prepared.event_features[:, 1] >= 0.8)\n )\n weights[organic_strong] = config.organic_strong_target_weight\n if config.target_weighting == \"discovery-graded\":\n like_outcomes = load_organic_like_outcomes(\n public, stage, prepared.item_ids\n )\n history_end = (\n base.TRAIN_END if stage == \"validation\" else base.VALIDATION_END\n )\n horizon = config.discovery_horizon_days * base.DAY\n grade_three = 0\n for user_row, uid in enumerate(prepared.user_ids):\n start = int(prepared.offsets[user_row])\n end = int(prepared.offsets[user_row + 1])\n if end - start < 2:\n continue\n last_training_target = start + int(training_counts[user_row])\n first_validation_target = last_training_target + 1\n split_timestamp = (\n int(prepared.sequence_timestamps[first_validation_target])\n if first_validation_target < end\n else history_end\n )\n user_likes = like_outcomes.get(int(uid), {})\n seen = {int(prepared.sequence_items[start])}\n for position in range(start + 1, end):\n item = int(prepared.sequence_items[position])\n novel = item not in seen\n seen.add(item)\n if not novel:\n # Repeats remain training examples, but they cannot be\n # eligible discovery outcomes at the corresponding\n # cutoff and therefore keep only the ordinary weight.\n weights[position] = config.ordinary_target_weight\n continue\n timestamps = user_likes.get(item)\n if timestamps is None:\n continue\n timestamp = int(prepared.sequence_timestamps[position])\n upper = min(\n timestamp + horizon,\n split_timestamp\n if position <= last_training_target\n else history_end,\n )\n outcome = int(np.searchsorted(timestamps, timestamp, side=\"left\"))\n if outcome < len(timestamps) and int(timestamps[outcome]) < upper:\n weights[position] = config.organic_like_target_weight\n grade_three += 1\n log(\n f\"grade-3 first-encounter targets={grade_three:,} \"\n f\"horizon_days={config.discovery_horizon_days}\"\n )\n log(\n f\"target weighting={config.target_weighting} \"\n f\"min={weights.min():.3f} mean={weights.mean():.3f} \"\n f\"max={weights.max():.3f}\"\n )\n return weights\n\n\ndef build_dense_block(\n prepared: SequenceData,\n target_weights: np.ndarray,\n users: np.ndarray,\n window_starts: np.ndarray,\n process_counts: np.ndarray,\n score_starts: np.ndarray,\n start: int,\n steps: int,\n config: Config,\n) -> tuple[np.ndarray, ...]:\n \"\"\"Build one causal sequence block and mark the target rows to score.\"\"\"\n batch_size = len(users)\n padding_item = len(prepared.item_ids)\n input_items = np.full((batch_size, steps), padding_item, dtype=np.uint16)\n event_features = np.zeros((batch_size, steps, 3), dtype=np.float16)\n activity_features = np.zeros((batch_size, steps, 4), dtype=np.float16)\n state_ids = np.zeros((batch_size, steps), dtype=np.uint32)\n targets = np.zeros((batch_size, steps), dtype=np.uint16)\n weights = np.ones((batch_size, steps), dtype=np.float32)\n scored = np.zeros((batch_size, steps), dtype=bool)\n histories = np.full(\n (batch_size, steps, config.bptt_steps),\n padding_item,\n dtype=np.uint16,\n )\n processed = np.zeros(batch_size, dtype=np.int32)\n\n lengths = prepared.lengths\n for row, user_row in enumerate(users):\n count = min(steps, max(int(process_counts[row]) - start, 0))\n if not count:\n continue\n offset = int(prepared.offsets[user_row])\n window_start = int(window_starts[row])\n absolute_start = window_start + start\n positions = slice(\n offset + absolute_start,\n offset + absolute_start + count,\n )\n target_positions = slice(\n offset + absolute_start + 1,\n offset + absolute_start + count + 1,\n )\n input_items[row, :count] = prepared.sequence_items[positions]\n event_features[row, :count] = prepared.event_features[positions]\n activity_features[row, :count] = prepared.activity_features[positions]\n state_ids[row, :count] = prepared.state_ids[positions]\n targets[row, :count] = prepared.sequence_items[target_positions]\n weights[row, :count] = target_weights[target_positions]\n score_from = min(count, max(int(score_starts[row]) - start, 0))\n"
}
}