Inverse Compton scattering (ICS) is a key mechanism for generating X-ray and gamma radiation. The scattering of high-energy electrons by intense laser pulses can significantly enhance photon yield. At high laser amplitudes, Compton scattering transitions to a nonlinear regime, necessitating consideration of the laser’s effect on the electron beam longitudinal component, particularly the light pressure. This study presents simulation methods for a Compton radiation source in both linear and nonlinear regimes. We investigate interactions between electron beams and laser radiation to predict unique radiation source parameters while accounting for various physical effects. A trajectory-based approach derived from the Liénard–Wiechert potentials is introduced. Additionally, we compare modeling results with the quadratic approximation approach for calculating total photon yield. Computations are performed on modern computing architectures (CPU and GPU) using high-performance computing technologies.