The recently introduced beam propagation method using complex Jacobi iteration adapted for modeling of non-paraxial beam propagation in nonlinear optical waveguides is presented in this paper. The beam propagation equation is based on our recently proposed modified Pad,(1,1) approximant operator. The resulting approach is very efficient and well-suited for large structures with long propagation paths
