This method is available for solving large matrix equations with (Finite Element method) or Finite Differential Element method.

Compared with the entire frequency equalizer, this algorithm has better performance because it avoid the large matrix inversion.

However, when spectral response function is calculated, the results are often not convergent, because large matrix is evolved, and special measure must be taken in order to obtain useful result.