OpenFOAM 中ldumatrix的spmv实现

OpenFOAM 中的 ldumatrix 类是一个稀疏矩阵类，包含了一些矩阵运算的功能，其中包括了向量-矩阵乘法（spmv）的实现。

ldumatrix 类中的 spmv 函数的代码如下：

template<class Type2>
void spmv(const scalarField& src, Field<Type2>& dst) const
{
    label nRows = rows().size();
    forAll(dst, i)
    {
        dst[i] = 0;
    }

    for (label i = 0; i < nRows; i++)
    {
        const label* rowPtr = rows()[i];
        const scalar* valPtr = vals()[i];
        const label rowSize = size(i);

        for (label j = 0; j < rowSize; j++)
        {
            const label colInd = rowPtr[j];
            const scalar val = valPtr[j];
            dst[colInd] += val*src[i];
        }
    }
}

该函数接受两个参数，一个是源向量（src），一个是目标向量（dst）。它首先将目标向量的所有元素设为0，然后按行遍历稀疏矩阵，对于每行中的非零元素，将该元素与源向量中的对应元素相乘，并加到目标向量的对应元素上。

这个实现比较简单，但是效率不是很高，因为它没有利用稀疏矩阵的特点。在实际使用中，可能需要对该函数进行优化，比如使用稀疏矩阵-向量乘法的算法，以提高计算效率。