工程电磁场数值计算4(有限元法1)剖析.pptVIP

2021/2/11 * 稀疏矩阵技术 没有稀疏矩阵技术就没有有限元的成功。 带状矩阵技术 通过合适的节点编码,可使系统矩阵的非零元素集中于主对角线附近的带形区域内。在使用LU分解法求解方程组的过程中,带形区域以外的0元素无需计算。 缺点:节点编码优化;带形区域内仍然有大量零元素 2021/2/11 * 非零元素存储技术 只存储矩阵的非零元素。一般用一个一维数组存储矩阵的非零元素,用另一个索引数组存储这些非零元素在原矩阵中的行和列值。例如: 2021/2/11 * 非零元素存储技术 在迭代法中,系统矩阵参与的运算只有矩阵左乘以某个相量。如果采用上面的存储方法,则实现 q = A p 的算法如下: C Compute the matrix–vector product DO k=1, not i=IA(k) j=JA(k) q(i)=q(i)+A(k)*p(j) ENDDO 2021/2/11 * 非零元素存储技术 下列的二维存储方案是我的发明:直接将矩阵挤扁。 数组A有n行m列,m是A各行非零元素个数的最大值。对角元素放在该行的第一个位置上,且不管是否为0,都要存贮。用一个整型数组JA存放各元素的列号。JA中第一列可以用来记录每行的非零元素个数。 2021/2/11 * 本节更多的参考文献: 金建