  This work presents a fast direct solver strategy allowing full-wave modeling
and dosimetry at terahertz (THz) frequencies. The novel scheme leverages a
preconditioned combined field integral equation together with a regularizer for
its elliptic spectrum to enable its compression into a non-hierarchical
skeleton, invertible in quasi-linear complexity. Numerical results will show
the effectiveness of the new scheme in a realistic skin modeling scenario.
