拡散項を有限要素法で定式化すると，係数行列として拡散行列が現れる．拡散行列は\[{L_{ab}} = \int {\nabla {N_a} \cdot } \nabla {N_b}d\it\Omega \]で定義される．この行列の計算にはかなりの時間が必要であるから，いろいろな高速化手法がある．