可列斯基分解
共轭梯度法
算法
解算器
迭代法
计算机科学
稀疏矩阵
计算复杂性理论
快速傅里叶变换
基质(化学分析)
数学
数学优化
计算化学
特征向量
复合材料
高斯分布
物理
化学
材料科学
量子力学
作者
Phani K. V. V. Nukala,Srđan Šimunović
出处
期刊:Journal of physics
[IOP Publishing]
日期:2003-10-28
卷期号:36 (45): 11403-11412
被引量:35
标识
DOI:10.1088/0305-4470/36/45/004
摘要
The high computational cost involved in modelling of the progressive fracture simulations using large discrete lattice networks stems from the requirement to solve a new large set of linear equations every time a new lattice bond is broken. To address this problem, we propose an algorithm that combines the multiple-rank sparse Cholesky downdating algorithm with the rank-p inverse updating algorithm based on the Sherman–Morrison–Woodbury formula for the simulation of progressive fracture in disordered quasi-brittle materials using discrete lattice networks. Using the present algorithm, the computational complexity of solving the new set of linear equations after breaking a bond reduces to the same order as that of a simple backsolve (forward elimination and backward substitution) using the already LU factored matrix. That is, the computational cost is O(nnz(L)), where nnz(L) denotes the number of non-zeros of the Cholesky factorization L of the stiffness matrix A. This algorithm using the direct sparse solver is faster than the Fourier accelerated preconditioned conjugate gradient (PCG) iterative solvers, and eliminates the critical slowing down associated with the iterative solvers that is especially severe close to the critical points. Numerical results using random resistor networks substantiate the efficiency of the present algorithm.
科研通智能强力驱动
Strongly Powered by AbleSci AI