Version 1.0.0

Code Download

  • Single
  • Double

MA46 Sparse unsymmetric finite-element system: multifrontal

To solve one or more set of sparse unsymmetric linear equations AX = B from finite-element applications, using a multifrontal elimination scheme. The matrix A must be input by elements and be of the form

A = k=1mA(k)

where A(k) is nonzero only in those rows and columns that correspond to variables of the nodes of the k-th element. Optionally, the user may pass an additional matrix Ad of coefficients for the diagonal. A is then of the form

A = k=1mA(k) + A d

The right-hand side B should be assembled through the summation

B = k=1mB(k),

before calling the solution routine.