Issues Related To The Numerical Implementation Of A Sparse Method For The Solution Of Volume Integral Equations At Low Frequencies