讲师中心 微信公众号
AI工具推荐 视频效率加速

c++ 稀疏矩阵存储 c++如何实现csr或csc格式的稀疏矩阵

夏芳姑娘_2953

夏芳姑娘_2953

发布时间:2026-03-23 19:24:00

|

789人浏览过

|

来源于php中文网

原创

CSR与CSC的核心区别在于indptr的语义:CSR中indptr[i]表示第i行首个非零元在values中的下标,适合按行访问;CSC中indptr[j]表示第j列首个非零元下标,适合按列访问。

c++ 稀疏矩阵存储 c++如何实现csr或csc格式的稀疏矩阵

CSR 和 CSC 的核心区别在哪?别只看存储顺序

CSR(Compressed Sparse Row)和 CSC(Compressed Sparse Column)本质不是“行优先 vs 列优先”的简单类比,而是**访问模式决定结构**。CSR 适合按行遍历(比如矩阵-向量乘 y = A * x),CSC 适合按列遍历(比如 A.transpose() * x 或求解器中 LU 的列主元选主)。误用会导致缓存失效、循环嵌套反直觉、甚至逻辑错误——比如用 CSR 写一个按列求和的函数,性能可能差 5 倍以上。

二者共享三数组结构:values(非零值)、indices(对应行号或列号)、indptr(行/列起始偏移)。关键差异在 indptr 含义:
– CSR 中 indptr[i] 是第 i 行第一个非零元在 values 中的下标;
– CSC 中 indptr[j] 是第 j 列第一个非零元在 values 中的下标。

手写 CSR 类时,indptr 容易越界或漏填 1 个元素

indptr 长度必须是 rows + 1(CSR)或 cols + 1(CSC),最后一个元素恒为非零元总数 nnz。新手常犯两个错:
– 把 indptr 初始化成长度 rows,导致访问 indptr[rows] 时越界;
– 构造时只填了前 rows 个位置,忘了设 indptr[rows] = nnz,后续遍历某行会读到错误范围。

正确做法:先统计每行非零个数,做前缀和:

std::vector<int> indptr(rows + 1, 0);
for (int i = 0; i < rows; ++i) {
    indptr[i + 1] = indptr[i] + row_nnz[i]; // row_nnz[i] 是第 i 行非零元数
}

注意:所有索引默认从 0 开始,别混用 1-based 的文献伪代码。

插入新元素?别直接 push_back 到 values 和 indices

C++ 稀疏矩阵一旦用 CSR/CSC 格式构建完成,就**不支持高效随机插入**。因为插入会破坏 indptr 的连续性,还可能触发整行数据搬移。常见错误场景:
– 读取 COO 格式(三元组)时边读边插,结果时间复杂度退化成 O(nnz²);
– 在已构造好的 CSR 上调用 insert(i, j, val),内部反复 resize 向量。

正确路径只有两条:
– 批量构建:先收集全部三元组,排序(按行→列对 CSR),再一次性构造三数组;
– 用中间格式过渡:COO → 排序 → CSR。标准库没提供现成排序,得自己写:
std::sort(coo_triplets.begin(), coo_triplets.end(), [](const auto& a, const auto& b) { return a.row < b.row || (a.row == b.row && a.col < b.col); });

用 Eigen 或 Intel MKL?先确认你真需要它们

Eigen 的 SparseMatrix<double, RowMajor> 默认就是 CSR,ColMajor 就是 CSC。但它封装太深:比如想直接访问 indptr,得调 innerIndexPtr()、outerIndexPtr()、valuePtr(),且这些指针只读;MKL 的 sparse_matrix_t 更是要求显式 convert,调试时看不到原始数组。

C++ Code Review Master
C++ Code Review Master

组合式C++代码评审方案,融合静态分析、AI推理、多轮迭代评审和C++专项检查,适用于PR审查、增量代码审查、全项目评审和代码质量评分,触发词包括review cpp、cpp代码评审、C++review、代码审查。

下载

立即学习“C++免费学习笔记(深入)”;

如果你只是做教学、小规模实验、或需深度定制(比如带符号压缩、混合精度索引),手写轻量 CSR 类更可控。示例骨架:

struct CSRMatrix {
    std::vector<double> values;
    std::vector<int> indices;   // 行索引(CSR)或列索引(CSC)
    std::vector<int> indptr;     // 长度 rows+1(CSR)或 cols+1(CSC)
    int rows, cols, nnz;
};

真正难的不是存数据,而是写对 matvec 和 transpose —— 这俩函数里 indptr 和 indices 的配合稍有偏差,结果就全错。建议先用小矩阵(3×4)手算一遍每步下标,再写代码验证。

相关文章

c++速学教程(入门到精通)
c++速学教程(入门到精通)

c++怎么学习?c++怎么入门?c++在哪学?c++怎么学才快?不用担心,这里为大家提供了c++速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

相关标签:

c++

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

热门AI工具

更多
咔片AIPPT

一款在线AI演示文稿制作工具,可根据主题和内容需求辅助生成PPT结构与页面,提高演示材料制作效率。

豆包大模型

豆包大模型是一款由字节跳动推出的企业级大语言模型服务平台。

AionClaw
AionClaw Hot

AionClaw是一款面向办公、创作和编程任务的AI桌面智能体。

蛙蛙写作

一款AI论文写作工具,主要用于超级AI智能写作助手,适合需要提升相关任务效率的用户。

UP简历
UP简历 Hot

一款AI办公效率工具,主要用于基于AI技术的免费在线简历制作工具,适合需要提升相关任务效率的用户。

Loomy
Loomy Hot

一款AI工具,主要用于科大讯飞发布的桌面级 AI 助理,比 OpenClaw 更易用、更安全!,适合需要提升相关任务效率的用户。

DeepSeek

DeepSeek是一款面向对话、写作、编程和推理场景的AI大模型工具。

超级简历WonderCV

一款AI办公效率工具,主要用于免费求职简历模版下载制作,应届生职场人必备简历制作神器,适合需要提升相关任务效率的用户。

WorkBuddy

一款AI办公效率工具,主要用于腾讯云推出的AI原生桌面智能体工作台,适合需要提升相关任务效率的用户。

相关专题

更多
sort排序函数用法
sort排序函数用法

sort排序函数的用法:1、对列表进行排序,默认情况下,sort函数按升序排序,因此最终输出的结果是按从小到大的顺序排列的;2、对元组进行排序,默认情况下,sort函数按元素的大小进行排序,因此最终输出的结果是按从小到大的顺序排列的;3、对字典进行排序,由于字典是无序的,因此排序后的结果仍然是原来的字典,使用一个lambda表达式作为key参数的值,用于指定排序的依据。

1138

2023.09.04

c语言const用法
c语言const用法

const是关键字,可以用于声明常量、函数参数中的const修饰符、const修饰函数返回值、const修饰指针。详细介绍:1、声明常量,const关键字可用于声明常量,常量的值在程序运行期间不可修改,常量可以是基本数据类型,如整数、浮点数、字符等,也可是自定义的数据类型;2、函数参数中的const修饰符,const关键字可用于函数的参数中,表示该参数在函数内部不可修改等等。

2078

2023.09.20

c++怎么把double转成int
c++怎么把double转成int

本专题整合了 c++ double相关教程,阅读专题下面的文章了解更多详细内容。

3648

2025.08.29

C++中int、float和double的区别
C++中int、float和double的区别

本专题整合了c++中int和double的区别,阅读专题下面的文章了解更多详细内容。

684

2025.10.23

C++ 智能指针与现代内存管理
C++ 智能指针与现代内存管理

深入讲解 C++ 现代内存管理的核心工具——智能指针,涵盖 unique_ptr 独占所有权语义、shared_ptr 引用计数机制与循环引用问题、weak_ptr 弱引用的应用场景、make_unique/make_shared 工厂函数的性能优势、自定义删除器的编写、RAII 资源管理思想的实践,以及从裸指针迁移到智能指针的重构策略,帮助开发者编写安全无泄漏的现代 C++ 代码。

339

2026.04.23

FrankenPHP集成Laravel详细教程
FrankenPHP集成Laravel详细教程

本专题提供FrankenPHP集成Laravel的详细配置指南,全面解析运行原理、开发环境搭建、Caddyfile配置、Octane工作模式、数据库连接、队列任务、定时任务和生产环境优化,解决部署过程中常见的报错与兼容性问题。

40

2026.10.08

LLVM自定义Pass怎么写
LLVM自定义Pass怎么写

本专题聚焦LLVM自定义Pass开发,整理Pass类结构、run()方法、PreservedAnalyses、CMake构建、插件注册、-load-pass-plugin加载和测试用例编写流程。

140

2026.09.30

LLVM RISC-V参数配置教程
LLVM RISC-V参数配置教程

本专题介绍LLVM对RISC-V基础ISA和扩展的支持方式,涵盖RV32、RV64、标准扩展、实验性扩展、厂商扩展、-menable-experimental-extensions和版本差异。

120

2026.09.30

LLVM IR中间表示入门指南
LLVM IR中间表示入门指南

本专题整理LLVM IR的核心概念,包括中间表示作用、模块结构、函数、基本块、SSA形式、类型系统和常见语法,帮助新手理解LLVM编译流程中的关键层。

100

2026.09.30

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
Conan 2 Essentials 免费课程
Conan 2 Essentials 免费课程

共0课时 | 0人学习

CMake 与 Conan 集成实践
CMake 与 Conan 集成实践

共0课时 | 0人学习

Conan 2 高级依赖模型介绍
Conan 2 高级依赖模型介绍

共0课时 | 0人学习

关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn