update

Author: ElektrikAkar

.claude/CLAUDE.md

@@ -26,3 +26,8 @@
 - [ ] CHANGELOG.md updated (Unreleased)
 - [ ] Docs updated (if user-facing)
 - [ ] Optional deps remain optional
+
+## Guidelines
+- Buffer > thread_local >> heap allocation: already enforced everywhere
+- No naked `new`/`delete` in core: already enforced
+- Contiguous arrays in hot paths: `Data::p_vec` as `vector<vector<data_t>>` is correct for variable-length series

.claude/TODO.md

@@ -45,16 +45,4 @@
 
 ### Platform
 - [ ] ARM Mac Studio investigation: test CPU path on Apple Silicon, evaluate Metal compute for GPU path
-- [ ] HIPify for AMD GPU — **accept community PRs only**, do not invest core developer time
-
-## Guidelines (not TODOs — current conventions)
-- Buffer > thread_local >> heap allocation: already enforced everywhere
-- No naked `new`/`delete` in core: already enforced
-- Contiguous arrays in hot paths: `Data::p_vec` as `vector<vector<data_t>>` is correct for variable-length series
-
-## Removed (completed or cut)
-- ~~NaN/-ffast-math robustification~~ — **DONE**: explicit fast-math sub-flags + `std::isnan()` wrapper
-- ~~Eigen 5.0.1 exploration~~ — **CUT**: Eigen used only as aligned allocator, no gap identified
-- ~~Condensed distance matrix~~ — **DONE**: `DenseDistanceMatrix` already uses packed triangular N*(N+1)/2
-- ~~Unnecessary memory allocations audit~~ — **CUT**: 30+ `thread_local` declarations already in place
-- ~~Device-side LB pruning (pair-level)~~ — **DONE**: `compact_active_pairs` in `cuda_dtw.cu`
+- [ ] HIPify for AMD GPU — **accept community PRs only**, do not invest core developer time

.claude/settings.json

@@ -1,7 +1,6 @@
 {
   "env": {
-    "CLAUDE_CODE_TASK_LIST_ID": "dtwcpp",
+    "CLAUDE_CODE_TASK_LIST_ID": "dtwcpp"
   },
-  "plansDirectory": "./.claude/reports",
-  "model": "opusplan"
+  "plansDirectory": "./.claude/reports"
 }

cmake/CPM.cmake

@@ -2,8 +2,8 @@
 #
 # SPDX-FileCopyrightText: Copyright (c) 2019-2023 Lars Melchior and contributors
 
-set(CPM_DOWNLOAD_VERSION 0.40.2)
-set(CPM_HASH_SUM "c8cdc32c03816538ce22781ed72964dc864b2a34a310d3b7104812a5ca2d835d")
+set(CPM_DOWNLOAD_VERSION 0.42.1)
+set(CPM_HASH_SUM "f3a6dcc6a04ce9e7f51a127307fa4f699fb2bade357a8eb4c5b45df76e1dc6a5")
 
 if(CPM_SOURCE_CACHE)
   set(CPM_DOWNLOAD_LOCATION "${CPM_SOURCE_CACHE}/cpm/CPM_${CPM_DOWNLOAD_VERSION}.cmake")

cmake/Dependencies.cmake

@@ -86,7 +86,7 @@ function(dtwc_setup_dependencies)
       CPMAddPackage(
         NAME benchmark
         GITHUB_REPOSITORY google/benchmark
-        VERSION 1.9.1
+        VERSION 1.9.5
         OPTIONS
           "BENCHMARK_ENABLE_TESTING OFF"
           "BENCHMARK_ENABLE_GTEST_TESTS OFF"
@@ -98,7 +98,7 @@ function(dtwc_setup_dependencies)
   if(DTWC_ENABLE_YAML AND NOT TARGET yaml-cpp)
     CPMAddPackage(
       NAME yaml-cpp
-      URL "https://github.com/jbeder/yaml-cpp/archive/refs/tags/0.8.0.tar.gz"
+      URL "https://github.com/jbeder/yaml-cpp/archive/refs/tags/0.9.0.tar.gz"
       SYSTEM
       EXCLUDE_FROM_ALL
       OPTIONS "YAML_CPP_BUILD_TESTS OFF" "YAML_CPP_BUILD_TOOLS OFF"

dtwc/missing_utils.hpp

@@ -7,8 +7,6 @@
  *     etc.) but NOT -ffinite-math-only, so std::isnan() is correct and available.
  *   - is_missing() is a thin semantic wrapper over std::isnan() — prefer this name
  *     in domain code for clarity ("missing value" not "not-a-number").
- *   - DenseDistanceMatrix tracks computed entries via std::vector<bool> (not NaN
- *     sentinels), which is safe regardless of floating-point mode.
  *
  * @author Volkan Kumtepeli
  * @date 02 Apr 2026

dtwc/scores.cpp

@@ -132,7 +132,7 @@ double daviesBouldinIndex(Problem &prob)
     double max_ratio = 0.0;
     for (int j = 0; j < Nc; ++j) {
       if (i != j) {
-        double d_ij = prob.distByInd(prob.centroids_ind[i], prob.centroids_ind[j]);
+        const double d_ij = prob.distByInd(prob.centroids_ind[i], prob.centroids_ind[j]);
         if (d_ij > 0) {
           double ratio = (scatter[i] + scatter[j]) / d_ij;
           max_ratio = std::max(max_ratio, ratio);
@@ -167,14 +167,11 @@ double dunnIndex(Problem &prob)
 
   for (int i = 0; i < N; ++i) {
     for (int j = i + 1; j < N; ++j) {
-      double d = prob.distByInd(i, j);
-      if (prob.clusters_ind[i] == prob.clusters_ind[j]) {
-        // Same cluster: contributes to intra-cluster diameter
-        max_intra = std::max(max_intra, d);
-      } else {
-        // Different clusters: contributes to inter-cluster distance
-        min_inter = std::min(min_inter, d);
-      }
+      const double d = prob.distByInd(i, j);
+      if (prob.clusters_ind[i] == prob.clusters_ind[j]) 
+        max_intra = std::max(max_intra, d); // Same cluster: contributes to intra-cluster diameter
+       else 
+        min_inter = std::min(min_inter, d); // Different clusters: contributes to inter-cluster distance
     }
   }