sublinear 0.2.0

# BMSSP Integration Implementation Plan

## 🎯 Overview

This plan outlines the integration of **@ruvnet/bmssp** (Bounded Multi-Source Shortest Path) with the existing sublinear-time-solver codebase. BMSSP provides WebAssembly-powered graph pathfinding that's 10-15x faster than JavaScript implementations.

## 📊 Integration Analysis

### Current Architecture
- **Core**: Sublinear solver algorithms (Neumann, random-walk, push methods)
- **Graph Tools**: PageRank, effective resistance, centrality measures
- **Matrix Operations**: Dense/sparse matrix support
- **MCP Interface**: Model Context Protocol server with solver tools

### BMSSP Capabilities
- **Performance**: 10-15x faster than JS implementations via WASM
- **Multi-source**: Simultaneous pathfinding from multiple sources
- **Bidirectional**: Optimized search from both ends
- **Neural Features**: WasmNeuralBMSSP for semantic pathfinding
- **Zero Dependencies**: Pure WASM with TypeScript support

## 🔗 Integration Points

### 1. Core Solver Enhancement
**Location**: `src/core/`

#### New BMSSP Solver Class
```typescript
// src/core/bmssp-solver.ts
import { BmsSpGraph } from '@ruvnet/bmssp';
import { WasmGraph, WasmNeuralBMSSP } from '@ruvnet/bmssp';

export class BMSSPSolver extends SublinearSolver {
  private bmsspGraph?: BmsSpGraph;
  private wasmGraph?: WasmGraph;
  private neuralBMSSP?: WasmNeuralBMSSP;
}
```

#### Integration Methods
- **Graph Construction**: Convert matrices to BMSSP graph format
- **Hybrid Solving**: Use BMSSP for shortest paths, sublinear for linear systems
- **Performance Switching**: Automatic method selection based on problem size

### 2. Graph Tools Enhancement
**Location**: `src/mcp/tools/graph.ts`

#### Enhanced Features
- **Fast PageRank**: Use BMSSP for graph traversal optimization
- **Multi-source Centrality**: Leverage BMSSP's multi-source capabilities
- **Semantic Pathfinding**: Neural BMSSP for embeddings-based paths

### 3. Matrix Operations Bridge
**Location**: `src/core/matrix.ts`

#### Conversion Utilities
- **Matrix to Graph**: Convert adjacency matrices to BMSSP format
- **Sparse Optimization**: Leverage BMSSP's efficient sparse handling
- **Memory Management**: WASM memory lifecycle integration

## 🛠 Implementation Strategy

### Phase 1: Core Integration (Week 1-2)
#### Deliverables
1. **BMSSP Wrapper Class**
   - `src/core/bmssp-wrapper.ts`
   - WASM lifecycle management
   - Memory safety patterns

2. **Matrix Conversion Utilities**
   - `src/core/bmssp-bridge.ts`
   - Adjacency matrix → BMSSP graph
   - Laplacian matrix → BMSSP format

3. **Hybrid Solver**
   - `src/core/hybrid-solver.ts`
   - Automatic method selection
   - Performance benchmarking

### Phase 2: Graph Algorithms (Week 3)
#### Deliverables
1. **Enhanced PageRank**
   ```typescript
   // Multi-source PageRank using BMSSP
   async pageRankBMSSP(adjacency: Matrix, sources?: number[])
   ```

2. **Fast Shortest Paths**
   ```typescript
   // Leverage BMSSP's O(m·log^(2/3) n) complexity
   async shortestPathsBMSSP(graph: Matrix, sources: number[], targets: number[])
   ```

3. **Centrality Measures**
   ```typescript
   // Betweenness centrality using BMSSP pathfinding
   async betweennessCentralityBMSSP(adjacency: Matrix)
   ```

### Phase 3: Neural Integration (Week 4)
#### Deliverables
1. **Semantic Pathfinding**
   ```typescript
   // Neural BMSSP for embedding-based paths
   class SemanticPathfinder {
     constructor(embeddings: Float64Array[], embeddingDim: number)
     async findSemanticPath(source: number, target: number, alpha: number)
   }
   ```

2. **Graph Embeddings**
   ```typescript
   // Update embeddings based on graph structure
   async updateGraphEmbeddings(gradients: Float64Array[], learningRate: number)
   ```

### Phase 4: MCP Tools Integration (Week 5)
#### Deliverables
1. **New MCP Tools**
   - `bmssp_shortest_path`
   - `bmssp_multi_source_pagerank`
   - `bmssp_semantic_pathfinding`

2. **Performance Tools**
   - `bmssp_benchmark`
   - `bmssp_memory_profile`

## 🏗 File Structure

```
src/
├── core/
│   ├── bmssp-wrapper.ts          # WASM lifecycle management
│   ├── bmssp-bridge.ts           # Matrix conversion utilities
│   ├── hybrid-solver.ts          # Hybrid BMSSP + sublinear solver
│   ├── semantic-pathfinder.ts    # Neural BMSSP integration
│   └── bmssp-benchmarks.ts       # Performance comparison
├── mcp/tools/
│   ├── bmssp-tools.ts            # BMSSP MCP tools
│   └── hybrid-graph-tools.ts    # Enhanced graph tools
└── integrations/
    ├── bmssp/
    │   ├── examples/             # Usage examples
    │   ├── benchmarks/           # Performance tests
    │   └── tests/                # Integration tests
```

## 📈 Performance Optimization Strategy

### 1. Automatic Method Selection
```typescript
class PerformanceOracle {
  selectOptimalMethod(
    problemSize: number,
    sparsity: number,
    queryType: 'single' | 'multi' | 'batch'
  ): 'bmssp' | 'sublinear' | 'hybrid' {
    // Intelligence-based selection
    if (queryType === 'multi' && problemSize > 1000) return 'bmssp';
    if (sparsity > 0.95 && problemSize > 10000) return 'bmssp';
    return 'hybrid';
  }
}
```

### 2. Memory Management
```typescript
class BMSSPMemoryManager {
  private wasmInstances: Map<string, any> = new Map();

  async getOrCreateInstance(graphId: string, config: any) {
    // Efficient WASM instance pooling
  }

  cleanup() {
    // Proper WASM memory cleanup
    this.wasmInstances.forEach(instance => instance.free());
  }
}
```

### 3. Batch Processing
```typescript
class BMSSPBatchProcessor {
  async processBatch(queries: PathQuery[]): Promise<PathResult[]> {
    // Leverage BMSSP's batch processing capabilities
    const graph = new BmsSpGraph();
    return graph.batch_shortest_paths(queries);
  }
}
```

## 🧪 Testing Strategy

### 1. Unit Tests
```typescript
// tests/bmssp-integration.test.ts
describe('BMSSP Integration', () => {
  test('Matrix conversion accuracy', () => {
    // Verify matrix → BMSSP graph conversion
  });

  test('Performance benchmarks', () => {
    // Compare BMSSP vs traditional methods
  });

  test('Memory safety', () => {
    // Ensure proper WASM cleanup
  });
});
```

### 2. Performance Tests
```typescript
// benchmarks/bmssp-vs-traditional.ts
const results = await benchmarkComparison({
  graphSizes: [1000, 10000, 100000],
  methods: ['javascript', 'bmssp', 'hybrid'],
  metrics: ['time', 'memory', 'accuracy']
});
```

### 3. Integration Tests
```typescript
// tests/hybrid-solver.test.ts
describe('Hybrid Solver', () => {
  test('Automatic method selection', () => {
    // Test intelligent algorithm switching
  });

  test('Cross-validation', () => {
    // Verify BMSSP and sublinear produce same results
  });
});
```

## 🎯 API Design

### 1. Enhanced Graph Tools
```typescript
interface BMSSPGraphTools extends GraphTools {
  // Multi-source pathfinding
  async multiSourceShortestPaths(
    adjacency: Matrix,
    sources: number[],
    targets?: number[]
  ): Promise<MultiPathResult>;

  // Semantic pathfinding
  async semanticPathfinding(
    graph: Matrix,
    embeddings: Float64Array[],
    source: number,
    target: number,
    alpha: number
  ): Promise<SemanticPathResult>;

  // Batch centrality computation
  async batchCentralityMeasures(
    adjacency: Matrix,
    measures: CentralityType[],
    nodes?: number[]
  ): Promise<BatchCentralityResult>;
}
```

### 2. MCP Tool Extensions
```typescript
// New MCP tools for BMSSP integration
const bmsspTools = [
  {
    name: 'bmssp_shortest_path',
    description: 'Ultra-fast shortest path using WASM',
    parameters: {
      adjacency: 'Matrix',
      source: 'number',
      target: 'number'
    }
  },
  {
    name: 'bmssp_multi_source_pagerank',
    description: 'Multi-source PageRank using BMSSP',
    parameters: {
      adjacency: 'Matrix',
      sources: 'number[]',
      damping: 'number?'
    }
  },
  {
    name: 'bmssp_semantic_pathfinding',
    description: 'Neural pathfinding with embeddings',
    parameters: {
      graph: 'Matrix',
      embeddings: 'Float64Array[]',
      source: 'number',
      target: 'number',
      alpha: 'number'
    }
  }
];
```

## 🚀 Usage Examples

### 1. Hybrid Pathfinding
```typescript
import { BMSSPHybridSolver } from './core/hybrid-solver.js';

const solver = new BMSSPHybridSolver({
  autoSelectMethod: true,
  bmsspEnabled: true
});

// Automatically selects optimal method
const result = await solver.shortestPath(adjacencyMatrix, source, target);
```

### 2. Multi-source Analysis
```typescript
import { BMSSPGraphTools } from './mcp/tools/bmssp-tools.js';

const sources = [0, 5, 10]; // Multiple starting points
const results = await BMSSPGraphTools.multiSourceShortestPaths(
  graph,
  sources
);

console.log(`Found ${results.paths.length} optimal paths`);
```

### 3. Semantic Pathfinding
```typescript
import { SemanticPathfinder } from './core/semantic-pathfinder.js';

const pathfinder = new SemanticPathfinder(embeddings, embeddingDim);
const semanticPath = await pathfinder.findSemanticPath(
  source,
  target,
  0.7 // alpha parameter for semantic weight
);
```

## 🎛 Configuration

### 1. Performance Tuning
```typescript
interface BMSSPConfig {
  // Automatic method selection
  autoSelect: boolean;

  // Performance thresholds
  bmsspThreshold: {
    minGraphSize: number;
    minSparsity: number;
    multiSourceMin: number;
  };

  // Memory management
  wasmPoolSize: number;
  memoryLimitMB: number;

  // Neural features
  enableSemanticPath: boolean;
  embeddingDim: number;
}
```

### 2. Integration Settings
```typescript
const config: BMSSPConfig = {
  autoSelect: true,
  bmsspThreshold: {
    minGraphSize: 1000,
    minSparsity: 0.9,
    multiSourceMin: 2
  },
  wasmPoolSize: 4,
  memoryLimitMB: 512,
  enableSemanticPath: true,
  embeddingDim: 128
};
```

## 📊 Expected Benefits

### 1. Performance Improvements
- **10-15x faster** shortest path computation
- **Sub-quadratic complexity** for large graphs
- **Batch processing** efficiency for multiple queries

### 2. New Capabilities
- **Multi-source pathfinding** - simultaneous computation
- **Semantic pathfinding** - embedding-based routes
- **Neural graph analysis** - learning-based optimization

### 3. Better Resource Utilization
- **WASM efficiency** - near-native performance
- **Memory optimization** - smart pooling and cleanup
- **Automatic scaling** - method selection based on problem size

## 🗓 Timeline

- **Week 1**: Core BMSSP wrapper and bridge utilities
- **Week 2**: Hybrid solver with automatic method selection
- **Week 3**: Enhanced graph algorithms integration
- **Week 4**: Neural BMSSP and semantic pathfinding
- **Week 5**: MCP tools integration and documentation

## 🔧 Dependencies

### Required Updates
```json
{
  "dependencies": {
    "@ruvnet/bmssp": "^1.0.0"
  },
  "devDependencies": {
    "@types/wasm": "^1.0.0"
  }
}
```

### TypeScript Configuration
```json
{
  "compilerOptions": {
    "experimentalDecorators": true,
    "allowSyntheticDefaultImports": true
  }
}
```

This integration will significantly enhance the sublinear-time-solver's graph processing capabilities while maintaining compatibility with existing APIs and adding powerful new features for semantic and multi-source pathfinding.