selen 0.15.5

Constraint Satisfaction Problem (CSP) solver
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66

use selen::prelude::*;
use std::time::Instant;

fn main() -> SolverResult<()> {
    println!("=== Solver Hanging Investigation ===\n");
    
    // Test 1: Single variable (this works)
    println!("Test 1: Single variable constraint");
    let start = Instant::now();
    let mut m = Model::default();
    let x = m.new_var(0.0.into(), 100.0.into());
    post!(m, x >= 10.0);
    let result = m.minimize(x)?;
    println!("  ✅ Single variable: {:?}", start.elapsed());
    if let Val::ValF(value) = result[x] {
        println!("     Value: {:.2}", value);
    }
    
    // Test 2: Two variables, no relation (should work)
    println!("\nTest 2: Two independent variables");
    let start = Instant::now();
    let mut m = Model::default();
    let x = m.new_var(0.0.into(), 100.0.into());
    let y = m.new_var(0.0.into(), 100.0.into());
    post!(m, x >= 10.0);
    post!(m, y >= 20.0);
    let result = m.minimize(x)?;
    println!("  ✅ Two independent: {:?}", start.elapsed());
    if let Val::ValF(value) = result[x] {
        println!("     Value: {:.2}", value);
    }
    
    // Test 3: Two variables with relation (this hangs)
    println!("\nTest 3: Two variables with x <= y relation");
    let start = Instant::now();
    let mut m = Model::default();
    let x = m.new_var(0.0.into(), 100.0.into());
    let y = m.new_var(0.0.into(), 100.0.into());
    post!(m, x >= 10.0);
    post!(m, y >= 20.0);
    
    println!("  About to add: x <= y");
    post!(m, x <= y);  // This is where it likely hangs
    println!("  ✅ Constraint added successfully");
    
    println!("  About to solve...");
    let result = m.minimize(x);
    let duration = start.elapsed();
    
    match result {
        Ok(solution) => {
            println!("  ✅ Solved in {:?}", duration);
            if let Val::ValF(value) = solution[x] {
                println!("     x = {:.2}", value);
            }
            if let Val::ValF(value) = solution[y] {
                println!("     y = {:.2}", value);
            }
        }
        Err(e) => {
            println!("  ❌ Error after {:?}: {:?}", duration, e);
        }
    }
    
    Ok(())
}