basin 0.8.0

Numerical optimization in pure Rust, with pluggable linear-algebra backends and WASM support.
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
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82

//! Backend-agnostic slice views over the parameter vector.
//!
//! L-BFGS-B's inner numerics (cauchy, subsm, formk, compact-form
//! helpers) all operate on `&[f64]` / `&mut [f64]`. To stay generic
//! over the user-chosen parameter backend (`Vec<f64>`, nalgebra
//! `DVector<f64>`, faer `Col<f64>`, ndarray `Array1<f64>`), the
//! top-level solver views each vector as a contiguous f64 slice at
//! iteration boundaries. The two
//! tiny traits below capture that view; impls are per-backend and
//! `pub(crate)` to keep them off the public API surface.
//!
//! This is the slice-extraction layer noted in the L-BFGS-B port plan
//! (Stage 6 in `~/.claude/plans/i-want-to-add-ticklish-breeze.md`).
//! It only matters at the solver / state boundary — once a slice has
//! been extracted, every cauchy / subsm / formk call already works on
//! `&[f64]`.

/// Read-only view of a parameter vector as a contiguous `&[f64]`.
pub(crate) trait AsFloatSlice {
    /// Borrow the underlying storage as a contiguous slice.
    fn as_float_slice(&self) -> &[f64];
}

/// Mutable companion to [`AsFloatSlice`].
pub(crate) trait AsFloatSliceMut: AsFloatSlice {
    /// Borrow the underlying storage as a contiguous mutable slice.
    fn as_float_slice_mut(&mut self) -> &mut [f64];
}

impl AsFloatSlice for Vec<f64> {
    fn as_float_slice(&self) -> &[f64] {
        self.as_slice()
    }
}
impl AsFloatSliceMut for Vec<f64> {
    fn as_float_slice_mut(&mut self) -> &mut [f64] {
        self.as_mut_slice()
    }
}

#[cfg(feature = "nalgebra")]
impl AsFloatSlice for nalgebra::DVector<f64> {
    fn as_float_slice(&self) -> &[f64] {
        self.as_slice()
    }
}
#[cfg(feature = "nalgebra")]
impl AsFloatSliceMut for nalgebra::DVector<f64> {
    fn as_float_slice_mut(&mut self) -> &mut [f64] {
        self.as_mut_slice()
    }
}

#[cfg(feature = "faer")]
impl AsFloatSlice for faer::Col<f64> {
    fn as_float_slice(&self) -> &[f64] {
        self.try_as_col_major()
            .expect("faer::Col<f64> backing storage must be col-major contiguous")
            .as_slice()
    }
}
#[cfg(feature = "faer")]
impl AsFloatSliceMut for faer::Col<f64> {
    fn as_float_slice_mut(&mut self) -> &mut [f64] {
        self.try_as_col_major_mut()
            .expect("faer::Col<f64> backing storage must be col-major contiguous")
            .as_slice_mut()
    }
}

#[cfg(feature = "ndarray")]
impl AsFloatSlice for ndarray::Array1<f64> {
    fn as_float_slice(&self) -> &[f64] {
        self.as_slice().expect("Array1 is contiguous")
    }
}
#[cfg(feature = "ndarray")]
impl AsFloatSliceMut for ndarray::Array1<f64> {
    fn as_float_slice_mut(&mut self) -> &mut [f64] {
        self.as_slice_mut().expect("Array1 is contiguous")
    }
}