Struct JacobianOperations 

pub struct JacobianOperations;

Jacobian matrix operations

Implementations

impl JacobianOperations

pub fn compute( functions: &[Expression], variables: &[Symbol], ) -> Vec<Vec<Expression>>

Compute Jacobian matrix for vector-valued function F: ℝⁿ → ℝᵐ

Examples

use mathhook_core::simplify::Simplify;
use mathhook_core::calculus::derivatives::Derivative;
use mathhook_core::calculus::derivatives::PartialUtils;
use mathhook_core::calculus::derivatives::MatrixUtils;
use mathhook_core::{Expression};
use mathhook_core::symbol;
use mathhook_core::calculus::derivatives::JacobianOperations;

let x = symbol!(x);
let y = symbol!(y);
let functions = vec![
    Expression::mul(vec![Expression::symbol(x.clone()), Expression::symbol(y.clone())]),
    Expression::add(vec![Expression::symbol(x.clone()), Expression::symbol(y.clone())])
];
let variables = vec![x, y];
let jacobian = JacobianOperations::compute(&functions, &variables);

pub fn compute_cached( functions: &[Expression], variables: &[Symbol], cache: &mut HashMap<String, Expression>, ) -> Vec<Vec<Expression>>

Compute Jacobian matrix with caching for repeated variable sets

Examples

use mathhook_core::{Expression};
use mathhook_core::symbol;
use mathhook_core::calculus::derivatives::JacobianOperations;
use std::collections::HashMap;

let x = symbol!(x);
let y = symbol!(y);
let functions = vec![
    Expression::symbol(x.clone()),
    Expression::symbol(y.clone())
];
let variables = vec![x, y];
let mut cache = HashMap::new();
let jacobian = JacobianOperations::compute_cached(&functions, &variables, &mut cache);

pub fn transpose( functions: &[Expression], variables: &[Symbol], ) -> Vec<Vec<Expression>>

Compute transpose of Jacobian matrix

Examples

use mathhook_core::{Expression};
use mathhook_core::symbol;
use mathhook_core::calculus::derivatives::JacobianOperations;

let x = symbol!(x);
let y = symbol!(y);
let functions = vec![
    Expression::symbol(x.clone()),
    Expression::symbol(y.clone())
];
let variables = vec![x, y];
let jacobian_t = JacobianOperations::transpose(&functions, &variables);