码农知识堂 - 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

  • 相关阅读:
    DMA+PWM驱动彩色RGB灯
    2023-11-18 Android Linux资源限制命令 ulimit,比如ulimit -d 是设置进程占用的最大数据段大小,默认是unlimited。
    Python 数据分析与挖掘(一)
    【面试专栏】java线程第二篇:Java线程/并发编程带来的问题(含线程题解答)
    面试官:TCC解决方案是什么?如何解决TCC幂等问题?
    Excelize 开源基础库 2.8.0 版本正式发布
    Android基础二:常见的几种组件
    Python解析MDX词典数据并保存到Excel
    Go语言入门【1】数据类型、变量、常量
    Aqwa 带你掌握船舶与海洋工程水动力理论与工程应用
  • 原文地址:https://blog.csdn.net/xhtchina/article/details/125472270
  • 最新文章
  • 沪漂五周年了:我越来越迷茫了
    Agentic Skill Routing 实战:别再把所有 Skill 塞进 AI Agent 上下文
    MySQL-Seconds_behind_master的精度误差
    [MAF预定义ChatClient中间件-03]CachingChatClient——利用缓存省钱省时间
    AI的至暗历史:从万众期待到被政府撤资,AI的两次死亡徘徊
    Agent OS :五种驯服不确定性的范式
    PortSwigger SQL注入LAB11
    数据库即时编译JIT
    [Begin]AI Learn Data Day 0
    深度学习进阶(二十七)现代 LLM 的核心架构设计其二:SwiGLU
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
小工具 小游戏
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1

京公网安备 11010502049817号