Formulations and efficient numerical solution techniques for transient 3D magneto- and electro-quasistatic field problems