ferrotherm 0.35.0

Thermodynamic computing in pure Rust: sparse energy-based models, chromatic block-Gibbs, parallel tempering, thermodynamic linear algebra, stochastic differentiable programs, a variational compiler onto device topologies, exact inference by variable elimination, planted instances with known optima, sampler certificates, and a first-class joules ledger. std-only, zero dependencies, wasm-clean, deterministic by seed.
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
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83

"""Solve the agreed model through the MCP server and print the answer.

Its own file rather than a heredoc inside `check-answers.sh`: the script already nests one, and a
second whose terminator is also `PY` silently closes the outer one. Separate files also mean the
driver can be run by hand when a surface disagrees, which is when you most want to.

Prints nothing on failure, so the gate reports PRODUCED NOTHING and shows stderr.
"""

import json
import os
import subprocess
import sys

ARGS = {
    "variables": [
        {"name": "a", "values": 3},
        {"name": "b", "values": 3},
        {"name": "t", "lo": 10, "hi": 13},
    ],
    "constraints": [
        {"type": "not_equal", "a": "a", "b": "b"},
        {"type": "at_most", "k": 1, "of": [{"var": "a", "value": 0}, {"var": "b", "value": 0}]},
        {"type": "fix", "var": "t", "value": 12},
    ],
    "objective": {
        "maximize": True,
        "terms": [
            {"var": "a", "value": 1, "weight": 3},
            {"var": "b", "value": 2, "weight": 4},
        ],
    },
    "tries": 64,
}

REQS = [
    {
        "jsonrpc": "2.0", "id": 1, "method": "initialize",
        "params": {
            "protocolVersion": "2025-06-18",
            "capabilities": {},
            "clientInfo": {"name": "check-answers", "version": "1"},
        },
    },
    {
        "jsonrpc": "2.0", "id": 2, "method": "tools/call",
        "params": {"name": "ferrotherm_solve", "arguments": ARGS},
    },
]


def main() -> int:
    out_dir = sys.argv[1]
    proc = subprocess.run(
        ["./target/release/ferrotherm-mcp"],
        input="".join(json.dumps(r) + "\n" for r in REQS),
        capture_output=True,
        text=True,
        timeout=180,
    )
    for line in proc.stdout.splitlines():
        try:
            msg = json.loads(line)
        except Exception:
            continue
        if msg.get("id") != 2:
            continue
        payload = json.loads(msg["result"]["content"][0]["text"])
        v = payload["values"]
        answer = (
            f"a={v['a']} b={v['b']} t={v['t']} "
            f"feasible={str(payload['feasible']).lower()}"
        )
        with open(os.path.join(out_dir, "mcp.txt"), "w") as fh:
            fh.write(answer)
        return 0
    print("no tools/call result in the MCP response", file=sys.stderr)
    print(proc.stderr[:400], file=sys.stderr)
    return 1


if __name__ == "__main__":
    raise SystemExit(main())