matten 0.20.1

A family-car multidimensional array (tensor) library for small numerical trials / PoCs.
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

//! Gram matrix: G = X · Xᵀ, used in kernel methods and feature covariance.
//!
//! Run: cargo run --example gram_matrix

use matten::Tensor;

fn main() {
    // 4 data points × 3 features
    let x = Tensor::new(
        vec![1.0, 0.0, 1.0, 0.0, 1.0, 1.0, 1.0, 1.0, 0.0, 0.0, 0.0, 1.0],
        &[4, 3],
    );

    // G = X · Xᵀ: shape [4,4]
    let g = x.matmul(&x.transpose());
    println!("Gram matrix shape = {:?}", g.shape()); // [4,4]

    // Diagonal: ||x_i||²
    for i in 0..4 {
        let diag = g.get(&[i, i]).unwrap();
        println!("G[{i},{i}] = {diag:.1}");
    }

    // Symmetric check
    for i in 0..4 {
        for j in 0..4 {
            assert!((g.get(&[i, j]).unwrap() - g.get(&[j, i]).unwrap()).abs() < 1e-10);
        }
    }
    println!("Gram matrix is symmetric: OK");
}