数值风云

OpenFOAM耦合ARPACK++计算软件包

我想利用ARPACK++中的一些函数处理矩阵的特征值和特征向量的问题,在open foam的求解器中使用ARPACK++库中相关函数,现在问题是不知道该如何在linux中将两周建立联系,感谢大家!

可能你需要了解fvMatrix,他是存离散方程的,比如

fvScalarMatrix TEqn
(
    fvm::ddt(T) - fvm::laplacian(DT,T) == fvOptions(T)
)

fvMatrix继承自lduMatrix,有各种各样的getter函数,比如diag,upper,lower,source等,可以提取线性方程组对应的矩阵元素
可能你需要了解一些lduMatrix的存储细节

矩阵有了之后,无非就是转成需要的格式然后用其他工具处理