Compare commits

..

2 Commits

Author SHA1 Message Date
5caccbaf50 improve var creation for SumProdCircuit 2026-04-28 17:29:11 +02:00
0c4bddf3b0 building quotient from ideal to make it more explicit 2026-04-27 23:29:03 +02:00
6 changed files with 20 additions and 30 deletions

View File

@@ -6,10 +6,10 @@ use circuit_cas::circuit::dag::CircuitExt;
fn main() {
let circuit: Rc<RefCell<ProbCircuit>> = ProbCircuit::new();
// Build (x + y) * (x + z)
// vars accept anything that implements Into<StaticVar>: &'static str, (&str, u32), (&str, u32, u32)
let x = circuit.var("x");
let y = circuit.var("y");
let z = circuit.var("z");
let y = circuit.var(("y", 1)); // indexed variable y_1
let z = circuit.var(("z", 0, 1)); // doubly-indexed variable z_{0,1}
let x_plus_y = circuit.var("x") + y;
let x_plus_z = circuit.var("x") + z;
@@ -19,7 +19,7 @@ fn main() {
let x2 = circuit.var("x");
assert_eq!(x.id, x2.id);
println!("(x + y) * (x + z) root node id: {:?}", expr.id);
println!("(x + y_1) * (x + z_{{0,1}}) root node id: {:?}", expr.id);
println!("x node id: {:?}", x.id);
println!("x deduplicated node id: {:?}", x2.id);
}

View File

@@ -3,21 +3,23 @@ use std::rc::Rc;
use circuit_cas::circuit::dag::CircuitExt;
use circuit_cas::circuit::quotient::QuotientCircuit;
use circuit_cas::circuit::traits::Circuit;
use circuit_cas::poly::ideal::{Generators, Ideal};
use circuit_cas::var;
fn main() {
let x = var!("x");
let nx = var!("x\u{0304}");
let idem = vec![
let idem: Ideal<_, Generators> = vec![
1 * (&x ^ 2) - 1 * (&x ^ 1),
1 * (&nx ^ 2) - 1 * (&nx ^ 1),
1 * ((&x ^ 1) * (&nx ^ 1)) - 1 * (&x ^ 1),
];
].into();
let quotient: Rc<RefCell<QuotientCircuit>> = idem.into_iter().collect();
let quotient: Rc<RefCell<QuotientCircuit>> = idem.into();
// Build x * x̄ + x in the DAG
// var accepts anything that implements Into<StaticVar>: &'static str, (&str, u32), (&str, u32, u32)
let xn = quotient.var("x");
let nxn = quotient.var("x\u{0304}");
let prod = xn * nxn;

View File

@@ -60,5 +60,5 @@ pub struct RefNode<C: Circuit> {
pub trait CircuitExt {
type C: Circuit;
type Var;
fn var(&self, v: impl Into<Self::Var>) -> RefNode<Self::C>;
fn var<T: Into<Self::Var>>(&self, v: T) -> RefNode<Self::C>;
}

View File

@@ -43,7 +43,7 @@ impl<V: Var> ProbCircuit<V> {
impl<V: Var> SumProdCircuit for ProbCircuit<V> {
type Var = V;
fn var(&mut self, v: V) -> NodeId { self.node(PNode::Var(v)) }
fn var<T: Into<V>>(&mut self, v: T) -> NodeId { self.node(PNode::Var(v.into())) }
fn add(&mut self, l: NodeId, r: NodeId) -> NodeId {
let (l, r) = if l <= r { (l, r) } else { (r, l) };

View File

@@ -3,7 +3,6 @@ use std::ops::{Deref, DerefMut};
use std::rc::Rc;
use crate::poly::var::{StaticVar, Var};
use crate::poly::flat::Poly;
use crate::poly::ideal::{Generators, GroebnerBasis, Ideal};
use super::dag::{Dag, NodeId};
use super::traits::{Node, SumProdCircuit};
@@ -34,26 +33,15 @@ pub struct QuotientCircuit<V: Var = StaticVar> {
dag: Dag<QNode<V>>,
}
impl<V: Var> QuotientCircuit<V> {
pub fn from_ideal(basis: Ideal<V, GroebnerBasis>) -> Rc<RefCell<Self>> {
Rc::new(RefCell::new(Self { basis, dag: Default::default() }))
}
pub fn from_polys(iter: impl IntoIterator<Item = Poly<V>>) -> Rc<RefCell<Self>> {
let ideal: Ideal<V, Generators> = iter.into_iter().collect();
Rc::new(RefCell::new(Self { basis: ideal.groebner_basis(), dag: Default::default() }))
impl<V: Var> From<Ideal<V, Generators>> for Rc<RefCell<QuotientCircuit<V>>> {
fn from(ideal: Ideal<V, Generators>) -> Self {
ideal.groebner_basis().into()
}
}
impl<V: Var> From<Ideal<V, GroebnerBasis>> for Rc<RefCell<QuotientCircuit<V>>> {
fn from(basis: Ideal<V, GroebnerBasis>) -> Self {
QuotientCircuit::from_ideal(basis)
}
}
impl<V: Var> FromIterator<Poly<V>> for Rc<RefCell<QuotientCircuit<V>>> {
fn from_iter<T: IntoIterator<Item = Poly<V>>>(iter: T) -> Self {
QuotientCircuit::from_polys(iter)
Rc::new(RefCell::new(QuotientCircuit { basis, dag: Default::default() }))
}
}
@@ -69,7 +57,7 @@ impl<V: Var> DerefMut for QuotientCircuit<V> {
impl<V: Var> SumProdCircuit for QuotientCircuit<V> {
type Var = V;
fn var(&mut self, v: V) -> NodeId { self.node(QNode::Var(v)) }
fn var<T: Into<V>>(&mut self, v: T) -> NodeId { self.node(QNode::Var(v.into())) }
fn add(&mut self, l: NodeId, r: NodeId) -> NodeId {
let (l, r) = if l <= r { (l, r) } else { (r, l) };

View File

@@ -29,7 +29,7 @@ where
pub trait SumProdCircuit: Circuit {
type Var;
fn var(&mut self, v: Self::Var) -> NodeId;
fn var<T: Into<Self::Var>>(&mut self, v: T) -> NodeId;
fn add(&mut self, l: NodeId, r: NodeId) -> NodeId;
fn mul(&mut self, l: NodeId, r: NodeId) -> NodeId;
}
@@ -56,8 +56,8 @@ impl<C: SumProdCircuit> CircuitExt for Rc<RefCell<C>> {
type C = C;
type Var = C::Var;
fn var(&self, v: impl Into<C::Var>) -> RefNode<C> {
let id = self.borrow_mut().var(v.into());
fn var<T: Into<C::Var>>(&self, v: T) -> RefNode<C> {
let id = self.borrow_mut().var(v);
RefNode { id, circuit: self.clone() }
}
}