As a first approximation beyond linearity, the nonlinear Schrödinger equation reliably describes a broad class of physical systems. Though numerical solutions of this model are well-established, these methods can be computationally complex, especially when system-specific details are incorporated. In this paper, we demonstrate how numerical computations that exploit the features of a graphics processing unit (GPU) result in 40-70× reduction in the time for solutions (depending on hardware details). As a specific case study, we investigate the Gross-Pitaevskii equation, a specific version of the nonlinear Schrödinger model, as it describes a trapped, interacting two-component Bose-Einstein condensate subject to spatially dependent interspin coupling, resulting in an analog to a spin-Hall system. This computational method allows us to probe high-resolution spatial features -revealing an interaction dependent phase transition -all in a reasonable amount of time with readily available hardware.