symplex 0.2.0

Exact symbolic mathematics for Rust: calculus, summation, solving, linear algebra, transforms, compile-time dimensional analysis, and Rust/C code generation
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
58
59
60
61
62
63

//! SymPy oracle, 0.2 surface — `Ex::compile` (native closure) and
//! `Ex::to_c_fn` (C source) for expressions with special functions, checked
//! against SymPy `N(f(x))` at seeded random points.
//!
//! The C source is not compiled here (no toolchain dependency in the test
//! suite); we check that generation succeeds and mentions the argument, and
//! that the compiled closure agrees with the oracle.

mod v02_oracle_common;

use v02_oracle_common::*;

const KNOWN_BUGS: &[KnownBug] = &[];

#[test]
fn codegen_compile_matches_sympy_n() {
    run_with_known_bugs("codegen", "compile", KNOWN_BUGS, |ctx, fx| {
        let e = match parse(ctx, fx.str("input").unwrap_or("")) {
            Ok(e) => e,
            Err(e) => return Status::NotImplemented(e),
        };
        let vars = fx.str_list("variables");
        let var_refs: Vec<&str> = vars.iter().map(String::as_str).collect();
        let f = match e.compile(&var_refs) {
            Ok(f) => f,
            Err(err) => return Status::NotImplemented(format!("compile: {err}")),
        };
        let pts = fx.eval_points("eval_points");
        if pts.is_empty() {
            return Status::SkippedOracle("no eval points".into());
        }
        let mut checked = 0;
        for pt in &pts {
            let Some(Num::Finite(re, im)) = pt.value else {
                continue;
            };
            if im.abs() > 1e-12 {
                continue; // real-valued closures only
            }
            let args: Vec<f64> = vars
                .iter()
                .map(|v| pt.subs.get(v).copied().unwrap_or(f64::NAN))
                .collect();
            let got = f.call(&args);
            checked += 1;
            // Compiled f64 special functions: 1e-10 relative.
            if !approx_eq_tol(got, re, 1e-10) {
                return Status::Fail(format!("compiled f({args:?}) = {got}, SymPy {re}"));
            }
        }
        if checked == 0 {
            return Status::SkippedOracle("no real-valued oracle points".into());
        }
        // C source generation must succeed for anything we can compile.
        match e.to_c_fn("f", &var_refs) {
            Ok(src) if src.contains("double f(") && vars.iter().all(|v| src.contains(v)) => {
                Status::Pass
            }
            Ok(src) => Status::Fail(format!("to_c_fn produced unexpected source: {src}")),
            Err(err) => Status::NotImplemented(format!("to_c_fn: {err}")),
        }
    });
}