use super::Constraint;
pub struct Rectangle<'a> {
xmin: Option<&'a [f64]>,
xmax: Option<&'a [f64]>,
}
impl<'a> Rectangle<'a> {
pub fn new(xmin: Option<&'a [f64]>, xmax: Option<&'a [f64]>) -> Rectangle<'a> {
assert!(xmin != None || xmax != None);
Rectangle { xmin, xmax }
}
}
impl<'a> Constraint for Rectangle<'a> {
fn project(&self, x: &mut [f64]) {
if let Some(xmin) = &self.xmin {
x.iter_mut().zip(xmin.iter()).for_each(|(x_, xmin_)| {
if *x_ < *xmin_ {
*x_ = *xmin_
};
});
}
if let Some(xmax) = &self.xmax {
x.iter_mut().zip(xmax.iter()).for_each(|(x_, xmax_)| {
if *x_ > *xmax_ {
*x_ = *xmax_
};
});
}
}
}