refactor!: completely overhaul the domain type

examples/rosenbrock.rs

@@ -1,7 +1,7 @@
 use gomez::nalgebra as na;
 use gomez::prelude::*;
 use gomez::solver::TrustRegion;
-use na::{Dim, DimName, IsContiguous};
+use na::{DimName, IsContiguous};
 
 // https://en.wikipedia.org/wiki/Rosenbrock_function
 struct Rosenbrock {
@@ -37,7 +37,7 @@ impl System for Rosenbrock {
 
 fn main() -> Result<(), String> {
     let f = Rosenbrock { a: 1.0, b: 1.0 };
-    let dom = Domain::with_dim(f.dim().value());
+    let dom = f.domain();
     let mut solver = TrustRegion::new(&f, &dom);
 
     // Initial guess.

src/analysis/initial.rs

@@ -55,8 +55,10 @@ where
         Sx: StorageMut<F::Scalar, F::Dim> + IsContiguous,
         Sfx: StorageMut<F::Scalar, F::Dim>,
     {
-        let scale_iter = dom.vars().iter().map(|var| var.scale());
-        let scale = OVector::from_iterator_generic(f.dim(), U1::name(), scale_iter);
+        let scale = dom
+            .scale()
+            .map(|scale| OVector::from_iterator_generic(f.dim(), U1::name(), scale.iter().copied()))
+            .unwrap_or_else(|| OVector::from_element_generic(f.dim(), U1::name(), convert(1.0)));
 
         // Compute F'(x) in the initial point.
         f.eval(x, fx)?;

src/core.rs

@@ -4,8 +4,8 @@
 //! traits, optionally specifying the [domain](Domain).
 //!
 //! Algorithms *developers* are interested in implementing the [`Solver`] trait
-//! (or possibly the [`Optimizer`] trait too) and using extension traits (e.g.,
-//! [`VectorDomainExt`]) as well as tools in [derivatives](crate::derivatives).
+//! (or possibly the [`Optimizer`] trait too) as well as tools in
+//! [derivatives](crate::derivatives).
 
 mod base;
 mod domain;

src/core/base.rs

@@ -21,7 +21,7 @@ pub trait Problem {
     /// Get the domain (bound constraints) of the system. If not overridden, the
     /// system is unconstrained.
     fn domain(&self) -> Domain<Self::Scalar> {
-        Domain::with_dim(self.dim().value())
+        Domain::unconstrained(self.dim().value())
     }
 }