oxieml 0.1.3

EML operator: all elementary functions from exp(x) - ln(y)
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

//! Python bindings for symbolic equation solving.

use pyo3::prelude::*;

fn map_err(e: impl std::fmt::Display) -> PyErr {
    pyo3::exceptions::PyValueError::new_err(e.to_string())
}

/// Find all real symbolic solutions of `expr_str = 0` for the given variable.
///
/// Parameters
/// ----------
/// expr_str : str
///     Expression string for the left-hand side (set equal to zero).
/// var : int
///     Variable index (0-based) to solve for.
///
/// Returns
/// -------
/// list of str
///     Each element is a LaTeX string for a symbolic root.
#[pyfunction]
pub fn solve_for_all_py(expr_str: &str, var: usize) -> PyResult<Vec<String>> {
    let tree = crate::parse(expr_str).map_err(map_err)?;
    let lowered = tree.lower().simplify();
    let zero = crate::LoweredOp::Const(0.0);
    let result = crate::solve_for_all(&lowered, &zero, var).map_err(map_err)?;
    Ok(result.roots.iter().map(|r| r.to_latex()).collect())
}

/// Find all complex roots of the polynomial `expr_str` in variable `var`.
///
/// Parameters
/// ----------
/// expr_str : str
///     Expression string for a polynomial expression.
/// var : int
///     Variable index (0-based) to extract the polynomial in.
///
/// Returns
/// -------
/// list of tuple[float, float]
///     Each element is `(re, im)` for a complex root.
#[pyfunction]
pub fn solve_polynomial_complex_py(expr_str: &str, var: usize) -> PyResult<Vec<(f64, f64)>> {
    let tree = crate::parse(expr_str).map_err(map_err)?;
    let lowered = tree.lower().simplify();
    let poly = crate::Poly::from_lowered(&lowered, var)
        .map_err(|e| pyo3::exceptions::PyValueError::new_err(e.to_string()))?;
    let complex_roots = crate::solve_poly::solve_polynomial_complex(&poly)
        .map_err(|e| pyo3::exceptions::PyValueError::new_err(e.to_string()))?;
    Ok(complex_roots
        .roots
        .into_iter()
        .map(|c| (c.re, c.im))
        .collect())
}