Zero modes in magnetic systems: General theory and an efficient computational scheme