buchberger 0.1.0

Provides polynomial operation and buchberger's algorithm implementation
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

use crate::{Field, Polynomial};
use std::fmt::{Debug, Display};

pub fn groebner<F: Field + Debug + Display>(polys: &mut Vec<Polynomial<F>>) {
    let mut changed = true;
    while changed {
        changed = false;
        'outer: for i in 0..polys.len() {
            for j in 0..i {
                let s = polys[i].s_polynomial(&polys[j]) / &*polys;
                if !s.is_zero() {
                    polys.push(s);
                    changed = true;
                    break 'outer;
                }
            }
        }
    }
}