An efficient finite-element formulation for the simulation of ultrashort optical pulses is thoroughly presented. This numerical propagator, capable to take into account spatial and temporal effects, incorporates a wide-angle-band Padé approach and perfectly matched layers at the truncated boundaries. Its effectiveness and higher performance over other similar approaches are clearly demonstrated through key simulations.