矩阵列拼接_线性矩阵求解器（一）：稀疏矩阵

扩散项离散的讨论大致就到上篇位置，这篇文章开始新的部分：线性方程组求解器。

CFD计算主要工作还是求解线性化后的动量、压力方程。方程求解的变量依然用

表示，那么根据之前介绍的扩散项离散结构，含有扩散项会引入相邻节点的

值。线性化的动量方程中含有扩散项、对流项、时间通量、表面力、体积力等等，但这些项大致都会用临近节点的值表达出来。这种差分方式叫差分模板（stencil）。对于一个节点的方程，差分模板所需要的信息只有临近的几个点而已。所有

个节点一起列出方程时，形成的方程组

中每一行实际上只有几个非零元素。比如，在100个网格上求解对流扩散方程，而对流项选取的是三点差分模板，那么100*100的矩阵中每行只有三个非零元素。

为了储存具有大量零元素的矩阵，诞生了稀疏矩阵这一技术。稀疏矩阵只存储非零元素及其相应位置信息，本文以为BLAS^[1]和Intel MKL开发者手册^[2]为参考，介绍通用COO/CSR/DIA格式。稀疏矩阵可能因为太基础了，反倒网上很少有人讨论，大部分都是用现成轮子。

1.COO

意为坐标表示。COO原理很简单，用一维数组把矩阵