oxiphysics-core 0.1.0

Core types, traits, and abstractions for the OxiPhysics engine
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

//! # LinearConstraint - Trait Implementations
//!
//! This module contains trait implementations for `LinearConstraint`.
//!
//! ## Implemented Traits
//!
//! - `IntervalConstraint`
//!
//! 🤖 Generated with [SplitRS](https://github.com/cool-japan/splitrs)

use super::functions::IntervalConstraint;
#[allow(unused_imports)]
use super::functions::*;
use super::types::{Interval, LinearConstraint};

impl IntervalConstraint for LinearConstraint {
    fn propagate(&self, box_: &mut [Interval]) -> bool {
        let mut changed = false;
        for &(target_idx, target_coeff) in &self.coeffs {
            if target_coeff.abs() < 1e-30 {
                continue;
            }
            let mut other_sum = Interval::point(0.0);
            for &(idx, coeff) in &self.coeffs {
                if idx != target_idx {
                    other_sum = other_sum + box_[idx] * coeff;
                }
            }
            let needed = (self.rhs - other_sum) / Interval::point(target_coeff);
            if let Some(narrowed) = Interval::intersection(box_[target_idx], needed)
                && (narrowed.lo > box_[target_idx].lo + 1e-15
                    || narrowed.hi < box_[target_idx].hi - 1e-15)
            {
                box_[target_idx] = narrowed;
                changed = true;
            }
        }
        changed
    }
}