码农知识堂 - 1000bd
  •   Python
  •   PHP
  •   JS/TS
  •   JAVA
  •   C/C++
  •   C#
  •   GO
  •   Kotlin
  •   Swift
  • Eigen稀疏矩阵操作


    在许多应用中(例如,有限元方法),通常处理非常大的矩阵,其中只有少数系数不为零。 在这种情况下,可以通过使用仅存储非零系数的专用表示来减少内存消耗并提高性能。 这样的矩阵称为稀疏矩阵。

    1. 稀疏矩阵模块

    模块

    头文件

    内容

    SparseCore

    #include <Eigen/SparseCore>

    SparseMatrix 和 SparseVector 类、矩阵组装、

    基本稀疏线性代数(包括稀疏三角求解器)

    SparseCholesky

    #include <Eigen/SparseCholesky>

    直接用稀疏 LLT 和 LDLT Cholesky 分解来解决稀疏自伴随正定问题

    SparseLU

    #include<Eigen/SparseLU>

    稀疏 LU 分解来解决一般平方稀疏系统

    SparseQR

    #include<Eigen/SparseQR

    用于解决稀疏线性最小二乘问题的稀疏 QR 分解

    IterativeLinearSolvers

    #include <Eigen/IterativeLinearSolvers>

    求解大型一般线性平方问题的迭代求解器(包括自伴随正定问题)

    Sparse

    #include <Eigen/Sparse>

    包含以上所有模块

    2. 稀疏矩阵格式

    SparseMatrix 类是 Eigen 稀疏模块的主要稀疏矩阵表示; 它提供高性能和低内存使用率。 它实现了广泛使用的压缩列(或行)存储方案的更通用的变体。 它由四个紧凑的数组组成:

    Values:存储非零系数值。
    InnerIndices:存储非零的行(或列)索引。
    OuterStarts:为每一列(或行)存储前两个数组中第一个非零的索引。
    InnerNNZs:存储每列(相应行)的非零数。 inner指的是一个内部向量,它是列主矩阵的列,或行主矩阵的行。 外部这个词指的是另一个方向。

    目前,给定内部向量的元素保证总是通过增加内部索引进行排序。 “_”表示可以快速插入新元素的可用空间。 假设不需要重新分配,随机元素的插入因此在 O(nnz_j) 中,其中 nnz_j 是相应内部向量的非零数。 另一方面,在给定的内部向量中插入具有增加内部索引的元素效率更高,因为这仅需要增加相应的 InnerNNZs 条目,即 O(1) 操作。

    没有可用空间的情况是一种特殊情况,称为压缩模式。 它对应于广泛使用的压缩列(或行)存储方案(CCS 或 CRS)。 通过调用 SparseMatrix::makeCompressed() 函数,可以将任何 SparseMatrix 转换为这种形式。 在这种情况下,可以注意到 InnerNNZs 数组与 OuterStarts 是冗余的,因为我们等式:InnerNNZs[j] = OuterStarts[j+1]-OuterStarts[j]。 因此,实际上调用 SparseMatrix::makeCompressed() 会释放此缓冲区。

    参考文献

    Eigen关于稀疏矩阵_爱数学的虹猫的博客-CSDN博客_eigen sparsematrix

    Eigen教程3 - 稀疏矩阵操作_xuezhisdc的博客-CSDN博客_eigen 稀疏矩阵

    Eigen教程5 - 求解稀疏线性方程组_xuezhisdc的博客-CSDN博客_eigen 线性方程组

    Eigen: Solving Sparse Linear Systems 

    Eigen: Quick reference guide for sparse matrices

  • 相关阅读:
    大模型解决方案:具体业务场景下的智能表单填充(附代码)
    用html编写的招聘简历
    QT--Qt的QScrollArea的使用核心
    正则校验的多种方式
    网络协议--动态选路协议
    一文讲透,商业智能BI未来发展趋势如何
    得了胆囊结石怎么办?
    Spring
    Windows【工具 04】WinSW官网使用说明及实例分享(将exe和jar注册成服务)实现服务器重启后的服务自动重启
    【1470. 重新排列数组】
  • 原文地址:https://blog.csdn.net/xhtchina/article/details/125472270
  • 最新文章
  • 【JVM】编译执行与解释执行的区别是什么?JVM 使用哪种方式?
    用 Hashids 优雅解决 C 端自增 ID 暴露问题
    V8引擎 精品漫游指南--Ignition篇(上) 指令 栈帧 槽位 调用约定 内存布局 基础内容
    LLVM Pass快速入门(四):代码插桩
    milkup:桌面端 markdown AI续写和即时渲染
    基于项目工程构建SBOM(软件物料清单)的研究
    鸿蒙应用开发UI基础第二节:鸿蒙应用程序框架核心解析与实操
    .NET 中如何快速实现 List 集合去重?
    扣子Coze实战:从0到1打造抖音+小红书热点监控智能体
    浅谈数据访问层
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
小工具 小游戏
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1

京公网安备 11010502049817号