The Mixed Finite-Element Method With Mass Lumping for Computing Optical Waveguide Modes