A Rust library for manipulation and evaluation of symbolic integer polynomials.
All of the symbolic expressions have three template types associated with them.
- I - the type that uniquely identifies a single symbolic variable
- C - the type of the free coefficient in every monomial
- P - the type of the power used in every monomial
The main class you most likely will be using is Polynomial<I, C, P>, which
represents any symbolic polynomial. The easiest way to create single variables
(e.g. like a, b, c ...) is by calling variable(I id). The will be a
unique identification of the variable. From there you can use standard
arithmetic operators with both other symbolic expressions and with constants.
If you want to evaluate a symbolic expression, you can call its eval method,
which requires you to specify a mapping from unique identifiers to their assignments.
You can also use automatic deduction to solve a system of equations.
The ordering of both the polynomials and monomials are based on
Graded reverse lexicographic order
derived from the ordering on I. Note that this requires the comparison operators
to be implemented for type I.
All of the sybolic variable implement Display to convert any expression to a humanly
readable format. Additionally the to_code method renders powers as repeated
multiplications, and the output string would look like code snippet.
Just add the dependency in your Cargo.toml file and then import the crate.
There are two features of the crate that can be specified when adding it as a dependency:
-
serialize - Implements serde's serialization/deserialization on all symbolic expressions.
-
repr_c - Adds the attirbute repr(C) to all symbolic expressions, specifically useful for crossing FFI boundaries.
Below is the code the example found in the examples folder.
use std::collections::HashMap;
extern crate symbolic_polynomials;
use symbolic_polynomials::*;
type SymInt = Polynomial<String, i64, u8>;
pub fn main() {
// Create symbolic variables
let a: &SymInt = &variable("a".into());
let b: &SymInt = &variable("b".into());
let c: &SymInt = &variable("c".into());
// Build polynomials
// 5b + 2
let poly1 = 5 * b + 2;
// ab
let poly2 = a * b;
// ab + ac + b + c
let poly3 = a * b + a * c + b + c;
// a^2 - ab + 12
let poly4 = a * a - a * b + 12;
// ac^2 + 3a + bc^2 + 3b + c^2 + 3
let poly5 = a * c * c + 3 * a + b * c * c + 3 * b + c * c + 3;
// floor(a^2, b^2)
let poly6 = floor(a * a, b * b);
// ceil(a^2, b^2)
let poly7 = ceil(a * a, b * b);
// min(ab + 12, ab + a)
let poly8 = min(a * b + 12, a * b + a);
// max (ab + 12, ab + a)
let poly9 = max(a * b + 12, a * b + a);
// max(floor(a^2, b) - 4, ceil(c, b) + 1)
let poly10 = max(floor(a * a, b) - 2, ceil(c, b) + 1);
// (5b + 2)^2
let poly11 = &poly1 * &poly1;
// floor((5b + 2)^2, 5b + 2) = 5b + 2
let poly12 = floor(&poly11, &poly1);
// ceil((5b + 2)^2, 5b + 2) = 5b + 2
let poly13 = ceil(&poly11, &poly1);
// Polynomial printing
let print_function = &|x: String| x;
println!("{}", (0..50).map(|_| "=").collect::<String>());
println!("Displaying polynomials (string representation = code representation):");
println!("{} = {}", poly1, poly1.to_code(print_function));
println!("{} = {}", poly2, poly2.to_code(print_function));
println!("{} = {}", poly3, poly3.to_code(print_function));
println!("{} = {}", poly4, poly4.to_code(print_function));
println!("{} = {}", poly5, poly5.to_code(print_function));
println!("{} = {}", poly6, poly6.to_code(print_function));
println!("{} = {}", poly7, poly7.to_code(print_function));
println!("{} = {}", poly8, poly8.to_code(print_function));
println!("{} = {}", poly9, poly9.to_code(print_function));
println!("{} = {}", poly10, poly10.to_code(print_function));
println!("(5b + 2)^2 = {}", poly11);
println!("floor((5b + 2)^2, 5b + 2) = {}", poly12);
println!("ceil((5b + 2)^2, 5b + 2) = {}", poly13);
println!("{}", (0..50).map(|_| "=").collect::<String>());
// Polynomial evaluation
let values = &mut HashMap::<String, i64>::new();
values.insert("a".into(), 3);
values.insert("b".into(), 2);
values.insert("c".into(), 5);
println!("Evaluating for a = 3, b = 2, c = 5.");
println!("{} = {} [Expected 12]", poly1, poly1.eval(values).unwrap());
println!("{} = {} [Expected 6]", poly2, poly2.eval(values).unwrap());
println!("{} = {} [Expected 28]", poly3, poly3.eval(values).unwrap());
println!("{} = {} [Expected 15]", poly4, poly4.eval(values).unwrap());
println!("{} = {} [Expected 168]", poly5, poly5.eval(values).unwrap());
println!("{} = {} [Expected 2]", poly6, poly6.eval(values).unwrap());
println!("{} = {} [Expected 3]", poly7, poly7.eval(values).unwrap());
println!("{} = {} [Expected 9]", poly8, poly8.eval(values).unwrap());
println!("{} = {} [Expected 18]", poly9, poly9.eval(values).unwrap());
println!("{} = {} [Expected 4]", poly10, poly10.eval(values).unwrap());
println!("{} = {} [Expected 144]", poly11, poly11.eval(values).unwrap());
println!("{} = {} [Expected 12]", poly12, poly12.eval(values).unwrap());
println!("{} = {} [Expected 12]", poly13, poly13.eval(values).unwrap());
println!("{}", (0..50).map(|_| "=").collect::<String>());
// Variable deduction
values.insert("a".into(), 5);
values.insert("b".into(), 3);
values.insert("c".into(), 8);
let implicit_values = &vec![
(&poly1, poly1.eval(values).unwrap()),
(&poly2, poly2.eval(values).unwrap()),
(&poly3, poly3.eval(values).unwrap()),
];
let deduced_values = deduce_values(implicit_values).unwrap();
println!("Deduced values:");
println!("a = {} [Expected 5]", deduced_values["a"]);
println!("b = {} [Expected 3]", deduced_values["b"]);
println!("c = {} [Expected 8]", deduced_values["c"]);
println!("{}", (0..50).map(|_| "=").collect::<String>());
}The output of the program:
==================================================
Displaying polynomials (string representation = code representation):
5b + 2 = 5 * b + 2
ab = a * b
ab + ac + b + c = a * b + a * c + b + c
a^2 - ab + 12 = a * a - a * b + 12
ac^2 + 3a + bc^2 + 3b + c^2 + 3 = a * c * c + 3 * a + b * c * c + 3 * b + c * c + 3
floor(a^2, b^2) = floor(a * a, b * b)
ceil(a^2, b^2) = ceil(a * a, b * b)
min(ab + 12, ab + a) = min(a * b + 12, a * b + a)
max(ab + 12, ab + a) = max(a * b + 12, a * b + a)
max(floor(a^2, b) - 2, ceil(c, b) + 1) = max(floor(a * a, b) - 2, ceil(c, b) + 1)
(5b + 2)^2 = 25b^2 + 20b + 4
floor((5b + 2)^2, 5b + 2) = 5b + 2
ceil((5b + 2)^2, 5b + 2) = 5b + 2
==================================================
Evaluating for a = 3, b = 2, c = 5.
5b + 2 = 12 [Expected 12]
ab = 6 [Expected 6]
ab + ac + b + c = 28 [Expected 28]
a^2 - ab + 12 = 15 [Expected 15]
ac^2 + 3a + bc^2 + 3b + c^2 + 3 = 168 [Expected 168]
floor(a^2, b^2) = 2 [Expected 2]
ceil(a^2, b^2) = 3 [Expected 3]
min(ab + 12, ab + a) = 9 [Expected 9]
max(ab + 12, ab + a) = 18 [Expected 18]
max(floor(a^2, b) - 2, ceil(c, b) + 1) = 4 [Expected 4]
25b^2 + 20b + 4 = 144 [Expected 144]
5b + 2 = 12 [Expected 12]
5b + 2 = 12 [Expected 12]
==================================================
Deduced values:
a = 5 [Expected 5]
b = 3 [Expected 3]
c = 8 [Expected 8]
==================================================You can check out the tests in the tests folder for more examples.
Currently, the automatic deduction for solving system of equations is pretty limited. The main reason is that for the purposes that the project has been developed it is sufficient. A more powerful and complete algorithm would probably use Grobner basis.
Symbolic Polynomials is distributed under the terms of both the MIT license and the Apache License (Version 2.0). See LICENSE-APACHE and LICENSE-MIT for details.