Abstract
In this paper, we present a scheme for efficient numerical modeling of operating at millimeter wavelength, e.g., 30 GHz. One of the key features of the method is the GPU acceleration adapted for the Characteristic Basis Function Method (CBFM) acceleration for problems involving layered media, which to the best of our knowledge has not been done in the past.