SUMMARYThe multigrid method widely used in computational fluid dynamics is applied to the problem of propagation of an optical wave in an optical waveguide. This method is a technique for accelerating the convergence speed of the relaxation method. The electromagnetic field components contain various frequency components. Higher-frequency components contain local effects of the shape. On the other hand, the lower-frequency components contain more global aspects. By using this phenomenon, the discretized equations are made to reduce errors corresponding to the grid size. As a numerical example, linear and nonlinear directional couplers are treated and the effects of the number of levels and the number of iterations in the multigrid method are studied.