Fast and accurate solutions of low-frequency electromagnetics problems are obtained with an iterative solver based on hierarchical matrices. Iterative solvers require matrix-vector multiplications (MVMs). The results show significant reductions both in CPU time and memory consumption compared to the O(N2) complexity of ordinary matrix filling and MVM.