A computational scheme is developed which permits tractable calculation of three-dimensional full-wave solutions to the Vlasov-Maxwell equations under typical ion cyclotron range of frequencies (ICRF) experimental conditions. The method is unique in that power deposition to the plasma is determined via the anti-Hermitian part of a truncated warm plasma dielectric operator, rather than as the result of an assumed phenomenological collision frequency. The resulting computer code allows arbitrary variation of density, temperature, magnetic field and minority concentration in the poloidal plane by performing a convolution of poloidal modes to produce a coupled system of differential equations in the radial variable. By assuming no inhomogeneity along the toroidal axis, an inverse transform over k|| is performed, yielding the global three-dimensional fast wave field solutions. The application of the code to TFTR-like plasmas shows a mild resonance structure in antenna loading related to the changing number of wavelengths between the antenna and the resonance layer.