Using the Penalty Method for Imposing Boundary Conditions in the Pseudospectral Optical Waveguide Mode Solver