1. A method of reformulating the finite-difference time-domain algorithm for hardware-based accelerators, comprising: splitting fields at every node in a mesh into two subfields; combining total-field and scattered-field formulations into a single, hybrid formulation; and determining sources by means of a lookup table, rather than through direct hardware computations.