@ohtaman/washi 0.2.0 → 0.4.0

package/dist/browser/washi.js

@@ -1,32 +1,37 @@
 import { a as e, c as t, i as n, l as r, n as i, o as a, r as o, s, t as c } from "./loader-CVt4CloI.js";
+//#region src/modeler/types.ts
+function l(e) {
+	return e && e.toLowerCase() === "maximize" ? "Maximize" : "Minimize";
+}
+//#endregion
 //#region src/modeler/boolean-expr.ts
-var l = 0, u = class {
+var u = 0, d = class {
 	and(e) {
-		return new f("and", [this, e]);
+		return new p("and", [this, e]);
 	}
 	or(e) {
-		return new f("or", [this, e]);
+		return new p("or", [this, e]);
 	}
 	implies(e) {
-		return new f("implies", [this, e]);
+		return new p("implies", [this, e]);
 	}
 	xor(e) {
-		return new f("xor", [this, e]);
+		return new p("xor", [this, e]);
 	}
 	not() {
-		return new f("not", [this]);
+		return new p("not", [this]);
 	}
-}, d = class extends u {
+}, f = class extends d {
 	id;
 	name;
 	value = null;
 	constructor(e) {
-		super(), this.id = ++l, this.name = e || `p_${this.id}`;
+		super(), this.id = ++u, this.name = e || `p_${this.id}`;
 	}
 	toString() {
 		return this.name;
 	}
-}, f = class extends u {
+}, p = class extends d {
 	op;
 	children;
 	constructor(e, t) {
@@ -36,27 +41,27 @@ var l = 0, u = class {
 		return this.op === "not" ? `not(${this.children[0].toString()})` : `(${this.children.map((e) => e.toString()).join(` ${this.op} `)})`;
 	}
 };
-function p(...e) {
+function m(...e) {
 	if (e.length === 0) throw Error("and() requires at least one expression");
-	return new f("and", e);
+	return new p("and", e);
 }
-function m(...e) {
+function h(...e) {
 	if (e.length === 0) throw Error("or() requires at least one expression");
-	return new f("or", e);
+	return new p("or", e);
 }
-function h(e) {
-	return new f("not", [e]);
+function g(e) {
+	return new p("not", [e]);
 }
-function g(e, t) {
-	return new f("implies", [e, t]);
+function _(e, t) {
+	return new p("implies", [e, t]);
 }
-function _(...e) {
+function v(...e) {
 	if (e.length === 0) throw Error("xor() requires at least one expression");
-	return new f("xor", e);
+	return new p("xor", e);
 }
 //#endregion
 //#region src/modeler/constraint.ts
-var ee = 0, v = class {
+var y = 0, b = class {
 	id;
 	lhs;
 	sense;
@@ -65,25 +70,22 @@ var ee = 0, v = class {
 	slack = null;
 	dual = null;
 	constructor(e, t, n, r) {
-		this.id = ++ee, this.lhs = e, this.sense = t, this.rhs = n, this.name = r || `c_${this.id}`;
-	}
-	get pi() {
-		return this.dual;
+		this.id = ++y, this.lhs = e, this.sense = t, this.rhs = n, this.name = r || `c_${this.id}`;
 	}
 	and(e) {
-		return new f("and", [this, e]);
+		return new p("and", [this, e]);
 	}
 	or(e) {
-		return new f("or", [this, e]);
+		return new p("or", [this, e]);
 	}
 	implies(e) {
-		return new f("implies", [this, e]);
+		return new p("implies", [this, e]);
 	}
 	xor(e) {
-		return new f("xor", [this, e]);
+		return new p("xor", [this, e]);
 	}
 	not() {
-		return new f("not", [this]);
+		return new p("not", [this]);
 	}
 	toString() {
 		let e = this.sense === "==" ? "=" : this.sense;
@@ -92,11 +94,11 @@ var ee = 0, v = class {
 };
 //#endregion
 //#region src/modeler/quad-expr.ts
-function y(e, t) {
+function x(e, t) {
 	let n = e.id, r = t.id;
 	return n <= r ? `${n}:${r}` : `${r}:${n}`;
 }
-var b = class e {
+var S = class e {
 	linearTerms;
 	quadTerms;
 	constant;
@@ -109,7 +111,7 @@ var b = class e {
 		return new e(this.linearTerms, this.quadTerms, this.constant);
 	}
 	getQuadCoeff(e, t) {
-		let n = y(e, t), r = this.quadTerms.get(n);
+		let n = x(e, t), r = this.quadTerms.get(n);
 		return r ? r.coeff : 0;
 	}
 	add(t) {
@@ -144,20 +146,45 @@ var b = class e {
 		return typeof t == "number" ? this.add(-t) : t instanceof e ? this.add(t.neg()) : ("terms" in t && t.terms instanceof Map, this.add(t.mul(-1)));
 	}
 	mul(t) {
-		let n = new e(/* @__PURE__ */ new Map(), /* @__PURE__ */ new Map(), this.constant * t);
-		for (let [e, r] of this.linearTerms) {
-			let i = r * t;
-			Math.abs(i) >= 1e-12 && n.linearTerms.set(e, i);
+		if (typeof t == "number") {
+			let n = new e(/* @__PURE__ */ new Map(), /* @__PURE__ */ new Map(), this.constant * t);
+			for (let [e, r] of this.linearTerms) {
+				let i = r * t;
+				Math.abs(i) >= 1e-12 && n.linearTerms.set(e, i);
+			}
+			for (let [e, r] of this.quadTerms) {
+				let i = r.coeff * t;
+				Math.abs(i) >= 1e-12 && n.quadTerms.set(e, {
+					var1: r.var1,
+					var2: r.var2,
+					coeff: i
+				});
+			}
+			return n;
 		}
-		for (let [e, r] of this.quadTerms) {
-			let i = r.coeff * t;
-			Math.abs(i) >= 1e-12 && n.quadTerms.set(e, {
-				var1: r.var1,
-				var2: r.var2,
-				coeff: i
-			});
+		if ("id" in t && "name" in t && !("terms" in t) && !("quadTerms" in t)) {
+			if (this.quadTerms.size > 0) throw Error("Multiplication creates a non-linear term of degree > 2 (not allowed in MILP/LP/QP)");
+			let n = new e();
+			Math.abs(this.constant) >= 1e-12 && n.linearTerms.set(t, this.constant);
+			for (let [r, i] of this.linearTerms) n = n.add(e.fromProduct(r, t, i));
+			return n;
 		}
-		return n;
+		let n = t, r = "quadTerms" in n ? n.quadTerms.size : 0, i = "quadTerms" in n ? n.linearTerms : n.terms;
+		if (this.quadTerms.size > 0 && (r > 0 || i.size > 0) || r > 0 && this.linearTerms.size > 0) throw Error("Multiplication creates a non-linear term of degree > 2 (not allowed in MILP/LP/QP)");
+		if (r === 0 && i.size === 0) return this.mul(n.constant);
+		if (this.quadTerms.size === 0 && this.linearTerms.size === 0) {
+			let t = n.mul(this.constant);
+			if (t instanceof e) return t;
+			let r = new e();
+			r.constant = t.constant;
+			for (let [e, n] of t.terms) r.linearTerms.set(e, n);
+			return r;
+		}
+		let a = new e();
+		if (a.constant = this.constant * n.constant, Math.abs(n.constant) >= 1e-12) for (let [e, t] of this.linearTerms) a.linearTerms.set(e, (a.linearTerms.get(e) || 0) + t * n.constant);
+		if (Math.abs(this.constant) >= 1e-12) for (let [e, t] of i) a.linearTerms.set(e, (a.linearTerms.get(e) || 0) + t * this.constant);
+		for (let [t, n] of this.linearTerms) for (let [r, o] of i) a = a.add(e.fromProduct(t, r, n * o));
+		return a;
 	}
 	div(e) {
 		if (Math.abs(e) < 1e-12) throw Error("Division by zero in QuadExpr.div()");
@@ -181,15 +208,31 @@ var b = class e {
 	getValue() {
 		return this.value;
 	}
+	toString() {
+		let e = [], t = !0;
+		for (let n of this.quadTerms.values()) {
+			let r = n.coeff, i = r >= 0 ? t ? "" : "+ " : "- ", a = Math.abs(r), o = a === 1 ? "" : `${a} `;
+			n.var1 === n.var2 ? e.push(`${i}${o}${n.var1.name}^2`) : e.push(`${i}${o}${n.var1.name} * ${n.var2.name}`), t = !1;
+		}
+		for (let [n, r] of this.linearTerms) {
+			let i = r >= 0 ? t ? "" : "+ " : "- ", a = Math.abs(r), o = a === 1 ? "" : `${a} `;
+			e.push(`${i}${o}${n.name}`), t = !1;
+		}
+		if (this.constant !== 0 || e.length === 0) {
+			let n = this.constant >= 0 ? t ? "" : "+ " : "- ";
+			e.push(`${n}${Math.abs(this.constant)}`);
+		}
+		return e.join(" ") || "0";
+	}
 	static fromProduct(t, n, r = 1) {
-		let i = new e(), a = y(t, n);
+		let i = new e(), a = x(t, n);
 		return i.quadTerms.set(a, {
 			var1: t.id <= n.id ? t : n,
 			var2: t.id <= n.id ? n : t,
 			coeff: r
 		}), i;
 	}
-}, x = class e {
+}, C = class e {
 	terms;
 	constant;
 	constructor(e = /* @__PURE__ */ new Map(), t = 0) {
@@ -199,7 +242,7 @@ var b = class e {
 		return new e(this.terms, this.constant);
 	}
 	add(e) {
-		if (e instanceof b) return e.add(this);
+		if (e instanceof S) return e.add(this);
 		let t = this.clone();
 		if (typeof e == "number") t.constant += e;
 		else if ("terms" in e) {
@@ -215,7 +258,7 @@ var b = class e {
 		return t;
 	}
 	sub(e) {
-		return e instanceof b ? e.neg().add(this) : typeof e == "number" ? this.add(-e) : ("terms" in e, this.add(e.mul(-1)));
+		return e instanceof S ? e.neg().add(this) : typeof e == "number" ? this.add(-e) : ("terms" in e, this.add(e.mul(-1)));
 	}
 	mul(t) {
 		if (typeof t == "number") {
@@ -226,9 +269,27 @@ var b = class e {
 			}
 			return n;
 		}
-		let n = new b();
-		Math.abs(this.constant) >= 1e-12 && n.linearTerms.set(t, this.constant);
-		for (let [e, r] of this.terms) n = n.add(b.fromProduct(e, t, r));
+		if (t instanceof D || typeof t == "object" && t && "id" in t && !("terms" in t) && !("quadTerms" in t)) {
+			let e = t, n = new S();
+			Math.abs(this.constant) >= 1e-12 && n.linearTerms.set(e, this.constant);
+			for (let [t, r] of this.terms) n = n.add(S.fromProduct(t, e, r));
+			return n;
+		}
+		if (t instanceof S) {
+			if (t.quadTerms.size > 0 && this.terms.size > 0) throw Error("Multiplication creates a non-linear term of degree > 2 (not allowed in MILP/LP/QP)");
+			if (this.terms.size === 0) return t.mul(this.constant);
+			let e = new S();
+			if (e.constant = this.constant * t.constant, Math.abs(t.constant) >= 1e-12) for (let [n, r] of this.terms) e.linearTerms.set(n, (e.linearTerms.get(n) || 0) + r * t.constant);
+			if (Math.abs(this.constant) >= 1e-12) for (let [n, r] of t.linearTerms) e.linearTerms.set(n, (e.linearTerms.get(n) || 0) + r * this.constant);
+			for (let [n, r] of this.terms) for (let [i, a] of t.linearTerms) e = e.add(S.fromProduct(n, i, r * a));
+			return e;
+		}
+		if (this.terms.size === 0) return t.mul(this.constant);
+		if (t.terms.size === 0) return this.mul(t.constant);
+		let n = new S();
+		if (n.constant = this.constant * t.constant, Math.abs(t.constant) >= 1e-12) for (let [e, r] of this.terms) n.linearTerms.set(e, (n.linearTerms.get(e) || 0) + r * t.constant);
+		if (Math.abs(this.constant) >= 1e-12) for (let [e, r] of t.terms) n.linearTerms.set(e, (n.linearTerms.get(e) || 0) + r * this.constant);
+		for (let [e, r] of this.terms) for (let [i, a] of t.terms) n = n.add(S.fromProduct(e, i, r * a));
 		return n;
 	}
 	div(e) {
@@ -271,7 +332,7 @@ var b = class e {
 		let i = this.clone(), a = 0;
 		typeof n == "number" ? a = n : "terms" in n ? (i = i.sub(new e(n.terms, 0)), a = n.constant) : i = i.sub(n);
 		let o = a - i.constant;
-		return i.constant = 0, new v(i, t, o, r);
+		return i.constant = 0, new b(i, t, o, r);
 	}
 	toString() {
 		let e = [], t = !0;
@@ -286,27 +347,113 @@ var b = class e {
 		return e.join(" ") || "0";
 	}
 };
-function S(e, t) {
-	let n = /* @__PURE__ */ new Map(), r = 0, i = 0;
-	for (let a of e) {
-		let e = t ? t(a, i++) : a;
-		if (typeof e == "number") r += e;
-		else if (e && typeof e == "object" && "terms" in e) {
-			let t = e;
-			for (let [e, r] of t.terms) {
-				let t = n.get(e) || 0;
-				n.set(e, t + r);
+function w(e, t) {
+	let n = /* @__PURE__ */ new Map(), r = 0;
+	function i(e) {
+		if (e != null) {
+			if (typeof e == "number") r += e;
+			else if (typeof e == "object" && "terms" in e) {
+				let t = e;
+				for (let [e, r] of t.terms) {
+					let t = n.get(e) || 0;
+					n.set(e, t + r);
+				}
+				r += t.constant;
+			} else {
+				let t = e, r = n.get(t) || 0;
+				n.set(t, r + 1);
 			}
-			r += t.constant;
-		} else if (e) {
-			let t = e, r = n.get(t) || 0;
-			n.set(t, r + 1);
 		}
 	}
+	if (typeof e == "number") {
+		let n = e, r = t;
+		for (let e = 0; e < n; e++) i(r(e, e));
+	} else {
+		let n = 0;
+		for (let r of e) i(t ? t(r, n++) : r);
+	}
 	for (let [e, t] of n) Math.abs(t) < 1e-12 && n.delete(e);
-	return new x(n, r);
+	return new C(n, r);
 }
-var C = S, w = S, T = 0, E = class e {
+function ee(e, t) {
+	let n = /* @__PURE__ */ new Map(), r = 0;
+	function i(e) {
+		return typeof e == "object" && !!e && typeof e.id == "number" && typeof e.name == "string" && typeof e.mul == "function" && typeof e.le == "function";
+	}
+	function a(e) {
+		return i(e) || e instanceof C || e instanceof S;
+	}
+	function o(e, t) {
+		if (typeof t == "number") r += e * t;
+		else if (t instanceof C) {
+			for (let [r, i] of t.terms) {
+				let t = n.get(r) || 0;
+				n.set(r, t + e * i);
+			}
+			r += e * t.constant;
+		} else if (i(t)) {
+			let r = t, i = n.get(r) || 0;
+			n.set(r, i + e);
+		} else throw Error("Cannot multiply number by non-variable/expression operand in dot().");
+	}
+	function s(e) {
+		return typeof e == "object" && !!e && !Array.isArray(e) && !a(e);
+	}
+	function c(e, t, n = "") {
+		if (e === void 0 || t === void 0 || e === null || t === null) throw Error(`Nullish or undefined operand encountered in dot() at '${n || "root"}'.`);
+		if (typeof e == "number") {
+			o(e, t);
+			return;
+		}
+		if (typeof t == "number") {
+			o(t, e);
+			return;
+		}
+		if (Array.isArray(e) || Array.isArray(t)) {
+			if (Array.isArray(e) && Array.isArray(t)) {
+				if (e.length !== t.length) throw Error(`Array length mismatch in dot() at '${n || "root"}': received lengths ${e.length} and ${t.length}.`);
+				for (let r = 0; r < e.length; r++) c(e[r], t[r], n ? `${n}[${r}]` : `[${r}]`);
+				return;
+			}
+			let r = Array.isArray(e) ? e : t, i = Array.isArray(e) ? t : e;
+			if (s(i)) {
+				let a = Object.keys(i);
+				if (a.length > 0 && a.every((e) => /^\d+$/.test(e))) {
+					if (a.length !== r.length) throw Error(`Array length mismatch in dot() at '${n || "root"}': array has length ${r.length}, but object has ${a.length} indices.`);
+					let o = [];
+					for (let e = 0; e < r.length; e++) !(e in i) && !(String(e) in i) && o.push(e);
+					if (o.length > 0) throw Error(`Key mismatch in dot() at '${n || "root"}': missing indices [${o.join(", ")}] in numeric record.`);
+					for (let i = 0; i < r.length; i++) c(e[i], t[i], n ? `${n}[${i}]` : `[${i}]`);
+					return;
+				}
+				throw Error(`Type mismatch in dot() at '${n || "root"}': one operand is an Array while the other is an Object with keys [${a.slice(0, 5).join(", ")}].`);
+			}
+			throw Error(`Type mismatch in dot() at '${n || "root"}': one operand is an Array while the other is not.`);
+		}
+		if (s(e) && s(t)) {
+			let r = Object.keys(e), i = Object.keys(t), a = r.filter((e) => !(e in t)), o = i.filter((t) => !(t in e));
+			if (a.length > 0 || o.length > 0) {
+				let e = [];
+				throw a.length > 0 && e.push(`missing in second operand: [${a.join(", ")}]`), o.length > 0 && e.push(`missing in first operand: [${o.join(", ")}]`), Error(`Key mismatch in dot() at '${n || "root"}': ${e.join("; ")}.`);
+			}
+			for (let i of r) c(e[i], t[i], n ? `${n}.${i}` : i);
+			return;
+		}
+		throw Error(`Unsupported operand types in dot() at '${n || "root"}'.`);
+	}
+	c(e, t);
+	for (let [e, t] of n) Math.abs(t) < 1e-12 && n.delete(e);
+	return new C(n, r);
+}
+//#endregion
+//#region src/modeler/variable.ts
+var T = 0;
+function E(e) {
+	if (!e) return "Continuous";
+	let t = e.toLowerCase();
+	return t === "binary" || t === "bin" ? "Binary" : t === "integer" || t === "int" ? "Integer" : "Continuous";
+}
+var D = class e {
 	id;
 	name;
 	type;
@@ -315,48 +462,49 @@ var C = S, w = S, T = 0, E = class e {
 	value = null;
 	reducedCost = null;
 	constructor(e = {}) {
-		this.id = ++T, this.name = e.name || `x_${this.id}`, this.type = e.type || "Continuous", this.type === "Binary" ? (this.lb = 0, this.ub = 1) : (this.lb = e.lb === void 0 ? 0 : e.lb, this.ub = e.ub === void 0 ? Infinity : e.ub);
-	}
-	get x() {
-		return this.value;
-	}
-	get varValue() {
-		return this.value;
+		this.id = ++T, this.name = e.name || `x_${this.id}`, this.type = E(e.type), this.type === "Binary" ? (this.lb = 0, this.ub = 1) : (this.lb = e.lb === void 0 ? 0 : e.lb, this.ub = e.ub === void 0 ? Infinity : e.ub);
 	}
 	mul(t) {
-		if (typeof t == "number") return new x(/* @__PURE__ */ new Map([[this, t]]), 0);
-		if (t instanceof e) return b.fromProduct(this, t, 1);
-		let n = new b();
+		if (typeof t == "number") return new C(/* @__PURE__ */ new Map([[this, t]]), 0);
+		if (t instanceof e) return S.fromProduct(this, t, 1);
+		if (t instanceof S) {
+			if (t.quadTerms.size > 0) throw Error("Multiplication creates a non-linear term of degree > 2 (not allowed in MILP/LP/QP)");
+			let e = new S();
+			Math.abs(t.constant) >= 1e-12 && e.linearTerms.set(this, t.constant);
+			for (let [n, r] of t.linearTerms) e = e.add(S.fromProduct(this, n, r));
+			return e;
+		}
+		let n = new S();
 		Math.abs(t.constant) >= 1e-12 && n.linearTerms.set(this, t.constant);
 		for (let [e, r] of t.terms) {
-			let t = b.fromProduct(this, e, r);
+			let t = S.fromProduct(this, e, r);
 			n = n.add(t);
 		}
 		return n;
 	}
 	add(e) {
-		return e instanceof b ? e.add(this) : new x(/* @__PURE__ */ new Map([[this, 1]]), 0).add(e);
+		return e instanceof S ? e.add(this) : new C(/* @__PURE__ */ new Map([[this, 1]]), 0).add(e);
 	}
 	sub(e) {
-		return e instanceof b ? e.neg().add(this) : new x(/* @__PURE__ */ new Map([[this, 1]]), 0).sub(e);
+		return e instanceof S ? e.neg().add(this) : new C(/* @__PURE__ */ new Map([[this, 1]]), 0).sub(e);
 	}
 	neg() {
 		return this.mul(-1);
 	}
 	le(e, t) {
-		return new x(/* @__PURE__ */ new Map([[this, 1]]), 0).le(e, t);
+		return new C(/* @__PURE__ */ new Map([[this, 1]]), 0).le(e, t);
 	}
 	lessThanOrEqual(e, t) {
 		return this.le(e, t);
 	}
 	ge(e, t) {
-		return new x(/* @__PURE__ */ new Map([[this, 1]]), 0).ge(e, t);
+		return new C(/* @__PURE__ */ new Map([[this, 1]]), 0).ge(e, t);
 	}
 	greaterThanOrEqual(e, t) {
 		return this.ge(e, t);
 	}
 	eq(e, t) {
-		return new x(/* @__PURE__ */ new Map([[this, 1]]), 0).eq(e, t);
+		return new C(/* @__PURE__ */ new Map([[this, 1]]), 0).eq(e, t);
 	}
 	equalTo(e, t) {
 		return this.eq(e, t);
@@ -367,7 +515,7 @@ var C = S, w = S, T = 0, E = class e {
 };
 //#endregion
 //#region src/modeler/parser.ts
-function te(e) {
+function O(e) {
 	let t = [], n = 0, r = e.length;
 	for (; n < r;) {
 		let i = e[n];
@@ -383,9 +531,16 @@ function te(e) {
 			type: "MINUS",
 			value: "-"
 		}), n++;
-		else if (i === "*") t.push({
+		else if (i === "*") n + 1 < r && e[n + 1] === "*" ? (t.push({
+			type: "CARET",
+			value: "**"
+		}), n += 2) : (t.push({
 			type: "STAR",
 			value: "*"
+		}), n++);
+		else if (i === "^") t.push({
+			type: "CARET",
+			value: "^"
 		}), n++;
 		else if (i === "/") t.push({
 			type: "SLASH",
@@ -423,7 +578,24 @@ function te(e) {
 	}
 	return t;
 }
-var ne = class {
+function k(e, t) {
+	return e instanceof S ? e.add(t) : t instanceof S ? t.add(e) : e.add(t);
+}
+function te(e, t) {
+	return e instanceof S ? e.sub(t) : t instanceof S ? t.neg().add(e) : e.sub(t);
+}
+function A(e, t) {
+	let n = e instanceof S ? e.quadTerms : /* @__PURE__ */ new Map(), r = t instanceof S ? t.quadTerms : /* @__PURE__ */ new Map(), i = e instanceof S ? e.linearTerms : e.terms, a = t instanceof S ? t.linearTerms : t.terms;
+	if (n.size > 0 && (r.size > 0 || a.size > 0) || r.size > 0 && (n.size > 0 || i.size > 0)) throw Error("Multiplication creates a non-linear term of degree > 2 (not allowed in MILP/LP/QP)");
+	if (n.size === 0 && i.size === 0) return t.mul(e.constant);
+	if (r.size === 0 && a.size === 0) return e.mul(t.constant);
+	let o = new S();
+	if (o.constant = e.constant * t.constant, Math.abs(t.constant) >= 1e-12) for (let [e, n] of i) o.linearTerms.set(e, (o.linearTerms.get(e) || 0) + n * t.constant);
+	if (Math.abs(e.constant) >= 1e-12) for (let [t, n] of a) o.linearTerms.set(t, (o.linearTerms.get(t) || 0) + n * e.constant);
+	for (let [e, t] of i) for (let [n, r] of a) o = o.add(S.fromProduct(e, n, t * r));
+	return o;
+}
+var j = class {
 	tokens;
 	pos = 0;
 	resolver;
@@ -444,11 +616,11 @@ var ne = class {
 		for (; this.pos < this.tokens.length;) {
 			let t = this.peek();
 			if (!t) break;
-			if (t.type === "PLUS") this.consume(), e = e.add(this.parseTerm());
-			else if (t.type === "MINUS") this.consume(), e = e.sub(this.parseTerm());
+			if (t.type === "PLUS") this.consume(), e = k(e, this.parseTerm());
+			else if (t.type === "MINUS") this.consume(), e = te(e, this.parseTerm());
 			else break;
 		}
-		return e;
+		return e instanceof S && e.quadTerms.size === 0 ? new C(e.linearTerms, e.constant) : e;
 	}
 	parseTerm() {
 		let e = this.parseFactor();
@@ -458,57 +630,83 @@ var ne = class {
 			if (t.type === "STAR") {
 				this.consume();
 				let t = this.parseFactor();
-				e = this.multiplyExpressions(e, t);
+				e = A(e, t);
 			} else if (t.type === "SLASH") {
 				this.consume();
-				let t = this.parseFactor();
-				if (t.terms.size > 0) throw Error("Cannot divide by variable in linear expression");
+				let t = this.parseFactor(), n = t instanceof S ? t.linearTerms : t.terms, r = t instanceof S ? t.quadTerms.size : 0;
+				if (n.size > 0 || r > 0) throw Error("Cannot divide by variable in expression");
 				if (Math.abs(t.constant) < 1e-12) throw Error("Division by zero in expression");
 				e = e.div(t.constant);
 			} else break;
 		}
 		return e;
 	}
-	multiplyExpressions(e, t) {
-		if (e.terms.size > 0 && t.terms.size > 0) throw Error("Multiplication of two variables creates a non-linear term (not allowed in MILP/LP)");
-		return e.terms.size === 0 ? t.mul(e.constant) : e.mul(t.constant);
-	}
 	parseFactor() {
 		let e = this.peek();
 		if (!e) throw Error("Unexpected end of expression");
 		if (e.type === "PLUS") return this.consume(), this.parseFactor();
 		if (e.type === "MINUS") return this.consume(), this.parseFactor().neg();
+		let t;
 		if (e.type === "NUMBER") {
 			this.consume();
-			let t = parseFloat(e.value);
-			if (this.peek()?.type === "IDENT") {
-				let e = this.consume();
-				return this.resolveVariable(e.value).mul(t);
+			let n = parseFloat(e.value);
+			if (this.peek()?.type === "IDENT" || this.peek()?.type === "LPAREN") {
+				let e;
+				if (this.peek()?.type === "IDENT") {
+					let t = this.consume(), n = this.resolve(t.value);
+					e = n instanceof D || "id" in n && !("terms" in n) && !("quadTerms" in n) ? new C(/* @__PURE__ */ new Map([[n, 1]]), 0) : n.clone();
+				} else this.consume("LPAREN"), e = this.parseExpression(), this.consume("RPAREN");
+				return this.applyPower(e).mul(n);
 			}
-			return new x(/* @__PURE__ */ new Map(), t);
-		}
-		if (e.type === "IDENT") {
+			t = new C(/* @__PURE__ */ new Map(), n);
+		} else if (e.type === "IDENT") {
+			this.consume();
+			let n = this.resolve(e.value);
+			t = n instanceof D || "id" in n && !("terms" in n) && !("quadTerms" in n) ? new C(/* @__PURE__ */ new Map([[n, 1]]), 0) : n.clone();
+		} else if (e.type === "LPAREN") this.consume("LPAREN"), t = this.parseExpression(), this.consume("RPAREN");
+		else throw Error(`Unexpected token '${e.value}' in expression`);
+		return this.applyPower(t);
+	}
+	applyPower(e) {
+		if (this.peek()?.type !== "CARET") return e;
+		let t = [];
+		for (; this.peek()?.type === "CARET";) {
+			this.consume();
+			let e = this.peek();
+			if (!e || e.type !== "NUMBER") throw Error("Exponent must be a non-negative integer number");
 			this.consume();
-			let t = this.resolveVariable(e.value);
-			return new x(/* @__PURE__ */ new Map([[t, 1]]), 0);
+			let n = parseFloat(e.value);
+			if (!Number.isInteger(n) || n < 0) throw Error("Exponent must be a non-negative integer");
+			t.push(n);
 		}
-		if (e.type === "LPAREN") {
-			this.consume("LPAREN");
-			let e = this.parseExpression();
-			return this.consume("RPAREN"), e;
+		let n = t[t.length - 1];
+		for (let e = t.length - 2; e >= 0; e--) n = t[e] ** +n;
+		if (!Number.isFinite(n) || n < 0) throw Error("Exponent creates an invalid or out-of-range power");
+		let r = e instanceof S, i = r ? e.quadTerms.size : 0, a = r ? e.linearTerms.size : e.terms.size;
+		if (i === 0 && a === 0) {
+			let t = e.constant ** +n;
+			return new C(/* @__PURE__ */ new Map(), t);
 		}
-		throw Error(`Unexpected token '${e.value}' in expression`);
+		if (n === 0) return new C(/* @__PURE__ */ new Map(), 1);
+		if (n === 1) return e;
+		if (n === 2) return A(e, e);
+		throw Error(`Exponent ${n} creates a non-linear term of degree > 2 (not allowed in MILP/LP/QP)`);
 	}
-	resolveVariable(e) {
+	getRemainingToken() {
+		return this.peek();
+	}
+	resolve(e) {
 		let t = this.resolver(e);
 		if (!t) throw Error(`Unknown variable '${e}' in expression`);
 		return t;
 	}
 };
-function D(e, t) {
-	return new ne(te(e), t).parseExpression();
+function M(e, t) {
+	let n = new j(O(e), t), r = n.parseExpression(), i = n.getRemainingToken();
+	if (i) throw Error(`Unexpected token '${i.value}' after expression in "${e}"`);
+	return r;
 }
-function O(e, t, n) {
+function N(e, t, n) {
 	let r = e.match(/(<=|>=|==|=<|=>|=)/);
 	if (!r) throw Error(`Constraint must contain a comparison operator (<=, >=, ==, =): "${e}"`);
 	let i = r[0], a = r.index, o = e.substring(0, a).trim(), s = e.substring(a + i.length).trim(), c;
@@ -516,20 +714,19 @@ function O(e, t, n) {
 	else if (i === ">=" || i === "=>") c = ">=";
 	else if (i === "==" || i === "=") c = "==";
 	else throw Error(`Unsupported operator "${i}"`);
-	let l = D(o, t), u = D(s, t);
-	return c === "<=" ? l.le(u, n) : c === ">=" ? l.ge(u, n) : l.eq(u, n);
+	let l = M(o, t), u = M(s, t);
+	if (l instanceof S && l.quadTerms.size > 0 || u instanceof S && u.quadTerms.size > 0) throw Error("Quadratic constraints are not supported (only linear constraints are allowed in MILP/LP/QP)");
+	let d = l instanceof S ? new C(l.linearTerms, l.constant) : l, f = u instanceof S ? new C(u.linearTerms, u.constant) : u;
+	return c === "<=" ? d.le(f, n) : c === ">=" ? d.ge(f, n) : d.eq(f, n);
 }
-function k(e, t) {
+function P(e, t) {
 	let n = /* @__PURE__ */ new Map(), r = "";
 	for (let i = 0; i < e.length; i++) if (r += e[i], i < t.length) {
 		let e = t[i];
 		if (typeof e == "number") r += `(${e})`;
-		else if (e instanceof E) {
-			let t = `__v_${i}__`;
+		else if (e instanceof D || e instanceof C || e instanceof S || typeof e == "object" && e && ("id" in e || "terms" in e || "quadTerms" in e)) {
+			let t = `__ph_${i}__`;
 			n.set(t, e), r += t;
-		} else if (e instanceof x) {
-			`${i}`, r += `(${e.toString()})`;
-			for (let [t] of e.terms) n.set(t.name, t);
 		}
 	}
 	return {
@@ -537,42 +734,42 @@ function k(e, t) {
 		resolver: (e) => n.get(e)
 	};
 }
-function re(e, ...t) {
-	let { fullText: n, resolver: r } = k(e, t);
-	return O(n, r);
+function F(e, ...t) {
+	let { fullText: n, resolver: r } = P(e, t);
+	return N(n, r);
 }
-function ie(e, ...t) {
-	let { fullText: n, resolver: r } = k(e, t);
-	return /(<=|>=|==|=<|=>|=)/.test(n) ? O(n, r) : D(n, r);
+function I(e, ...t) {
+	let { fullText: n, resolver: r } = P(e, t);
+	return /(<=|>=|==|=<|=>|=)/.test(n) ? N(n, r) : M(n, r);
 }
 //#endregion
 //#region src/modeler/lp-writer.ts
-function A(e, t, n) {
+function L(e, t, n) {
 	let r = e.replace(/[^a-zA-Z0-9_]/g, "_");
 	return /^[0-9]/.test(r) && (r = `${t}_${r}`), (!r || r.length === 0) && (r = `${t}_${n}`), r;
 }
-function j(e) {
+function R(e) {
 	let { sense: t, objective: n, constraints: r, variables: i } = e, a = /* @__PURE__ */ new Map(), o = /* @__PURE__ */ new Map(), s = /* @__PURE__ */ new Map(), c = /* @__PURE__ */ new Set();
 	for (let e of i) {
-		let t = A(e.name, "x", e.id), n = t, r = 1;
+		let t = L(e.name, "x", e.id), n = t, r = 1;
 		for (; c.has(n);) n = `${t}_${r++}`;
 		c.add(n), o.set(e, n), a.set(n, e);
 	}
 	let l = /* @__PURE__ */ new Set(), u = /* @__PURE__ */ new Map();
 	for (let e of r) {
-		let t = A(e.name, "c", e.id), n = t, r = 1;
+		let t = L(e.name, "c", e.id), n = t, r = 1;
 		for (; l.has(n);) n = `${t}_${r++}`;
 		l.add(n), u.set(e, n), s.set(n, e);
 	}
 	let d = [];
-	d.push(t === "Maximize" ? "Maximize" : "Minimize");
-	let f = [], p = !0, m = n instanceof b ? n.linearTerms : n.terms;
+	d.push(t && t.toLowerCase() === "maximize" ? "Maximize" : "Minimize");
+	let f = [], p = !0, m = n instanceof S ? n.linearTerms : n.terms;
 	if (m) for (let [e, t] of m) {
 		if (Math.abs(t) < 1e-12) continue;
 		let n = o.get(e), r = t >= 0 ? p ? "" : "+ " : "- ", i = Math.abs(t), a = i === 1 ? "" : `${i} `;
 		f.push(`${r}${a}${n}`), p = !1;
 	}
-	if (n instanceof b && n.quadTerms.size > 0) {
+	if (n instanceof S && n.quadTerms.size > 0) {
 		let e = [], t = !0;
 		for (let r of n.quadTerms.values()) {
 			if (Math.abs(r.coeff) < 1e-12) continue;
@@ -624,38 +821,38 @@ function j(e) {
 }
 //#endregion
 //#region src/modeler/solvers/registry.ts
-var M = Symbol.for("washi.solverRegistry"), N = Symbol.for("washi.defaultSolverFactory"), P = globalThis;
-P[M] || (P[M] = /* @__PURE__ */ new Map());
-var F = P[M];
-function I(e, t) {
-	F.set(e.toLowerCase(), t);
+var z = Symbol.for("washi.solverRegistry"), B = Symbol.for("washi.defaultSolverFactory"), V = globalThis;
+V[z] || (V[z] = /* @__PURE__ */ new Map());
+var H = V[z];
+function U(e, t) {
+	H.set(e.toLowerCase(), t);
 }
-function L(e) {
+function W(e) {
 	if (typeof e == "string") {
-		let t = F.get(e.toLowerCase());
+		let t = H.get(e.toLowerCase());
 		if (!t) throw Error(`Solver '${e}' is not registered.`);
-		P[N] = t;
-	} else P[N] = typeof e == "function" ? e : () => e;
+		V[B] = t;
+	} else V[B] = typeof e == "function" ? e : () => e;
 }
-function R(e, t) {
+function G(e, t) {
 	if (!e) {
-		let e = P[N] || F.get("highs");
-		if (!e) throw Error(`No default solver available. Available solvers: ${z().join(", ")}`);
+		let e = V[B] || H.get("highs");
+		if (!e) throw Error(`No default solver available. Available solvers: ${K().join(", ")}`);
 		return e(t);
 	}
 	if (typeof e == "string") {
-		let n = e.toLowerCase(), r = F.get(n);
-		if (!r) throw Error(n === "glpk" ? `Solver 'glpk' not found. Available solvers: ${z().join(", ")}. Please import '@ohtaman/washi/solvers/glpk' before selecting this solver.` : n === "z3" ? `Solver 'z3' not found. Available solvers: ${z().join(", ")}. Please import '@ohtaman/washi/solvers/z3' before selecting this solver.` : `Solver '${e}' not found. Available solvers: ${z().join(", ")}`);
+		let n = e.toLowerCase(), r = H.get(n);
+		if (!r) throw Error(n === "glpk" ? `Solver 'glpk' not found. Available solvers: ${K().join(", ")}. Please import '@ohtaman/washi/solvers/glpk' before selecting this solver.` : n === "z3" ? `Solver 'z3' not found. Available solvers: ${K().join(", ")}. Please import '@ohtaman/washi/solvers/z3' before selecting this solver.` : `Solver '${e}' not found. Available solvers: ${K().join(", ")}`);
 		return r(t);
 	}
 	return e;
 }
-function z() {
-	return Array.from(F.keys());
+function K() {
+	return Array.from(H.keys());
 }
 //#endregion
 //#region \0washi-browser-worker-factory
-function ae() {
+function q() {
 	let e = r("highs-worker.js");
 	if (new URL(e).origin === window.location.origin) return new Worker(e, { type: "module" });
 	let t = URL.createObjectURL(new Blob(["import " + JSON.stringify(e) + ";"], { type: "text/javascript" }));
@@ -668,41 +865,41 @@ function ae() {
 }
 //#endregion
 //#region src/modeler/solvers/highs/worker-client.ts
-var B = null, V = 0, H = /* @__PURE__ */ new Map();
-function U() {
+var J = null, Y = 0, X = /* @__PURE__ */ new Map();
+function Z() {
 	if (typeof window > "u" || typeof Worker > "u") return null;
-	if (!B) try {
-		B = ae();
-		let e = B;
-		B.onmessage = (e) => {
-			let { type: t, id: n, result: r, runTime: i, error: a } = e.data, o = H.get(n);
-			o && (H.delete(n), t === "SOLVE_SUCCESS" ? o.resolve({
+	if (!J) try {
+		J = q();
+		let e = J;
+		J.onmessage = (e) => {
+			let { type: t, id: n, result: r, runTime: i, error: a } = e.data, o = X.get(n);
+			o && (X.delete(n), t === "SOLVE_SUCCESS" ? o.resolve({
 				result: r,
 				runTime: i
 			}) : o.reject(Error(a || "Worker solve failed")));
 		};
 		let t = (t) => {
-			if (B === e) {
-				B = null, e.terminate();
-				for (let e of H.values()) e.reject(t);
-				H.clear();
+			if (J === e) {
+				J = null, e.terminate();
+				for (let e of X.values()) e.reject(t);
+				X.clear();
 			}
 		};
-		B.onerror = (e) => t(/* @__PURE__ */ Error(`[WaSHI] HiGHS Worker failed to load or execute. Check Worker assets, CORS and CSP. ${e?.message || ""}`)), B.onmessageerror = () => t(/* @__PURE__ */ Error("[WaSHI] Could not deserialize the HiGHS Worker response."));
+		J.onerror = (e) => t(/* @__PURE__ */ Error(`[WaSHI] HiGHS Worker failed to load or execute. Check Worker assets, CORS and CSP. ${e?.message || ""}`)), J.onmessageerror = () => t(/* @__PURE__ */ Error("[WaSHI] Could not deserialize the HiGHS Worker response."));
 	} catch {
-		B = null;
+		J = null;
 	}
-	return B;
+	return J;
 }
-async function W(e, n = {}) {
-	let r = await a(n.wasmUrl), i = n.useWorker ?? typeof window < "u" ? U() : null;
+async function ne(e, n = {}) {
+	let r = await a(n.wasmUrl), i = n.useWorker ?? typeof window < "u" ? Z() : null;
 	if (!i) return t(e, n);
 	let o = {
 		...n,
 		wasmUrl: r
-	}, s = `req_${++V}`;
+	}, s = `req_${++Y}`;
 	return new Promise((t, n) => {
-		H.set(s, {
+		X.set(s, {
 			resolve: t,
 			reject: n
 		});
@@ -714,13 +911,13 @@ async function W(e, n = {}) {
 				options: o
 			});
 		} catch (e) {
-			H.delete(s), n(e);
+			X.delete(s), n(e);
 		}
 	});
 }
 //#endregion
 //#region src/modeler/solvers/highs/solver.ts
-function G(e) {
+function re(e) {
 	if (typeof e != "string") return "Unknown";
 	let t = e.trim();
 	if (!t) return "Unknown";
@@ -749,7 +946,7 @@ function G(e) {
 	let n = t.toLowerCase();
 	return n.includes("infeasible") && n.includes("unbounded") ? "Unknown" : n.includes("infeasible") ? "Infeasible" : n.includes("unbounded") ? "Unbounded" : n.includes("time limit") || n.includes("timelimit") ? "TimeLimit" : n.includes("iteration limit") || n.includes("iterationlimit") ? "IterationLimit" : n.includes("error") ? "Error" : n === "optimal" ? "Optimal" : "Unknown";
 }
-var K = class e {
+var Q = class e {
 	static capabilities = {
 		supportsMIP: !0,
 		supportsQP: !0
@@ -765,7 +962,18 @@ var K = class e {
 		let r = {
 			...this.options,
 			...n
-		}, { lpString: i, nameToVar: a, nameToConstraint: o } = t, { result: s, runTime: c } = await W(i, r), l = G(s?.Status), u = s.ObjectiveValue, d = typeof u == "number" ? u + e.objective.constant : null, f = {}, p = {}, m = {}, h = {};
+		}, { lpString: i, nameToVar: a, nameToConstraint: o } = t, { result: s, runTime: c } = await ne(i, r), l = re(s?.Status);
+		if (l === "Infeasible" || l === "Unbounded") return {
+			status: l,
+			objectiveValue: null,
+			values: {},
+			reducedCosts: {},
+			duals: {},
+			slacks: {},
+			runTime: c,
+			raw: s
+		};
+		let u = s.ObjectiveValue, d = typeof u == "number" ? u + e.objective.constant : null, f = {}, p = {}, m = {}, h = {};
 		if (s.Columns) for (let [e, t] of Object.entries(s.Columns)) {
 			let n = a.get(e);
 			if (n) {
@@ -794,14 +1002,14 @@ var K = class e {
 			raw: s
 		};
 	}
-}, q = K;
-function oe() {
-	I("highs", (e) => new K(e)), L((e) => new K(e));
+}, ie = Q;
+function ae() {
+	U("highs", (e) => new Q(e)), W((e) => new Q(e));
 }
-oe();
+ae();
 //#endregion
 //#region src/modeler/solvers/mock.ts
-var J = class {
+var $ = class {
 	name;
 	capabilities;
 	options;
@@ -828,13 +1036,13 @@ var J = class {
 		};
 	}
 };
-function se() {
-	I("mock", (e) => new J(e));
+function oe() {
+	U("mock", (e) => new $(e));
 }
-se();
+oe();
 //#endregion
 //#region src/modeler/model.ts
-var Y = class {
+var se = class {
 	name;
 	sense;
 	objective;
@@ -845,6 +1053,8 @@ var Y = class {
 	unsatCore = [];
 	values = /* @__PURE__ */ new Map();
 	varMapByName = /* @__PURE__ */ new Map();
+	boolVarMapByName = /* @__PURE__ */ new Map();
+	boolVarPrefixCounters = /* @__PURE__ */ new Map();
 	constrMapByName = /* @__PURE__ */ new Map();
 	status = "NotSolved";
 	objectiveValue = null;
@@ -856,7 +1066,7 @@ var Y = class {
 		n = typeof e == "string" ? {
 			...t,
 			name: e
-		} : e, this.name = n.name || "OptimizationModel", this.sense = n.sense || "Minimize", this.objective = new x(), n.solver && (this.solver = R(n.solver));
+		} : e, this.name = n.name || "OptimizationModel", this.sense = l(n.sense), this.objective = new C(), n.solver && (this.solver = G(n.solver));
 	}
 	get vars() {
 		return this.variables;
@@ -870,70 +1080,120 @@ var Y = class {
 	get numConstrs() {
 		return this.constraints.length;
 	}
-	addVar(e, t, n, r) {
-		let i;
-		i = typeof e == "string" ? {
-			name: e,
-			type: t,
-			lb: n,
-			ub: r
-		} : e || {};
-		let a = new E(i);
-		if (this.varMapByName.has(a.name)) throw Error(`Variable with name '${a.name}' already exists in model`);
-		return this.variables.push(a), this.varMapByName.set(a.name, a), a;
-	}
-	addVars(e, t, n) {
+	createDomainVars(e, t, n, r) {
+		if (typeof n == "number" || Array.isArray(n)) {
+			let i = typeof t == "number" ? Array.from({ length: t }, (e, t) => t) : t, a = typeof n == "number" ? Array.from({ length: n }, (e, t) => t) : n, o = r || {}, s = o.prefix || o.name || "x", c = {};
+			for (let t of i) {
+				c[t] = {};
+				for (let n of a) c[t][n] = this.addVar({
+					...o,
+					type: e,
+					name: `${s}_${t}_${n}`
+				});
+			}
+			return c;
+		}
+		if (typeof t == "string") {
+			let r = n || {};
+			return this.addVar({
+				...r,
+				type: e,
+				name: t
+			});
+		}
+		let i = n || {}, a = i.prefix || i.name;
+		if (typeof t == "number") {
+			let n = [], r = a || "x";
+			for (let a = 0; a < t; a++) n.push(this.addVar({
+				...i,
+				type: e,
+				name: `${r}_${a}`
+			}));
+			return n;
+		}
+		if (Array.isArray(t) || t && typeof t[Symbol.iterator] == "function") {
+			let n = {};
+			for (let r of t) {
+				let t = a ? `${a}_${r}` : `${r}`;
+				n[r] = this.addVar({
+					...i,
+					type: e,
+					name: t
+				});
+			}
+			return n;
+		}
+		throw Error(`Invalid arguments provided for ${e.toLowerCase()} variable creation`);
+	}
+	binary(e, t, n) {
+		return this.createDomainVars("Binary", e, t, n);
+	}
+	integer(e, t, n) {
+		return this.createDomainVars("Integer", e, t, n);
+	}
+	continuous(e, t, n) {
+		return this.createDomainVars("Continuous", e, t, n);
+	}
+	bool(e, t, n) {
 		if (typeof t == "number" || Array.isArray(t)) {
-			let r = typeof e == "number" ? Array.from({ length: e }, (e, t) => t) : e, i = typeof t == "number" ? Array.from({ length: t }, (e, t) => t) : t, a = n || {}, o = a.prefix || a.name || "x", s = {};
+			let r = typeof e == "number" ? Array.from({ length: e }, (e, t) => t) : e, i = typeof t == "number" ? Array.from({ length: t }, (e, t) => t) : t, a = n || {}, o = a.prefix || a.name || "p", s = {};
 			for (let e of r) {
 				s[e] = {};
-				for (let t of i) s[e][t] = this.addVar({
-					...a,
-					name: `${o}_${e}_${t}`
-				});
+				for (let t of i) s[e][t] = this.createBoolVar(`${o}_${e}_${t}`);
 			}
 			return s;
 		}
+		if (typeof e == "string" || e === void 0) return this.createBoolVar(e);
 		let r = t || {}, i = r.prefix || r.name;
 		if (typeof e == "number") {
-			let t = [], n = i || "x";
-			for (let i = 0; i < e; i++) t.push(this.addVar({
-				...r,
-				name: `${n}_${i}`
-			}));
+			let t = [], n = i || "p";
+			for (let r = 0; r < e; r++) t.push(this.createBoolVar(this.getNextBoolVarName(n)));
 			return t;
 		}
-		{
+		if (Array.isArray(e) || e && typeof e[Symbol.iterator] == "function") {
 			let t = {};
 			for (let n of e) {
 				let e = i ? `${i}_${n}` : `${n}`;
-				t[n] = this.addVar({
-					...r,
-					name: e
-				});
+				t[n] = this.createBoolVar(e);
 			}
 			return t;
 		}
+		throw Error("Invalid arguments provided for boolean variable creation");
 	}
-	addVars2D(e, t, n) {
-		let r = {}, i = n?.name || "x";
-		for (let a of e) {
-			r[a] = {};
-			for (let e of t) r[a][e] = this.addVar({
-				...n,
-				name: `${i}_${a}_${e}`
-			});
-		}
-		return r;
+	hasVarName(e) {
+		return this.varMapByName.has(e) || this.boolVarMapByName.has(e);
+	}
+	getNextBoolVarName(e) {
+		let t = this.boolVarPrefixCounters.get(e) ?? 0, n = `${e}_${t}`;
+		for (; this.hasVarName(n);) t++, n = `${e}_${t}`;
+		return this.boolVarPrefixCounters.set(e, t + 1), n;
+	}
+	createBoolVar(e) {
+		let t = typeof e == "string" ? e.trim() : void 0, n = t && t.length > 0 ? t : this.getNextBoolVarName("p");
+		if (this.hasVarName(n)) throw Error(`Variable with name '${n}' already exists in model`);
+		let r = new f(n);
+		return this.boolVars.push(r), this.boolVarMapByName.set(r.name, r), r;
+	}
+	addVar(e) {
+		let t = e || {}, n = typeof t.name == "string" ? t.name.trim() : void 0, r = n && n.length > 0 ? n : void 0, i = new D({
+			...t,
+			name: r
+		});
+		if (this.hasVarName(i.name)) throw Error(`Variable with name '${i.name}' already exists in model`);
+		return this.variables.push(i), this.varMapByName.set(i.name, i), i;
 	}
 	getVarByName(e) {
 		return this.varMapByName.get(e);
 	}
+	getBoolVarByName(e) {
+		return this.boolVarMapByName.get(e);
+	}
 	getConstraintByName(e) {
 		return this.constrMapByName.get(e);
 	}
 	setObjective(e, t) {
-		t && (this.sense = t), this.objective = typeof e == "string" ? D(e, (e) => this.getVarByName(e)) : typeof e == "number" ? new x(/* @__PURE__ */ new Map(), e) : e instanceof b || "terms" in e ? e.clone() : new x(/* @__PURE__ */ new Map([[e, 1]]), 0);
+		if (e instanceof b) throw Error("Cannot set objective to a constraint comparison (<=, >=, ==). Did you mean to use m.constrain(expr`... <= ...`)? For objectives, supply an unconstrained expression like expr`45 * ${x} + 80 * ${y}`.");
+		t && (this.sense = l(t)), this.objective = typeof e == "string" ? M(e, (e) => this.getVarByName(e)) : typeof e == "number" ? new C(/* @__PURE__ */ new Map(), e) : e instanceof S || "terms" in e ? e.clone() : new C(/* @__PURE__ */ new Map([[e, 1]]), 0);
 	}
 	maximize(e) {
 		this.setObjective(e, "Maximize");
@@ -941,12 +1201,30 @@ var Y = class {
 	minimize(e) {
 		this.setObjective(e, "Minimize");
 	}
-	addBoolVar(e) {
-		let t = new d(e);
-		return this.boolVars.push(t), t;
+	constrain(e, t) {
+		if (Array.isArray(e)) {
+			for (let t of e) this.constrain(t);
+			return this;
+		}
+		if (e instanceof S) throw Error("Cannot add a quadratic expression without a comparison operator (<=, >=, ==) as a constraint. Quadratic constraints are not supported (only linear constraints are allowed in MILP/LP/QP).");
+		if (e instanceof C) throw Error("Cannot add a linear expression without a comparison operator (<=, >=, ==) as a constraint. Use expr`... <= ...` or .le(), .ge(), .eq().");
+		return this.addConstraint(e, t), this;
+	}
+	getValues(e) {
+		if (e == null) return null;
+		if (e instanceof D || e instanceof f) return e.value;
+		if (Array.isArray(e)) return e.map((e) => this.getValues(e));
+		if (typeof e == "object") {
+			let t = {};
+			for (let [n, r] of Object.entries(e)) t[n] = this.getValues(r);
+			return t;
+		}
+		return null;
 	}
 	addConstraint(e, t) {
-		if (e instanceof u) {
+		if (e instanceof S) throw Error("Cannot add a quadratic expression without a comparison operator (<=, >=, ==) as a constraint. Quadratic constraints are not supported (only linear constraints are allowed in MILP/LP/QP).");
+		if (e instanceof C) throw Error("Cannot add a linear expression without a comparison operator (<=, >=, ==) as a constraint. Use expr`... <= ...` or .le(), .ge(), .eq().");
+		if (e instanceof d) {
 			let n = t || `lc_${this.logicalConstraints.length + 1}`;
 			return this.logicalConstraints.push({
 				name: n,
@@ -954,47 +1232,43 @@ var Y = class {
 			}), e;
 		}
 		let n;
-		if (typeof e == "string" ? n = O(e, (e) => this.getVarByName(e), t) : (n = e, t && (n.name = t)), this.constrMapByName.has(n.name)) throw Error(`Constraint with name '${n.name}' already exists in model`);
+		if (typeof e == "string" ? n = N(e, (e) => this.getVarByName(e), t) : (n = e, t && (n.name = t)), this.constrMapByName.has(n.name)) throw Error(`Constraint with name '${n.name}' already exists in model`);
 		return this.constraints.push(n), this.constrMapByName.set(n.name, n), n;
 	}
-	addConstr(e, t) {
-		return this.addConstraint(e, t);
-	}
-	addConstrs(e) {
-		return e.map((e) => this.addConstraint(e));
-	}
 	add(e, t) {
 		if (Array.isArray(e)) {
 			for (let t of e) this.add(t);
 			return this;
 		}
-		return e instanceof E ? (this.varMapByName.has(e.name) || (this.variables.push(e), this.varMapByName.set(e.name, e)), this) : (this.addConstraint(e, t), this);
-	}
-	addIndicator(e, t, n = {}) {
-		let r = n.activeValue === void 0 ? 1 : n.activeValue, i = n.M === void 0 ? 1e5 : n.M, a = n.name || `ind_${e.name}_${t.name}`, o = [], { lhs: s, sense: c, rhs: l } = t;
-		if (r === 1) {
-			if (c === "<=" || c === "==") {
-				let t = this.addConstraint(s.add(e.mul(i)).le(l + i), c === "==" ? `${a}_le` : a);
-				o.push(t);
+		return e instanceof D ? (this.varMapByName.has(e.name) || (this.variables.push(e), this.varMapByName.set(e.name, e)), this) : (this.addConstraint(e, t), this);
+	}
+	indicator(e, t, n = {}) {
+		if (t instanceof S) throw Error("Quadratic expressions or constraints are not supported in indicator constraints (only linear constraints are allowed).");
+		if (!("sense" in t) || !("lhs" in t)) throw Error("Invalid constraint expression passed to indicator(). Expression must include a comparison operator (<=, >=, ==).");
+		let r = t, i = n.activeValue === void 0 ? 1 : n.activeValue, a = n.M === void 0 ? 1e5 : n.M, o = n.name || `ind_${e.name}_${r.name}`, s = [], { lhs: c, sense: l, rhs: u } = r;
+		if (i === 1) {
+			if (l === "<=" || l === "==") {
+				let t = this.addConstraint(c.add(e.mul(a)).le(u + a), l === "==" ? `${o}_le` : o);
+				s.push(t);
 			}
-			if (c === ">=" || c === "==") {
-				let t = this.addConstraint(s.sub(e.mul(i)).ge(l - i), c === "==" ? `${a}_ge` : a);
-				o.push(t);
+			if (l === ">=" || l === "==") {
+				let t = this.addConstraint(c.sub(e.mul(a)).ge(u - a), l === "==" ? `${o}_ge` : o);
+				s.push(t);
 			}
 		} else {
-			if (c === "<=" || c === "==") {
-				let t = this.addConstraint(s.sub(e.mul(i)).le(l), c === "==" ? `${a}_le` : a);
-				o.push(t);
+			if (l === "<=" || l === "==") {
+				let t = this.addConstraint(c.sub(e.mul(a)).le(u), l === "==" ? `${o}_le` : o);
+				s.push(t);
 			}
-			if (c === ">=" || c === "==") {
-				let t = this.addConstraint(s.add(e.mul(i)).ge(l), c === "==" ? `${a}_ge` : a);
-				o.push(t);
+			if (l === ">=" || l === "==") {
+				let t = this.addConstraint(c.add(e.mul(a)).ge(u), l === "==" ? `${o}_ge` : o);
+				s.push(t);
 			}
 		}
-		return o;
+		return s;
 	}
 	toLpString() {
-		let { lpString: e } = j({
+		let { lpString: e } = R({
 			name: this.name,
 			sense: this.sense,
 			objective: this.objective,
@@ -1005,7 +1279,8 @@ var Y = class {
 	}
 	async optimize(e, t) {
 		let n, r;
-		if (e && typeof e == "object" && "solve" in e && typeof e.solve == "function" ? (n = e, r = t || {}) : typeof e == "function" ? (n = new e(t), r = t || {}) : typeof e == "string" ? (n = R(e, t), r = t || {}) : (r = { ...e }, r.solver ? (n = R(r.solver, r), delete r.solver) : n = this.solver || R()), this.objective instanceof b && this.objective.quadTerms.size > 0 && n.capabilities?.supportsQP === !1) throw Error(`Solver "${n.name}" does not support quadratic objectives (QP). Please use a QP-capable solver such as HiGHS or Clarabel.`);
+		if (e && typeof e == "object" && "solve" in e && typeof e.solve == "function" ? (n = e, r = t || {}) : typeof e == "function" ? (n = new e(t), r = t || {}) : typeof e == "string" ? (n = G(e, t), r = t || {}) : (r = { ...e }, r.solver ? (n = G(r.solver, r), delete r.solver) : n = this.solver || G()), this.objective instanceof S && this.objective.quadTerms.size > 0 && n.capabilities?.supportsQP === !1) throw Error(`Solver "${n.name}" does not support quadratic objectives (QP). Please use a QP-capable solver such as HiGHS or Clarabel.`);
+		if (this.logicalConstraints.length > 0 && n.capabilities?.supportsSMT !== !0) throw Error(`Solver "${n.name}" does not support SMT / logical constraints (or, implies, not, xor). Please use an SMT-capable solver such as Z3.`);
 		let i = {
 			name: this.name,
 			sense: this.sense,
@@ -1014,13 +1289,23 @@ var Y = class {
 			variables: this.variables,
 			logicalConstraints: this.logicalConstraints,
 			boolVars: this.boolVars
-		}, a = j(i), o = await n.solve(i, a, r);
+		}, a = R(i), o = await n.solve(i, a, r);
 		this.runTime = o.runTime ?? 0, this.status = o.status, this.objectiveValue = o.objectiveValue ?? null;
 		let s = {}, c = {}, l = {}, u = {};
-		if (this.values.clear(), o.values) for (let [e, t] of Object.entries(o.values)) {
-			this.values.set(e, t);
-			let n = this.getVarByName(e) || a.nameToVar.get(e);
-			n && (n.value = t, s[n.name] = t);
+		this.values.clear();
+		for (let e of this.variables) e.value = null, e.reducedCost = null;
+		for (let e of this.boolVars) e.value = null;
+		for (let e of this.constraints) e.dual = null, e.slack = null;
+		if (a?.nameToVar) for (let e of a.nameToVar.values()) e.value = null, e.reducedCost = null;
+		if (a?.nameToConstraint) for (let e of a.nameToConstraint.values()) e.dual = null, e.slack = null;
+		if (o.values && o.status !== "Infeasible" && o.status !== "Unbounded") {
+			for (let [e, t] of Object.entries(o.values)) if (typeof t == "number" && !Number.isNaN(t)) {
+				this.values.set(e, t);
+				let n = this.getVarByName(e) || a.nameToVar.get(e);
+				n && (n.value = t, s[n.name] = t);
+				let r = this.boolVars.find((t) => t.name === e);
+				r && (r.value = !!t);
+			}
 		}
 		if (o.reducedCosts) for (let [e, t] of Object.entries(o.reducedCosts)) {
 			let n = this.getVarByName(e) || a.nameToVar.get(e);
@@ -1050,7 +1335,7 @@ var Y = class {
 	}
 	async check(e, t) {
 		let n, r;
-		if (e && typeof e == "object" && "solve" in e && typeof e.solve == "function" ? (n = e, r = t || {}) : typeof e == "function" ? (n = new e(t), r = t || {}) : typeof e == "string" ? (n = R(e, t), r = t || {}) : (r = { ...e }, r.solver ? (n = R(r.solver, r), delete r.solver) : n = this.solver || R()), this.logicalConstraints.length > 0 && n.capabilities?.supportsSMT !== !0) throw Error(`Solver "${n.name}" does not support SMT / logical constraints (or, implies, not, xor). Please use an SMT-capable solver such as Z3.`);
+		if (e && typeof e == "object" && "solve" in e && typeof e.solve == "function" ? (n = e, r = t || {}) : typeof e == "function" ? (n = new e(t), r = t || {}) : typeof e == "string" ? (n = G(e, t), r = t || {}) : (r = { ...e }, r.solver ? (n = G(r.solver, r), delete r.solver) : n = this.solver || G()), this.logicalConstraints.length > 0 && n.capabilities?.supportsSMT !== !0) throw Error(`Solver "${n.name}" does not support SMT / logical constraints (or, implies, not, xor). Please use an SMT-capable solver such as Z3.`);
 		if (typeof n.check == "function") {
 			let e = await n.check({
 				name: this.name,
@@ -1059,21 +1344,26 @@ var Y = class {
 				constraints: this.constraints,
 				logicalConstraints: this.logicalConstraints
 			}, r);
-			if (this.status = e.status, this.runTime = e.runTime ?? 0, this.unsatCore = e.unsatCore || [], e.values) {
-				this.values.clear();
-				for (let [t, n] of Object.entries(e.values)) {
-					this.values.set(t, n);
-					let e = this.getVarByName(t);
-					e && typeof n == "number" && (e.value = n);
-				}
+			this.status = e.status, this.runTime = e.runTime ?? 0, this.unsatCore = e.unsatCore || [], this.values.clear();
+			for (let e of this.variables) e.value = null;
+			for (let e of this.boolVars) e.value = null;
+			if (e.values && e.status === "Sat") for (let [t, n] of Object.entries(e.values)) {
+				this.values.set(t, n);
+				let e = this.getVarByName(t);
+				e && typeof n == "number" && (e.value = n);
+				let r = this.boolVars.find((e) => e.name === t);
+				r && (r.value = !!n);
 			}
 			return e;
 		}
 		let i = this.objective, a = this.sense;
 		try {
-			this.objective = new x(), this.sense = "Minimize";
+			this.objective = new C(), this.sense = "Minimize";
 			let e = await this.optimize(n, r), t = e === "Optimal" || e === "Feasible", i = e === "Infeasible", a = {};
-			for (let e of this.variables) a[e.name] = e.value ?? 0;
+			if (t) {
+				for (let e of this.variables) e.value !== null && (a[e.name] = e.value);
+				for (let e of this.boolVars) e.value !== null && (a[e.name] = e.value);
+			}
 			return {
 				status: t ? "Sat" : i ? "Unsat" : "Unknown",
 				values: a,
@@ -1083,131 +1373,6 @@ var Y = class {
 			this.objective = i, this.sense = a;
 		}
 	}
-}, X = "Minimize", ce = "Maximize", le = "Continuous", ue = "Integer", de = "Binary", Z = {
-	Optimal: "Optimal",
-	Feasible: "Feasible",
-	Infeasible: "Infeasible",
-	Unbounded: "Unbounded",
-	TimeLimit: "Not Solved (Time Limit)",
-	IterationLimit: "Not Solved (Iteration Limit)",
-	Error: "Error",
-	NotSolved: "Not Solved",
-	Sat: "Sat",
-	Unsat: "Unsat",
-	Unknown: "Unknown"
-}, fe = class e extends E {
-	constructor(e, t, n, r = "Continuous") {
-		super({
-			name: e,
-			lb: t === void 0 ? 0 : t,
-			ub: n === void 0 ? Infinity : n,
-			type: r
-		});
-	}
-	static dicts(t, n, r, i, a = "Continuous") {
-		let o = {};
-		for (let s of n) {
-			let n = `${t}_${s}`;
-			o[s] = new e(n, r, i, a);
-		}
-		return o;
-	}
-};
-function pe(e) {
-	return e == null ? null : typeof e == "number" ? e : e.value;
-}
-var me = class {
-	model;
-	constructor(e = "NoName", t = X) {
-		this.model = new Y({
-			name: e,
-			sense: t
-		});
-	}
-	get name() {
-		return this.model.name;
-	}
-	get sense() {
-		return this.model.sense;
-	}
-	set sense(e) {
-		this.model.sense = e;
-	}
-	get objective() {
-		return this.model.objective;
-	}
-	get variables() {
-		return this.model.variables;
-	}
-	get constraints() {
-		return this.model.constraints;
-	}
-	get status() {
-		return this.model.status;
-	}
-	get objectiveValue() {
-		return this.model.objectiveValue;
-	}
-	addVariable(e) {
-		return this.model.getVarByName(e.name) || (this.model.variables.push(e), this.model.varMapByName.set(e.name, e)), e;
-	}
-	setObjective(e) {
-		if (typeof e == "object" && e) {
-			if ("terms" in e) for (let [t] of e.terms) this.addVariable(t);
-			else e instanceof E && this.addVariable(e);
-		}
-		this.model.setObjective(e);
-	}
-	addConstraint(e, t) {
-		if (typeof e == "object" && e && "lhs" in e) for (let [t] of e.lhs.terms) this.addVariable(t);
-		return this.model.addConstraint(e, t);
-	}
-	add(e, t) {
-		typeof e == "string" ? /(<=|>=|==|=<|=>|=)/.test(e) ? this.addConstraint(e, t) : this.setObjective(e) : e instanceof v ? this.addConstraint(e, t) : this.setObjective(e);
-	}
-	async solve(e, t) {
-		return this.model.optimize(e, t);
-	}
-	toLpString() {
-		return this.model.toLpString();
-	}
-}, Q = "Minimize", $ = "Maximize", he = "Binary", ge = "Integer", _e = "Continuous";
-function ve(e) {
-	return {
-		expr: e,
-		sense: $
-	};
-}
-function ye(e) {
-	return {
-		expr: e,
-		sense: Q
-	};
-}
-var be = class extends Y {
-	constructor(e = "MipModel", t = Q) {
-		super({
-			name: e,
-			sense: t
-		});
-	}
-	add_var(e = {}) {
-		return this.addVar({
-			name: e.name,
-			type: e.var_type,
-			lb: e.lb,
-			ub: e.ub
-		});
-	}
-	set objective_spec(e) {
-		this.setObjective(e.expr, e.sense);
-	}
-	get objective_value() {
-		return this.objectiveValue;
-	}
-	add_constr(e, t) {
-		return this.addConstraint(e, t);
-	}
 };
 //#endregion
-export { he as BINARY, u as BoolExpr, f as BoolOp, d as BoolVar, _e as CONTINUOUS, v as Constraint, q as HiGHS, K as HighsSolver, ge as INTEGER, x as LinExpr, de as LpBinary, le as LpContinuous, ue as LpInteger, ce as LpMaximize, X as LpMinimize, me as LpProblem, Z as LpStatus, fe as LpVariable, $ as MAXIMIZE, Q as MINIMIZE, be as MipModel, J as MockSolver, Y as Model, c as OFFICIAL_HIGHS_CDN_WASM_URL, b as QuadExpr, E as Var, p as and, re as c, ie as expr, j as generateLpString, i as getHighs, o as getHighsWasmUrl, R as getSolver, g as implies, n as isValidWasmHeader, z as listAvailableSolvers, w as lpSum, ve as maximize, ye as minimize, h as not, m as or, O as parseConstraint, D as parseExpr, e as probeWasmUrl, I as registerSolver, a as resolveBrowserWasmUrl, L as setDefaultSolver, s as setHighsWasmUrl, t as solveLpWithHighs, S as sum, pe as value, _ as xor, C as xsum };
+export { d as BoolExpr, p as BoolOp, f as BoolVar, b as Constraint, ie as HiGHS, Q as HighsSolver, C as LinExpr, $ as MockSolver, se as Model, c as OFFICIAL_HIGHS_CDN_WASM_URL, S as QuadExpr, D as Var, m as and, F as c, ee as dot, I as expr, R as generateLpString, i as getHighs, o as getHighsWasmUrl, G as getSolver, _ as implies, n as isValidWasmHeader, K as listAvailableSolvers, l as normalizeSense, g as not, h as or, N as parseConstraint, M as parseExpr, e as probeWasmUrl, U as registerSolver, a as resolveBrowserWasmUrl, W as setDefaultSolver, s as setHighsWasmUrl, t as solveLpWithHighs, w as sum, v as xor };