RustedSciThe 0.3.29

Rust framework for symbolic and numerical computing: symbolic calculations, nonlinear systems, IVP and BVP for ODE, optimization, fitting and more
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

// Simple optimization: Pre-allocate with capacity
pub fn calc_jacobian_with_capacity(&mut self) {
    let variable_string_vec = self.variable_string.clone();
    let new_jac: Vec<Vec<Expr>> = self
        .vector_of_functions
        .par_iter()
        .enumerate()
        .map(|(i, function)| {
            // Pre-allocate with exact capacity
            let mut vector_of_partial_derivatives = Vec::with_capacity(self.vector_of_variables.len());
            
            for j in 0..self.vector_of_variables.len() {
                let variable = &variable_string_vec[j];
                let list_of_variables_for_this_eq = &self.variables_for_all_disrete[i];
                
                if list_of_variables_for_this_eq.contains(variable) {
                    let mut partial = Expr::diff(function, variable);
                    partial = partial.symplify();
                    vector_of_partial_derivatives.push(partial);
                } else {
                    vector_of_partial_derivatives.push(Expr::Const(0.0));
                }
            }
            vector_of_partial_derivatives
        })
        .collect();
    
    self.symbolic_jacobian = new_jac;
}

// Even simpler: Just avoid the clone
pub fn calc_jacobian_no_clone(&mut self) {
    let new_jac: Vec<Vec<Expr>> = self
        .vector_of_functions
        .par_iter()
        .enumerate()
        .map(|(i, function)| {
            let mut vector_of_partial_derivatives = Vec::new();
            
            for j in 0..self.vector_of_variables.len() {
                let variable = &self.variable_string[j]; // Direct reference
                let list_of_variables_for_this_eq = &self.variables_for_all_disrete[i];
                
                if list_of_variables_for_this_eq.contains(variable) {
                    let mut partial = Expr::diff(function, variable);
                    partial = partial.symplify();
                    vector_of_partial_derivatives.push(partial);
                } else {
                    vector_of_partial_derivatives.push(Expr::Const(0.0));
                }
            }
            vector_of_partial_derivatives
        })
        .collect();
    
    self.symbolic_jacobian = new_jac;
}