solvr 0.2.0

Advanced computing library for real-world problem solving - optimization, differential equations, interpolation, statistics, and more
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

//! Generic distance computation implementations.
//!
//! These are thin wrappers around numr's DistanceOps trait, providing
//! a consistent interface within solvr's spatial module.
use crate::DType;

use crate::spatial::{validate_matching_dims, validate_points_2d, validate_points_dtype};
use numr::error::Result;
use numr::ops::DistanceOps;
use numr::runtime::{Runtime, RuntimeClient};
use numr::tensor::Tensor;

use crate::spatial::traits::distance::DistanceMetric;

/// Generic implementation of cdist (pairwise distances between two point sets).
pub fn cdist_impl<R, C>(
    client: &C,
    x: &Tensor<R>,
    y: &Tensor<R>,
    metric: DistanceMetric,
) -> Result<Tensor<R>>
where
    R: Runtime<DType = DType>,
    C: DistanceOps<R> + RuntimeClient<R>,
{
    validate_points_dtype(x.dtype(), "cdist")?;
    validate_points_dtype(y.dtype(), "cdist")?;
    validate_points_2d(x.shape(), "cdist")?;
    validate_points_2d(y.shape(), "cdist")?;
    validate_matching_dims(x.shape(), y.shape(), "cdist")?;

    // Delegate to numr's DistanceOps
    client.cdist(x, y, metric)
}

/// Generic implementation of pdist (pairwise distances within a point set).
pub fn pdist_impl<R, C>(client: &C, x: &Tensor<R>, metric: DistanceMetric) -> Result<Tensor<R>>
where
    R: Runtime<DType = DType>,
    C: DistanceOps<R> + RuntimeClient<R>,
{
    validate_points_dtype(x.dtype(), "pdist")?;
    validate_points_2d(x.shape(), "pdist")?;

    // Delegate to numr's DistanceOps
    client.pdist(x, metric)
}

/// Generic implementation of squareform (condensed to square).
pub fn squareform_impl<R, C>(client: &C, condensed: &Tensor<R>, n: usize) -> Result<Tensor<R>>
where
    R: Runtime<DType = DType>,
    C: DistanceOps<R> + RuntimeClient<R>,
{
    // Delegate to numr's DistanceOps
    client.squareform(condensed, n)
}

/// Generic implementation of squareform_inverse (square to condensed).
pub fn squareform_inverse_impl<R, C>(client: &C, square: &Tensor<R>) -> Result<Tensor<R>>
where
    R: Runtime<DType = DType>,
    C: DistanceOps<R> + RuntimeClient<R>,
{
    // Delegate to numr's DistanceOps
    client.squareform_inverse(square)
}