当一些边界条件值已知时，如何存储矩阵向量的稀疏矩阵？

Question

我有一个表示3D矩形空间的稀疏矩阵。沿着一些边界，我知道价值将是什么（它是一个常数）。其他边界可能是反射的，差异的等等。

我是否应该将问题设置为好像所有边界都说，差异，然后返回并将解决方案向量b中的节点设置为常量？

谢谢！

Answer 1

在有限元方法中，您对Dirchelet（值约束）和Neumann（衍生约束）的处理方式不同。通常先组装矩阵而不考虑边界条件，然后应用边界条件，然后进行LU分解求解。

通过修改汇编矩阵和RHS矢量来应用边界条件。我必须知道更多细节，准确地告诉你你需要做什么。