Text this: Calculation of normal mode spectra in laterally heterogeneous earth models using an iterative direct solution method