4个版本 (破坏性)

0.4.0 2022年6月20日
0.3.0 2021年3月18日
0.2.0 2020年9月29日
0.1.0 2020年7月27日

#2196 in 算法

MIT/Apache

475KB
12K SLoC

sprs, Rust的稀疏矩阵

Build status crates.io

sprs在纯Rust中实现了某些稀疏矩阵数据结构和线性代数算法。

API仍在开发中,对其粗糙边缘的反馈非常受欢迎 :)

功能

结构

  • CSR/CSC矩阵
  • 三元组矩阵
  • 稀疏向量

操作

  • 稀疏矩阵/稀疏向量乘积
  • 稀疏矩阵/稀疏矩阵乘积
  • 稀疏矩阵/稀疏矩阵加法、减法
  • 稀疏向量/稀疏向量加法、减法、点积
  • 稀疏/稠密矩阵操作

算法

  • 压缩稀疏矩阵的外迭代器
  • 稀疏向量迭代
  • 稀疏向量联合非零迭代
  • 简单的稀疏Cholesky分解(需要选择LGPL许可)
  • 稀疏三角求解,右侧为稠密矩阵

示例

矩阵构建

  use sprs::{CsMat, CsMatOwned, CsVec};
  let eye : CsMatOwned<f64> = CsMat::eye(3);
  let a = CsMat::new_csc((3, 3),
                         vec![0, 2, 4, 5],
                         vec![0, 1, 0, 2, 2],
                         vec![1., 2., 3., 4., 5.]);

矩阵向量乘法

  use sprs::{CsMat, CsVec};
  let eye = CsMat::eye(5);
  let x = CsVec::new(5, vec![0, 2, 4], vec![1., 2., 3.]);
  let y = &eye * &x;
  assert_eq!(x, y);

矩阵矩阵乘法、加法

  use sprs::{CsMat, CsVec};
  let eye = CsMat::eye(3);
  let a = CsMat::new_csc((3, 3),
                         vec![0, 2, 4, 5],
                         vec![0, 1, 0, 2, 2],
                         vec![1., 2., 3., 4., 5.]);
  let b = &eye * &a;
  assert_eq!(a, b.to_csr());

更多示例,请务必查看热扩散示例。

文档

文档可在docs.rs找到。

变更日志

请参阅变更日志

最低支持的Rust版本

目前最低支持的Rust版本是1.64。在1.0版本之前,提高MSRV将不会被视为破坏性更改,但将尽力避免破坏。

许可证

根据您的选择,许可如下

贡献

除非您明确声明,否则根据Apache-2.0许可证定义,您有意提交给作品中的任何贡献都将双重许可,如上所述,没有任何附加条款或条件。

请参阅贡献指南,了解有关贡献的更多信息。


lib.rs:

随机稀疏矩阵生成

依赖关系

~3MB
~53K SLoC