8.4 Symbolic Factorization of Sparse Matrices

The goal of symbolic factorization is to define the sparsity pattern of the LU decomposition of a sparse matrix A. Recall that

LU = PAQ

where P and Q are row and column permutations that reflect the pivot strategy associated with the factorization process.

8.4.1 Symbolic Factorization with Minimum Degree Pivot

The goal of the current section is somewhat more specific. It describes a symbolic factorization algorithm that simulates the decomposition of A when a minimum degree pivot strategy is applied. The algorithm operates on an undirected graph G = (V, E) whose vertices V are labeled from 1 to |V|\left|V\right|. G is defined by its adjacency list A. The algorithm arrives at a symbolic factorization by eliminating vertices from G. Each stage of the process creates a reduced graph G ′ = (V ′, E ′). The vertex v with the minimum degree in G ′ is chosen as the next elimination candidate. The algorithm describes the structure of L, U, P, and Q in the following manner:

  • The adjacency list A is augmented to account for fill-ups that will occur during numeric factorization.

  • A list L is created. The list orders the vertices in V according to a minimum degree pivot strategy. If the vertices of V are labeled according to their position in LU, a minimum degree ordering of V is generated.

  • A mapping ψ is created. The domain of ψ is the initial label of vertex v. The range of ψ is the minimum degree label of v. That is, ψ is the minimum degree label of v.

  • A mapping τ is created. The domain of τ is the the minimum degree label of vertex v. The range of τ is the initial label of v. That is, if v¯\overline{v} is the minimum degree label of a vertex, τ(v¯)\tau(\overline{v}) is the initial label of the vertex.

One field associated with edge eij in the adjacency list contains a binary variable indicating whether edge (i,j) is a fill-up or not. The fill-up indicator is communicated through the buffer e in the normal manner.

Algorithm 22 computes a symbolic factorization of G that consists of these four data structures.

Algorithm 22: Symbolic Factorization of a Sparse Matrix
for i=1,,|V|i=1,\cdots,|V|
    v=v= minimum_degree(VV^{\prime})
    while[w=w=iterate(AA,vertex v,1,|V|,ev,1,|V|,e)] \neq finished
       if in_graph(V,wV^{\prime},w)
           decrease_degree(V,wV^{\prime},w)
           push_iteration
           while[z=z=iterate(AA,vertex v,w+1,|V|,ev,w+1,|V|,e)] \neq finished
               e.fillup=e.fillup=true
               if insert(A,z,w,eA,z,w,e)
                  increase_degree(V,zV^{\prime},z)
                  insert(A,w,z,eA,w,z,e)
                  increase_degree(V,wV^{\prime},w)
           pop_iteration
    remove(V,vV^{\prime},v)
    map(ψ,v,i\psi,v,i)
    map(τ,i,v\tau,i,v)
    link(L,vL,v)
† Operations and data types are defined in Section 8.2.

A few observations concerning the factorization procedure:

  • The reduced graph G ′ is tracked in data structures designed to model vertex elimination–see Section 8.2.3 of this monograph and Section 8 (Vertex Elimination) of its companion Graph Algorithms for implementation details.

  • The adjacency list A is augmented to log fill-ups but it is never diminished to reflect graph reduction.

  • If the algorithm is efficiently implemented, ψ, τ, and l can occupy space that is vacated as G ′ shrinks.

8.4.2 Computational Complexity of Symbolic Factorization

The complexity of the symbolic factorization is determined by the size of the adjacency list. One iteration of the main loop requires e2-e2\frac{{e}^{2}-e}{2} adjacency list accesses, where e is the number of edges incident upon vertex v. If all vertices in V have the same number of incident edges, the total operation count is |V|e2-e2\left|V\right|\frac{{e}^{2}-e}{2}. In this case, the time complexity of symbolic factorization is O(|V|e2)O\left(\left|V\right|{e}^{2}\right) if all operations on the data structures are O(1).

Consider the following comparison of symbolic factorization operation counts to those associated with the LU decomposition of a dense matrix. If a graph has 100 vertices each of which has 5 incident edges, computing the symbolic factorization requires 1000 operations. Computing the dense factorization requires 661,650 operations. If the number of vertices increases to 1000 while the incident edges per vertex remains constant (the typical scenario in network analysis problems associated with electric power systems), symbolic factorization requires 1 × 104 operations and dense matrix decomposition requires 6 × 108 operations.